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 343 records · Page 19Linked to original sources
PubMed2002

Docking unbound proteins using shape complementarity, desolvation, and electrostatics.

A comprehensive docking study was performed on 27 distinct protein-protein complexes. For 13 test systems, docking was performed with the unbound X-ray structures of both the receptor and the ligand. For the remaining systems, the unbound X-ray structure of only molecule was available; therefore the bound structure for the other molecule was used. Our method optimizes desolvation, shape complementarity, and electrostatics using a Fast Fourier Transform algorithm. A global search in the rotational and translational space without any knowledge of the binding sites was performed for all proteins except nine antibodies recognizing antigens. For these antibodies, we docked their well-characterized binding site-the complementarity-determining region defined without information of the antigen-to the entire surface of the antigen. For 24 systems, we were able to find near-native ligand orientations (interface C(alpha) root mean square deviation less than 2.5 A from the crystal complex) among the top 2,000 choices. For three systems, our algorithm could identify the correct complex structure unambiguously. For 13 other complexes, we either ranked a near-native structure in the top 20 or obtained 20 or more near-native structures in the top 2,000 or both. The key feature of our algorithm is the use of target functions that are highly tolerant to conformational changes upon binding. If combined with a post-processing method, our algorithm may provide a general solution to the unbound docking problem. Our program, called ZDOCK, is freely available to academic users (http://zlab.bu.edu/~rong/dock/).

Algorithms↗
PubMed2002

Understanding protein structure-function relationships in Family 47 alpha-1,2-mannosidases through computational docking of ligands.

Family 47 alpha-1,2-mannosidases are crucial enzymes involved in N-glycan maturation in the endoplasmic reticula and Golgi apparati of eukaryotic cells. High-resolution crystal structures of the human and yeast endoplasmic reticulum alpha-1,2-mannosidases have been recently determined, the former complexed with the inhibitors 1-deoxymannojirimycin and kifunensine, both of which bind in its active site in the unusual 1C4 conformation. However, unambiguous identification of the catalytic proton donor and nucleophile involved in glycoside bond hydrolysis was not possible from this structural information. In this work, alpha-D-galactose, alpha-D-glucose, and alpha-D-mannose were computationally docked in the active site in the energetically stable 4C1 conformation as well as in the 1C4 conformation to compare their interaction energetics. From these docked structures, a model for substrate and conformer selectivity based on the dimensions of the active site was proposed. Alpha-D-galactopyranosyl-(1-->2)-alpha-D-mannopyranose, alpha-D-glucopyranosyl-(1-->2)-alpha-D-mannopyranose, and alpha-D-mannopyranosyl-(1-->2)-alpha-D-mannopyranose were also docked into the active site with their nonreducing-end residues in the 1C4 and E4 (representing the transition state) conformations. Based on the docked structure of alpha-D-mannopyranosyl-E4-(1-->2)-alpha-D-mannopyranose, the catalytic acid and base are Glu132 and Glu435, respectively.

1-Deoxynojirimycin↗
PubMed2003

Pharmacophore-based molecular docking to account for ligand flexibility.

Rapid computational mining of large 3D molecular databases is central to generating new drug leads. Accurate virtual screening of large 3D molecular databases requires consideration of the conformational flexibility of the ligand molecules. Ligand flexibility can be included without prohibitively increasing the search time by docking ensembles of precomputed conformers from a conformationally expanded database. A pharmacophore-based docking method whereby conformers of the same or different molecules are overlaid by their largest 3D pharmacophore and simultaneously docked by partial matches to that pharmacophore is presented. The method is implemented in DOCK 4.0.

Algorithms↗
PubMed2003

ICM-DISCO docking by global energy optimization with fully flexible side-chains.

The ICM-DISCO (Docking and Interface Side-Chain Optimization) protein-protein-docking method is a direct stochastic global energy optimization from multiple starting positions of the ligand. The first step is performed by docking of a rigid all-atom ligand molecule to a set of soft receptor potentials precalculated on a 0.5 A grid from realistic solvent-corrected force-field energies. This step finds the correct solution as the lowest energy conformation in almost 100% of the cases in which interfaces do not change on binding. The second step is needed to deal with the induced changes and includes the global optimization of the interface side-chains of up to 400 best solutions. The CAPRI predictions were performed fully automatically with this method. Available experimental information was included as a filtering step to favor expected docking surfaces. In three of the seven proposed targets, the ICM-DISCO method found a good solution (>50% of correct contacts) within the five submitted models. The procedure is global and fully automated. We demonstrate that the algorithm handles the induced changes of surface side-chains but is less successful if the backbone undergoes large-scale rearrangements.

Algorithms↗
PubMed2003

Weighted geometric docking: incorporating external information in the rotation-translation scan.

Weighted geometric docking is a prediction algorithm that matches weighted molecular surfaces. Each molecule is represented by a grid of complex numbers, storing information about the shape of the molecule in the real part and weight information in the imaginary part. The weights are based on experimental biochemical and biophysical data or on theoretical analyses of amino acid conservation or correlation patterns in multiple-sequence alignments of homologous proteins. Only a few surface residues on either one or both molecules are weighted. In contrast to methods that use postscan filtering based on biochemical information, our method incorporates the external data in the rotation-translation search, producing a different set of docking solutions biased toward solutions in which the up-weighted residues are at the interface. Similarly, interactions involving specified residues can be impeded. The weighted geometric algorithm was applied to five systems for which regular geometric docking of the unbound molecules gave poor results. We obtained much better ranking of the nearly correct prediction and higher statistical significance when weighted geometric docking was used. The method was successful even when the weighted portion of the surface corresponded only partially and approximately to the binding site.

Algorithms↗
PubMed2003

The tissue factor/factor VIIa/factor Xa complex: a model built by docking and site-directed mutagenesis.

Factor X is activated to factor Xa (fXa) in the extrinsic coagulation pathway by the tissue factor (TF)/factor VIIa (fVIIa) complex. Upon activation, the fXa molecule remains associated with the TF/fVIIa complex, and this ternary complex is known to activate protease-activated receptors (PARs) 1 and 2. Activation of fVII in the TF complex by fXa is also seen at physiologic concentrations. The ternary complexes TF/fVII/fXa, TF/fVIIa/fX, and TF/fVIIa/fXa are therefore all physiologically relevant and of interest as targets for inhibition of both coagulation and cell-signaling pathways that are important in cardiovascular disease and inflammation. We therefore present a model of the TF/fVIIa/fXa complex, built with the use of the available structures of the TF/fVIIa complex and fXa by protein-protein docking calculations with the program Surfdock. The fXa model has an extended conformation, similar to that of fVIIa in the TF/fVIIa complex, with extensive interactions with TF and the protease domain of fVIIa. All four domains of fXa are involved in the interaction. The gamma-carboxyglutamate (Gla) and epithelial growth factor (EGF1 and EGF2) domains of fVIIa are not significantly involved in the interaction. Docking of the Gla domain of fXa to TF/fVIIa has been reported previously. The docking results identify potential interface residues, allowing rational selection of target residues for site-directed mutagenesis. This combination of docking and mutagenesis confirms that residues Glu51 and Asn57 in the EGF1 domain, Asp92 and Asp95 in the EGF2 domain, and Asp 185a, Lys 186, and Lys134 in the protease domain of factor Xa are involved in the interaction with TF/fVIIa. Other fX protease domain residues predicted to be involved in the interaction come from the 160s loop and the N-terminus of the fX protease domain, which is oriented in such a way that activation of both fVII by fXa, and the reciprocal fX activation by fVIIa, is possible.

Binding Sites↗
PubMed2004

Type I collagen N-telopeptides adopt an ordered structure when docked to their helix receptor during fibrillogenesis.

The in vitro rate and specificity of fibrillogenesis in type I collagen depends on the integrity of the amino (N)-telopeptide domain. In vivo an intact N-telopeptide domain is also required for normal fibril assembly. Although Chou-Fasman predictions and NMR studies suggested that a type I beta-turn could be induced in alpha1(I) N-telopeptide chains, computer modeling did not identify ordered structures. Nevertheless, X-ray analysis and electron tomography studies have shown that the N-telopeptide is in one of the most highly ordered fibril domains. This study was undertaken to determine if the docking of the N-telopeptide to its helix receptor domain could induce the telopeptides to take up a specific conformation. With use of molecular modeling suite of programs, a (Gly-Pro-Pro)(n) triple-helical structure was built on the basis of high-resolution X-ray crystallographic coordinates and then replaced with the actual bovine collagen residues 924-938, the triple-helical alpha1(I)-N-telopeptide-receptor sequences. Energy minimization produced a modified triple-helical conformation. The bovine alpha1(I) N-telopeptide sequence was similarly minimized and docked to this receptor. The docking induced an ordered conformation with a stabilizing hydrogen bond in the N-telopeptide and, importantly, a reciprocal reordering of the triple-helical conformation in the binding domain. This docked structure placed Lys residues in both telopeptide and helix in the correct locations for cross-link formation. The modeling has been extended to the three-chain N-telopeptide domain and finally to the construction of the Hulmes-Miller quasi-hexagonal packing structure. Each N-telopeptide domain can form linkages with two adjacent, aligned helix receptor domains. The telopeptides and the order of staggering of the three chains in the helix play crucial roles in the packing and intrafibrillar cross-linking patterns and the relative azimuthal orientations of adjacent molecules in the fibril. The models confirm the high order in the N-telopeptide 4D overlap zone.

Amino Acid Sequence↗
PubMed2001

Model for the three-dimensional structure of vitronectin: predictions for the multi-domain protein from threading and docking.

The structure of vitronectin, an adhesive protein that circulates in high concentrations in human plasma, was predicted through a combination of computational methods and experimental approaches. Fold recognition and sequence-structure alignment were performed using the threading program PROSPECT for each of three structural domains, i.e., the N-terminal somatomedin B domain (residues 1-53), the central region that folds into a four-bladed beta-propeller domain (residues 131-342), and the C-terminal heparin-binding domain (residues 347-459). The atomic structure of each domain was generated using MODELLER, based on the alignment obtained from threading. Docking experiments between the central and C-terminal domains were conducted using the program GRAMM, with limits on the degrees of freedom from a known inter-domain disulfide bridge. The docked structure has a large inter-domain contact surface and defines a putative heparin-binding groove at the inter-domain interface. We also docked heparin together with the combined structure of the central and C-terminal domains, using GRAMM. The predictions from the threading and docking experiments are consistent with experimental data on purified plasma vitronectin pertaining to protease sensitivity, ligand-binding sites, and buried cysteines.

Amino Acid Sequence↗
PubMed2004

Hydrophobic complementarity in protein-protein docking.

Formation of hydrophobic contacts across a newly formed interface is energetically favorable. Based on this observation we developed a geometric-hydrophobic docking algorithm that estimates quantitatively the hydrophobic complementarity at protein-protein interfaces. Each molecule to be docked is represented as a grid of complex numbers, storing information regarding the shape of the molecule in the real part and information regarding the hydropathy of the surface in the imaginary part. The grid representations are correlated using fast Fourier transformations. The algorithm is used to compare the extent of hydrophobic complementarity in oligomers (represented by D2 tetramers) and in hetero-dimers of soluble proteins (complexes). We also test the implication of hydrophobic complementarity in distinguishing correct from false docking solutions. We find that hydrophobic complementarity at the interface exists in oligomers and in complexes, and in both groups the extent of such complementarity depends on the size of the interface. Thus, the non-polar portions of large interfaces are more often juxtaposed than non-polar portions of small interfaces. Next we find that hydrophobic complementarity helps to point out correct docking solutions. In oligomers it significantly improves the ranks of nearly correct reassembled and modeled tetramers. Combining geometric, electrostatic and hydrophobic complementarity for complexes gives excellent results, ranking a nearly correct solution < 10 for 5 of 23 tested systems, < 100 for 8 systems and < 1000 for 19 systems.

Algorithms↗
PubMed2005

The effect of resolution-dependent global shape modifications on rigid-body protein-protein docking.

Docking unbound molecules presents a challenge in the case where no prior biological or bioinformatic knowledge exists. This is mainly due to differences between the structures of the molecules when in a complex and in the free state. Presumably, these differences interfere with the ability of protein-protein docking algorithms, which rely on a dominant shape descriptor, to identify the correct solution and rank it higher than false solutions. In this study we verify the notion that small discords in the molecular fit can be eliminated by using appropriately designed low-resolution shape descriptors, thereby improving the docking results. We exploit the inherent gradual resolution dependency of Fourier transforms and formulate a resolution-dependent shape descriptor by truncating selected Fourier transform terms. Thus, different levels of shape modification are attained, affecting the degree of detail in the depiction of the molecular surface. We applied the modified descriptor to a selection of 23 protein-protein systems, using the unbound structures where possible. The docking results obtained with the new geometric descriptor were considerably superior to former results, improving the ranks of nearly correct solutions for 17 systems. Unification of the results of scans in which different resolutions were employed further improved the ranks of nearly correct solutions to less than 100 for 12 of the 23 systems and less than 300 for 20 systems. The new geometric descriptor can be combined with other descriptors, which typify the electrostatic or hydrophobic character of the molecular surface, and with external experimental or bioinformatic data.

Binding Sites↗
PubMed2005

Assessment of CAPRI predictions in rounds 3-5 shows progress in docking procedures.

The current status of docking procedures for predicting protein-protein interactions starting from their three-dimensional (3D) structure is reassessed by evaluating blind predictions, performed during 2003-2004 as part of Rounds 3-5 of the community-wide experiment on Critical Assessment of PRedicted Interactions (CAPRI). Ten newly determined structures of protein-protein complexes were used as targets for these rounds. They comprised 2 enzyme-inhibitor complexes, 2 antigen-antibody complexes, 2 complexes involved in cellular signaling, 2 homo-oligomers, and a complex between 2 components of the bacterial cellulosome. For most targets, the predictors were given the experimental structures of 1 unbound and 1 bound component, with the latter in a random orientation. For some, the structure of the free component was derived from that of a related protein, requiring the use of homology modeling. In some of the targets, significant differences in conformation were displayed between the bound and unbound components, representing a major challenge for the docking procedures. For 1 target, predictions could not go to completion. In total, 1866 predictions submitted by 30 groups were evaluated. Over one-third of these groups applied completely novel docking algorithms and scoring functions, with several of them specifically addressing the challenge of dealing with side-chain and backbone flexibility. The quality of the predicted interactions was evaluated by comparison to the experimental structures of the targets, made available for the evaluation, using the well-agreed-upon criteria used previously. Twenty-four groups, which for the first time included an automatic Web server, produced predictions ranking from acceptable to highly accurate for all targets, including those where the structures of the bound and unbound forms differed substantially. These results and a brief survey of the methods used by participants of CAPRI Rounds 3-5 suggest that genuine progress in the performance of docking methods is being achieved, with CAPRI acting as the catalyst.

Algorithms↗
PubMed2005

Geometry-based flexible and symmetric protein docking.

We present a set of geometric docking algorithms for rigid, flexible, and cyclic symmetry docking. The algorithms are highly efficient and have demonstrated very good performance in CAPRI Rounds 3-5. The flexible docking algorithm, FlexDock, is unique in its ability to handle any number of hinges in the flexible molecule, without degradation in run-time performance, as compared to rigid docking. The algorithm for reconstruction of cyclically symmetric complexes successfully assembles multimolecular complexes satisfying C(n) symmetry for any n in a matter of minutes on a desktop PC. Most of the algorithms presented here are available at the Tel Aviv University Structural Bioinformatics Web server (http://bioinfo3d.cs.tau.ac.il/).

Algorithms↗
PubMed2005

Data-driven docking: HADDOCK's adventures in CAPRI.

We have shown previously that given high-resolution structures of the unbound molecules, structure determination of protein complexes is possible by including biochemical and/or biophysical data as highly ambiguous distance restraints in a docking approach. We applied this method, implemented in the HADDOCK (High Ambiguity Driven DOCKing) package (Dominguez et al., J Am Chem Soc 2003;125:1731-1737), to the targets in the fourth and fifth rounds of CAPRI. Here we describe our results and analyze them in detail. Special attention is given to the role of flexibility in our docking method and the way in which this improves the docking results. We describe extensions to our approach that were developed as a direct result of our participation in CAPRI. In addition to experimental information, we also included interface residue predictions from PPISP (Protein-Protein Interaction Site Predictor; Zhou and Shan, Proteins 2001;44:336-343), a neural network method. Using HADDOCK we were able to generate acceptable structures for 6 of the 8 targets, and to submit at least 1 acceptable structure for 5 of them. Of these 5 submissions, 3 were of medium quality (Targets 10, 11, and 15) and 2 of high quality (Targets 13 and 14). In all cases, predictions were obtained containing at least 40% of the correct epitope at the interface for both ligand and receptor simultaneously.

Algorithms↗
PubMed2005

ATTRACT: protein-protein docking in CAPRI using a reduced protein model.

Protein-protein complex structures have been predicted for CAPRI Rounds 3 and 5 using a reduced protein model. Proteins are represented by up to 3 pseudoatoms per amino acid. The docking approach termed ATTRACT is based on energy minimization in translational and rotational degrees of freedom of one protein with respect to another protein. The reduced protein model allows one to perform systematic docking minimization of many thousand start structures in reasonable computer time. Flexibility of critical surface side-chains can be accounted for by a multiple conformational copy approach. The multicopy approach allows simultaneous adjustment of side-chain conformations and optimization of translational and rotational degrees of freedom of one protein with respect to the partner during docking. For 3 (Targets 8, 14, and 19) out of 5 CAPRI targets, the approach resulted in predictions in close agreement with experiment [root-mean-square deviation (RMSD) of backbone atoms within 10 A of the protein-protein interface < 1.8 A]. The comparison of predicted and experimental structures of the CAPRI targets indicates that besides local conformational changes (e.g., changes in side-chain conformations), global conformational changes of the protein backbone can be critical for complex formation. These conformational changes not accounted for during docking are a likely reason for the unrealistic predictions in 2 cases (Targets 9 and 18).

Algorithms↗
PubMed2005

Docking essential dynamics eigenstructures.

This article describes our attempts to dock the targets in CAPRI Rounds 3-5 using Hex 4.2, and it introduces a novel essential dynamics approach to generate multiple feasible conformations for docking. In the blind trial, the basic Hex algorithm found 1 high-accuracy solution for CAPRI Target 12, and several further medium- and low-accuracy solutions for Targets 11, 12, 13, and 14. Subsequent a posteriori docking of the targets using essential dynamics "eigenstructures" was found to give consistently better predictions than rigidly docking only the unbound or model-built starting structures. Some suggestions to improve this promising new approach are presented.

Algorithms↗
PubMed2006

Scoring a diverse set of high-quality docked conformations: a metascore based on electrostatic and desolvation interactions.

Predicting protein-protein interactions involves sampling and scoring docked conformations. Barring some large structural rearrangement, rapidly sampling the space of docked conformations is now a real possibility, and the limiting step for the successful prediction of protein interactions is the scoring function used to reduce the space of conformations from billions to a few, and eventually one high affinity complex. An atomic level free-energy scoring function that estimates in units of kcal/mol both electrostatic and desolvation interactions (plus van der Waals if appropriate) of protein-protein docked conformations is used to rerank the blind predictions (860 in total) submitted for six targets to the community-wide Critical Assessment of PRediction of Interactions (CAPRI; http://capri.ebi.ac.uk). We found that native-like models often have varying intermolecular contacts and atom clashes, making unlikely that one can construct a universal function that would rank all these models as native-like. Nevertheless, our scoring function is able to consistently identify the native-like complexes as those with the lowest free energy for the individual models of 16 (out of 17) human predictors for five of the targets, while at the same time the modelers failed to do so in more than half of the cases. The scoring of high-quality models developed by a wide variety of methods and force fields confirms that electrostatic and desolvation forces are the dominant interactions determining the bound structure. The CAPRI experiment has shown that modelers can predict valuable models of protein-protein complexes, and improvements in scoring functions should soon solve the docking problem for complexes whose backbones do not change much upon binding. A scoring server and programs are available at http://structure.pitt.edu.

Models, Biological↗
PubMed2006

Docking of human interleukin-15 to its specific receptor alpha chain: correlation between molecular modeling and mutagenesis experimental data.

A structural model of the sushi domain of IL-15Ralpha was first obtained by homology modeling to study its interactions with IL-15 by means of molecular modeling, peptide scanning, and site-directed mutagenesis. From these experimental data, a putative interacting surface of IL-15Ralpha with a previously published IL-15 model was inferred: Leu25, Leu44, and Glu46 of IL-15 and Arg35 of IL-15Ralpha were found to be key interfacial residues and were subsequently used as filters for the construction of docking solutions. Human IL-15/IL-15Ralpha complexes were constructed in two stages, with a preliminary docking procedure, treating the two partners as rigid bodies and using these filters. In this first stage, two classes of docking solutions were characterized. From a topological point of view, each solution could be derived from the other by reverse orientation of one partner in relation to the other. In a second stage, several further energy refinements clearly favored one solution. Moreover, this unique docking solution was confirmed by molecular modeling of IL-15 mutants previously built and tested in our laboratory. Finally, this complex model, which is a useful tool to study the IL-15/IL-15Ralpha interface, was topologically compared to IL-2/IL-2Ralpha complexes (previous model in the literature and recent crystal structure).

Amino Acid Sequence↗
PubMed1996

A fast flexible docking method using an incremental construction algorithm.

We present an automatic method for docking organic ligands into protein binding sites. The method can be used in the design process of specific protein ligands. It combines an appropriate model of the physico-chemical properties of the docked molecules with efficient methods for sampling the conformational space of the ligand. If the ligand is flexible, it can adopt a large variety of different conformations. Each such minimum in conformational space presents a potential candidate for the conformation of the ligand in the complexed state. Our docking method samples the conformation space of the ligand on the basis of a discrete model and uses a tree-search technique for placing the ligand incrementally into the active site. For placing the first fragment of the ligand into the protein, we use hashing techniques adapted from computer vision. The incremental construction algorithm is based on a greedy strategy combined with efficient methods for overlap detection and for the search of new interactions. We present results on 19 complexes of which the binding geometry has been crystallographically determined. All considered ligands are docked in at most three minutes on a current workstation. The experimentally observed binding mode of the ligand is reproduced with 0.5 to 1.2 A rms deviation. It is almost always found among the highest-ranking conformations computed.

Algorithms↗