Search PubMedSearch

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

Identification of the basic residues of cytochrome f responsible for electrostatic docking interactions with plastocyanin in vitro: relevance to the electron transfer reaction in vivo.

The prominent basic patch seen in the atomic structure of the lumen-side domain of turnip cytochrome f, consisting of Arg209 and Lys187, 58, 65, and 66, was proposed to be an electrostatically complementary docking site for its physiological electron acceptor, plastocyanin [Martinez, S. E., Huang, D., Szczepaniak, A., Cramer, W. A., and Smith, J. L. (1994) Structure 2, 95-105]. This proposal agrees with solution studies on the cytochrome f/plastocyanin electron-transfer reaction that showed a major contribution of electrostatic interactions to the docking, but not with studies on the oxidation rate of cyt f in vivo using mutants in which the basic patch of cyt f was neutralized. The apparent contradiction might be explained by an unknown electron acceptor protein for cyt f. However, (i) flash-induced oxidation of cyt f is absent in a PC-deficient mutant. (ii) Lys58, 65, and 66 in the large domain and Lys188 and 189 in the small domain are major contributors to the ionic strength dependence of the electron-transfer reaction in solution. Replacement of Lys58 and 65 by neutral residues and of Lys66 by the acidic residue Glu66 resulted in a >10-fold decrease in the rate of electron transfer in solution and complete loss of its ionic strength dependence. Replacement of Lys188 and Lys189 in the small domain of cyt f resulted in a 3-4-fold decrease in the second-order rate constant and a smaller dependence of the overall rate of electron transfer on ionic strength, corresponding to a loss of two positive charges. (iii) Acidification of the thylakoid lumen cannot explain the absence of electrostatic interactions. (iv) Changing the five lysines to acidic residues did not result in any significant retardation of the rate of cyt f oxidation in vivo. If the docking of cyt f and plastocyanin in vivo is mediated by basic residues of cyt f, they are different from those that mediate electron transfer in vitro or that are implicated by simulations of electrostatic interactions of the docking. Alternatively, docking of cyt f/PC in vivo is limited by spatial constraints or release of PC from P700 that precludes a rate-limiting mediation of the cyt f/PC reaction by specific electrostatic interactions. The cyt f/PC system in Chlamydomonas reinhardtii is the first electron-transfer couple for which the role of electrostatics in mediating the docking reaction has been studied both in vitro and in vivo.

Animals

Previously identified protein of uncertain function is karyopherin alpha and together with karyopherin beta docks import substrate at nuclear pore complexes.

Previously, we had purified a cytosolic protein complex, termed karyopherin, that functions in docking import substrate at the nuclear envelope in digitonin-permeabilized cells and also had molecularly cloned and sequenced its 97-kDa beta subunit. We now report that the karyopherin alpha subunit is the previously identified protein NPI-1/SRP-1 of hitherto uncertain function. Using purified recombinant karyopherin alpha or beta subunit, we showed that neither karyopherin alpha nor karyopherin beta alone was sufficient for docking of import substrate at the nuclear envelope. Docking occurred only when both subunits were present. Moreover, docking of import substrate by the two recombinant karyopherin subunits was productive, as it led to nuclear internalization of the docked substrate in the presence of additional, previously characterized cytosolic factors. In a binding assay using immobilized karyopherin alpha and beta subunits and import substrate as a ligand, we found that only karyopherin alpha bound ligand. We suggest that karyopherin beta functions as an adaptor that binds both to karyopherin alpha and to any of a large number of docking sites that are represented by a repetitive peptide motif containing nucleoporins on both the cytoplasmic and nucleoplasmic side of the nuclear pore complex (NPC), bidirectionally ferrying a complex of karyopherin alpha-substrate across the NPC.

Animals

A method for biomolecular structural recognition and docking allowing conformational flexibility.

In this work, we present an algorithm developed to handle biomolecular structural recognition problems, as part of an interdisciplinary research endeavor of the Computer Vision and Molecular Biology fields. A key problem in rational drug design and in biomolecular structural recognition is the generation of binding modes between two molecules, also known as molecular docking. Geometrical fitness is a necessary condition for molecular interaction. Hence, docking a ligand (e.g., a drug molecule or a protein molecule), to a protein receptor (e.g., enzyme), involves recognition of molecular surfaces. Conformational transitions by "hinge-bending" involves rotational movements of relatively rigid parts with respect to each other. The generation of docked binding modes between two associating molecules depends on their three dimensional structures (3-D) and their conformational flexibility. In comparison to the particular case of rigid-body docking, the computational difficulty grows considerably when taking into account the additional degrees of freedom intrinsic to the flexible molecular docking problem. Previous docking techniques have enabled hinge movements only within small ligands. Partial flexibility in the receptor molecule is enabled by a few techniques. Hinge-bending motions of protein receptors domains are not addressed by these methods, although these types of transitions are significant, e.g., in enzymes activity. Our approach allows hinge induced motions to exist in either the receptor or the ligand molecules of diverse sizes. We allow domains/subdomains/group of atoms movements in either of the associating molecules. We achieve this by adapting a technique developed in Computer Vision and Robotics for the efficient recognition of partially occluded articulated objects. These types of objects consist of rigid parts which are connected by rotary joints (hinges). Our method is based on an extension and generalization of the Hough transform and the Geometric Hashing paradigms for rigid object recognition. We show experimental results obtained by the successful application of the algorithm to cases of bound and unbound molecular complexes, yielding fast matching times. While the "correct" molecular conformations of the known complexes are obtained with small RMS distances, additional, predictive good-fitting binding modes are generated as well. We conclude by discussing the algorithm's implications and extensions, as well as its application to investigations of protein structures in Molecular Biology and recognition problems in Computer Vision.

Algorithms

Initial docking of ER-derived vesicles requires Uso1p and Ypt1p but is independent of SNARE proteins.

ER-to-Golgi transport in yeast may be reproduced in vitro with washed membranes, purified proteins (COPII, Uso1p and LMA1) and energy. COPII coated vesicles that have budded from the ER are freely diffusible but then dock to Golgi membranes upon the addition of Uso1p. LMA1 and Sec18p are required for vesicle fusion after Uso1p function. Here, we report that the docking reaction is sensitive to excess levels of Sec19p (GDI), a treatment that removes the GTPase, Ypt1p. Once docked, however, vesicle fusion is no longer sensitive to GDI. In vitro binding experiments demonstrate that the amount of Uso1p associated with membranes is reduced when incubated with GDI and correlates with the level of membrane-bound Ypt1p, suggesting that this GTPase regulates Uso1p binding to membranes. To determine the influence of SNARE proteins on the vesicle docking step, thermosensitive mutations in Sed5p, Bet1p, Bos1p and Sly1p that prevent ER-to-Golgi transport in vitro at restrictive temperatures were employed. These mutations do not interfere with Uso1p-mediated docking, but block membrane fusion. We propose that an initial vesicle docking event of ER-derived vesicles, termed tethering, depends on Uso1p and Ypt1p but is independent of SNARE proteins.

Animals

Acute cortisol responses of lambs castrated and docked using rubber rings with or without a castration clamp.

OBJECTIVE: To measure the plasma cortisol response in lambs castrated and docked by three different methods (ring castration and tail docking, ring and clamp castration and tail docking; ring and clamp castration and ring tail docking) for 4 h after treatment. DESIGN: A physiological study with controls. PROCEDURE: A rubber ring was applied to the scrotum. Then the castration clamp was placed, distal to the ring, on each spermatic cord such that there was no overlap of the crush lines. RESULTS: There were no differences in the cortisol secretion of the lambs castrated and tail docked by the ring only and by the ring plus castration clamp methods. Using the castration clamp in addition to the ring on the tail had no effect on cortisol secretion. CONCLUSIONS: Using the castration clamp did not appear to reduce the cortisol response to ring castration and docking, as demonstrated elsewhere, because the areas of uncrushed tissue between the two castration clamp crush lines allowed nociception from ischemic scrotal tissue to be transmitted cranially via undamaged nerve fibers.

Analysis of Variance

Lamb tail docking: a controlled field study of the effects of tail amputation on health and productivity.

A detailed study of the effects of tail docking on lamb health and productivity was prompted by current concern about the welfare aspects of lamb tail amputation. Using a controlled field trial, comparing over 3000 docked and undocked lambs on seven farms, the effects of tail docking on mortality, blowfly strike and production variables were examined. The incidence of blowfly strike was strongly and consistently higher in undocked than docked lambs (rate ratio 6.03, 95 per cent confidence interval (CI) 2.99 to 12.19 for male lambs and 4.25, 95 per cent CI 2.25 to 8.01 for female lambs). The incidence of faecal soiling of the breech was slightly higher in undocked lambs and was identified as an important independent risk factor for blowfly strike. Both the mortality and production parameters were similar for docked and undocked lambs. It was concluded that tail docking protected against blowfly strike, with little evidence of any detrimental effect on lamb mortality and production.

Amputation, Surgical

Flexible ligand docking: a multistep strategy approach.

A flexible ligand docking protocol based on a divide-and-conquer strategy is investigated. This approach first separates total search space into conformation and orientation space. It uses a grid-based method to sample the conformation of an unbound ligand and to select the low-energy conformers. Rigid docking is then carried out to locate the low-energy binding orientations for these conformers. These docking structures are subsequently subjected to structure refinement including molecular mechanics minimization, conformational scanning at the binding site and a short period of molecular dynamics-based simulated annealing. This approach has been applied to twelve ligand-protein complexes with three to sixteen rotatable bonds. The docked lowest-energy structures have root mean square deviations ranging from 0.64 A to 2.01 A with respect to the corresponding crystal structures. The effect of atomic charges and van der Waals parameters on the docking results, and the role of the dielectric constant in the conformation sampling are discussed in detail. A fragment-based docking approach that takes advantages of the divide-and-conquer strategy has also been explored and the results are compared with those produced by a whole molecule-based approach.

Hydrogen Bonding

A structural and functional analysis of the docking protein. Characterization of active domains by proteolysis and specific antibodies.

Docking protein is a 73-kDa integral membrane protein of the rough endoplasmic reticulum. It is essential for translocation of nascent secretory proteins into the lumen of the endoplasmic reticulum. Monoclonal and polyclonal antibodies have been generated which, in conjunction with limited proteolysis, have been used to characterize several subspecies of docking protein. These proteolytic fragments have been analyzed with respect to the various functions ascribed to docking protein which can be assayed in vitro. Proteolytic digestion of membrane-associated or of affinity-purified intact docking protein showed that: elastase cleavage generates a 59-kDa soluble fragment and one of 14 kDa which contains the membrane anchoring domain; trypsin as well as endogenous proteolysis generates a 46-kDa fragment, leaving a 27-kDa domain containing the membrane anchor. This 27-kDa fragment can be reduced to a 13- and a 14-kDa piece by elastase digestion. The characteristics of these various subspecies were examined. The 59-kDa soluble fragment, which can reconstitute full translocation activity to docking protein-depleted microsomes (Meyer, D. I., and Dobberstein, B. (1980) J. Cell Biol. 87, 503-508) was capable of releasing a signal recognition particle-mediated translation arrest. The 46-kDa fragment was neither able to reassociate with nor to reconstitute the activity of docking protein-depleted microsomes. Moreover this fragment was unable to release a signal recognition particle-mediated arrest. This suggests that the 13-kDa fragment (the difference between 46 and 59 kDa) is both essential for association with the membrane, and for the release of translation arrests.

Animals

Surface solid angle-based site points for molecular docking.

We are developing a new site descriptor for the DOCK molecular modeling program suite. Sphgen, the current site description program for the DOCK suite, describes the pockets of a macromolecule by filling a volume with intersecting spheres. DOCK then identifies possible ligand orientations in the pocket by overlapping the atoms of proposed ligands with the sphere centers. Sphgen limits use of the DOCK program to concave binding regions, but macromolecular binding regions can be solvent-exposed rather than buried pockets. We present a more general site descriptor, based on the surface solid angle, which generates site points by determining the solid angle of exposure for points on the surface of the molecule, then identifying patches of surface with similar solid angle values which are then built into site points. We find possible ligand orientations by matching shape-based site points on the ligand and protein and demanding complementary solid angle values. Orientations are evaluated using the DOCK's force field-based score, which evaluates the Coulombic and van der Waals energy. The surface solid angle descriptor displays the complementary characteristics of the interfaces of our test systems: trypsin/trypsin inhibitor, chymotrypsin/turkey ovomucoid third domain, and subtilisin/chymotrypsin inhibitor. The solid angle site points can be used by DOCK to generate orientations within 1.5 A r.m.s.d. of the crystal structure orientation.

Algorithms

Flexible docking allowing induced fit in proteins: insights from an open to closed conformational isomers.

Here we dock a ligand onto a receptor surface allowing hinge-bending domain/substructural movements. Our approach mimics and manifests induced fit in molecular recognition. All angular rotations are allowed on the one hand, while a conformational space search is avoided on the other. Rather than dock each of the molecular parts separately with subsequent reconstruction of the consistently docked molecules, all parts are docked simultaneously while still utilizing the position of the hinge from the start. Like pliers closing on a screw, the receptor automatically closes on its ligand in the best surface-matching way. Movements are allowed either in the ligand or in the larger receptor, hence reproducing induced molecular fit. Hinge bending movements are frequently observed when molecules associate. There are numerous examples of open versus closed conformations taking place upon binding. Such movements are observed when the substrate binds to its respective enzyme. In particular, such movements are of interest in allosteric enzymes. The movements can involve entire domains, subdomains, loops, (other) secondary structure elements, or between any groups of atoms connected by flexible joints. We have implemented the hinges at points and at bonds. By allowing 3-dimensional (3-D) rotation at the hinge, several rotations about (consecutive or nearby) bonds are implicitly taken into account. Alternatively, if required, the point rotation can be restricted to bond rotation. Here we illustrate this hinge-bending docking approach and the insight into flexibility it provides on a complex of the calmodulin with its M13 ligand, positioning the hinges either in the ligand or in the larger receptor. This automated and efficient method is adapted from computer vision and robotics. It enables utilizing entire molecular surfaces rather than focusing a priori on active sites. Hence, allows attaining the overall optimally matching surfaces, the extent and type of motions which are involved. Here we do not treat the conformational flexibility of side-chains or of very small pieces of the molecules. Therefore, currently available methods addressing these issues and the method presented here, are complementary to each other, expanding the repertoire of computational docking tools foreseen to aid in studies of recognition, conformational flexibility and drug design.

Algorithms

Flexible docking using Tabu search and an empirical estimate of binding affinity.

This article describes the implementation of a new docking approach. The method uses a Tabu search methodology to dock flexibly ligand molecules into rigid receptor structures. It uses an empirical objective function with a small number of physically based terms derived from fitting experimental binding affinities for crystallographic complexes. This means that docking energies produced by the searching algorithm provide direct estimates of the binding affinities of the ligands. The method has been tested on 50 ligand-receptor complexes for which the experimental binding affinity and binding geometry are known. All water molecules are removed from the structures and ligand molecules are minimized in vacuo before docking. The lowest energy geometry produced by the docking protocol is within 1.5 A root-mean square of the experimental binding mode for 86% of the complexes. The lowest energies produced by the docking are in fair agreement with the known free energies of binding for the ligands.

Algorithms

Critical evaluation of the research docking program for the CASP2 challenge.

The binding positions of six small-molecule ligands in their complexes with target proteins were predicted using our Research docking program for the CASP2 challenge. Research uses a Monte Carlo procedure with pairwise energies and allows for the conformational searching of ligand torsional space. We were able to predict 2 of the 5 noncovalent complexes within 2 A root-mean-square (RMS) of the experimental structures as ranked by interaction energy or by a score calculated using our interaction evaluation program, Out-rank. In addition, for 4 of the 5 noncovalent structures we found a docking within 2 A RMS of the experimental structure within the top 20 dockings as ranked by energy. The main limitation in our approach is in the ability of the energy function and Outrank to discriminate among the lowest energy dockings. On the other hand, our success in exploring the multi-dimensional docking space of position, orientation and conformation is encouraging.

Arabinose

Prediction of IgE(Lb4)-ligand complex structures by automated docking.

A mouse monoclonal anti-2,4,6-trinitrophenyl IgE (clone Lb4) was screened with a random set of over 2000 compounds, and several ligands were found to bind with affinities comparable to that of the immunizing hapten (KD in the microM range). An automated docking algorithm was used for the prediction of complex structures formed by 2,4-dinitrophenyl (DNP) and non-DNP ligands in the fragment variable region of IgE(Lb4). All ligands were found to dock in an L-shaped cavity of 15 x 16 x 10 A, surrounded by complementary-determining regions L1, L3, H2 and H3. The ligands were found to occupy the same binding site in different orientations. For rigid ligands the most stable orientation could be predicted with high probability, based on the calculated energy of binding and the occurrence frequencies of identical complexes obtained by repeated simulations. The localization of a flexible ligand (cycrimine-R) was more ambiguous, but it still docked in the same site. The results support a model for heteroligating antibody (Ab) binding sites, where different ligands utilize the total set of available contacts in different combinations. It is suggested that although pseudoenergies calculated by the docking algorithm do not correlate with experimentally measured binding energies, the screening-and-docking procedure can be useful for the mapping of Ab and other receptor binding sites ligating small molecules.

2,4-Dinitrophenol

Flexible ligand docking using conformational ensembles.

Molecular docking algorithms suggest possible structures for molecular complexes. They are used to model biological function and to discover potential ligands. A present challenge for docking algorithms is the treatment of molecular flexibility. Here, the rigid body program, DOCK, is modified to allow it to rapidly fit multiple conformations of ligands. Conformations of a given molecule are pre-calculated in the same frame of reference, so that each conformer shares a common rigid fragment with all other conformations. The ligand conformers are then docked together, as an ensemble, into a receptor binding site. This takes advantage of the redundancy present in differing conformers of the same molecule. The algorithm was tested using three organic ligand protein systems and two protein-protein systems. Both the bound and unbound conformations of the receptors were used. The ligand ensemble method found conformations that resembled those determined in X-ray crystal structures (RMS values typically less than 1.5 A). To test the method's usefulness for inhibitor discovery, multi-compound and multi-conformer databases were screened for compounds known to bind to dihydrofolate reductase and compounds known to bind to thymidylate synthase. In both cases, known inhibitors and substrates were identified in conformations resembling those observed experimentally. The ligand ensemble method was 100-fold faster than docking a single conformation at a time and was able to screen a database of over 34 million conformations from 117,000 molecules in one to four CPU days on a workstation.

Algorithms

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

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