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 721 records · Page 40Linked to original sources

Docking prediction using biological information, ZDOCK sampling technique, and clustering guided by the DFIRE statistical energy function.

We entered the CAPRI experiment during the middle of Round 4 and have submitted predictions for all 6 targets released since then. We used the following procedures for docking prediction: (1) the identification of possible binding region(s) of a target based on known biological information, (2) rigid-body sampling around the binding region(s) by using the docking program ZDOCK, (3) ranking of the sampled complex conformations by employing the DFIRE-based statistical energy function, (4) clustering based on pairwise root-mean-square distance and the DFIRE energy, and (5) manual inspection and relaxation of the side-chain conformations of the top-ranked structures by geometric constraint. Reasonable predictions were made for 4 of the 6 targets. The best fraction of native contacts within the top 10 models are 89.1% for Target 12, 54.3% for Target 13, 29.3% for Target 14, and 94.1% for Target 18. The origin of successes and failures is discussed. .

Algorithms↗

Structural basis for SUMO-E2 interaction revealed by a complex model using docking approach in combination with NMR data.

The interaction between small ubiquitin-related modifier SUMO and its conjugating-enzyme Ubc9 (E2) is an essential step in SUMO conjugation cascade. However, an experimental structure of such a transient complex is still unavailable. Here, a structural model of SUMO-3-Ubc9 complex was obtained with HADDOCK, combining NMR chemical shift mapping information. Docking calculations were performed using SUMO-3 and Ubc9 structures as input. The resulting complex reveals that the complementary surface electrostatic potentials contribute dominantly to the specific interaction. At the interface, similar numbers of oppositely-charged conserved residues are identified on the respective binding partners. Hydrogen bonds are formed in the vicinity of the interface to stabilize the complex. Comparison of the structure of SUMO-3-Ubc9 complex generated by HADDOCK and the experimental structures in free form indicates that SUMO-3 and Ubc9 maintain their respective fold as a whole after docking. However, the N-terminal helix alpha1 and its subsequent L1 loop of Ubc9 experience sizeable changes upon complex formation. They cooperatively move towards the hydrophilic side of the beta-sheet of SUMO-3. Our observations are consistent with the data from previous Ubc9 mutational analysis and conformational flexibility studies. Together, we have proposed that the SUMO-3-Ubc9 interaction is strongly electrostatically driven and the N terminus of Ubc9 shifts to SUMO-3 to facilitate the interaction. The NMR-based structural model, which provides considerable insights into the molecular basis of the specific SUMO-E2 recognition and interaction, implicates the general interaction mode between SUMO-3 and Ubc9 homologues from yeast to humans.

Binding Sites↗

Interaction of (-)-epigallocatechin-3-gallate with human serum albumin: fluorescence, fourier transform infrared, circular dichroism, and docking studies.

(-)-Epigallocatechin-3-gallate (EGCG), the major constituent of green tea has been reported to prevent many diseases by virtue of its antioxidant properties. The binding of EGCG with human serum albumin (HSA) has been investigated for the first time by using fluorescence, circular dichroism (CD), Fourier transform infrared (FTIR) spectroscopy, and protein-ligand docking. We observed a quenching of fluorescence of HSA in the presence of EGCG. The binding parameters were determined by a Scatchard plot and the results were found to be consistent with those obtained from a modified Stern-Volmer equation. From the thermodynamic parameters calculated according to the van't Hoff equation, the enthalpy change deltaH degrees and entropy change deltaS degrees were found to be -22.59 and 16.23 J/mol K, respectively. These values suggest that apart from an initial hydrophobic association, the complex is held together by van der Waals interactions and hydrogen bonding. Data obtained by fluorescence spectroscopy, CD, and FTIR experiments along with the docking studies suggest that EGCG binds to residues located in subdomains IIa and IIIa of HSA. Specific interactions are observed with residues Trp 214, Arg 218, Gln 221, Asn 295 and Asp 451. We have also looked at changes in the accessible surface area of the interacting residues on binding EGCG for a better understanding of the interaction.

Catechin↗

Protein-protein recognition analyzed by docking simulation.

Antibody-lysozyme and protease-inhibitor complexes are reconstituted by docking lysozyme as a rigid body onto the combining site of the antibodies and the inhibitors onto the active site of the proteases. Simplified protein models with one sphere per residue are subjected to simulated annealing using a crude energy function where the attractive component is proportional to the interface area. The procedure finds clusters of orientations in which a steric fit between the two protein components is achieved over a large contact surface. With five out of six complexes, the native structure of the complexes determined by X-ray crystallography is among those retained. Docked complexes are then subjected to conformational energy refinement with full atomic detail. With Fab HyHEL 5 and lysozyme, a native-like complex has the lowest refined energy. It can also be retrieved when starting with the X-ray structure of free lysozyme. However, some non-native complexes cannot be rejected: they form large interfaces, have a large number of H-bonds, and few unpaired polar groups. While these are necessary features of protein-protein recognition, they are not sufficient in determining specificity.

Algorithms↗

Orientational sampling and rigid-body minimization in molecular docking.

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

Carboxypeptidases↗

Rigid-body docking with mutant constraints of influenza hemagglutinin with antibody HC19.

An automatic docking algorithm has been applied to the modeling of the complex between hemagglutinin from influenza virus and the Fab fragment of a monoclonal antibody raised against this antigen. We have introduced here the use of biochemical information provided by mutants of hemagglutinin. The docking procedure finds a small number of candidate solutions where three sites of escape mutations are buried and form hydrogen bonds in the interface. The localization of the epitope is improved by additional biochemical data about mutants that do not affect antibody binding. Five candidate solutions with low energy, reasonably well-packed interfaces, and six to ten hydrogen bonds are compatible with mutant information. One of the five stands out as generally better than the others from these points of views.

Algorithms↗

Docking of a human rhinovirus neutralizing antibody onto the viral capsid.

The structure of the complex between the Fab fragment of a human rhinovirus serotype 2 (HRV2) neutralizing antibody (8F5) and a cross-reactive synthetic peptide derived from the viral capsid protein VP2 has been recently determined by crystallographic methods. The conformation adopted by the peptide was very similar to and could be superimposed onto the corresponding region of the viral protein VP2 of human rhinovirus 1A (HRV1A) whose three-dimensional structure is known. The structure of the Fab fragment determined in the complex was docked onto the viral capsid using the superimposition transformation found for the peptide. In the resulting model the Fab protrudes almost radially to about 60 A from the surface of the virion without any major steric problem. The Fab fragment was then placed on each one of the 60 equivalent epitopes using the T = 1 icosahedral symmetry of the virus. The closest pairs of Fab fragments are related by viral 2-fold axes and run almost parallel to each other without clashing. These axes of symmetry from the viral particle could thus be coincident with the dyad axes of the antibodies. Furthermore, comparison of the three-dimensional structure of the Fab/peptide complex with the structure of the Fab fragment alone indicates that the flexibility of the antibody's elbow would facilitate bivalent attachment to the same viral particle. In accordance with the docking results, experimental determination of the stoichiometry of binding yielded a ratio of 30 IgG molecules per virion also suggesting bivalent attachment of antibody 8F5 onto the viral particle. The neutralization of viral infectivity, being neither aggregation (this paper) nor inhibition of receptor binding, might be mainly achieved by reducing viral spread from cell to cell and/or inhibition of uncoating.

Amino Acid Sequence↗

Structure-function relations of antithrombin III-heparin interactions as assessed by biophysical and biological assays and molecular modeling of peptide-pentasaccharide-docked complexes.

The serine proteinase inhibitor antithrombin III (ATIII) is a key regulatory protein of intrinsic blood coagulation. ATIII attains its full biological activity only upon binding polysulfated glycosaminoglycans, such as heparin. Peptide K121-A134, based on the sequence of ATIII in the D helix region, was previously shown by us (Tyler-Cross et al., Protein Sci. 3, 620-627, 1994) to encompass part (or all) of the purported high-affinity heparin binding region of ATIII. A series of peptide analogs has now been prepared whose sequences are identical to K121-A134 except that single cationic residues of K121-A134 have been successively replaced with Ala. In one case, the Arg residue of the reference peptide corresponding to R129 of ATIII has been replaced by Gln (R129deltaQ peptide), thus mimicking the naturally occurring mutant protein, ATIII Geneva. The heparin affinity of all peptides was quantitated by isothermal titration calorimetry and by peptide/ATIII competition binding assays. Replacement of any single cationic residue with Ala had a deleterious effect on heparin binding. The greatest reduction in heparin affinity (more than 30-fold) was observed with the R129deltaQ peptide (KD = 1.5 +/- 0.06 microM vs 51 +/- 2 nM for the reference peptide, K121-A134). Furthermore, each of the Ala-replacement peptides was a less-effective inhibitor of ATIII-heparin complex formation than the reference peptide. The poorest inhibitor was the R129deltaQ peptide which showed nearly 30% decrease in inhibition potency (60% inhibition at 100 microM peptide vs 90% inhibition with the reference peptide). The relative heparin affinities of the peptides measured by biological assay were the same as determined by titration calorimetry. Consequently, we modeled the complexes formed between the pentasaccharide unit structure and the R129deltaQ peptide or the reference peptide, K121-A134. In the "docked" complex, the assumed conformation of K121-A134 permitted juxtaposition of the cationic residues of the peptide with functional anionic groups of the pentasaccharide known to be involved in binding. A docked complex could also be formed between the R129deltaQ peptide and the pentasaccharide, but misalignment of critical peptide and saccharide functional groups was observed. The structure of the R129deltaQ-pentasaccharide complex was highly irregular because F123 and Y131 were completely surface exposed, likely yielding an unfavorable structure in aqueous solution. The observations from molecular modeling allow us to suggest that ATIII Geneva displays decreased heparin binding affinity due to its inability to form a productive binding complex in which essential electrostatic contacts are made between suitably juxtaposed saccharide anionic functional groups and cationic amino acid side chains.

Amino Acid Sequence↗

The unitary evoked potential at the frog nerve-muscle junction results from synchronous gating of fusion pores at docked vesicles.

Exocytosis of a single vesicle has been proposed as the mechanism which determines quantal size by releasing a prepackaged and standard amount of acetylcholine. As first described by del Castillo and Katz (1954) the endplate potential is composed of 100 unitary events and the small variance suggests a binomial release from 100 "discrete patches of membrane". However, exocytosis of 100 vesicles selected randomly from 5000 docked vesicles would yield a variance that is 7 times greater than observed values. We propose that the presynaptic ridge with its compliment of docked vesicles functions as the "discrete patch of membrane" such that arrays of calcium activated fusion pores meter transmitter to form the unit of release. A model based on the synchronous flicker of a large number of fusion pores produces the small variance of both miniature end plate potentials and unitary end plate potentials. Release from a single locus (fusion pore) would generate the sub-MEPP. This model permits vesicle trafficking and vesicular content depletion during tetanic stimulation and explains the frequency dependency of MEPP amplitudes and changes in sub-MEPP to bell-MEPP class ratios.

Animals↗

Sequential docking, molecular differentiation, and positioning of T-Tubule/SR junctions in developing mouse skeletal muscle.

Skeletal muscle Ca(2+) release units (CRUs) are junctions of the surface membrane/T-tubule system and the sarcoplasmic reticulum (SR) that function in excitation-contraction coupling. They contain high concentrations of dihydropyridine receptors (DHPRs) in the T-tubules and of ryanodine receptors (RyR) in the SR and they are positioned at specific locations in the sarcomere. In order to characterize the sequence of developmental steps leading to the specific molecular and structural organization of CRUs, we applied a range of imaging techniques that allowed us to follow the differentiation of the membrane compartments and the expression of junctional proteins in developing mouse diaphragm muscle. We find that docking of the two membrane systems precedes the incorporation of the RyRs into the junctions, and that T-tubule/SR junctions are formed and positioned at the I-A interface at a stage when the orientation of T-tubule is predominantly longitudinal. Thus, the sequence of developmental events is first the docking of T-tubules and SR, secondly the incorporation of RyR in the junctions, thirdly the positioning of the junctions in the sarcomere, and only much later the transverse orientation of the T-tubules. These sequential stages suggests an order of inductive processes for the molecular differentiation and structural organization of the CRUs in skeletal muscle development.

Animals↗

PUZZLE: a new method for automated protein docking based on surface shape complementarity.

We describe here a novel procedure for automated protein docking, based only on geometric criteria. In our algorithm we project protein surfaces into bi-dimensional matrices; the search for complementary regions is performed by detecting matching sub-matrices. An exhaustive sampling of the rotation space is made in order to analyse all the possible relative orientations of the two proteins, but nevertheless this procedure requires a relatively short processing time (3 h to 24 h cpu time on a SG4D320, depending on the complexity of the input information). When tested with co-crystallized, free components and models of components of known protein-protein complexes, the method gave very satisfactory results. The procedure selects no more than four relative orientations of the molecular components, but the correct orientation is always present among them, ranking either first or second. In more than half the cases the "wrong" solutions nevertheless correctly identify most of the residues involved in the interaction. This is remarkable also in view of the fact that the chosen test complexes (trypsin-trypsin inhibitor and antibody-lysozyme) have a very different geometry of surface complementarity: trypsin inhibitor inserts a long side-chain into the deep specificity pocket of the protease, while the interface between antibody and lysozyme is rather flat and contains buried water molecules (not included in the calculation). In order to simulate a more realistic protein docking problem, we also used a trypsin inhibitor and an anti-lysozyme antibody model in our simulations, again with satisfying results.

Algorithms↗

Hydrogen bonding and molecular surface shape complementarity as a basis for protein docking.

A geometric docking algorithm based upon correlation analysis for quantification of geometric complementarity between protein molecular surfaces in close interfacial contact has been developed by a detailed optimization of the conformational search of the algorithm. In order to reduce the entire conformation space search required by the method a physico-chemical pre-filter of conformation space has been developed based upon the a priori assumption that two or more intermolecular hydrogen bonds are intrinsic to the mechanism of binding within protein complexes. Donor sites are defined spatially and directionally by the positions of explicitly calculated donor hydrogen atoms, and the vector space within a defined range about the donor atom-hydrogen atom bond vector. Acceptor sites are represented spatially and directionally by the van der Waals molecular surface points having normal vectors within a predefined range of vector space about the acceptor atom covalent bond vector(s). Geometric conditions necessary for the simultaneous hydrogen bonding interaction between both sites of functionally congruent hydrogen bonding site pairs, located on the individual proteins, are then tested on the basis of a transformation invariant parameterization of the site pair spatial and directional properties. Sterically acceptable conformations defined by interaction of functionally, spatially, and directionally compatible site pairs are then refined to a maximum contact of complementary contact surfaces using the simplex method for the angular search and correlation techniques for the translational search. The utility of the spatial and directional properties of hydrogen bonding donor and acceptor sites for the identification of candidate docking conformations is demonstrated by the reliable preliminary reduction of conformation space, the improved geometric ranking of the minimum RMS conformations of some complexes and the overall reduction of CPU time obtained.

Algorithms↗

An example of a protein ligand found by database mining: description of the docking method and its verification by a 2.3 A X-ray structure of a thrombin-ligand complex.

A computer program (SANDOCK) has been developed for the automated docking of small ligands to a target protein. It uses a guided matching algorithm to fit ligand atoms into the protein binding pocket. The protein is described by a modified Lee-Richard's dotted surface with each dot coded by chemical property and accessibility. Orientations of the ligand in the active site are generated such that a chemical and a shape complementary between the ligand and the active site cavity have to be fulfilled. The generated fits are evaluated with scoring functions which account for van der Waals, hydrophobic and hydrogen bonding interactions. This newly developed docking program can efficiently screen very large databases in a reasonable time and has been used to successfully identify novel ligands. The X-ray structure of a thrombin-ligand complex predicted by SANDOCK is described. The ligand binds to thrombin with a Kd of 65 microM and has an rmsd of 0.7 A for all ligand atoms from the predicted binding mode by SANDOCK.

Algorithms↗

The discovery of steroids and other novel FKBP inhibitors using a molecular docking program.

The molecular docking computer program SANDOCK was used to screen small molecule three-dimensional databases in the hunt for novel FKBP inhibitors. Spectroscopic measurements confirmed binding of over 20 compounds to the target protein, some with dissociation constants in the low micromolar range. The discovery that FK506 binding protein is a steroid binding protein may be of wider biological significance. Two-dimensional NMR was used to determine the steroid binding mode and confirmed the interactions predicted by the docking program.

Binding Sites↗

Guiding a docking mode by phage display: selection of correlated mutations at the staphylokinase-plasmin interface.

During co-evolution of interacting proteins, functionally disruptive mutations on one side of the interface may be compensated by local amino acid changes on the other to restore binding affinity. This information can be useful for geometry-based docking approaches by reducing the translational and rotational space available to the proteins. Here, we demonstrate that correlated mutations at a protein-protein interface can be rapidly identified by selecting a phage-displayed library of a randomly mutated component of the complex for complementation of mutations that decreased binding in the interacting partner. This approach was used to deduce the binding mode of staphylokinase (Sak), a 15.5 kDa "indirect" plasminogen activator on microplasmin (microPli), the 28 kDa serine protease domain of plasmin. Biopanning indicated that residues Arg94 and Gly174 in microPli are located in close proximity to Glu75 and the Glu88:Ile128 pair in Sak, respectively. The coupled mutations Glu94<-->Lys75 reversed and Gly174<-->Lys88:Val128 introduced a salt bridge, whereby the binding affinities (with coupling energies of 1.8 to 2.3 kcal mol-1, respectively) and the plasminogen activation ability of the mutated complexes were partially restored. These findings suggested a unique docking mode of Sak at the western rim of the active-site cleft of microPli, that is in agreement with the structure of the Sak-microPli complex as recently derived by other methods.

Amino Acid Substitution↗

Automated structure-based prediction of functional sites in proteins: applications to assessing the validity of inheriting protein function from homology in genome annotation and to protein docking.

A major problem in genome annotation is whether it is valid to transfer the function from a characterised protein to a homologue of unknown activity. Here, we show that one can employ a strategy that uses a structure-based prediction of protein functional sites to assess the reliability of functional inheritance. We have automated and benchmarked a method based on the evolutionary trace approach. Using a multiple sequence alignment, we identified invariant polar residues, which were then mapped onto the protein structure. Spatial clusters of these invariant residues formed the predicted functional site. For 68 of 86 proteins examined, the method yielded information about the observed functional site. This algorithm for functional site prediction was then used to assess the validity of transferring the function between homologues. This procedure was tested on 18 pairs of homologous proteins with unrelated function and 70 pairs of proteins with related function, and was shown to be 94 % accurate. This automated method could be linked to schemes for genome annotation. Finally, we examined the use of functional site prediction in protein-protein and protein-DNA docking. The use of predicted functional sites was shown to filter putative docked complexes with a discrimination similar to that obtained by manually including biological information about active sites or DNA-binding residues.

Algorithms↗

Rethinking "shape space": evidence from simulated docking suggests that steric shape complementarity is not limiting for antibody-antigen recognition and idiotypic interactions.

The concept of "shape space" is based on the assumption that the relevant properties of individual molecules can be adequately specified by a finite list of N parameters; and that cij, the affinity between molecules i and j, can be specified by an equation of the form: cij = f(xi, xj), where xi and xj are N-dimensional vectors representing the absolute positions of molecules i and j in an objective, referential "shape space", and f is an appropriate function. We have performed simulated docking of the combining sites of immunoglobulin molecules, based on their crystallographic structures. The results suggest that shape complementarity cannot account for the specificity of idiotypic interactions, since in the simulations each pair of docked proteins had a buried surface area as great as that occurring in known complexes. It therefore seems likely that the atomic interactions accounting for the specificity of immunoglobulin recognition are highly relational. This casts doubt on the basic assumptions underlying the shape-space concept, at least in the simple form hitherto used in theoretical modelling of the immune system. In order to be realistic, the dimensionality N would have to be high (more than 20), and the function f would be irregular and discontinuous. Alternatively, if the equation cij = f(xi, xj) is interpreted as a purely formal construction in an abstract "inversion space", its validity is entirely relative to the empirical affinity matrix on which the construction is based. We conclude that at present there is no sure way of adequately characterizing the internal structure of idiotypic affinity matrices; and that models of the immune system should therefore aim at being generic and robust with respect to the structure of the idiotypic affinity matrices of unselected immunoglobulins.

Algorithms↗

Flexibases: a way to enhance the use of molecular docking methods.

Specially expanded databases containing three-dimensional structures are created to enhance the utility of docking methods to find new leads, i.e., active compounds of pharmacological interest. The expansion is based on the automatic generation of a set of maximally dissimilar conformations. The ligand receptor system of methotrexate and dihydrofolate reductase is used to demonstrate the feasibility of creating flexibases and their utility in docking studies.

Algorithms↗