Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Docking”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

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

At least 163 records · Page 9Linked to original sources

Hierarchical docking of databases of multiple ligand conformations.

Ligand flexibility is an important problem in molecular docking and virtual screening. To address this challenge, we investigate a hierarchical pre-organization of multiple conformations of small molecules. Such organization of pre-calculated conformations removes the exploration of ligand conformational space from the docking calculation and allows for concise representation of what can be thousands of conformations. The hierarchy also recognizes and prunes incompatible conformations early in the calculation, eliminating redundant calculations of fit. We investigate the method by docking the MDL Drug Data Report (MDDR), an annotated database of 100,000 molecules, into apo and holo forms of seven unrelated targets. This annotated database allows us to track the ranking of tens to hundreds of annotated ligands in each of the docking systems. The binding sites and database are prepared in an automated fashion in an attempt to remove some human bias from the calculations. Many thousands of explicit and implicit ligand conformations may be docked in calculations not much longer than required for single conformer docking. As long as internal energies are not considered, recombination with the hierarchy is additive as the number of degrees of freedom is increased. Molecules with even millions of conformations can be docked in a few minutes on a single desktop computer.

Computer Simulation↗

The FlexX database docking environment--rational extraction of receptor based pharmacophores.

We present an integrated docking environment that allows for iterative and interactive detailed analysis of many docking solutions. All docking information is stored in an ORACLE database. New scoring schemes (e.g. target-specific scoring functions) as well as various types of filters can be easily defined and tested within this environment. As an example application we investigated the validity of the following hypothesis: If a docking procedure can lead to enrichments significantly better than random then a bias towards (partially) correct placements should be detectable. Such bias in terms of a preference for certain interacting groups within the active site can be used to select a set of receptor-based pharmacophore constraints, which in turn might be used to enhance the docking procedure. As a proof of concept for this approach we performed docking studies on three targets: thrombin, the cyclin-dependent kinase 2 (CDK2) and the angiotensin converting enzyme (ACE). We docked a set of known active compounds with standard FlexX and derived three sets of target-specific receptor-based pharmacophore constraints by statistical analysis of the predicted placements. Applying these receptor-based constraints in a virtual screening protocol utilizing FlexX-Pharm led to significantly improved enrichments.

Calibration↗

Effect of docking on health, blood cells and metabolites and growth of Dorset lambs.

Docking by knife or elastrator had no sustained effect on growth of lambs 0 to 8 wk postpartum. However, during the 2 wk period postdocking, growth of lambs was limited or depressed, regardless of whether lambs were docked by knife at 14, 28 or 42 d of age. Blood hematocrit and hemoglobin were decreased (P less than .01) 1 d after docking by knife and 14 d after docking by elastrator. Neither method of docking proved superior. Regardless of whether lambs were docked by knife or elastrator, tail tissue was destroyed and became necrotic. Numbers of white blood cells increased (P less than .01) in response to docking. This necessitates that lambs be observed often through 2 wk postdocking and that medical treatment be administered when needed.

Age Factors↗

The development of a CDK2-docking site peptide that inhibits p53 and sensitizes cells to death.

Cyclin-dependent protein kinases play important roles in cell cycle progression and are attractive targets for the design of anti-proliferative drugs. Two distinct synthetic CDK1/2 inhibitors, Roscovitine and NU2058, are pharmacologically distinct in their ability to modify p53-dependent transcription and perturb cell cycle progression. Although such active-site CDK1/2 inhibitors comprise the most standard type of enzyme inhibitor, many protein kinases are proving to harbour high affinity docking sites that may provide a potentially novel interface for the design of kinase-inhibitors. We examined whether CDK2 has a docking site for its oligomeric substrate p53, whether small-peptide leads can be developed that inhibit CDK2 function, and whether such peptide-inhibitors are pharmacologically distinct from Roscovitine or NU2058. A docking site for CDK2 was identified in the tetramerization domain of p53 at a site that is distinct from the phospho-acceptor site. Peptides derived from the tetramerization domain of p53 block CDK2 phosphorylation and identification of critical CDK2 contacts in the tetramerization domain of p53 suggest that kinase docking does not require tetramerization of the substrate. Transient transfection assays were developed to show that the GFP-CDK2 docking site fusion protein (GFP-CIP) attenuates p53 activity in vivo and suppresses p21WAF1 induction which is similar to NU2058 but distinct from Roscovitine. A stable cell line with an inducible GFP-CIP gene attenuates p53 activity and induces significant cell death in a drug-resistant melanoma cell line, sensitizes cells to death induced by Doxorubicin, and suppresses cell growth in a colony formation assay. These data indicate that CDK2, in addition to cyclin A, can have a high affinity docking site for a substrate and highlights the possibility that CDK2 docking sites may represent effective targets for inhibitor design.

Antineoplastic Agents↗

Morphologically docked synaptic vesicles are reduced in synaptotagmin mutants of Drosophila.

Nerve terminal specializations include mechanisms for maintaining a subpopulation of vesicles in a docked, fusion-ready state. We have investigated the relationship between synaptotagmin and the number of morphologically docked vesicles by an electron microscopic analysis of Drosophila synaptotagmin (syt) mutants. The overall number of synaptic vesicles in a terminal was reduced, although each active zone continued to have a cluster of vesicles in its vicinity. In addition, there was an increase in the number of large vesicles near synapses. Examining the clusters, we found that the pool of synaptic vesicles immediately adjacent to the presynaptic membrane, the pool that includes the docked population, was reduced to 24 +/- 5% (means +/- SEM) of control in the sytnull mutation. To separate contributions of overall vesicle depletion and increased spontaneous release from direct effects of synaptotagmin on morphological docking, we examined syt mutants in an altered genetic background. Recombining syt alleles onto a second chromosome bearing an as yet uncharacterized mutation resulted in the expected decrease in evoked release but suppressed the increase in spontaneous release frequency. Motor nerve terminals in this genotype contained more synaptic vesicles than control, yet the number of vesicles immediately adjacent to the presynaptic membrane near active zones was still reduced (33 +/- 4% of control). Our findings demonstrate that there is a decrease in the number of morphologically docked vesicles seen in syt mutants. The decreases in docking and evoked release are independent of the increase in spontaneous release. These results support the hypothesis that synaptotagmin stabilizes the docked state.

Animals↗

CASP2 molecular docking predictions with the LIGIN software.

Seven docking predictions were made with the LIGIN program. In six cases the location of the binding pocket was identified correctly by systematically docking everywhere within the protein structure. In two cases the ligand was docked to within 1.8 A RMSD of the experimentally determined structure. LIGIN has not been optimized to deal with highly flexible ligands that dock at the surface of proteins. Consequently, in three cases the exposed part of the ligand was docked poorly, although the buried parts were docked well, and made similar atomic contacts with the protein as in the experimentally determined structure.

Amiloride↗

Evaluation of the FLEXX incremental construction algorithm for protein-ligand docking.

We report on a test of FLEXX, a fully automatic docking tool for flexible ligands, on a highly diverse data set of 200 protein-ligand complexes from the Protein Data Bank. In total 46.5% of the complexes of the data set can be reproduced by a FLEXX docking solution at rank 1 with an rms deviation (RMSD) from the observed structure of less than 2 A. This rate rises to 70% if one looks at the entire generated solution set. FLEXX produces reliable results for ligands with up to 15 components which can be docked in 80% of the cases with acceptable accuracy. Ligands with more than 15 components tend to generate wrong solutions more often. The average runtime of FLEXX on this test set is 93 seconds per complex on a SUN Ultra-30 workstation. In addition, we report on "cross-docking" experiments, in which several receptor structures of complexes with identical proteins have been used for docking all cocrystallized ligands of these complexes. In most cases, these experiments show that FLEXX can acceptably dock a ligand into a foreign receptor structure. Finally we report on screening runs of ligands out of a library with 556 entries against ten different proteins. In eight cases FLEXX is able to find the original inhibitor within the top 7% of the total library.

Algorithms↗

Automated docking of alpha-(1-->4)- and alpha-(1-->6)-linked glucosyl trisaccharides and maltopentaose into the soybean beta-amylase active site.

The Lamarckian genetic algorithm of AutoDock 3.0 was used to dock alpha-maltotriose, methyl alpha-panoside, methyl alpha-isopanoside, methyl alpha-isomaltotrioside, methyl alpha-(6(1)-alpha-glucopyranosyl)-maltoside, and alpha-maltopentaose into the closed and, except for alpha-maltopentaose, into the open conformation of the soybean beta-amylase active site. In the closed conformation, the hinged flap at the mouth of the active site closes over the substrate. The nonreducing end of alpha-maltotriose docks preferentially to subsites -2 or +1, the latter yielding nonproductive binding. Some ligands dock into less optimal conformations with the nonreducing end at subsite -1. The reducing-end glucosyl residue of nonproductively-bound alpha-maltotriose is close to residue Gln194, which likely contributes to binding to subsite +3. In the open conformation, the substrate hydrogen-bonds with several residues of the open flap. When the flap closes, the substrate productively docks if the nonreducing end is near subsites -2 or -1. Trisaccharides with alpha-(1-->6) bonds do not successfully dock except for methyl alpha-isopanoside, whose first and second glucosyl rings dock exceptionally well into subsites -2 and -1. The alpha-(1-->6) bond between the second and third glucosyl units causes the latter to be improperly positioned into subsite +1; the fact that isopanose is not a substrate of beta-amylase indicates that binding to this subsite is critical for hydrolysis.

Algorithms↗

Molecular docking towards drug discovery.

Fueled by advances in molecular structure determination, tools for structure-based drug design are proliferating rapidly. Lead discovery through searching of ligand databases with molecular docking techniques represents an attractive alternative to high-throughout random screening. The size of commercial databases imposes severe computational constraints on molecular docking, compromising the level of calculational detail permitted for each putative ligand. We describe alternative philosophies for docking which effectively address this challenge. With respect to the dynamic aspects of molecular recognition, these strategies lie along a spectrum of models bounded by the Lock-and-Key and Induced-Fit theories for ligand binding. We explore the potential of a rigid model in exploiting species specificity and of a tolerant model in predicting absolute ligand binding affinity. Current molecular docking methods are limited primarily by their ability to rank docked complexes; we therefore place particular emphasis on this aspect of the problem throughout our validation of docking strategies.

Animals↗

Docking molecules by families to increase the diversity of hits in database screens: computational strategy and experimental evaluation.

Molecular docking programs screen chemical databases for novel ligands that fit protein binding sites. When one compound fits the site well, close analogs typically do the same. Therefore, many of the compounds that are found in such screens resemble one another. This reduces the variety and novelty of the compounds suggested. In an attempt to increase the diversity of docking hit lists, the Available Chemicals Directory was grouped into families of related structures. All members of every family were docked and scored, but only the best scoring molecule of a high-ranking family was allowed in the hit list. The identity and scores of the other members of these families were recorded as annotations to the best family member, but they were not independently ranked. This family-based docking method was compared with molecule-by-molecule docking in screens against the structures of thymidylate synthase, dihydrofolate reductase (DHFR), and the cavity site of the mutant T4 lysozyme Leu99 --> Ala (L99A). In each case, the diversity of the hit list increased, and more families of known ligands were found. To investigate whether the newly identified hits were sensible, we tested representative examples experimentally for binding to L99A and DHFR. Of the six compounds tested against L99A, five bound to the internal cavity. Of the seven compounds tested against DHFR, six inhibited the enzyme with apparent K(i) values between 0.26 and 100 microM. The segregation of potential ligands into families of related molecules is a simple technique to increase the diversity of candidates suggested by database screens. The general approach should be applicable to most docking methods. Proteins 2001;42:279-293.

Algorithms↗

Ligand-protein inverse docking and its potential use in the computer search of protein targets of a small molecule.

Ligand-protein docking has been developed and used in facilitating new drug discoveries. In this approach, docking single or multiple small molecules to a receptor site is attempted to find putative ligands. A number of studies have shown that docking algorithms are capable of finding ligands and binding conformations at a receptor site close to experimentally determined structures. These algorithms are expected to be equally applicable to the identification of multiple proteins to which a small molecule can bind or weakly bind. We introduce a ligand-protein inverse-docking approach for finding potential protein targets of a small molecule by the computer-automated docking search of a protein cavity database. This database is developed from protein structures in the Protein Data Bank (PDB). Docking is conducted with a procedure involving multiple-conformer shape-matching alignment of a molecule to a cavity followed by molecular-mechanics torsion optimization and energy minimization on both the molecule and the protein residues at the binding region. Scoring is conducted by the evaluation of molecular-mechanics energy and, when applicable, by the further analysis of binding competitiveness against other ligands that bind to the same receptor site in at least one PDB entry. Testing results on two therapeutic agents, 4H-tamoxifen and vitamin E, showed that 50% of the computer-identified potential protein targets were implicated or confirmed by experiments. The application of this approach may facilitate the prediction of unknown and secondary therapeutic target proteins and those related to the side effects and toxicity of a drug or drug candidate. Proteins 2001;43:217-226.

Computer Simulation↗

Gaussian docking functions.

A shape-based Gaussian docking function is constructed which uses Gaussian functions to represent the shapes of individual atoms. A set of 20 trypsin ligand-protein complexes are drawn from the Protein Data Bank (PDB), the ligands are separated from the proteins, and then are docked back into the active sites using numerical optimization of this function. It is found that by employing this docking function, quasi-Newton optimization is capable of moving ligands great distances [on average 7 A root mean square distance (RMSD)] to locate the correctly docked structure. It is also found that a ligand drawn from one PDB file can be docked into a trypsin structure drawn from any of the trypsin PDB files. This implies that this scoring function is not limited to more accurate x-ray structures, as is the case for many of the conventional docking methods, but could be extended to homology models.

Binding Sites↗

The structural basis for docking in modular polyketide biosynthesis.

Polyketide natural products such as erythromycin and rapamycin are assembled on polyketide synthases (PKSs), which consist of modular sets of catalytic activities distributed across multiple protein subunits. Correct protein-protein interactions among the PKS subunits which are critical to the fidelity of biosynthesis are mediated in part by "docking domains" at the termini of the proteins. The NMR solution structure of a representative docking domain complex from the erythromycin PKS (DEBS) was recently solved, and on this basis it has been proposed that PKS docking is mediated by the formation of an intermolecular four-alpha-helix bundle. Herein, we report the genetic engineering of such a docking domain complex by replacement of specific helical segments and analysis of triketide synthesis by mutant PKSs in vivo. The results of these helix swaps are fully consistent with the model and highlight residues in the docking domains that may be targeted to alter the efficiency or specificity of subunit-subunit docking in hybrid PKSs.

Amino Acid Sequence↗

An efficient molecular docking using conformational space annealing.

Molecular docking falls into the general category of global optimization problems because its main purpose is to find the most stable complex consisting of a receptor and its ligand. Conformational space annealing (CSA), a powerful global optimization method, is incorporated with the Tinker molecular modeling package to perform molecular docking simulations of six receptor-ligand complexes (3PTB, 1ULB, 2CPP, 1STP, 3CPA, and 1PPH) from the Protein Data Bank. In parallel, Monte Carlo with the minimization (MCM) method is also incorporated into the Tinker package for comparison. The energy function, consisting of electrostatic interactions, van der Waals interactions, and torsional energy terms, is calculated using the AMBER94 all-atom empirical force field. Rigid docking simulations for all six complexes and flexible docking simulations for three complexes (1STP, 3CPA, and 1PPH) are carried out using the CSA and the MCM methods. The simulation results show that the docking procedures using the CSA method generally find the most stable complexes as well as the native-like complexes more efficiently and accurately than those using the MCM, demonstrating that CSA is a promising search method for molecular docking problems.

Algorithms↗

New and fast statistical-thermodynamic method for computation of protein-ligand binding entropy substantially improves docking accuracy.

We present a novel method to estimate the contributions of translational and rotational entropy to protein-ligand binding affinity. The method is based on estimates of the configurational integral through the sizes of clusters obtained from multiple docking positions. Cluster sizes are defined as the intervals of variation of center of ligand mass and Euler angles in the cluster. Then we suggest a method to consider the entropy of torsional motions. We validate the suggested methods on a set of 135 PDB protein-ligand complexes by comparing the averaged root-mean square deviations (RMSD) of the top-scored ligand docked positions, accounting and not accounting for entropy contributions, relative to the experimentally determined positions. We demonstrate that the method increases docking accuracy by 10-21% when used in conjunction with the AutoDock docking program, thus reducing the percent of incorrectly docked ligands by 1.4-fold to four-fold, so that in some cases the percent of ligands correctly docked to within an RMSD of 2 A is above 90%. We show that the suggested method to account for entropy of relative motions is identical to the method based on the Monte Carlo integration over intervals of variation of center of ligand mass and Euler angles in the cluster.

Algorithms↗

Adaptation of a fast Fourier transform-based docking algorithm for protein design.

Designing proteins with novel protein/protein binding properties can be achieved by combining the tools that have been developed independently for protein docking and protein design. We describe here the sequence-independent generation of protein dimer orientations by protein docking for use as scaffolds in protein sequence design algorithms. To dock monomers into sequence-independent dimer conformations, we use a reduced representation in which the side chains are approximated by spheres with atomic radii derived from known C2 symmetry-related homodimers. The interfaces of C2-related homodimers are usually more hydrophobic and protein core-like than the interfaces of heterodimers; we parameterize the radii for docking against this feature to capture and recreate the spatial characteristics of a hydrophobic interface. A fast Fourier transform-based geometric recognition algorithm is used for docking the reduced representation protein models. The resulting docking algorithm successfully predicted the wild-type homodimer orientations in 65 out of 121 dimer test cases. The success rate increases to approximately 70% for the subset of molecules with large surface area burial in the interface relative to their chain length. Forty-five of the predictions exhibited less than 1 A C(alpha) RMSD compared to the native X-ray structures. The reduced protein representation therefore appears to be a reasonable approximation and can be used to position protein backbones in plausible orientations for homodimer design.

Algorithms↗

Atomic solvation parameters in the analysis of protein-protein docking results.

Several sets of amino acid surface areas and transfer free energies were used to derive a total of nine sets of atomic solvation parameters (ASPs). We tested the accuracy of each of these sets of parameters in predicting the experimentally determined transfer free energies of the amino acid derivatives from which the parameters were derived. In all cases, the calculated and experimental values correlated well. We then chose three parameter sets and examined the effect of adding an energetic correction for desolvation based on these three parameter sets to the simple potential function used in our multiple start Monte Carlo docking method. A variety of protein-protein interactions and docking results were examined. In the docking simulations studied, the desolvation correction was only applied during the final energy calculation of each simulation. For most of the docking results we analyzed, the use of an octanol-water-based ASP set marginally improved the energetic ranking of the low-energy dockings, whereas the other ASP sets we tested disturbed the ranking of the low-energy dockings in many of the same systems. We also examined the correlation between the experimental free energies of association and our calculated interaction energies for a series of proteinase-inhibitor complexes. Again, the octanol-water-based ASP set was compatible with our standard potential function, whereas ASP sets derived from other solvent systems were not.

Amino Acids↗

Automated generation of MCSS-derived pharmacophoric DOCK site points for searching multiconformation databases.

All docking methods employ some sort of heuristic to orient the ligand molecules into the binding site of the target structure. An automated method, MCSS2SPTS, for generating chemically labeled site points for docking is presented. MCSS2SPTS employs the program Multiple Copy Simultaneous Search (MCSS) to determine target-based theoretical pharmacophores. More specifically, chemically labeled site points are automatically extracted from selected low-energy functional-group minima and clustered together. These pharmacophoric site points can then be directly matched to the pharmacophoric features of database molecules with the use of either DOCK or PhDOCK to place the small molecules into the binding site. Several examples of the ability of MCSS2SPTS to reproduce the three-dimensional pharmacophoric features of ligands from known ligand-protein complex structures are discussed. In addition, a site-point set calculated for one human immunodeficiency virus 1 (HIV1) protease structure is used with PhDOCK to dock a set of HIV1 protease ligands; the docked poses are compared to the corresponding complex structures of the ligands. Finally, the use of an MCSS2SPTS-derived site-point set for acyl carrier protein synthase is compared to the use of atomic positions from a bound ligand as site points for a large-scale DOCK search. In general, MCSS2SPTS-generated site points focus the search on the more relevant areas and thereby allow for more effective sampling of the target site.

Binding Sites↗