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 181 records · Page 10Linked to original sources

Molecular determinants of thapsigargin binding by SERCA Ca2+-ATPase: a computational docking study.

Thapsigargin (TG) is a potent and commonly used inhibitor of the ion transport activity of sarco/endoplasmic reticulum Ca2+-ATPases (SERCA). Based on the recently published crystal structures of rabbit muscle SERCA1a in the Ca2+/E1 (E1) and TG/E2 (E2) conformations, we performed computational docking studies to characterize the molecular interactions that govern binding of TG and TG-analogs by the enzyme. Using the program GOLD (genetic optimization for ligand docking) in combination with the scoring function ChemScore, TG was docked into the binding site of the E1 and E2 conformations of SERCA1a. The docking results revealed a consensus ligand-binding mode consistent with the crystal structure and showed that hydrophobic interactions are the primary driving force of TG binding by SERCA. Moreover, it was shown that the conformational changes accompanying the E2 to E1 transition in the enzyme likely displace TG from its favored orientation in the binding site, thereby substantially reducing its binding affinity. This finding illustrates on the molecular level how TG may exert its inhibitory effect in binding tightly to the E2 form and preventing it from converting into its E1 form, a requirement for catalytic function. We also docked 9 TG analogs into the E2 conformation of the enzyme. Eight of the analogs adopted a binding mode very similar to that of TG, whereas one compound preferred a different orientation in the binding site. Analysis of the predicted binding affinities showed a good correlation with the experimentally observed inhibitory potencies of the analogs. Docking was also performed with several modeled mutants of SERCA1a, whose phenylalanine residue in position 256 (Phe256) had been modified. The experimentally observed declines in TG sensitivity in most of the Phe256 mutants was qualitatively accounted for and appears, at least in part, be due to a slightly altered TG-binding mode.

Algorithms↗

The impact of protein flexibility on protein-protein docking.

Accounting for protein flexibility in protein-protein docking algorithms is challenging, and most algorithms therefore treat proteins as rigid bodies or permit side-chain motion only. While the consequences are obvious when there are large conformational changes upon binding, the situation is less clear for the modest conformational changes that occur upon formation of most protein-protein complexes. We have therefore studied the impact of local protein flexibility on protein-protein association by means of rigid body and torsion angle dynamics simulation. The binding of barnase and barstar was chosen as a model system for this study, because the complexation of these 2 proteins is well-characterized experimentally, and the conformational changes accompanying binding are modest. On the side-chain level, we show that the orientation of particular residues at the interface (so-called hotspot residues) have a crucial influence on the way contacts are established during docking from short protein separations of approximately 5 A. However, side-chain torsion angle dynamics simulations did not result in satisfactory docking of the proteins when using the unbound protein structures. This can be explained by our observations that, on the backbone level, even small (2 A) local loop deformations affect the dynamics of contact formation upon docking. Complementary shape-based docking calculations confirm this result, which indicates that both side-chain and backbone levels of flexibility influence short-range protein-protein association and should be treated simultaneously for atomic-detail computational docking of proteins.

Computer Simulation↗

Classification of protein complexes based on docking difficulty.

Based on the results of several groups using different docking methods, the key properties that determine the expected success rate in protein-protein docking calculations are measures of conformational change, interface area, and hydrophobicity. A classification of protein complexes in terms of these measures provides a prediction of docking difficulty. This classification is used to study the targets of the CAPRI docking experiment. Results show that targets with a moderate expected difficulty were indeed predicted well by a number of groups, whereas the use of additional a priori information was necessary to obtain good results for some very difficult targets. The analysis indicates that CAPRI and other relatively large-scale docking studies represent very important steps toward understanding the capabilities and limitations of current protein-protein docking methods.

Algorithms↗

Docking to single-domain and multiple-domain proteins: old and new challenges.

The diverse selection of targets in the CAPRI experiments provides grounds for determining the limits of our rigid-body docking program MolFit, and for extending it. We find that the sensitivity of MolFit is high, enabling it to produce reasonably accurate docking solutions when the structures undergo moderate local conformation changes upon complex formation or when the docked molecules are modeled. Yet the ranks of these solutions are sometimes too low to meet the requirements of CAPRI assessment. This indicates that the selectivity of MolFit, which was optimized for docking of unbound X-ray structures, and which relies on the availability of external data from biochemical and bioinformatic sources, needs readjustment in order to meet the challenges presented by NMR or modeled structures. A different challenge is presented by large global conformation changes such as movements of domains. We show that such changes can be accommodated within the rigid-body approximation by employing rigid multibody multistage docking procedures. We also address the difficulty of ranking results from 2-body and multibody docking scans in cases in which there are no external data favoring one option over the other.

Algorithms↗

Force calculations in automated docking: enzyme-substrate interactions in Fusarium oxysporum Cel7B.

AutoDock is a small-molecule docking program that uses an energy function to score docked ligands. Here AutoDock's grid-based method for energy evaluation was exploited to evaluate the force exerted by Fusarium oxysporum Cel7B on the atoms of docked cellooligosaccharides and a thiooligosaccharide substrate analog. Coupled with the interaction energies evaluated for each docked ligand, these forces give insight into the dynamics of the ligand in the active site, and help to elucidate the relative importance of specific enzyme-substrate interactions in stabilizing the substrate transition-state conformation. The processive force on the docked substrate in the F. oxysporum Cel7B active site is less than half of that on the docked substrate in the Hypocrea jecorina Cel7A active site. Hydrogen bonding interactions of the enzyme with the C2 hydroxyl group of the glucosyl residue in subsite -2 and with the C3 hydroxyl group of the glucosyl residue in subsite +1 are the most significant in stabilizing the distorted14B transition-state conformation of the glucosyl residue in subsite -1. The force calculations also help to elucidate the mechanism that prevents the active site from fouling.

Cellulase↗

Accounting for loop flexibility during protein-protein docking.

Although reliable docking can now be achieved for systems that do not undergo important induced conformational change upon association, the presence of flexible surface loops, which must adapt to the steric and electrostatic properties of a partner, generally presents a major obstacle. We report here the first docking method that allows large loop movements during a systematic exploration of the possible arrangements of the two partners in terms of position and rotation. Our strategy consists in taking into account an ensemble of possible loop conformations by a multi-copy representation within a reduced protein model. The docking process starts from regularly distributed positions and orientations of the ligand around the whole receptor. Each starting configuration is submitted to energy minimization during which the best-fitting loop conformation is selected based on the mean-field theory. Trials were carried out on proteins with significant differences in the main-chain conformation of the binding loop between isolated form and complexed form, which were docked to their partner considered in their bound form. The method is able to predict complexes very close to the crystal complex both in terms of relative position of the two partners and of the geometry of the flexible loop. We also show that introducing loop flexibility on the isolated protein form during systematic docking largely improves the predictions of relative position of the partners in comparison with rigid-body docking.

Binding Sites↗

Protein-ligand docking: current status and future challenges.

Understanding the ruling principles whereby protein receptors recognize, interact, and associate with molecular substrates and inhibitors is of paramount importance in drug discovery efforts. Protein-ligand docking aims to predict and rank the structure(s) arising from the association between a given ligand and a target protein of known 3D structure. Despite the breathtaking advances in the field over the last decades and the widespread application of docking methods, several downsides still exist. In particular, protein flexibility-a critical aspect for a thorough understanding of the principles that guide ligand binding in proteins-is a major hurdle in current protein-ligand docking efforts that needs to be more efficiently accounted for. In this review the key concepts of protein-ligand docking methods are outlined, with major emphasis being given to the general strengths and weaknesses that presently characterize this methodology. Despite the size of the field, the principal types of search algorithms and scoring functions are reviewed and the most popular docking tools are briefly depicted. Recent advances that aim to address some of the traditional limitations associated with molecular docking are also described. A selection of hand-picked examples is used to illustrate these features.

Algorithms↗

ROSETTALIGAND: protein-small molecule docking with full side-chain flexibility.

Protein-small molecule docking algorithms provide a means to model the structure of protein-small molecule complexes in structural detail and play an important role in drug development. In recent years the necessity of simulating protein side-chain flexibility for an accurate prediction of the protein-small molecule interfaces has become apparent, and an increasing number of docking algorithms probe different approaches to include protein flexibility. Here we describe a new method for docking small molecules into protein binding sites employing a Monte Carlo minimization procedure in which the rigid body position and orientation of the small molecule and the protein side-chain conformations are optimized simultaneously. The energy function comprises van der Waals (VDW) interactions, an implicit solvation model, an explicit orientation hydrogen bonding potential, and an electrostatics model. In an evaluation of the scoring function the computed energy correlated with experimental small molecule binding energy with a correlation coefficient of 0.63 across a diverse set of 229 protein- small molecule complexes. The docking method produced lowest energy models with a root mean square deviation (RMSD) smaller than 2 A in 71 out of 100 protein-small molecule crystal structure complexes (self-docking). In cross-docking calculations in which both protein side-chain and small molecule internal degrees of freedom were varied the lowest energy predictions had RMSDs less than 2 A in 14 of 20 test cases.

Algorithms↗

Complementarity of hydrophobic properties in ATP-protein binding: a new criterion to rank docking solutions.

ATP is an important substrate of numerous biochemical reactions in living cells. Molecular recognition of this ligand by proteins is very important for understanding enzymatic mechanisms. Considerable insight into the problem may be gained via molecular docking simulations. At the same time, standard docking protocols are often insufficient to predict correct conformations for protein-ATP complexes. Thus, in most cases the native-like solutions can be found among the docking poses, but current scoring functions have only limited ability to discriminate them from false positives. To improve the selection of correct docking solutions obtained with the GOLD software, we developed a new ranking criterion specific for ATP-protein binding. The method is based on detailed analysis of the intermolecular interactions in 40 high-resolution 3D structures of ATP-protein complexes (the training set). We found that the most important factors governing this recognition are hydrogen-bonding, stacking between adenine and aromatic protein residues, and hydrophobic contacts between adenine and protein residues. To address the latter, we applied the formalism of 3D molecular hydrophobicity potential. The results obtained were used to construct an ATP-oriented scoring criterion as a linear combination of the terms describing these intermolecular interactions. The criterion was then validated using the test set of 10 additional ATP-protein complexes. As compared with the standard scoring functions, the new ranking criterion significantly improved the selection of correct docking solutions in both sets and allowed considerable enrichment at the top of the list containing docking poses with correct solutions.

Adenine↗

Automated docking in crystallography: analysis of the substrates of aconitase.

Automated docking of substrates to proteins of known structure aids the process of crystallographic analysis in two ways. First, automated docking can be used to generate a small number of starting models for substrates using only protein coordinates from an early stage of refinement. Second, automated docking provides a method for exploring aspects of catalysis that are inaccessible to crystallography by postulating binding modes of catalytic intermediates. This paper describes the use of automated docking to explore the binding of substrates to aconitase. The technique starts with a substrate molecule in an arbitrary configuration and position and finds favorable docked configurations in a (static) protein active site based on a molecular mechanics type force field. Using protein coordinates from an early stage of refinement of an aconitase-isocitrate complex, we successfully predicted the binding configuration of isocitrate. Four configurations were found, the energetically most favorable of which fit the observed electron density well and was used as a starting model for further refinement. Two configurations were found in citrate docking experiments, the second of which approximates the mode of substrate binding in an aconitase-nitrocitrate complex. We were also able to propose two binding modes of the catalytic intermediate cis-aconitate. These correspond closely to the isocitrate and the citrate binding modes. The relation of these new results to the proposed reaction mechanism is discussed.

Aconitate Hydratase↗

The role and application of in silico docking in chemical genomics research.

In silico docking techniques are being used to investigate the complementarity at the molecular level of a ligand and a protein target. As such, docking studies can be used to identify the structural features that are important for binding and for in silico screening efforts in which suitable binding partners can be identified. Here we describe a practical approach for setting up docking simulations using different docking programs. We also cover the analysis and rescoring of the obtained docking poses. Possible pitfalls in the docking studies are discussed and hints are provided to resolve commonly occurring problems.

Binding Sites↗

Orientational sampling and rigid-body minimization in molecular docking revisited: on-the-fly optimization and degeneracy removal.

Strategies for computational association of molecular components entail a compromise between configurational exploration and accurate evaluation. Following the work of Meng et al. [Proteins, 17 (1993) 266], we investigate issues related to sampling and optimization in molecular docking within the context of the DOCK program. An extensive analysis of diverse sampling conditions for six receptor-ligand complexes has enabled us to evaluate the tractability and utility of on-the-fly force-field score minimization, as well as the method for configurational exploration. We find that the sampling scheme in DOCK is extremely robust in its ability to produce configurations near to those experimentally observed. Furthermore, despite the heavy resource demands of refinement, the incorporation of a rigid-body, grid-based simplex minimizer directly into the docking process results in a docking strategy that is more efficient at retrieving experimentally observed configurations than docking in the absence of optimization. We investigate the capacity for further performance enhancement by implementing a degeneracy checking protocol aimed at circumventing redundant optimizations of geometrically similar orientations. Finally, we present methods that assist in the selection of sampling levels appropriate to desired result quality and available computational resources.

Algorithms↗

Correlation of syntaxin-1 and SNAP-25 clusters with docking and fusion of insulin granules analysed by total internal reflection fluorescence microscopy.

AIMS/HYPOTHESIS: The interaction of syntaxin-1 and SNAP-25 with insulin exocytosis was examined using the diabetic Goto-Kakizaki (GK) rat and a total internal reflection fluorescence (TIRF) imaging system. METHODS: Primary rat pancreatic beta cells were immunostained with anti-syntaxin-1A, anti-SNAP-25 and anti-insulin antibodies, and then observed by TIRF microscopy. The real-time image of GFP-labelled insulin granules motion was monitored by TIRF. RESULTS: The number of syntaxin-1A and SNAP-25 clusters, and the number of docked insulin granules on the plasma membrane were reduced in GK beta cells. When GK rats were treated with daily insulin injection for 2 weeks, the number of syntaxin-1 and SNAP-25 clusters was restored, along with the number of docked insulin granules. The infection of GK beta cells with Adex1CA SNAP-25 increased the number of docked insulin granules. TIRF imaging analysis demonstrated that the decreased number of fusion events from previously docked insulin granules in GK beta cells was restored when the number of docked insulin granules increased by insulin treatment or Adex1CA SNAP-25 infection. CONCLUSIONS/INTERPRETATION: There was a close correlation between the number of syntaxin-1 and SNAP-25 clusters and the number of docked insulin granules, which is associated with the fusion of insulin granules.

Animals↗

Secretory organelle docking at the cell membrane of Paramecium cells: dedocking and synchronized redocking of trichocysts.

We present the first evidence that secretory organelle docking at the cell membrane can be reversed in vivo. In nondischarge (nd) mutants of Paramecium tetraurelia all trichocysts can be detached from the cell surface within 2-3 h by different means, including cytochalasin B (but not D), high cell density, or Ca2+ ionophores. Considering the well-established ultrastructural differences between nd and wild-type (wt) cells, one can conclude that trichocyst docking at the cell periphery involves two docking sites (I, II): Site I ties the organelles to the epiplasm, and site II is the connection to the cell membrane at the fusogenic zone (expressed only in wt cells); both sites are close to the cell surface and only 150 nm apart. When the trigger for detachment of cortically docked trichocysts (high cell density, cytochalasin B) is relieved, trichocysts are synchronously reattached at the cell membrane, within 40-50 min, with a rate of 20-40 organelles/min, which far exceeds spontaneous docking rates. This is therefore also the first report on synchronization of secretory organelle docking. It is shown by radioactive leucine labeling that the same organelles are redocked, because trichocyst biogenesis is minimal under the conditions of de/redocking used. Surprisingly not only redocking but also detachment of trichocysts from the cell surface can be abolished by inhibitors of protein synthesis. Since Ca2+ ionophores mimic the effects of other conditions sufficient to detach trichocysts from the cell surface, we assume that a protein-dependent mechanism sensitive to Ca2+ (or other ions in exchange) may operate in trichocyst detachment. The precise mechanism involved in attachment or detachment of trichocysts remains to be elucidated.

Animals↗

Effects of local anesthesia and intravenous naloxone on the changes in behaviour and plasma concentrations of cortisol produced by castration and tail docking with tight rubber rings in young lambs.

The reliability of some behavioral and physiological indices used for the recognition and assessment of acute pain in lambs after castration and tail docking has been examined. Changes in the indices were measured after blocking neural activity with local anaesthetic (lignocaine) and after an opioid antagonist (naloxone) was administered. Six lambs, aged less than one week, were allocated randomly to each of six treatments. (i) control handling and blood sampling; (ii) castration plus tail docking with tight rubber rings; (iii) local anaesthesia; (iv) local anaesthesia followed by castration and tail docking; (v) intravenous naloxone only (0.2 mg kg-1); and (vi) intravenous naloxone followed by castration and tail docking. Local anaesthesia eliminated the behavioural and plasma cortisol changes which usually follow castration and tail docking. Naloxone had a limited effect on the increase in cortisol but altered the behaviour. The results support the view that such indices are useful for assessment of the response to acute pain and that, although endogenous opioids do reduce pain in young lambs after castration and tail docking, the effect is small.

Anesthesia, Local↗

A rab protein is required for the assembly of SNARE complexes in the docking of transport vesicles.

Rab proteins are generally required for transport vesicle docking. We have exploited yeast secretion mutants to demonstrate that a rab protein is required for v-SNAREs and t-SNAREs to assemble. The absence of the rab protein in the docking complex suggests that, in a broad sense, rab proteins participate in a reaction catalyzing SNARE complex assembly. In so doing, rab proteins could help impart an additional layer of specificity to vesicle docking. This mechanism likely involves the Sec1 homolog Sly1, which we identified in isolated docking complexes. We also report the identification of a novel v-SNARE (Ykt6p) component of the yeast ER-Golgi docking complex that has a CAAX box and is predicted to be lipid anchored. The surprising finding that docking complexes can contain many distinct species of SNAREs (Sed5p, Bos1p, Sec22p, Ykt6p, and likely Bet1p, p28, and p14) suggests that multimeric interactions are features of the fusion machinery, and may also improve the fidelity of vesicle targeting.

Amino Acid Sequence↗

The peptide repeat domain of nucleoporin Nup98 functions as a docking site in transport across the nuclear pore complex.

We report the cDNA deduced primary structure of a wheat germ agglutinin-reactive nuclear pore complex (NPC) protein of rat. The protein, termed Nup98 (for nucleoporin of 98 kDa), contains numerous GLFG and FG repeats and some FXFG repeats and is thus a vertebrate member of a family of GLFG nucleoporins that were previously discovered in yeast. Immunoelectron microscopy showed Nup98 to be asymmetrically located at the nucleoplasmic side of the NPC. Nup98 functions as one of several docking site nucleoporins in a cytosolic docking activity-mediated binding of a model transport substrate. The docking site of Nup98 was mapped to its N-terminal half, which contains all of the peptide repeats. A recombinant segment of this region depleted the docking activity of cytosol. We suggest that the peptide repeat domain of Nup98, together with peptide repeat domains of other nucleoporins, forms an array of sites for mediated docking of transport substrate, and that bidirectional transport across the NPC proceeds by repeated docking and undocking reactions.

Amino Acid Sequence↗

Docking simulation with a purine nucleoside specific homology model of deoxycytidine kinase, a target enzyme for anticancer and antiviral therapy.

5'-Phosphorylation, catalyzed by human deoxycytidine kinase (dCK), is a crucial step in the metabolic activation of anticancer and antiviral nucleoside antimetabolites, such as cytarabine (AraC), gemcitabine, cladribine (CdA), and lamivudine. Recently, crystal structures of dCK (dCKc) with various pyrimidine nucleosides as substrates have been reported. However, there is no crystal structure of dCK with a bound purine nucleoside, although purines are good substrates for dCK. We have developed a model of dCK (dCKm) specific for purine nucleosides based on the crystal structure of purine nucleoside bound deoxyguanosine kinase (dGKc) as the template. dCKm is essential for computer aided molecular design (CAMD) of novel anticancer and antiviral drugs that are based on purine nucleosides since these did not bind to dCKc in our docking experiments. The active site of dCKm was larger than that of dCKc and the amino acid (aa) residues of dCKm and dCKc, in particular Y86, Q97, D133, R104, R128, and E197, were not in identical positions. Comparative docking simulations of deoxycytidine (dC), cytidine (Cyd), AraC, CdA, deoxyadenosine (dA), and deoxyguanosine (dG) with dCKm and dCKc were carried out using the FlexX docking program. Only dC (pyrimidine nucleoside) docked into the active site of dCKc but not the purine nucleosides dG and dA. As expected, the active site of dCKm appeared to be more adapted to bind purine nucleosides than the pyrimidine nucleosides. While water molecules were essential for docking experiments using dCKc, the absence of water molecules in dCKm did not affect the ability to correctly dock various purine nucleosides.

Amino Acid Sequence↗