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

Computational docking of L-arginine and its structural analogues to C-terminal domain of Escherichia coli arginine repressor protein (ArgRc).

The arginine repressor (ArgR) of Escherichia coli binds to six L-arginine molecules that act as its co-repressor in order to bind to DNA. The binding of L-arginine molecules as well as its structural analogues is compared by means of computational docking. A grid-based energy evaluation method combined with a Monte Carlo simulated annealing process was used in the automated docking. For all ligands, the docking procedure proposed more than one binding site in the C-terminal domain of ArgR (ArgRc). Interaction patterns of ArgRc with L-arginine were also observed for L-canavanine and L-citrulline. L-lysine and L-homoarginine, on the other hand, were shown to bind poorly at the binding site. Figure A general overview of the sites found from docking the various ligands into ArgRc ( grey ribbons). Red coloured sticks: residues in binding site H that was selected for docking

Arginine↗

New algorithm to model protein-protein recognition based on surface complementarity. Applications to antibody-antigen docking.

A novel algorithm is presented which models protein-protein interactions using surface complementarity. The method is applied to antibody-antigen docking. A steric scoring scheme, based upon a soft potential, is used to assess complementarity, and a simple electrostatic model is then used to remove infeasible interactions. The soft potential allows for structural changes that occur during docking. Biochemical knowledge is necessary to reduce the number of docking orientations produced by the method to a manageable size. The information used includes the known epitope residues and a single loose distance constraint. The method is applied to all three crystallographically determined antibody-lysozyme complexes, HyHEL-10, D1.3 and HyHEL-5. For the first time, a predicted antibody structure (that of D1.3) is used as a docking target. In the four systems modelled, the method identifies between 15 and 40 possible docking orientations. The root-mean-square (r.m.s.) deviation between these orientations and the relevant crystallographic complex is measured in the interface region. For all four complexes an orientation is found with r.m.s. deviation in the range 1.9 A and 4.8 A. The algorithm is implemented on a single instruction/multiple datastream (SI/MD) architecture computer. The use of a parallel architecture computer ensures detailed coverage of the search space, whilst still maintaining a search time of two days.

Algorithms↗

c-Jun can recruit JNK to phosphorylate dimerization partners via specific docking interactions.

Structurally related serine/threonine kinases recognize similar phosphoacceptor peptides in vitro yet in vivo, they phosphorylate distinct substrates. To understand the basis for this specificity, we studied the interaction between the Jun kinases (JNKs) and Jun proteins. JNKs phosphorylate c-Jun very efficiently, JunD less efficiently, but they do not phosphorylate JunB. Effective JNK substrates require a separate docking site and specificity-conferring residues flanking the phosphoacceptor. The docking site increases the efficiency and specificity of the phosphorylation reaction. JunB has a functional JNK docking site but lacks specificity-conferring residues. Insertion of such residues brings JunB under JNK control. JunD, by contrast, lacks a JNK docking site, but its phosphoacceptor peptide is identical to that of c-Jun. Substrates such as JunD can be phosphorylated by JNK through heterodimerization with docking competent partners. Therefore, heterodimerization can affect the recognition of transcription factors by signal-regulated protein kinases.

Amino Acid Sequence↗

Molecular docking studies on interaction of diverse retinol structures with human alcohol dehydrogenases predict a broad role in retinoid ligand synthesis.

Some members of the human alcohol dehydrogenase (ADH) family possess retinol dehydrogenase activity and may thus function in production of the active nuclear receptor ligand retinoic acid. Many diverse natural forms of retinol exist including all-trans-retinol (vitamin A(1)), 9-cis-retinol, 3,4-didehydroretinol (vitamin A(2)), 4-oxo-retinol, and 4-hydroxy-retinol as well as their respective carboxylic acid derivatives which are active ligands for retinoid receptors. This raises the question of whether ADHs can accommodate all these different retinols and thus participate in the activation of several retinoid ligands. The crystal structures of human ADH1B and ADH4 provide the opportunity to examine their active sites for potential binding to many diverse retinol structures using molecular docking algorithms. The criteria used to score successful docking included achievement of distances of 1.9-2.4 A between the catalytic zinc and the hydroxyl oxygen of retinol and 3.2-3.6 A between C-4 of the coenzyme NAD and C-15 of retinol. These distances are sufficient to enable hydride transfer during the oxidation of an alcohol to an aldehyde. By these criteria, all-trans-retinol, 4-oxo-retinol, and 4-hydroxy-retinol were successfully docked to both ADH1B and ADH4. However, 9-cis-retinol and 3,4-didehydroretinol, which have more restrictive conformations, were successfully docked to only ADH4 which possesses a wider active site than ADH1B and more easily accommodates the C-19 methyl group. Furthermore, docking of all retinols was more favorable in the active site of ADH4 rather than ADH1B as measured by force field and contact scores. These findings suggest that ADH1B has a limited capacity to metabolize retinols, but that ADH4 is well suited to function in the metabolism of many diverse retinols and is predicted to participate in the synthesis of the active ligands all-trans-retinoic acid, 9-cis-retinoic acid, 3, 4-didehydroretinoic acid, 4-oxo-retinoic acid, and 4-hydroxy-retinoic acid.

Alcohol Dehydrogenase↗

The SH2/SH3 adaptor protein dock interacts with the Ste20-like kinase misshapen in controlling growth cone motility.

Recent studies suggest that the SH2/SH3 adaptor Dock/Nck transduces tyrosine phosphorylation signals to the actin cytoskeleton in regulating growth cone motility. The signaling cascade linking the action of Dock/Nck to the reorganization of cytoskeleton is poorly understood. We now demonstrate that Dock interacts with the Ste20-like kinase Misshapen (Msn) in the Drosophila photoreceptor (R cell) growth cones. Loss of msn causes a failure of growth cones to stop at the target, a phenotype similar to loss of dock, whereas overexpression of msn induces pretarget growth cone termination. Physical and genetic interactions between Msn and Dock indicate a role for Msn in the Dock signaling pathway. We propose that Msn functions as a key controller of growth cone cytoskeleton in response to Dock-mediated signals.

Adaptor Proteins, Signal Transducing↗

Docking of sulfonamides to carbonic anhydrase II and IV.

Starting with a known active site of a protein and a database of compounds, one would like to quickly identify a few compounds that "dock" into the active site and obtain "good" binding free energies. The main goal of current automated docking procedures is to predict the "best" substrate-enzyme complex while other programs such as UHBD and DelPhi can be used to compute binding free energies. In this paper, we will focus on the application of docking methods and parameters to study substrate-enzyme interactions of a metalloenzyme system. Specifically, we report on the docking of sulfonamides to carbonic anhydrase II and IV, which are of interest due to their application in glaucoma therapy. Using a standard docking protocol, it is possible to correctly predict not only the orientation of inhibitors to a specific isozyme, but also determine the qualitative affinity for a group of inhibitor for an isozyme.

Algorithms↗

Evaluation of docking strategies for virtual screening of compound databases: cAMP-dependent serine/threonine kinase as an example.

In an effort to establish efficient docking routines for computational screening of compound databases on protein structures, cAMP-dependent protein kinase has been selected as a test case and a variety of docking options and scoring functions were compared. These included rigid-body and flexible docking and scoring based on surface complementarity and/or force field energy. Inhibitors were removed from complex crystal structures and added to compound libraries in their binding conformations and, in addition, deliberately modified conformations. Rigid-body docking and contact scoring well reproduced two of three experimental enzyme-inhibitor complexes. Ligand docking with flexible torsional angles failed to do so but anchored search of some inhibitors converged near to experimental structures, however, only when energy scoring was applied.

1-(5-Isoquinolinesulfonyl)-2-Methylpiperazine↗

Thermodynamic and activation parameters for binding of a pyrene-labeled substrate by the Tetrahymena ribozyme: docking is not diffusion-controlled and is driven by a favorable entropy change.

Association and dissociation rates for the pyrene-(pyr)-labeled oligoribonucleotide substrate pyrCUCU binding to the L-21 ScaI group I ribozyme are reported as a function of temperature. Combined with thermodynamic parameters for binding of pyrCUCU to rGGAGAA, the results allow calculation of the activation and thermodynamic parameters for docking of pyrCUCU into the catalytic core of the ribozyme. The activation enthalpy for docking is 22 kcal/mol, much larger than the approximately 4 kcal/mol expected for a diffusion-controlled process. Thus, docking is not diffusion-controlled. The activation and equilibrium entropies for docking are favorable at 21 and 37 eu, respectively. The results suggest the rate-limiting step and the driving force for docking may involve desolvation of RNA functional groups or of Mg2+ ions.

Animals↗

Mutations in the Tetrahymena ribozyme internal guide sequence: effects on docking of the P1 helix into the catalytic core and correlation with catalytic activity.

Binding of substrate by the ribozyme derived from the self-splicing intron of Tetrahymena thermophila involves at least two steps. In the first step, base pairing between the ribozyme internal guide sequence (IGS) and the substrate forms a helical duplex (P1). Through specific tertiary interactions between P1 and the ribozyme core, P1 is then docked into the ribozyme active site. We have investigated the effects of compensatory mutations in positions 2-6 of the P1 helix on docking of P1 into the ribozyme core. Equilibrium binding of matching oligonucleotides by catalytically active IGS mutant ribozymes was evaluated by gel-shift analysis. While the strength of base pairing changed with base composition as expected, the strength of tertiary interactions between P1 and the ribozyme core was not affected by the P1 mutations. These results support a model in which efficient docking of P1 is determined by P1 structure and the presence of a conserved G-U pair. Determination of the rate of dissociation of matching oligonucleotides from each ribozyme revealed that mutations in the IGS change the tightness of binding by increasing or decreasing the dissociation rate. Surprisingly, dissociation rates determined in this fashion were 20-900-fold less than the values of the multiple-turnover rate constant for these ribozymes, initially suggesting that turnover did not require product dissociation. A more detailed analysis for the wild-type ribozyme defined two distinct product dissociation rates. The slower rate equaled that determined under the conditions used for the equilibrium binding studies. The weighted average of the two dissociation rates equaled the multiple-turnover rate constant. These results are explained by a model in which ribozyme preparations consist of two ribozyme conformers: one with tight docking of P1 and another with weaker docking of P1.

Animals↗

New approach to molecular docking and its application to virtual screening of chemical databases

This paper describes the validation of a molecular docking method and its application to virtual database screening. The code flexibly docks ligand molecules into rigid receptor structures using a tabu search methodology driven by an empirically derived function for estimating the binding affinity of a protein-ligand complex. The docking method has been tested on 70 ligand-receptor complexes for which the experimental binding affinity and binding geometry are known. The lowest energy geometry produced by the docking protocol is within 2.0 A root mean square of the experimental binding mode for 79% of the complexes. The method has been applied to the problem of virtual database screening to identify known ligands for thrombin, factor Xa, and the estrogen receptor. A database of 10,000 randomly chosen "druglike" molecules has been docked into the three receptor structures. In each case known receptor ligands were included in the study. The results showed good separation between the predicted binding affinities of the known ligand set and the database subset.

Journal Article↗

Surflex: fully automatic flexible molecular docking using a molecular similarity-based search engine.

Surflex is a fully automatic flexible molecular docking algorithm that combines the scoring function from the Hammerhead docking system with a search engine that relies on a surface-based molecular similarity method as a means to rapidly generate suitable putative poses for molecular fragments. Results are presented evaluating reliability and accuracy of dockings compared with crystallographic experimental results on 81 protein/ligand pairs of substantial structural diversity. In over 80% of the complexes, Surflex's highest scoring docked pose was within 2.5 A root-mean-square deviation (rmsd), with over 90% of the complexes having one of the top ranked poses within 2.5 A rmsd. Results are also presented assessing Surflex's utility as a screening tool on two protein targets (thymidine kinase and estrogen receptor) using data sets on which competing methods were run. Performance of Surflex was significantly better, with true positive rates of greater than 80% at false positive rates of less than 1%. Docking time was roughly linear in number of rotatable bonds, beginning with a few seconds for rigid molecules and adding approximately 10 s per rotatable bond.

Algorithms↗

Multiple automatic base selection: protein-ligand docking based on incremental construction without manual intervention.

A possible way of tackling the molecular docking problem arising in computer-aided drug design is the use of the incremental construction method. This method consists of three steps: the selection of a part of a molecule, a so-called base fragment, the placement of the base fragment into the active site of a protein, and the subsequent reconstruction of the complete drug molecule. Assuming that a part of a drug molecule is known, which is specific enough to be a good base fragment, the method is proven to be successful for a large set of docking examples. In addition, it leads to the fastest algorithms for flexible docking published so far. In most real-world applications of docking, large sets of ligands have to be tested for affinity to a given protein. Thus, manual selection of a base fragment is not practical. On the other hand, the selection of a base fragment is critical in that only few selections lead to a low-energy structure. We overcome this limitation by selecting a representative set of base fragments instead of a single one. In this paper, we present a set of rules and algorithms to automate this selection. In addition, we extend the incremental construction method to deal with multiple fragmentations of the drug molecule. Our results show that with multiple automated base selection, the quality of the docking predictions is almost as good as with one manually preselected base fragment. In addition, the set of solutions is more diverse and alternative binding modes with low scores are found. Although the run time of the overall algorithm increases, the method remains fast enough to search through large ligand data sets.

Automation↗

MCDOCK: a Monte Carlo simulation approach to the molecular docking problem.

Prediction of the binding mode of a ligand (a drug molecule) to its macromolecular receptor, or molecular docking, is an important problem in rational drug design. We have developed a new docking method in which a non-conventional Monte Carlo (MC) simulation technique is employed. A computer program, MCDOCK, was developed to carry out the molecular docking operation automatically. The current version of the MCDOCK program (version 1.0) allows for the full flexibility of ligands in the docking calculations. The scoring function used in MCDOCK is the sum of the interaction energy between the ligand and its receptor, and the conformational energy of the ligand. To validate the MCDOCK method, 19 small ligands, the binding modes of which had been determined experimentally using X-ray diffraction, were docked into their receptor binding sites. To produce statistically significant results, 20 MCDOCK runs were performed for each protein-ligand complex. It was found that a significant percentage of these MCDOCK runs converge to the experimentally observed binding mode. The root-mean-square (rms) of all non-hydrogen atoms of the ligand between the predicted and experimental binding modes ranges from 0.25 to 1.84 A for these 19 cases. The computational time for each run on an SGI Indigo2/R10000 varies from less than 1 min to 15 min, depending upon the size and the flexibility of the ligands. Thus MCDOCK may be used to predict the precise binding mode of ligands in lead optimization and to discover novel lead compounds through structure-based database searching.

Binding Sites↗

DREAM++: flexible docking program for virtual combinatorial libraries.

We present a set of programs, DREAM+2 (Docking and Reaction programs using Efficient seArch Methods written in C++), for docking computationally generated ligands into macromolecular binding sites. DREAM++ is composed of three programs: ORIENT++, REACT++ and SEARCH++. The program ORIENT++ positions molecules in a binding site with the DOCK algorithm. Its output can be used as input to REACT++ and SEARCH+2. The program REACT++ performs user-specific chemical reactions on a docked molecule, so that reaction products can be evaluated for three dimensional complementarity with the macromolecular site. The program SEARCH++ performs an efficient conformation search on the reaction products using a hybrid backtrack and incremental construction algorithm. We have applied the programs to HIV protease-inhibitor complexes as test systems. We found that we can differentiate high-affinity ligands based on several measures: interaction energies, occupancy of protein subsites and the number of successfully docked conformations for each product. Encouraged by the results in the test case, we applied the programs to propose novel inhibitors of HIV protease. These inhibitors can be generated by organic reactions using commercially available reagents. They are alternatives to the inhibitors synthesized by Glaxo.

Algorithms↗

Protein ligand docking based on empirical method for binding affinity estimation.

An empirical protein-ligand binding affinity estimation method, SCORE, was incorporated into a popular docking program, DOCK4. The combined program, ScoreDock, was used to reconstruct the 200 protein-ligand complex structures and found to give good results for the complexes with high binding affinities. A quality assessment method for docking results from ScoreDock was developed based on the whole test set and tested by additionally selected complexes. The method significantly improves the docking accuracy and was shown to be reliable in docking quality assessment. As a docking tool in structural based drug design, ScoreDock can screen out final hits directly based on the predicted negative logarithms of dissociation equilibrium constants of protein-ligand complexes, and can explicitly deal with structure water molecules, as well as metal atoms.

Computer Simulation↗

SNAP-25 and synaptotagmin 1 function in Ca2+-dependent reversible docking of granules to the plasma membrane.

In neuroendocrine cells, Ca2+ triggers fusion of granules with the plasma membrane and functions at earlier steps by increasing the size of the readily releasable pool of vesicles. The effect of Ca2+ at early steps of secretion may be due to the recruitment at the plasma membrane of granules localized in the cytoplasm. To study the mechanism of granule docking, a new in vitro assay is designed using membrane fractions from mouse pituitary AtT-20 cells. By using this assay, it is found that granule docking to the plasma membrane is controlled by Ca2+ concentrations in the micromolar range, is reversible and requires intact SNAP-25, but not VAMP-2. In the docking assay, addition of Ca2+ induces the formation of a SNAP-25-Synaptotagmin 1 complex. The cytosolic domain C2AB of Synaptotagmin 1 and anti-Synaptotagmin 1 antibodies block granule docking. These results show that Ca2+ modulates dynamic docking of granules to the plasma membrane and that this process is due to a Ca2+-dependent interaction between SNAP-25 and Synaptotagmin 1.

Animals↗

Homology of 54K protein of signal-recognition particle, docking protein and two E. coli proteins with putative GTP-binding domains.

Most proteins exported from mammalian cells contain a signal sequence which mediates targeting to and insertion into the membrane of the endoplasmic reticulum (ER). Involved in this process are the signal-recognition particle (SRP) and docking protein (DP), the receptor for SRP in the ER membrane. SRP interacts with the signal sequence on nascent polypeptide chains and retards their further elongation, which resumes only after interaction of the arrested ribosomal complex with the docking protein. SRP is a ribonucleoprotein particle comprising a 7S RNA and six polypeptides with relative molecular masses (Mr) of 9,000 (9K) 14K, 19K, 54K, 68K and 72K (ref. 1). The 9K and 14K proteins are essential for elongation arrest and the 68K-72K heterodimer is required for docking to the ER membrane. The 54K protein binds to the signal sequence when it emerges from the ribosome. Docking protein consists of two polypeptides, a 72K alpha-subunit (DP alpha) and a 30K beta-subunit (DP beta). No components structurally homologous to SRP and docking protein have yet been found in yeast or Escherichia coli. To understand the molecular nature of the interaction between the signal sequence and its receptor(s) we have characterized a complementary DNA coding for the 54K protein of SRP. Significant sequence homology was found to part of DP alpha and two E. coli proteins of unknown function. The homologous region includes a putative GTP-binding domain.

Amino Acid Sequence↗

Docking sites on mitogen-activated protein kinase (MAPK) kinases, MAPK phosphatases and the Elk-1 transcription factor compete for MAPK binding and are crucial for enzymic activity.

Mitogen-activated protein kinase (MAPK) cascades control gene expression patterns in response to extracellular stimuli. MAPK/ERK (extracellular-signal-regulated kinase) kinases (MEKs) activate MAPKs by phosphorylating them; activated MAPKs, in turn, phosphorylate target transcription factors, and are deactivated by phosphatases. One mechanism for maintaining signal specificity and efficiency is the interaction of MAPKs with their substrates and regulators through high-affinity docking sites. In the present study, we show that peptides corresponding to the MAPK-docking sites of MEK1, MEK2, Ste7, Elk-1 and MAPK phosphatase (MKP)-2 potently inhibit MEK2 phosphorylation of ERK2, ERK2 phosphorylation of Elk-1, and MKP-1 dephosphorylation of ERK2. Each peptide inhibited multiple reactions; for example, the MEK2 peptide inhibited not only MEK2, but also ERK2 and MKP-1. In addition, these docking-site peptides inhibited MEK2-ERK2 binding. The MAPK-docking site of MEK1 also potently stimulated ERK2-mediated phosphorylation of a target site on the same peptide. Control peptides with mutations of conserved basic and hydrophobic residues of the MAPK-docking site consensus lacked biological activity. We conclude that MEKs, MKPs and the Elk-1 transcription factor compete for binding to the same region of ERK2 via protein-protein interactions that are crucial for kinase/phosphatase activity.

Amino Acid Sequence↗