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 793 records · Page 44Linked to original sources

PubMedDate not supplied

Molecular simulation of ligand-binding with DNA: implications for 125I-labeled pharmaceutical design.

PURPOSE: By using computer-assisted molecular modeling software, to assess the effects of structural modification on the interaction of 125I-labeled iodoHoechst ligands and DNA and to design new analogs with specified distances between the Auger-electron-emitting 125I atom and the DNA central axis. MATERIALS AND METHODS: The Lamarckian genetic algorithm (AutoDock 3.0) was used to model the interaction between DNA and m-iodo-p-methoxyHoechst (IMH), a ligand whose binding to the minor groove of DNA has been demonstrated (crystal structure) and which is available in the Protein Data Bank. m-Iodo-p-ethoxyHoechst (IEH), a radioligand we had previously synthesized and characterized, was then docked onto DNA, the IEH-DNA complex minimized, and the binding free energy and inhibition constant (Ki) estimated and compared with those for IMH-DNA. Using the protocol, several novel iodoHoechst analogs (IH-A and IH-B) were designed. Finally, Insight II was used to measure the distances between the iodine atom (e.g. 125I) of these Hoechst analogs and of 5-iodo-2'-deoxyuridine (IdUrd) and the central axis of the targeted DNA, and these values were correlated with the expected/measured DSB yield following 125I decay. RESULTS: The docking of IMH and DNA leads to a ligand-DNA complex approximately 1 A RMSD (root mean square deviation) from the crystal-structure position, and the IEH-DNA complex also has a small RMSD (1.27 A). The distances between the 125I atom and the DNA central axis are estimated as 8.61 A for IMH, 9.20 A for IEH and 5.7A for IdUrd. For the newly designed analogs, the distances from the 125I atom to the central DNA axis are 10.92 A for IH-A and 16.39 A for IH-B. CONCLUSION: These software programs can predict the reactivity of newly designed radiolabelled molecules with their targeted DNA molecules by molecular modeling prior to their chemical synthesis.

Binding Sites↗

PubMed2004

Molecular recognition in purinergic receptors. 2. Diastereoselectivity of the h-P2Y1-receptor.

In the companion paper, part 1, we described the construction of an improved molecular model for the h-P2Y1 receptor (h-P2Y1-R) and proposed a rational for the stereoelectronic selectivity of the receptor. Here, we extend our studies on the molecular recognition of the h-P2Y1-R to the exploration of the diastereoselectivity of this receptor. For this purpose, we implemented an integrative approach combining synthesis, spectral analysis, biochemical assays, and computational analysis. Specifically, we selected and synthesized novel ATP analogues bearing a chiral center on the phosphate chain. We analyzed the conformation of the chiral ATP analogues in solution by 1H/13C NMR and assigned the absolute configuration of the diastereoisomers. The coordination mode of these analogues with a Mg2+ ion was evaluated by 31P NMR. These chiral analogues were biochemically evaluated and found to be potent h-P2Y1-R ligands. An EC50 difference of ca. 20-fold was observed between the diastereoisomers. Their spectral absolute configuration assignment was confirmed by comparison of the biochemical results to those of ATP-alpha-S diastereoisomers whose chirality is known. Finally, a computational analysis was performed for the elucidation of molecular recognition employing molecular mechanics (docking) studies on the receptor:ligands complexes. On the basis of the current results, we hypothesize that h-P2Y1-R's chiral discrimination originates from the requirement that the nucleotide analogue interacts with a Mg2+ ion within the receptor binding site. This Mg2+ ion is possibly coordinated with both Asp204 and the ATP's alpha, beta, gamma-phosphates in a Lambda configuration.

Adenosine Triphosphate↗

PubMed2004

The selection and design of GPCR ligands: from concept to the clinic.

Virtual screening methods using structure-based, pharmacophore-based and descriptor based protocols may be used to identify ligands for the G-protein coupled receptor target family. A complementary approach is the synthesis and screening of compound libraries designed using privileged motifs and/or based on validated hit molecules. A virtual screening approach based on molecular docking performed with GOLD using a templated homology model and a consensus scoring procedure can identify vasopressin 1a receptor antagonists. In a separate project a library design and synthesis approach based around validated hit GPCR ligands led to the identification of potent oxytocin antagonists. Subsequent optimisation of the initial library compounds has provided compounds that are now being evaluated in the clinic for the treatment of preterm labour.

Animals↗

PubMed1993

Orientational sampling and rigid-body minimization in molecular docking.

The biological activities of proteins depend on specific molecular recognition and binding. Computational methods for predicting binding modes can facilitate the discovery and design of ligands and yield information on the factors governing complementarity. The DOCK suite of programs has been applied to several systems; here, the degree of orientational sampling required to reproduce and identify known binding modes, with and without rigid-body energy minimization, is investigated for four complexes. There is a tradeoff between sampling and minimization. The known binding modes can be identified with intensive sampling alone (10,000 to 20,000 orientations generated per system) or with moderate sampling combined with minimization. Optimization improves energies significantly, particularly when steric clashes are present, and brings many orientations closer to the experimentally observed position. Whether or not minimization is performed, however, sampling must be sufficient to find at least one structure in the vicinity of the presumed true binding mode. Hybrid approaches combining docking and minimization are promising and will become more viable with the use of faster algorithms and the judicious selection of fewer orientations for minimization.

Carboxypeptidases↗

PubMed1997

Spectroscopic and binding studies on the interaction of inorganic anions with lactoperoxidase.

The interaction of several inorganic species (SCN-, I-, Br-, Cl-, F-, NO2-, N3-, CN-) with bovine lactoperoxidase was investigated through kinetic and binding studies by using UV-Vis spectroscopy. The above ligands form 1:1 complexes with the protein and can be assigned to three different groups, on the basis of the dissociation constant values (KD) of the adducts: (1) SCN-, I-, Br-, and Cl- (KD increases along the series); (2) F- (which shows a singular behavior); (3) NO2-, N3-, and CN- (that bind at the iron site). KD values for the LPO/SCN- adduct appeared to be modified in the presence of other inorganic species; a strong competition between this substrate and all other anions (with the exception of F-) was evidentiated. Binding investigations on the natural substrates SCN- and I-, at varying pH and temperature, showed that their interaction with lactoperoxidase involves the protonation of a common site in proximity of the iron (possibly distal histidine). Michaelis-Menten constants for SCN-, I-, and Br- followed roughly the same trend as KD; KM for hydrogen peroxide is strongly dependent on the cosubstrate. Computer-assisted docking simulations showed that all ligands can penetrate inside the heme pocket.

Animals↗

PubMed1998

Modelling a 3D structure for EgDf1 from Echinococcus granulosus: putative epitopes, phosphorylation motifs and ligand.

EgDf1 is a developmentally regulated protein from the parasite Echinococcus granulosus related to a family of hydrophobic ligand binding proteins. This protein could play a crucial role during the parasite life cycle development since this organism is unable to synthetize most of their own lipids de novo. Furthermore, it has been shown that two related protein from other parasitic platyhelminths (Fh15 from Fasciola hepatica and Sm14 from Schistosoma mansoni) are able to confer protective inmunity against experimental infection in animal models. A three-dimensional structure would help establishing structure/function relationships on a knowledge based manner. 3D structures for EgDf1 protein were modelled by using myelin P2 (mP2) and intestine fatty acid binding protein (I-FABP) as templates. Molecular dynamics techniques were used to validate the models. Template mP2 yielded the best 3D structure for EgDf1. Palmitic and oleic acids were docked inside EgDf1. The present theoretical results suggest definite location in the secondary structure of the epitopic regions, consensus phosphorylation motifs and oleic acid as a good ligand candidate to EgDf1. This protein might well be involved in the process of supplying hydrophobic metabolites for membrane biosynthesis and for signaling pathways.

Amino Acid Sequence↗

PubMed2004

NMR solution structure of a potent cyclic nonapeptide inhibitor of ICAM-1-mediated leukocyte adhesion produced by homologous amino acid substitution.

We have previously described a disulfide-linked cyclic nonapeptide (inhibitory peptide-01, IP01), with the sequence CLLRMRSIC, which binds to intercellular adhesion molecule-1 (ICAM-1), and blocks binding to its counter-structure, the integrin alphaLbeta2 (leukocyte functional antigen-1, LFA-1) (Sillerud et al., J. Peptide Res. 62, 2003: 97). We now report the optimization of this peptide by means of single homologous amino acid substitutions to yield a new peptide (IP02-K6; CLLRMKSAC) which shows an approximately sixfold improvement in inhibitory activity of multivalent leukocyte binding (inhibition constant for 50% inhibition, IC50 = 90 microm) compared with IP01 (IC50 = 580 microm). This improvement in activity gives IP02-K6 potent in vivo activity in a murine model of ischemia reperfusion injury (Merchant et al., Am. J. Physiol. Heart Circ. 284, 2003: H1260). In order to determine the structural features relevant to ICAM-1-binding, we have determined the structure of IP02-K6 using proton nuclear magnetic resonance (NMR) spectroscopy and restrained molecular modeling. In our previously reported study of solution models of IP01, we observed three interconverting conformations during low-temperature molecular dynamics simulation. In the present study, we find a single conformation of IP02-K6 similar to one of the previously found conformations of IP01 (family C). In particular, an R4-S7 beta-turn is present in similar proportions in both conformation C of IP01 and in IP02-K6; this motif is important in binding to ICAM-1 because this turn enables the IP02-K6 backbone to drape over proline-36 on ICAM-1. The NMR-derived solution model of IP02-K6 was found to dock at the alphaLbeta2-binding site on ICAM-1 with no changes in peptide backbone conformation. This docking model displaced five of the 15 alphaLbeta2 residues at the ICAM-1-binding site and provided a rationale for understanding the quantitative relationship between IP02-K6 structure and biologic activity.

Amino Acid Motifs↗

PubMed2001

DOCK 4.0: search strategies for automated molecular docking of flexible molecule databases.

In this paper we describe the search strategies developed for docking flexible molecules to macomolecular sites that are incorporated into the widely distributed DOCK software, version 4.0. The search strategies include incremental construction and random conformation search and utilize the existing Coulombic and Lennard-Jones grid-based scoring function. The incremental construction strategy is tested with a panel of 15 crystallographic testcases, created from 12 unique complexes whose ligands vary in size and flexibility. For all testcases, at least one docked position is generated within 2 A of the crystallographic position. For 7 of 15 testcases, the top scoring position is also within 2 A of the crystallographic position. The algorithm is fast enough to successfully dock a few testcases within seconds and most within 100 s. The incremental construction and the random search strategy are evaluated as database docking techniques with a database of 51 molecules docked to two of the crystallographic testcases. Incremental construction outperforms random search and is fast enough to reliably rank the database of compounds within 15 s per molecule on an SGI R10000 cpu.

Algorithms↗

PubMed2006

Inherent limitations in protein-protein docking procedures.

MOTIVATION: The limited success rate of protein-protein docking procedures is generally attributed to structure differences between the bound and unbound states of the molecules. Herein we analyze a large dataset of protein-protein docking results and identify additional parameters that affect the performance of docking procedures. RESULTS: We find that the distinction between nearly correct models (NCMs) and decoys depends on the size of the interface to be predicted thus setting a limit to the prediction ability of docking procedures, particularly those in which the geometric complementarity descriptor is dominant. The geometric complementarity score in grid-based docking carries a large statistical error which further reduces the distinction between NCMs and decoys. We propose a method for correcting the statistical error and show that the distinction is improved when the docking models are ranked by statistically equivalent scores. AVAILABILITY: MolFit can be downloaded from our website http://www.weizmann.ac.il/Chemical_Research_Support/molfit. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Algorithms↗

PubMed2006

Molecular docking-based study of vasopressin analogues modified at positions 2 and 3 with N-methylphenylalanine: influence on receptor-bound conformations and interactions with vasopressin and oxytocin receptors.

In this study, four cyclic vasopressin (CYFQNCPRG-NH(2), AVP) analogues substituted at positions 2 and 3 with four combinations of enantiomers of N-methylphenylalanine have been investigated. Three-dimensional structures of analogues have been formerly determined using NMR spectroscopy in dimethyl sulfoxide. Three-dimensional models of the vasopressin and oxytocin receptors were constructed by combining the multiple sequence alignment and the RD crystal structure as a template. The analogues have been docked into the receptor using the AutoDock program. The relaxation of the receptor-ligand complexes using energy minimization, followed by the constrained simulated annealing protocols (CSA), has been performed. The receptor-bound conformations of the investigated analogues have been proposed. We concluded that the N-methylated residues at positions 2 and 3 act as a structural restraint, determining the conformation of analogues, their location inside the receptor cavity, and mutual arrangement of the aromatic side chains. The conserved polar residues constitute the handles keeping the biologically active analogues inside the binding cavity. The Arg(8)-D(2.65) salt bridge might be responsible for analogue-selective binding in OTR and V1aR versus V2R, where the positively charged K(2.65) 100 is present at the equivalent position.

Amino Acid Sequence↗

PubMed1996

New approaches in molecular structure prediction.

In the past years, much effort has been put on the development of new methodologies and algorithms for the prediction of protein secondary and tertiary structures from (sequence) data; this is reviewed in detail. New approaches for these predictions such as neural network methods, genetic algorithms, machine learning, and graph theoretical methods are discussed. Secondary structure prediction algorithms were improved mostly by considering families of related proteins; however, for the reliable tertiary structure modeling of proteins, knowledge-based techniques are still preferred. Methods and examples with more or less successful results are described. Also, programs and parameterizations for energy minimisations, molecular dynamics, and electrostatic interactions have been improved, especially with respect to their former limits of applicability. Other topics discussed in this review include the use of traditional and on-line databases, the docking problem and surface properties of biomolecules, packing of protein cores, de novo design and protein engineering, prediction of membrane protein structures, the verification and reliability of model structures, and progress made with currently available software and computer hardware. In summary, the prediction of the structure, function, and other properties of a protein is still possible only within limits, but these limits continue to be moved.

Chemical Phenomena↗

PubMed1995

Cimetidine inhibition of human gastric and liver alcohol dehydrogenase isoenzymes: identification of inhibitor complexes by kinetics and molecular modeling.

Cimetidine, an H2-receptor antagonist, is one of the most commonly prescribed drugs in the world. It has been reported to increase blood alcohol concentrations in drinking individuals. To determine if this increase could be due to inhibition of alcohol dehydrogenase activity, the effect of the drug on ethanol oxidation by gastric sigma sigma alcohol dehydrogenase and liver beta 2 beta 2, pi pi, and chi chi alcohol dehydrogenase isoenzymes was observed. Cimetidine inhibited all isoenzymes studied except chi chi; the chi chi isoenzyme showed no inhibition up to 5 mM cimetidine. Inhibition of the alcohol dehydrogenase isoenzymes by the H2-receptor antagonists nizatidine, ranitidine, and famotidine was negligible. Docking simulations with the beta 2.NAD+.4-iodopyrazole X-ray structure indicated that cimetidine fit well into the substrate binding site. The substitution on the thiazole ring of nizatidine, however, prevented docking into the binding site. Cimetidine inhibition of ethanol oxidation by sigma sigma and beta 2 beta 2 was competitive with varied ethanol, exhibiting Ki values of 2.8 +/- 0.4 mM and 0.77 +/- 0.07 mM, respectively. Cimetidine inhibition of ethanol oxidation by pi pi was noncompetitive with varied ethanol (Ki = 0.50 +/- 0.03 mM). Inhibition of ethanol oxidation by sigma sigma and beta 2 beta 2 with varied NAD+ was competitive. These results, together with the cimetidine inhibition kinetics of acetaldehyde reduction by sigma sigma and beta 2 beta 2, with either varied NADH or varied acetaldehyde, are consistent with cimetidine binding to two enzyme species. These species are free enzyme and the productive enzyme.NAD+ complex.

Alcohol Dehydrogenase↗

PubMed1996

Conformational analysis of retinoids and restriction of their dynamics by retinoid-binding proteins.

An exhaustive sampling of the configurational space of all-trans retinol using a 0.1 microsecond molecular-dynamics simulation is presented. The essential dynamics technique is used to describe the conformational changes in retinol using only three degrees of freedom. The different conformational states of retinol are analysed, and differences in free energy are calculated. The essential dynamics description allows a detailed comparison of free retinol and retinoids bound to retinoid-binding proteins and opens new possibilities in the small-molecule docking field. The dynamics of retinoids when complexed with their binding proteins are restricted, and they are forced into strained conformations. A "spring' model for retinoid binding is proposed. This model is extended to a hypothesis for retinoid binding to visual pigments and bacteriorhodopsin.

Computer Simulation↗

PubMed2005

Terphenyl-Based Bak BH3 alpha-helical proteomimetics as low-molecular-weight antagonists of Bcl-xL.

We describe a general method for the mimicry of one face of an alpha-helix based on a terphenyl scaffold that spatially projects functionality in a manner similar to that of two turns of an alpha-helix. The synthetic scaffold reduces the flexibility and molecular weight of the mimicked protein secondary structure. We have applied this design to the development of antagonists of the alpha-helix binding protein Bcl-x(L). Using a sequential synthetic strategy, we have prepared a library of terphenyl derivatives to mimic the helical region of the Bak BH3 domain that binds Bcl-x(L). Fluorescence polarization assays were carried out to evaluate the ability of terphenyl derivatives to displace the Bcl-x(L)-bound Bak peptide. Terphenyl 14 exhibited good in vitro affinity with a K(i) value of 0.114 muM. These terphenyl derivatives were more selective at disrupting the Bcl-x(L)/Bak over the HDM2/p53 interaction, which involves binding of the N-terminal alpha-helix of p53 to HDM2. Structural studies using NMR spectroscopy and computer-aided docking simulations suggested that the helix binding area on the surface of Bcl-x(L) is the target for the synthetic ligands. Treatment of human embryonic kidney 293 (HEK293) cells with terphenyl derivatives resulted in the disruption of the binding of Bcl-x(L) to Bax in intact cells.

Biomimetic Materials↗

PubMed1996

Low-resolution docking: prediction of complexes for underdetermined structures.

One of the most fundamental questions concerning ligand-receptor interaction is whether such a process of intermolecular association is generally determined by local structural elements of the participating molecules, or whether there are also large-scale motifs in molecule structures that facilitate complex formation. From the point of view of practical docking computations, the elaborate character of local structural details in ligand-receptor interaction creates a large number of false-positive matches, which interfere with determination of the best fit. Another significant obstacle in protein docking is the problem of structural data inaccuracy (poor structure resolution, conformational changes upon complex formation, etc.). Our study [Vakser (1995) Protein Eng., 8, 371-377], based on ultralow (approximately 7 A resolution) representation of molecular structures, allowes to average all high-resolution structural details, and still predict most of the structural features of the ligand-receptor complex. The approach dramatically improves the signal-to-noise ratio in determination of the best fit, and moves the structure inaccuracy tolerance to the range of the macrostructure. In the present paper, we describe a further validation of the main principles of this approach and a detailed analysis of the low-resolution docking results. This includes clustering of ligand positions around the receptor molecule and cross-validation of ligands and receptors from different complexes. We also discuss the important implications of the approach to the multiple-minima problem and a possible role of different structural elements in the recognition mechanism.

Animals↗

PubMed2004

Solution structure of BmKK2, a new potassium channel blocker from the venom of chinese scorpion Buthus martensi Karsch.

A natural K+ channel blocker, BmKK2 (a member of scorpion toxin subfamily alpha-KTx 14), which is composed of 31 amino acid residues and purified from the venom of the Chinese scorpion Buthus martensi Karsch, was characterized using whole-cell patch-clamp recording in rat hippocampal neurons. The three dimensional structure of BmKK2 was determined with two-dimensional NMR spectroscopy and molecular modelling techniques. In solution this toxin adopted a common alpha/beta-motif, but showed distinct local conformation in the loop between alpha-helix and beta-sheet in comparison with typical short-chain scorpion toxins (e.g., CTX and NTX). Also, the alpha helix is shorter and the beta-sheet element is smaller (each strand consisted only two residues). The unusual structural feature of BmKK2 was attributed to the shorter loop between the alpha-helix and beta-sheet and the presence of two consecutive Pro residues at position 21 and 22 in the loop. Moreover, two models of BmKK2/hKv1.3 channel and BmKK2/rSK2 channel complexes were simulated with docking calculations. The results demonstrated the existence of a alpha-mode binding between the toxin and the channels. The model of BmKK2/rSK2 channel complex exhibited favorable contacts both in electrostatic and hydrophobic, including a network of five hydrogen bonds and bigger interface containing seven pairs of inter-residue interactions. In contrast, the model of BmKK2/hKv1.3 channel complex, containing only three pairs of inter-residue interactions, exhibited poor contacts and smaller interface. The results well explained its lower activity towards Kv channel, and predicted that it may prefer a type of SK channel with a narrower entryway as its specific receptor.

Amino Acid Sequence↗

PubMed2005

Molecular docking study on the "back door" hypothesis for product clearance in acetylcholinesterase.

Acetylcholinesterase (AChE) is one of the fastest enzymes known, even though the active site is buried inside the protein at the end of a 20-A deep narrow gorge. Among the great variety of crystal structures of this enzyme, both in the absence and presence of various ligands and proteins, the structure of a complex of AChE with the pseudo-irreversible inhibitor Mf268 is of particular interest, as it assists in the proposal of a back door for product clearance from the active site. Binding of Mf268 to AChE results in the carbamoylation of Ser200 and liberation of an eseroline-fragment as the leaving group. The crystal structure of the AChE-Mf268 complex, however, proves that eseroline has escaped from the enzyme, despite the fact that the Ser-bound inhibitor fragment blocks the gorge entrance. The existence of alternative routes other than through the gorge for product clearance has been postulated but is still controversially discussed in the literature, as an experimental proof for such a back door is still missing. We have used Monte Carlo-based molecular docking methods in order to examine possible alternative pathways that could allow eseroline to be released from the protein after being cleaved from the substrate by Ser200. Based on our results, a short channel at the bottom of the gorge seems to be the most probable back-door site, which begins at amino acid Trp84 and ends at the enzyme surface in a cavity close to amino acid Glu445. [Figure: see text].

Acetylcholinesterase↗

PubMed2005

Pseudoazurin-nitrite reductase interactions.

The nitrite reductase-binding site on pseudoazurin has been determined by using NMR chemical-shift perturbations. It comprises residues in the hydrophobic patch surrounding the exposed copper ligand His81 as well as several positively charged residues. The binding site is similar for both redox states of pseudoazurin, despite differences in the binding mode. The results suggest that pseudoazurin binds in a well-defined orientation. Docking simulations provide a putative structure of the complex with a binding site on nitrite reductase that has several hydrophobic and polar residues as well as a ridge of negatively charged side chains and a copper-to-copper distance of 14 A.

Azurin↗