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

PubMed2006

Structural and functional characterization of human microsomal prostaglandin E synthase-1 by computational modeling and site-directed mutagenesis.

Microsomal prostaglandin (PG) E synthase-1 (mPGES-1) has recently been recognized as a novel, promising drug target for inflammation-related diseases. Functional and pathological studies on this enzyme further stimulate to understand its structure and the structure-function relationships. Using an approach of the combined structure prediction, molecular docking, site-directed mutagenesis, and enzymatic activity assay, we have developed the first three-dimensional (3D) model of the substrate-binding domain (SBD) of mPGES-1 and its binding with substrates prostaglandin H2 (PGH2) and glutathione (GSH). In light of the 3D model, key amino acid residues have been identified for the substrate binding and the obtained experimental activity data have confirmed the computationally determined substrate-enzyme binding mode. Both the computational and experimental results show that Y130 plays a vital role in the binding with PGH2 and, probably, in the catalytic reaction process. R110 and T114 interact intensively with the carboxyl tail of PGH2, whereas Q36 and Q134 only enhance the PGH2-binding affinity. The modeled binding structure indicates that substrate PGH2 interacts with GSH through hydrogen binding between the peroxy group of PGH2 and the -SH group of GSH. The -SH group of GSH is expected to attack the peroxy group of PGH2, initializing the catalytic reaction transforming PGH2 to prostaglandin E2 (PGE2). The overall agreement between the calculated and experimental results demonstrates that the predicted 3D model could be valuable in future rational design of potent inhibitors of mPGES-1 as the next-generation inflammation-related therapeutic.

Algorithms↗

PubMedDate not supplied

Considerations in compound database preparation--"hidden" impact on virtual screening results.

Structure-based virtual screening (SBVS) utilizing docking algorithms has become an essential tool in the drug discovery process, and significant progress has been made in successfully applying the technique to a wide range of receptor targets. In silico validation of virtual screening protocols before application to a receptor target using a corporate or commercially available compound collection is key to establishing a successful process. Ultimately, retrieval of a set of active compounds from a database of inactives is required, and the metric of enrichment (E) is habitually used to discern the quality of separation of the two. Numerous reports have addressed the performance of docking algorithms with regard to the quality of binding mode prediction and the issue of postprocessing "hit lists" of docked ligands. However, the impact of ligand database preprocessing has yet to be examined in the context of virtual screening and prioritization of compounds for biological evaluation. We provide an insight into the implications of cheminformatic preprocessing of a validation database of compounds where multiple protonated, tautomeric, stereochemical, and conformational states have been enumerated. Several commonly used methods for the generation of ligand conformations and conformational ensembles are examined, paired with an exhaustive rigid-body algorithm for the docking of different "multimeric" compound representations to the ligand binding site of the human estrogen receptor alpha. Chemgauss, a shapegaussian scoring function with intrinsic chemical knowledge, was combined with PLP as a consensus-scoring scheme to rank output from the docking protocol and enrichment rates calculated for each screen. The overheads of CPU consumption and the effect on relative database size (disk requirement) for each of the protocols employed are considered. Assessment of these parameters indicates that SBVS enrichments are highly dependent on the initial cheminformatic treatment(s) used in database construction. The interplay of SMILES representations, stereochemical information, protonation state enumeration, and ligand conformation ensembles are critical in achieving optimum enrichment rates in such screening.

Computer Simulation↗

PubMed2006

In-silico pharmacodynamics: correlation of adverse effects of H2-antihistamines with histamine N-methyl transferase binding potential.

Adverse effects are exhibited by most drugs in current clinical practice, the causes for which are often not known. In this post genomic era, bioinformatics has the potential to address several issues in understanding the mechanism of drug action and in designing improved drugs. This study describes the analysis of the possible pharmacodynamic behaviour of antihistamines blocking the histamine H(2) receptor (H(2)-antihistamines), by adopting the basic tenets of a systems biology approach. The different components that could form an appropriate sub-system are identified, thus providing a system landscape. Docking and analysis of the chosen antihistamines into each of these components resulted in identifying histamine N-methyl transferase (HNMT) as a potential unintended target for H(2)-antihistamines. Correlation with experimental data available from the literature indicates the inhibition of HNMT to be a possible cause for the adverse effects exhibited by these drugs. Implications for design of safer H(2)-antihistamines are discussed. The method reported here has the potential for application as a general strategy in understanding drug effects.

Binding Sites↗

PubMed2004

Inhibitory mode of indole-2-carboxamide derivatives against HLGPa: molecular docking and 3D-QSAR analyses.

The interaction of a series of indole-2-carboxamide compounds with human liver glycogen phosphorylase a (HLGPa) have been studied employing molecular docking and 3D-QSAR approaches. The Lamarckian Genetic Algorithm (LGA) of AutoDock 3.0 was employed to locate the binding orientations and conformations of the inhibitors interacting with HLGPa. The binding models were demonstrated in the aspects of inhibitor's conformation, subsite interaction, and hydrogen bonding. The very similar binding conformations of these inhibitors show that they interact with HLGPa in a very similar way. Good correlations between the calculated interaction free energies and experimental inhibitory activities suggest that the binding conformations of these inhibitors are reasonable. The structural and energetic differences in inhibitory potencies of indole-2-carboxamide compounds were reasonably explored. Using the binding conformations of indole-2-carboxamides, consistent and highly predictive 3D-QSAR models were developed by CoMFA and CoMSIA analyses. The q2 values are 0.697 and 0.622 for CoMFA and CoMSIA models, respectively. The predictive ability of these models was validated by four compounds that were not included in the training set. Mapping these models back to the topology of the active site of HLGPa leads to a better understanding of the vital indole-2-carboxamide-HLGPa interactions. Structure-based investigations and the final 3D-QSAR results provide clear guidelines and accurate activity predictions for novel inhibitor design.

Binding Sites↗

PubMed2005

Ligand docking in the gastric H+/K+-ATPase: homology modeling of reversible inhibitor binding sites.

Using the recent high-resolution X-ray structures determined for the Ca2+-ATPase, we have generated two homology models of the gastric H+/K+-ATPase reflecting the E1 and E2 conformations adopted by P-type ATPases in their catalytic cycle. In regimes where the in situ solid-state NMR-determined structure for 1,2,3-trimethyl-8-(pentafluorophenylmethoxy)imidazo[1,2-a]pyridinium iodide (TMPFPIP), a reversible inhibitor of the gastric H+/K+-ATPase, was retained in its predefined conformation and was allowed full torsional flexibility in docking, the ligands localized to discrete binding volumes in the E1 model and to a single central binding space, together with secondary peripheral locations, in the E2 conformation. The results of these binding studies are in good agreement with current site-directed mutagenesis data and support the suggestion that the binding site is proximal to the loop between TM5 and TM6 and TM8, the transmembrane (TM) region considered important for cation translocation. Furthermore, the results of the simulation with the flexible ligand complement the solid-state NMR structural constraints of this inhibitor when bound in situ to the protein.

Binding Sites↗

PubMed2004

Visualizing complexes of phospholipids with Streptomyces phospholipase D by automated docking.

The automated docking program AutoDock was used to dock nine phosphatidic acids (PAs), six phosphatidylcholines, five phosphatidylethanolamines, four phosphatidylglycerols, one phosphatidylinositol and two phosphatidylserines, which have two identical saturated fatty acid residues with an even numbers of carbon atoms, onto the active site of Streptomyces sp. PMF phospholipase D (PLD). Two PAs with one double bond on the fatty acid chain linked to the C2 of the glycerol residue were also docked. In general, binding energies become progressively more negative as fatty acid residues become longer. When these residues are of sufficient length, one is coiled against a hydrophobic cliff in a well that also holds the glycerol and phosphate residues and the head group, while the other generally is bound by a hydrophobic surface outside the well. Phosphatidylcholines have the only head group that is firmly bound by the active site, giving a possible structural explanation for the low selectivity of Streptomyces PLD for other phospholipid substrates.

Binding Sites↗

PubMed2004

Antagonist binding in the rat muscarinic receptor A study by docking and X-ray crystallography.

A series of agonists to the rat muscarinic receptor have been docked computationally to the active site of a homology model of rat M1 muscarinic receptor. The agonists were modelled on the X-ray crystal structure of atropine, which is reported here and the docking studies are shown to reproduce correctly the order of experimental binding affinities for the agonists as well as indicate where there appear to be inconsistencies in the experimental data. The crystal and molecular structure of atropine (tropine tropate; alpha-[hydroxymethyl]benzeneacetic acid 8-methyl[3.2.1]oct-3-yl ester C17H23NO3) has been determined by X-ray crystallography using an automated Patterson search method, and refined by full-matrix least-squares to a final R of 0.0452 for 2701 independent observed reflections and 192 parameters using Mo Kalpha radiation, lambda=0.71073A at 150K. The compound crystallises in space group Fdd2 with Z=16 molecules per unit cell.

Animals↗

PubMed2000

Synthetic peptide mimicking of binding sites on olfactory receptor protein for use in 'electronic nose'.

A putative tertiary structure model of the dog's olfactory receptor (olfd canfa) is established in this study. By using a target odorous compound (trimethylamine), it is possible to locate the most plausible binding sites between the receptor model structure and the target odorous molecules through computer docking simulations. The two short oligo-peptide sequences (orp61 and orp188) for trimethylamine sensing were identified, synthesized, purified and coated onto the surface of the separate piezoelectric gold electrodes. These two peptides show a high binding capability for trimethylamine. To further enhance the sensitivity of the polypeptides towards the target compound, the polarity and the degree of docking were changed by a site-specific modification technique. The orp61 sequence was modified by substituting two amino acids in the binding pocket resulting in 33% increase in sensitivity towards trimethylamine and reduced noises from other non-target chemicals. The techniques used in the present study offer a unique approach for synthesizing peptides in mimicking binding domain of olfactory receptors. The approach can be easily applied to further development of recognized molecules for gas sensing, especially for use in 'electronic noses'.

Amino Acid Sequence↗

PubMed2004

Novel technologies for virtual screening.

There are several methods for virtual screening of databases of small organic compounds to find tight binders to a given protein target. Recent reviews in Drug Discovery Today have concentrated on screening by docking and by pharmacophore searching. Here, we complement these reviews by focusing on virtual screening methods that are based on analyzing ligand similarity on a structural level. Specifically, we concentrate on methods that exploit structural properties of the complete ligand molecules, as opposed to using just partial structural templates, such as pharmacophores. The in silico procedure of virtual screening (VS) and its relationship to the experimental procedure, HTS, is discussed, new developments in the field are summarized and perspectives on future research are offered.

Combinatorial Chemistry Techniques↗

PubMed1998

Prediction of an anti-IgE binding site on IgE.

A combination of computation techniques have been used to determine the binding site and amino acid residues on IgE that are critical for binding to a therapeutic anti-IgE. Homology modeling was used to model parts of IgE and of an antibody that binds to IgE. Docking simulations using shape descriptions were then carried out using the models to determine which residues in the IgE are involved in the binding interaction. The anti-IgE has been determined to bind close to some of the residues that are believed to be in the Fc(epsilon)RI receptor site on IgE, therefore preventing IgE from binding to Fc(epsilon)RI on mast cells and basophils and causing the release of pharmacologic mediators from these cells. Experiments have been suggested to verify the binding site and the residues involved in binding to anti-IgE.

Amino Acid Sequence↗

PubMed2002

Electron transfer between cytochrome b(5) and some oxidised haemoglobins: the role of ionic strength.

We have compared experimental measurements and Brownian dynamic calculations for the reduction of oxidised adult human haemoglobin by reduced bovine cytochrome b(5) over a range of ionic strengths. Our calculations suggest that the presence of molecular electrostatic fields have a significant role to play in the formation of the electron transfer complexes. These results predict that electron transfer occurs within an ensemble of similarly weakly docked complexes, the formation of which is strongly ionic strength dependent. Application of electron tunneling analysis to the complexes allows us to predict the rates of electron transfer within each ensemble of complexes as a function of ionic strength. The outcome of this theoretical study is compared with experimental rate measurements. A comparison of the results obtained from adult and embryonic haemoglobins, at a fixed ionic strength, indicates a significant difference in the characteristics of complex formation. These data emphasise the role played by electrostatic interactions in this important physiological reaction.

Animals↗

PubMed2004

Identification of protein-protein interaction sites from docking energy landscapes.

Protein recognition is one of the most challenging and intriguing problems in structural biology. Despite all the available structural, sequence and biophysical information about protein-protein complexes, the physico-chemical patterns, if any, that make a protein surface likely to be involved in protein-protein interactions, remain elusive. Here, we apply protein docking simulations and analysis of the interaction energy landscapes to identify protein-protein interaction sites. The new protocol for global docking based on multi-start global energy optimization of an all-atom model of the ligand, with detailed receptor potentials and atomic solvation parameters optimized in a training set of 24 complexes, explores the conformational space around the whole receptor without restrictions. The ensembles of the rigid-body docking solutions generated by the simulations were subsequently used to project the docking energy landscapes onto the protein surfaces. We found that highly populated low-energy regions consistently corresponded to actual binding sites. The procedure was validated on a test set of 21 known protein-protein complexes not used in the training set. As much as 81% of the predicted high-propensity patch residues were located correctly in the native interfaces. This approach can guide the design of mutations on the surfaces of proteins, provide geometrical details of a possible interaction, and help to annotate protein surfaces in structural proteomics.

Binding Sites↗

PubMed1996

Crystal structure implies that cyclophilin predominantly catalyzes the trans to cis isomerization.

The crystal structure of human recombinant cyclophilin A complexed with a substrate of succinyl-Ala-Ala-Pro-Phe-p-nitroanilide (AAPF) has been determined and refined to an R-factor of 0.189 at 2.4 A resolution. The structure revealed only the cis form of the substrate bound to cyclophilin A in a stoichiometry of 1:1. This binding ratio is different from the structure of cyclophilin A complexed with the tetrapeptide N-acetyl-Ala-Ala-Pro-Ala-amidomethylcourmarin. Model docking revealed that the trans form of AAPF does not fit into the active site. The observation that only the trans cis form of AAPF binds to cyclophilin A implies that cyclophilin A predominantly catalyzes the trans to cis isomerization of a peptidylprolyl amide bond. On the basis of the structure, it is proposed that Arg55 hydrogen-bonds to the nitrogen to deconjugate the resonance of the prolyl amide bond and thus facilitates the cis-trans rotation.

Amino Acid Isomerases↗

PubMed1999

The conformations of a substrate and a product bound to the active site of S-adenosylmethionine synthetase.

S-Adenosylmethionine (AdoMet) is the most widely used alkyl group donor in biological systems. The formation of AdoMet from ATP and L-methionine is catalyzed by S-adenosylmethionine synthetase (AdoMet synthetase). Elucidation of the conformations of enzyme-bound substrates, product, and inhibitors is important for the understanding of the catalytic mechanism of the enzyme and the design of new inhibitors. To obtain structural data for enzyme-bound substrates and product, we have used two-dimensional transferred nuclear Overhauser effect spectroscopy to determine the conformation of enzyme-bound AdoMet and 5'-adenylyl imidodiphosphate (AMPPNP). AMPPNP, an analogue of ATP, is resistant to the ATP hydrolysis activity of AdoMet synthetase because of the presence of a nonhydrolyzable NH-link between the beta- and gamma-phosphates but is a substrate for AdoMet formation during which tripolyphosphate is produced. AdoMet and AMPPNP both bind in an anti conformation about the glycosidic bond. The ribose rings are in C3'-exo and C4'-exo conformations in AdoMet and AMPPNP, respectively. The differences in ribose ring conformations presumably reflect the different steric requirements of the C5' substituents in AMPPNP and AdoMet. The NMR-determined conformations of AdoMet and AMPPNP were docked into the E. coli AdoMet synthetase active site taken from the enzyme.ADP. Pi crystal structure. Since there are no nonexchangeable protons either in the carboxy-terminal end of the methionine segment of AdoMet or in the tripolyphosphate segment of AMPPNP, these portions of the molecules were modeled into the enzyme active site. The interactions of AdoMet and AMPPNP with the enzyme predict the location of the methionine binding site and suggest how the positive charge formed on the sulfur during AdoMet synthesis is stabilized.

Adenylyl Imidodiphosphate↗

PubMedDate not supplied

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↗

PubMed2000

Fast prediction and visualization of protein binding pockets with PASS.

PASS (Putative Active Sites with Spheres) is a simple computational tool that uses geometry to characterize regions of buried volume in proteins and to identify positions likely to represent binding sites based upon the size, shape, and burial extent of these volumes. Its utility as a predictive tool for binding site identification is tested by predicting known binding sites of proteins in the PDB using both complexed macromolecules and their corresponding apoprotein structures. The results indicate that PASS can serve as a front-end to fast docking. The main utility of PASS lies in the fact that it can analyze a moderate-size protein (approximately 30 kDa) in under 20 s, which makes it suitable for interactive molecular modeling, protein database analysis, and aggressive virtual screening efforts. As a modeling tool, PASS (i) rapidly identifies favorable regions of the protein surface, (ii) simplifies visualization of residues modulating binding in these regions, and (iii) provides a means of directly visualizing buried volume, which is often inferred indirectly from curvature in a surface representation. PASS produces output in the form of standard PDB files, which are suitable for any modeling package, and provides script files to simplify visualization in Cerius2, InsightII, MOE, Quanta, RasMol, and Sybyl. PASS is freely available to all.

Algorithms↗

PubMedDate not supplied

Determination of alpha-conotoxin binding modes on neuronal nicotinic acetylcholine receptors.

alpha-Conotoxins, from cone snails, and alpha-neurotoxins, from snakes, are competitive inhibitors of nicotinic acetylcholine receptors (nAChRs) that have overlapping binding sites in the ACh binding pocket. These disulphide-rich peptides are used extensively as tools to localize and pharmacologically characterize specific nAChRs subtypes. Recently, a homology model based on the high-resolution structure of an ACh binding protein (AChBP) allowed the three-fingered alpha-neurotoxins to be docked onto the alpha7 nAChR. To investigate if alpha-conotoxins interact with the nAChR in a similar manner, we built homology models of human alpha7 and alpha3beta2 nAChRs, and performed docking simulations of alpha-conotoxins ImI, PnIB, PnIA and MII using the program GOLD. Docking revealed that alpha-conotoxins have a different mode of interaction compared with alpha-neurotoxins, with surprisingly few nAChR residues in common between their overlapping binding sites. These docking experiments show that ImI and PnIB bind to the ACh binding pocket via a small cavity located above the beta9/beta10 hairpin of the (+)alpha7 nAChR subunit. Interestingly, PnIB, PnIA and MII were found to bind in a similar location on alpha7 or alpha3beta2 receptors mostly through hydrophobic interactions, while ImI bound further from the ACh binding pocket, mostly through electrostatic interactions. These findings, which distinguish alpha-conotoxin and alpha-neurotoxin binding modes, have implications for the rational design of selective nAChR antagonists.

Amino Acid Sequence↗

PubMed2002

The ternary complex of cytochrome f and cytochrome c: identification of a second binding site and competition for plastocyanin binding.

The complex of yeast cytochrome c and cytochrome f from the cyanobacterium Phormidium laminosum was investigated by NMR spectroscopy. Chemical shift perturbation analysis reveals that residues around the haem edge of cytochrome c are involved in the complex interface. Binding curves derived from an NMR spectroscopy titration at 10 mM ionic strength indicate that there are two sites for cytochrome c with binding constants of approximately 2 x 10(4) M(-1) and 4 x 10(3) M(-1). A protein docking simulation with NMR-derived constraints identifies two sites, at the front (Site I) and back faces (Site II) of the haem region of cytochrome f. Site I is homologous to the binding site previously determined for the natural cytochrome f partner plastocyanin. Site II may represent the binding site for the Rieske protein in the cytochrome bf complex. Cytochrome c and plastocyanin are shown to compete for binding at Site I. The competition appears to involve electrostatic screening rather than simple steric occlusion of the binding site.

Amides↗