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,081 records · Page 60Linked to original sources

PubMed2002

Molecular modeling as a powerful technique for understanding small-large molecules interactions.

In the present review we summarize recent work, aimed at a better understanding of the interactions in macromolecule ligand complexes, performed by means of computational tools such as pseudoreceptor generation, molecular docking, conformational search and energy minimization. While the first approach has been applied when the three-dimensional structural properties of the biological target were unknown, the remaining protocols exploited the knowledge of the overall structure of the involved macromolecules and their active sites. Molecular modeling techniques were used in the cases reported to study and propose macromolecular binding sites and to predict their interactions with bioactive conformers of the ligands.

Anti-HIV Agents↗

PubMed2001

A putative three-dimensional arrangement of the human serotonin transporter transmembrane helices: a tool to aid experimental studies.

The human serotonin transporter is the molecular target for selective serotonin reuptake inhibitor drugs which are being used for treatment of depression. A three-dimensional model of the membrane spanning parts of the transporter was constructed. The transporter was assumed to consist of 12 transmembrane alpha-helices. The model was based on published experimental data of cocaine binding to mutant transporters, amino acid sequence analysis, and interactive molecular graphics. The model suggests that a high affinity cocaine binding site is situated in a region of the model where Asp98 acts like an anchor, while a putative low affinity site is situated in another region with Glu508 as the anchoring amino acid. A series of docking experiments with various reuptake inhibitors were conducted, using interactive molecular graphics techniques combined with energy calculations and analysis of the transporter-ligand complexes. Experiments involving molecular mapping of ligand binding areas may benefit from using the current model in experimental design. From the current model, several amino acids were proposed as prime candidates for mutagenesis and subsequent ligand binding studies. Also for evaluation of results from site directed mutagenesis experiments with SERT and similar transporters we assume the model will be helpful.

Binding Sites↗

PubMed1999

Modulation of fibroblast growth factor-2 receptor binding, signaling, and mitogenic activity by heparin-mimicking polysulfonated compounds.

Basic fibroblast growth factor (FGF-2) interacts with high-affinity tyrosine-kinase fibroblast growth factor receptors (FGFRs) and low-affinity heparan sulfate proteoglycans (HSPGs) in target cells. Both interactions are required for FGF-2-mediated biological responses. Here we report the FGF-2 antagonist activity of novel synthetic sulfonic acid polymers with distinct chemical structures and molecular masses (MMs). PAMPS [poly(2-acrylamido-2-methyl-1-propanesulfonic acid)], (MM approximately 7,000-10,000), PAS [poly(anetholesulfonic acid)], (MM approximately 9,000-11,000), PSS [poly(4-styrenesulfonic acid)], (MM = 70,000), and poly(vinylsulfonic acid) (MM = 2,000), inhibited FGF-2 binding to HSPGs and FGFRs in fetal bovine aortic endothelial GM 7373 cells. They also abrogated the formation of the HSPG/FGF-2/FGFR ternary complex, as evidenced by their capacity to prevent FGF-2-mediated cell-cell attachment of FGFR-1-overexpressing, HSPG-deficient Chinese hamster ovary cells to wild-type HSPG-bearing cells. Direct interaction of the polysulfonates with FGF-2 was demonstrated by their ability to protect the growth factor from proteolytic cleavage. Accordingly, molecular modeling, based on the crystal structure of the interaction of FGF-2 with a heparin hexamer, showed the feasibility of docking PAMPS into the heparin-binding domain of FGF-2. In agreement with their FGF-2-binding capacity, PSS, PAS, and PAMPS inhibited FGF-2-induced cell proliferation in GM 7373 cells and murine brain microvascular endothelial cells. The antiproliferative activity of these compounds was associated with the abrogation of FGF-2-induced tyrosine phosphorylation of FGFR-1. Moreover, the polysulfonates PSS and PAS inhibited FGF-2-induced activation of mitogen-activated protein kinase-1/2, involved in FGF-2 signal transduction. In conclusion, sulfonic acid polymers bind FGF-2 by mimicking heparin interaction. These compounds may provide a tool to inhibit FGF-2-induced endothelial cell proliferation in angiogenesis and tumor growth.

Animals↗

PubMed2005

Evaluation of aldose reductase inhibition and docking studies of some secondary metabolites, isolated from Origanum vulgare L. ssp. hirtum.

Five polar constituents of Origanum vulgare L. ssp. hirtum were investigated for their ability to inhibit aldose reductase (ALR2), the first enzyme of the polyol pathway implicated in the secondary complications of diabetes. The most active compound was found to be lithospermic acid B. Caffeic acid was inactive as it showed no inhibitory activity against the enzyme. The order of the inhibitory activity of the remaining compounds was: rosmarinic acid >12-hydroxyjasmonic acid 12-O-beta-glucopyranoside > p-menth-3-ene-1,2-diol 1-O-beta-glucopyranoside. Docking studies have been undertaken to gain insight into the binding mode of the investigated compounds at the active site of ALR2. The predicted hydrogen bonding and hydrophobic interactions may explain the observed inhibitory activity.

Aldehyde Reductase↗

PubMed2006

Phytoestrogens as inhibitors of the human progesterone metabolizing enzyme AKR1C1.

Phytoestrogens are plant-derived, non-steroidal constituents of our diets. They can act as agonists or antagonists of estrogen receptors, and they can modulate the activities of the key enzymes in estrogen biosynthesis. Much less is known about their actions on the androgen and progesterone metabolizing enzymes. We have examined the inhibitory action of phytoestrogens on the key human progesterone-metabolizing enzyme, 20alpha-hydroxysteroid dehydrogenase (AKR1C1). This enzyme inactivates progesterone and the neuroactive 3alpha,5alpha-tetrahydroprogesterone, to form their less active counterparts, 20alpha-hydroxyprogesterone and 5alpha-pregnane-3alpha,20alpha-diol, respectively. We overexpressed recombinant human AKR1C1 in Escherichia coli, purified it to homogeneity, and examined the selected phytoestrogens as inhibitors of NADPH-dependent reduction of a common AKR substrate, 9,10-phenantrenequinone, and progesterone. The most potent inhibitors were 7-hydroxyflavone, 3,7-dihydroxyflavone and flavanone naringenin with IC(50) values in the low microM range. Docking of the flavones in the active site of AKR1C1 revealed their possible binding modes, in which they are sandwiched between the Leu308 and Trp227 of AKR1C1.

20-Hydroxysteroid Dehydrogenases↗

PubMed2003

Automated docking of ligands to an artificial active site: augmenting crystallographic analysis with computer modeling.

The W191G cavity of cytochrome c peroxidase is useful as a model system for introducing small molecule oxidation in an artificially created cavity. A set of small, cyclic, organic cations was previously shown to bind in the buried, solvent-filled pocket created by the W191G mutation. We docked these ligands and a set of non-binders in the W191G cavity using AutoDock 3.0. For the ligands, we compared docking predictions with experimentally determined binding energies and X-ray crystal structure complexes. For the ligands, predicted binding energies differed from measured values by +/- 0.8 kcal/mol. For most ligands, the docking simulation clearly predicted a single binding mode that matched the crystallographic binding mode within 1.0 A RMSD. For 2 ligands, where the docking procedure yielded an ambiguous result, solutions matching the crystallographic result could be obtained by including an additional crystallographically observed water molecule in the protein model. For the remaining 2 ligands, docking indicated multiple binding modes, consistent with the original electron density, suggesting disordered binding of these ligands. Visual inspection of the atomic affinity grid maps used in docking calculations revealed two patches of high affinity for hydrogen bond donating groups. Multiple solutions are predicted as these two sites compete for polar hydrogens in the ligand during the docking simulation. Ligands could be distinguished, to some extent, from non-binders using a combination of two trends: predicted binding energy and level of clustering. In summary, AutoDock 3.0 appears to be useful in predicting key structural and energetic features of ligand binding in the W191G cavity.

Automation↗

PubMed1999

Structural dynamics of ligand diffusion in the protein matrix: A study on a new myoglobin mutant Y(B10) Q(E7) R(E10).

A triple mutant of sperm whale myoglobin (Mb) [Leu(B10) --> Tyr, His(E7) --> Gln, and Thr(E10) --> Arg, called Mb-YQR], investigated by stopped-flow, laser photolysis, crystallography, and molecular dynamics (MD) simulations, proved to be quite unusual. Rebinding of photodissociated NO, O2, and CO from within the protein (in a "geminate" mode) allows us to reach general conclusions about dynamics and cavities in proteins. The 3D structure of oxy Mb-YQR shows that bound O2 makes two H-bonds with Tyr(B10)29 and Gln(E7)64; on deoxygenation, these two residues move toward the space occupied by O2. The bimolecular rate constant for NO binding is the same as for wild-type, but those for CO and O2 binding are reduced 10-fold. While there is no geminate recombination with O2 and CO, geminate rebinding of NO displays an unusually large and very slow component, which is pretty much abolished in the presence of xenon. These results and MD simulations suggest that the ligand migrates in the protein matrix to a major "secondary site," located beneath Tyr(B10)29 and accessible via the motion of Ile(G8)107; this site is different from the "primary site" identified by others who investigated the photolyzed state of wild-type Mb by crystallography. Our hypothesis may rationalize the O2 binding properties of Mb-YQR, and more generally to propose a mechanism of control of ligand binding and dissociation in hemeproteins based on the dynamics of side chains that may (or may not) allow access to and direct temporary sequestration of the dissociated ligand in a docking site within the protein. This interpretation suggests that very fast (picosecond) fluctuations of amino acid side chains may play a crucial role in controlling O2 delivery to tissue at a rate compatible with physiology.

Animals↗

PubMed2006

The mid-region of parathyroid hormone (1-34) serves as a functional docking domain in receptor activation.

Elucidating the bimolecular interface between parathyroid hormone (PTH) and its cognate G protein-coupled receptor (PTHR1) should yield insights into the basis of molecular recognition and the mechanism of ligand-mediated intracellular signaling for a system that is critically important in regulating calcium levels in blood. We used photoaffinity scanning (PAS) to identify key ligand-receptor interactions for residues from the unstructured mid-region domain of PTH-(1-34). Four PTH analogues, containing a single photoreactive p-benzoylphenylalanine (Bpa) residue in position 11, 15, 18, or 21, were found to photo-cross-link within receptor regions [165-176], [183-189], [190-298], and [165-176], respectively. Addition of these mid-region contacts as constraints to our previously proposed model of the PTH-PTHR1 complex and extensive molecular simulation experiments enables substantial refinement of the model. Specifically, (1) the overall receptor-bound conformation of the hormone is not extended, but bent; (2) helix [169-176] of the N-terminal extracellular domain (N-ECD) of the receptor is redirected toward the heptahelical bundle; and (3) the hormone traverses between the top of transmembrane (TM) helices 1 and 2, rather than between TM-7 and TM-1. This significantly alters the model of both the receptor-bound tertiary structure of the hormone and the topological orientation of the C-terminus of the N-ECD in the hormone-receptor bimolecular complex. We propose that the mid-region of PTH-(1-34) has a role in fixing, by extensive contacts with the receptor, the entry of the N-terminal helix of the hormone into the heptahelical bundle between TM-1 and TM-2. This anchorage would orient the amino terminus into position to activate the receptor.

Animals↗

PubMed2000

A possible role of N-cadherin in thalidomide teratogenicity.

Several experimental findings indicate that the adhesion molecule N-cadherin participates in distinct processes of embryogenesis that spatiotemporarily correlate with high sensitivity to thalidomide. Therefore, we suppose that thalidomide might interfere with N-cadherin-mediated interactions. This hypothesis is supported by protein-ligand docking studies simulating and characterizing the binding of thalidomide to N-cadherin molecules. Thalidomide was found to bind at the N-terminal domain of N-cadherin mimicking a tryptophan residue which is critical for the homodimerization of the adhesion molecule. Based on these results, we suggest that thalidomide might disturb cellular recognition and migration processes in morphogenesis by interaction with N-cadherin.

Abnormalities, Drug-Induced↗

PubMed2002

Modified AutoDock for accurate docking of protein kinase inhibitors.

Protein kinases are an important class of enzymes controlling virtually all cellular signaling pathways. Consequently, selective inhibitors of protein kinases have attracted significant interest as potential new drugs for many diseases. Computational methods, including molecular docking, have increasingly been used in the inhibitor design process [1]. We have considered several docking packages in order to strengthen our kinase inhibitor work with computational capabilities. In our experience, AutoDock offered a reasonable combination of accuracy and speed, as opposed to methods that specialize either in fast database searches or detailed and computationally intensive calculations. However, AutoDock did not perform well in cases where extensive hydrophobic contacts were involved, such as docking of SB203580 to its target protein kinase p38. Another shortcoming was a hydrogen bonding energy function, which underestimated the attraction component and, thus, did not allow for sufficiently accurate modeling of the key hydrogen bonds in the kinase-inhibitor complexes. We have modified the parameter set used to model hydrogen bonds, which increased the accuracy of AutoDock and appeared to be generally applicable to many kinase-inhibitor pairs without customization. Binding to largely hydrophobic sites, such as the active site of p38, was significantly improved by introducing a correction factor selectively affecting only carbon and hydrogen energy grids, thus, providing an effective, although approximate, treatment of solvation.

Adenosine Diphosphate↗

PubMed2005

General and targeted statistical potentials for protein-ligand interactions.

We present a novel atom-atom potential derived from a database of protein-ligand complexes. First, we clarify the similarities and differences between two statistical potentials described in the literature, PMF and Drugscore. We highlight shortcomings caused by an important factor unaccounted for in their reference states, and describe a new potential, which we name the Astex Statistical Potential (ASP). ASP's reference state considers the difference in exposure of protein atom types towards ligand binding sites. We show that this new potential predicts binding affinities with an accuracy similar to that of Goldscore and Chemscore. We investigate the influence of the choice of reference state by constructing two additional statistical potentials that differ from ASP only in this respect. The reference states in these two potentials are defined along the lines of Drugscore and PMF. In docking experiments, the potential using the new reference state proposed for ASP gives better success rates than when these literature reference states were used; a success rate similar to the established scoring functions Goldscore and Chemscore is achieved with ASP. This is the case both for a large, general validation set of protein-ligand structures and for small test sets of actives against four pharmaceutically relevant targets. Virtual screening experiments for these targets show less discrimination between the different reference states in terms of enrichment. In addition, we describe how statistical potentials can be used in the construction of targeted scoring functions. Examples are given for cdk2, using four different targeted scoring functions, biased towards increasingly large target-specific databases. Using these targeted scoring functions, docking success rates as well as enrichments are significantly better than for the general ASP scoring function. Results improve with the number of structures used in the construction of the target scoring functions, thus illustrating that these targeted ASP potentials can be continuously improved as new structural data become available.

Binding Sites↗

PubMed2006

Minimizing false positives in kinase virtual screens.

In spite of recent improvements in docking and scoring methods, high false-positive rates remain a common issue in structure-based virtual screening. In this study, the distinctive features of false positives in kinase virtual screens were investigated. A series of retrospective virtual screens on kinase targets was performed on specifically designed test sets, each combining true ligands and experimentally confirmed inactive compounds. A systematic analysis of the docking poses generated for the top-ranking compounds highlighted key aspects differentiating true hits from false positives. The most recurring feature in the poses of false positives was the absence of certain key interactions known to be required for kinase binding. A systematic analysis of 444 crystal structures of ligand-bound kinases showed that at least two hydrogen bonds between the ligand and the backbone protein atoms in the kinase hinge region are present in 90% of the complexes, with very little variability across targets. Closer inspection showed that when the two hydrogen bonds are present, one of three preferred hinge-binding motifs is involved in 96.5% of the cases. Less than 10% of the false positives satisfied these two criteria in the minimized docking poses generated by our standard protocol. Ligand conformational artifacts were also shown to contribute to the occurrence of false positives in a number of cases. Application of this knowledge in the form of docking constraints and post-processing filters provided consistent improvements in virtual screening performance on all systems. The false-positive rates were significantly reduced and the enrichment factors increased by an average of twofold. On the basis of these results, a generalized two-step protocol for virtual screening on kinase targets is suggested.

Animals↗

PubMed2006

Structure basis for the inhibitory mechanism of a novel DNase gamma-specific inhibitor, DR396.

DNase gamma, a member of the DNase I family, has been suggested to cause DNA fragmentation during apoptosis. We recently identified 4-(4,6-dichloro-[1,3,5]-triazine-2-ylamino)-2-(6-hydroxy-3-oxo-3H-xanthen-9-yl)-benzoic acid (DR396) as a novel specific inhibitor for human DNase gamma [Sunaga, S.; Kobayashi, T.; Yoshimori, A.; Shiokawa, D.; Tanuma, S. Biochem. Biophys. Res. Commun.2004, 325, 1292]. However, the binding mode (coordinate) of DR396 to DNase gamma has not yet been defined. Here, we examined the molecular basis for the inhibitory activity of DR396 to DNase gamma by structure-based computational docking studies. In the blind-docking study using a human DNase gamma homology model, a unique binding site of DR396 was predicted, which is tentatively named the 'DNA trapping site' because of the binding domain of the unhydrolyzed DNA strand, but not the active site. Targeting the DNA trapping site as a hot spot, new human DNase gamma inhibitors were obtained from our diverse chemical library in silico. These inhibitors showed high correlations between their predicted binding-free energies (DeltaGs) and observed IC50 values in the DNA trapping site but not the active site. The IC50 of a regioisomer of DR396, 5-(4,6-dichloro-[1,3,5]-triazine-2-ylamino)-2-(6-hydroxy-3-oxo-3H-xanthen-9-yl)-benzoic acid (DF365), was 73 microM (DeltaG=-9.75 kcal/mol), a 20-fold weaker inhibitory ability than that of DR396 (IC50=3.2 microM, DeltaG=-11.22 kcal/mol). Fluorescein and triazine derivatives, partial structures of DR396, had little inhibitory activity for DNase gamma. Docking analyses of the interaction between DR396 and DNase gamma revealed that DR396 binds tightly to three subsites (S1, S2, and S3) in the trapping site of DNase gamma by forming six hydrogen bonds, whereas DF365 and the partial structures are unable to form hydrogen bonds at all three subsites. These findings suggest that the specificity and potency of the inhibitory activity of DR396 for DNase gamma is due to the specific interaction of DR396 with three subsites in the DNA trapping site of DNase gamma.

Amino Acid Sequence↗

PubMed2006

Arylthiopyrrole (AThP) derivatives as non-nucleoside HIV-1 reverse transcriptase inhibitors: synthesis, structure-activity relationships, and docking studies (part 2).

Arylthio isopropyl pyridinylmethylpyrrolemethanols (AThPs) have been recently reported as a new class of human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT) inhibitors acting at the non-nucleoside binding site (NNBS) of this enzyme. Docking experiments of the potent inhibitors 4k (IC(50) = 0.24 microM, SI = 167) and 5e (IC(50) = 0.11 microM, SI > 1667) of wild-type RT prompted the synthesis and biological evaluation of novel AThP derivatives featuring a number of polar groups in position 3 of the pyrrole ring and larger and more hydrophobic alicyclic substituents in place of the isopropyl group at position 4. Among the compounds synthesized and tested in cell-based assays against HIV-1 infected cells, 19b was the most active, with EC(50) = 0.007 microM, CC(50) = 114.5 microm, and SI = 16357. This compound and its precursor 18b retained interesting activities against clinically relevant drug-resistant RT forms carrying K103N, Y181I, and L100I mutations. Docking calculations of 10, 14, 18b, and 19b were also performed to investigate their binding mode into the RT NNBS and to rationalize both structure-activity relationship and resistance data.

Binding Sites↗

PubMed2005

The solution structure of BmTx3B, a member of the scorpion toxin subfamily alpha-KTx 16.

This article reports the solution structure of BmTx3B (alpha-KTx16.2), a potassium channel blocker belonging to the subfamily alpha-KTx16, purified from the venom of the Chinese scorpion Buthus martensi Karsch. In solution, BmTx3B assumes a typical CSalphabeta motif, with an alpha-helix connected to a triple-stranded beta-sheet by 3 disulfide bridges, which belongs to the first structural group of short-chain scorpion toxins. On the other hand, BmTx3B is quite different from other toxins (such as ChTx and AgTx2) of this group in terms of the electrostatic and hydrophobic surface distribution. The functional surface (beta-face) of the molecule is characterized by less basic residues (only 2: Lys28 and Arg35) and extra aromatic residues (Phe1, Phe9, Trp15, and Tyr37). The peptide shows a great preference for the Kca1.1 channel over the Kv channel (about a 10(3)-fold difference). The model of BmTx3B/Kca1.1 channel complex generated by docking and dynamic simulation reveals that the stable binding between the BmTx3B and Kca1.1 channel is favored by a number of aromatic pi-pi stacking interactions. The influences of these structural features on the kinetic behavior of the toxin binding to Kca1.1 channel are also discussed.

Amino Acid Sequence↗

PubMed2003

Identification of small-molecule inhibitors of human angiogenin and characterization of their binding interactions guided by computational docking.

Angiogenin (ANG) is a potent inducer of angiogenesis and an RNase A homologue whose ribonucleolytic activity is essential for its biological action. Recently, we reported the identification of small non-nucleotide inhibitors of the enzymatic activity of ANG by high-throughput screening (HTS) [Kao, R. Y. T., et al. (2002) Proc. Natl. Acad. Sci. U.S.A. 99, 10066-10071]. Two of the inhibitors that were obtained, National Cancer Institute compound NSC-65828 [8-amino-5-(4'-hydroxybiphenyl-4-ylazo)naphthalene-2-sulfonate] and ChemBridge compound C-181431 [4,4'-dicarboxy-3,3'-bis(naphthylamido)diphenylmethanone], were judged to be suitable for further development, and one of these (NSC-65828) was shown to possess antitumor activity in mice. Here we have used computational docking as a guide for the identification of available NSC-65828 and C-181431 analogues that bind more tightly to ANG, and for the characterization of inhibitor binding modes. Numerous analogues were found to have greater avidity than the HTS compounds or any small nucleotide inhibitors; four were considered to be of interest as potential leads (K(i) = 5-25 microM). Two of these analogues bind more tightly to ANG than to RNase A, and are the first small molecules shown to exhibit this selectivity. The predicted binding orientations of the HTS compounds and the new lead inhibitors were evaluated by determining the effects of ANG active site mutations on inhibitory potency. The results with ANG variants R5A, H8A, N68A, and des(121-123) are highly consistent with the docking models. Affinity changes observed with Q12A and Q117G reveal aspects of active site function that are not apparent from the free ANG crystal structure or from the modeled complexes. These findings should prove to be useful in the design of more effective and specific ANG antagonists.

Amino Acid Substitution↗

PubMed2006

C-(4,5,6-trimethoxyindan-1-yl)methanamine: a mescaline analogue designed using a homology model of the 5-HT2A receptor.

A conformationally restricted analogue of mescaline, C-(4,5,6-trimethoxyindan-1-yl)-methanamine, was designed using a 5-HT(2A) receptor homology model. The compound possessed 3-fold higher affinity and potency than and efficacy equal to that of mescaline at the 5-HT(2A) receptor. The new analogue substituted fully for LSD in drug discrimination studies and was 5-fold more potent than mescaline. Resolution of this analogue into its enantiomers corroborated the docking experiments, showing the R-(+) isomer to have higher affinity and potency and to have efficacy similar to that of mescaline at the 5-HT(2A) receptor.

Animals↗

PubMed1996

Three-dimensional models for agonist and antagonist complexes with beta 2 adrenergic receptor.

Computer-modeling techniques have been used to generate docked complexes for a series of beta adrenergic agonists and antagonists with a three-dimensional model of the beta 2 adrenergic receptor. For all ligands tested, it proved possible to dock low-energy conformers in the receptor model, with sensible electrostatic, steric, and hydrogen-bonding interactions, many of which are supported by experimental studies of the beta 2 receptor. Our results illustrate the power of molecular modeling techniques, when coupled with appropriate experimental methods and data, to investigate structure-function properties of integral membrane receptor proteins that cannot yet be studied by direct structural methods.

Adrenergic beta-Agonists↗