SEARCH · PubMed Health

Results for “Molecular Docking Simulation”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 811 records · Page 45Linked to original sources

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↗

PubMed2000

Possible influence of intramolecular hydrogen bonds on the three-dimensional structure of polyamine amides and their interaction with ionotropic glutamate receptors.

Polyamine amides are potent antagonists of many classes of ionotropic receptor. Here, calculations of the conformations of 26 polyamine amides using molecular mechanics methodology have shown that intramolecular hydrogen bonds strongly influence the in vacuo three-dimensional structure of a polyamine amide. Although these bonds are less stable in an aqueous environment, they may occur more when a polyamine amide interacts with a binding site. The estimated three-dimensional structures of polyamine amides provide an explanation for the differences in their antagonist potency at quisqualate-sensitive ionotropic glutamate receptors (qGluR) observed in experimental studies. Relative antagonist potency at qGluR is correlated with the number of free amino groups on a polyamine amide, i.e. those not involved in intramolecular hydrogen bonds. Also, intramolecular hydrogen bonds significantly restrict the conformational freedom of the uncharged moiety of a polyamine amide. Docking of polyamine amides to a molecular model of a mammalian AMPA receptor (GluR1) channel shows that intramolecular H-bonds may also provide a good structural explanation for the action of these compounds at this site.

Amides↗

PubMed2004

3D-QSAR CoMFA, CoMSIA studies on substituted ureas as Raf-1 kinase inhibitors and its confirmation with structure-based studies.

Three-dimensional quantitative structure activity relationship (3D-QSAR) analyses were carried out on 91 substituted ureas in order to understand their Raf-1 kinase inhibitory activities. The studies include Comparative Molecular Field Analysis (CoMFA) and Comparative Molecular Similarity Indices Analysis (CoMSIA). Models with good predictive abilities were generated with the cross validated r2 (r2cv) values for CoMFA and CoMSIA being 0.53 and 0.44, respectively. The conventional r2 values are 0.93 and 0.87 for CoMFA and CoMSIA, respectively. In addition, a homology model of Raf-1 was also constructed using the crystal structure of the kinase domain of B-Raf isoform with one of the most active Raf-1 inhibitors (48) inside the active site. The ATP binding pocket of Raf-1 is virtually similar to that of B-Raf. Selected ligands were docked in the active site of Raf-1. Molecule 48 adopts an orientation similar to that inside the B-Raf active site. The 4-pyridyl group bearing amide substituent is located in the adenosine binding pocket, and anchored to the protein through a pair of hydrogen bonds with Cys424 involving ring N-atom and amide NH group. The results of best 3D-QSAR model were compared with structure-based studies using the Raf-1 homology model. The results of 3D-QSAR and docking studies validate each other and provided insight into the structural requirements for activity of this class of molecules as Raf-1 inhibitors. Based on these results, novel molecules with improved activity can be designed.

Binding Sites↗

PubMed2004

Putative implication of alpha-amylase loop 7 in the mechanism of substrate binding and reaction products release.

Alpha-amylases are widespread endo-enzymes involved in the hydrolysis of internal alpha-(1,4) glycosidic linkages of starch polymers. Molecular modeling of amylose-amylase interactions is a step toward enzymatic mechanism understanding and rational design of new enzymes. From the crystallographic complex of barley alpha-amylase AMY2-acarbose, the static aspects of amylose-amylase docking have been characterized with a model of maltododecaose (DP12) (G. André, A. Buléon, R. Haser, and V. Tran, Biopolymers 1999, Vol. 50, pp. 751-762; G. André and V. Tran, Special Publication no. 246 1999, The Royal Society of Chemistry, H. J. Gilbert, G. J. Davies, B. Henrissat, and B. Svensson, Eds., Cambridge, pp. 165-174). These studies, consistent with the experimental subsite mapping (K. Bak-Jensen, G. André, V. Tran, and B. Svensson, Journal of Biological Chemistry, to be published), propose a propagation scheme for an amylose chain in the active cleft of AMY2. The topographical overview of alpha-amylases identified loop 7 as a conserved segment flanking the active site. Since some crystallographic experiments suspected its high flexibility, its putative motion was explored through a robotic scheme, an alternate route to dynamics simulations that consume CPU time. The present article describes the characteristics of the flexibility of loop 7: location and motion in AMY2. A back-and-forth motion with a large amplitude of more than 0.6 nm was evaluated. This movement could be triggered by two hinge residues. It results in the loop flipping over the active site to enhance the docking of the native helical substrate through specific interactions, it positions the catalytic residues, it distorts the substrate towards its transition state geometry, and finally monitors the release of the products after hydrolysis. The residues involved in the process are now rational mutation points in the hands of molecular biologists.

Animals↗

PubMed2001

Biophysical characterization of recombinant human Bcl-2 and its interactions with an inhibitory ligand, antimycin A.

Apoptosis is an essential physiological process, regulated by the family of Bcl-2-related proteins. However, the molecular mechanism by which Bcl-2 regulates apoptosis still remains elusive. Here we report the functional studies of recombinant human Bcl-2 with the deletion of 22 residues at the C-terminal membrane-anchoring region (rhBcl-2Delta22). Characterization of rhBcl-2Delta22 showed that the recombinant protein is homogeneous and monodisperse in nondenaturing solutions, stable at room temperature in the presence of a metal chelator, and an alpha-helical protein with unfolding of secondary structure at a T(m) of 62.8 degrees C. Optimal membrane pore formation by rhBcl-2Delta22 required negatively charged phospholipids. The existence of a hydrophobic groove in rhBcl-2Delta22 was demonstrated by the fluorescence enhancement of the hydrophobic ANS probe with which a pro-apoptotic Bak BH3 peptide competed. The respiratory inhibitor antimycin A also bound to the hydrophobic groove of rhBcl-2Delta22 with a K(d) of 0.82 microM. The optimal binding conformation of antimycin A was predicted from molecular docking of antimycin A with the hBcl-2 model created by homology modeling. Antimycin A selectively induces apoptosis in cells overexpressing Bcl-2, suggesting that hydrophobic groove-binding compounds may act as selective apoptotic triggers in tumor cells.

Amino Acid Sequence↗

PubMed2002

Sensitivity of molecular docking to induced fit effects in influenza virus neuraminidase.

Many proteins undergo small side chain or even backbone movements on binding of different ligands into the same protein structure. This is known as induced fit and is potentially problematic for virtual screening of databases against protein targets. In this report we investigate the limits of the rigid protein approximation used by the docking program, GOLD, through cross-docking using protein structures of influenza neuraminidase. Neuraminidase is known to exhibit small but significant induced fit effects on ligand binding. Some neuraminidase crystal structures caused concern due to the bound ligand conformation and GOLD performed poorly on these complexes. A 'clean' set, which contained unique, unambiguous complexes, was defined. For this set, the lowest energy structure was correctly docked (i.e. RMSD < 1.5 A away from the crystal reference structure) in 84% of proteins, and the most promiscuous protein (1mwe) was able to dock all 15 ligands accurately including those that normally required an induced fit movement. This is considerably better than the 70% success rate seen with GOLD against general validation sets. Inclusion of specific water molecules involved in water-mediated hydrogen bonds did not significantly improve the docking performance for ligands that formed water-mediated contacts but it did prevent docking of ligands that displaced these waters. Our data supports the use of a single protein structure for virtual screening with GOLD in some applications involving induced fit effects, although care must be taken to identify the protein structure that performs best against a wide variety of ligands. The performance of GOLD was significantly better than the GOLD implementation of ChemScore and the reasons for this are discussed. Overall, GOLD has shown itself to be an extremely good, robust docking program for this system.

Algorithms↗

PubMed2002

Kinesin: a molecular motor with a spring in its step.

A key step in the processive motion of two-headed kinesin along a microtubule is the 'docking' of the neck linker that joins each kinesin head to the motor's dimerized coiled-coil neck. This process is similar to the folding of a protein beta-hairpin, which starts in a highly mobile unfolded state that has significant entropic elasticity and finishes in a more rigid folded state. We therefore suggest that neck-linker docking is mechanically equivalent to the thermally activated shortening of a spring that has been stretched by an applied load. This critical tension-dependent step utilizes Brownian motion and it immediately follows the binding of ATP, the hydrolysis of which provides the free energy that drives the kinesin cycle. A simple three-state model incorporating neck-linker docking can account quantitatively for both the kinesin force-velocity relation and the unusual tension-dependence of its Michaelis constant. However, we find that the observed randomness of the kinesin motor requires a more detailed four-state model. Monte Carlo simulations of single-molecule stepping with this model illustrate the possibility of sub-8 nm steps, the size of which is predicted to vary linearly with the applied load.

Adenosine Triphosphate↗

PubMed2001

Protein conformational relaxation and ligand migration in myoglobin: a nanosecond to millisecond molecular movie from time-resolved Laue X-ray diffraction.

A time-resolved Laue X-ray diffraction technique has been used to explore protein relaxation and ligand migration at room temperature following photolysis of a single crystal of carbon monoxymyoglobin. The CO ligand is photodissociated by a 7.5 ns laser pulse, and the subsequent structural changes are probed by 150 ps or 1 micros X-ray pulses at 14 laser/X-ray delay times, ranging from 1 ns to 1.9 ms. Very fast heme and protein relaxation involving the E and F helices is evident from the data at a 1 ns time delay. The photodissociated CO molecules are detected at two locations: at a distal pocket docking site and at the Xe 1 binding site in the proximal pocket. The population by CO of the primary, distal site peaks at a 1 ns time delay and decays to half the peak value in 70 ns. The secondary, proximal docking site reaches its highest occupancy of 20% at approximately 100 ns and has a half-life of approximately 10 micros. At approximately 100 ns, all CO molecules are accounted for within the protein: in one of these two docking sites or bound to the heme. Thereafter, the CO molecules migrate to the solvent from which they rebind to deoxymyoglobin in a bimolecular process with a second-order rate coefficient of 4.5 x 10(5) M(-1) s(-1). Our results also demonstrate that structural changes as small as 0.2 A and populations of CO docking sites of 10% can be detected by time-resolved X-ray diffraction.

Animals↗

PubMed1994

Comparison of the three-dimensional molecular models of bovine submicellar caseins with small-angle X-ray scattering. Influence of protein hydration.

To test the applicability of two energy-minimized, three-dimensional structures of the bovine casein submicelle, theoretical small-angle X-ray scattering curves in the presence and absence of water were compared to experimental data. The published method simulates molecular dynamics of proteins in solution by employing adjustable Debye-Waller temperature factors (B factors) for the protein, for the solvent, and for protein-bound water. The programs were first tested upon bovine pancreatic trypsin inhibitor beginning with its known X-ray crystal structure. To approximate the degree of protein hydration previously determined by NMR relaxation experiments (0.01 g water/g protein), 120 water molecules were docked into the large void of the kappa-casein portion of the structure for both the symmetric and asymmetric casein submicelle models. To approximate hydrodynamic hydration (0.244 g water/g protein), 2703 water molecules were added to each of the above structures using the "droplet" algorithm in the Sybyl molecular modeling package. All structures were then energy-minimized and their solvation energies calculated. Theoretical small-angle X-ray scattering curves were calculated for all unhydrated and hydrated structures and compared with experimentally determined scattering profiles for submicellar casein. Best results were achieved with the 120-bound-water structure for both the symmetric and asymmetric submicelle models. Comparison of results for the protein submicelle models with those for the theoretical and literature values of bovine pancreatic trypsin inhibitor demonstrates the applicability of the methodology.

Animals↗

PubMed2004

Consequences of molecular-level Ca2+ channel and synaptic vesicle colocalization for the Ca2+ microdomain and neurotransmitter exocytosis: a monte carlo study.

Morphological and biochemical studies indicate association between voltage-gated Ca2+ channels and the vesicle docking complex at vertebrate presynaptic active zones, which constrain the separation between some Ca2+ channels and vesicles to 20 nm or less. To address the effect of the precise geometrical relationship among the vesicles, the Ca2+ channel, and the proteins of the release machinery on neurotransmitter release, we developed a Monte Carlo simulation of Ca2+ diffusion and buffering with nanometer resolution. We find that the presence of a vesicle as a diffusion barrier alters the shape of the Ca2+ microdomain of a single Ca2+ channel around the vesicle. This effect is maximal in the vicinity of the vesicle and depends critically on the vesicle's distance from the plasmalemma. Ca2+-sensor(s) for release would be exposed to markedly different [Ca2+], varying by up to 13-fold, depending on their position around the vesicle. As a result, the precise position of Ca2+-sensor(s) with respect to the vesicle and the channel can be critical to determining the release probability. Variation in the position of Ca2+-sensor molecule(s) and their accessibility could be an important source of heterogeneity in vesicle release probability.

Calcium↗

PubMed2004

The 'functional' dyad of scorpion toxin Pi1 is not itself a prerequisite for toxin binding to the voltage-gated Kv1.2 potassium channels.

Pi1 is a 35-residue scorpion toxin cross-linked by four disulphide bridges that acts potently on both small-conductance Ca2+-activated (SK) and voltage-gated (Kv) K+ channel subtypes. Two approaches were used to investigate the relative contribution of the Pi1 functional dyad (Tyr-33 and Lys-24) to the toxin action: (i) the chemical synthesis of a [A24,A33]-Pi1 analogue, lacking the functional dyad, and (ii) the production of a Pi1 analogue that is phosphorylated on Tyr-33 (P-Pi1). According to molecular modelling, this phosphorylation is expected to selectively impact the two amino acid residues belonging to the functional dyad without altering the nature and three-dimensional positioning of other residues. P-Pi1 was directly produced by peptide synthesis to rule out any possibility of trace contamination by the unphosphorylated product. Both Pi1 analogues were compared with synthetic Pi1 for bioactivity. In vivo, [A24,A33]-Pi1 and P-Pi1 are lethal by intracerebroventricular injection in mice (LD50 values of 100 and 40 microg/mouse, respectively). In vitro, [A24,A33]-Pi1 and P-Pi1 compete with 125I-apamin for binding to SK channels of rat brain synaptosomes (IC50 values of 30 and 10 nM, respectively) and block rat voltage-gated Kv1.2 channels expressed in Xenopus laevis oocytes (IC50 values of 22 microM and 75 nM, respectively), whereas they are inactive on Kv1.1 or Kv1.3 channels at micromolar concentrations. Therefore, although both analogues are less active than Pi1 both in vivo and in vitro, the integrity of the Pi1 functional dyad does not appear to be a prerequisite for the recognition and binding of the toxin to the Kv1.2 channels, thereby highlighting the crucial role of other toxin residues with regard to Pi1 action on these channels. The computed simulations detailing the docking of Pi1 peptides on to the Kv1.2 channels support an unexpected key role of specific basic amino acid residues, which form a basic ring (Arg-5, Arg-12, Arg-28 and Lys-31 residues), in toxin binding.

Amino Acid Sequence↗

PubMed1991

Protein docking and complementarity.

Predicting the structures of protein-protein complexes is a difficult problem owing to the topographical and thermodynamic complexity of these structures. Past efforts in this area have focussed on fitting the interacting proteins together using rigid body searches, usually with the conformations of the proteins as they occur in crystal structure complexes. Here we present work which uses a rigid body docking method to generate the structures of three known protein complexes, using both the bound and unbound conformations of the interacting molecules. In all cases we can regenerate the geometry of the crystal complexes to high accuracy. We also are able to find geometries that do not resemble the crystal structure but nevertheless are surprisingly reasonable both mechanistically and by some simple physical criteria. In contrast to previous work in this area, we find that simple methods for evaluating the complementarity at the protein-protein interface cannot distinguish between the configurations that resemble the crystal structure complex and those that do not. Methods that could not distinguish between such similar and dissimilar configurations include surface area burial, solvation free energy, packing and mechanism-based filtering. Evaluations of the total interaction energy and the electrostatic interaction energy of the complexes were somewhat better. Of the techniques that we tried, energy minimization distinguished most clearly between the "true" and "false" positives, though even here the energy differences were surprisingly small. We found the lowest total interaction energy from amongst all of the putative complexes generated by docking was always within 5 A root-mean-square of the crystallographic structure. There were, however, several putative complexes that were very dissimilar to the crystallographic structure but had energies that were close to that of the low energy structure. The magnitude of the error in energy calculations has not been established in macromolecular systems, and thus the reliability of the small differences in energy remains to be determined. The ability of this docking method to regenerate the crystallographic configurations of the interacting proteins using their unbound conformations suggests that it will be a useful tool in predicting the structures of unsolved complexes.

Algorithms↗

PubMed2006

eHiTS: a new fast, exhaustive flexible ligand docking system.

The flexible ligand docking problem is divided into two subproblems: pose/conformation search and scoring function. For successful virtual screening the search algorithm must be fast and able to find the optimal binding pose and conformation of the ligand. Statistical analysis of experimental data of bound ligand conformations is presented with conclusions about the sampling requirements for docking algorithms. eHiTS is an exhaustive flexible-docking method that systematically covers the part of the conformational and positional search space that avoids severe steric clashes, producing highly accurate docking poses at a speed practical for virtual high-throughput screening. The customizable scoring function of eHiTS combines novel terms (based on local surface point contact evaluation) with traditional empirical and statistical approaches. Validation results of eHiTS are presented and compared to three other docking software on a set of 91 PDB structures that are common to the validation sets published for the other programs.

Algorithms↗

PubMed2004

Protein-protein docking: is the glass half-full or half-empty?

Are current docking methods capable of building complexes from putative component protein structures? Results of recent computational studies, including those of the CAPRI (Critical Assessment of Protein Interactions) competition, were used to determine the key properties for successful docking and introduce a classification of protein complexes based on docking difficulty. Enzyme-inhibitor complexes could be determined with reasonable accuracy - possibly to within a few alternative structures. Results for antigen-antibody pairs are less predictable, and data for small signaling complexes are generally poor. However, moderate amounts of experimental data can remove uncertainty and the methodology is rapidly improving. Transient complexes with large interface areas undergo substantial conformational change and are beyond the reach of current docking methods. The docking of such complexes might therefore require fundamentally new approaches.

Algorithms↗

PubMed2002

Characterization of the vitamin E-binding properties of human plasma afamin.

Human plasma afamin, the fourth member of the albumin gene family, is shown to be a specific binding protein for vitamin E. A radio ligand-binding assay followed by Scatchard and Hill analysis are used to show that afamin has a binding affinity for both alpha-tocopherol and gamma-tocopherol, two of the most important forms of vitamin E, in vitro. The binding-dissociation constant was determined to be 18 microM, indicating that afamin plays a role as vitamin E carrier in body fluids such as human plasma and follicular fluid under physiological conditions. Additionally, we demonstrate that afamin has multiple binding sites for both alpha- and gamma-tocopherol. Due to the large binding capacity of afamin for vitamin E, it might take over the role of vitamin E transport in body fluids under conditions where the lipoprotein system is not sufficient for vitamin E transport. To confirm the experimental results, we performed homology modeling and docking calculations on the predicted tertiary structure, which showed coincidence between calculated and in vitro results.

Binding Sites↗

PubMed2006

Targeting protein-protein interactions with small molecules: challenges and perspectives for computational binding epitope detection and ligand finding.

A promising way to interfere with biological processes is through the control of protein-protein interactions by means of small molecules that modulate the formation of protein-protein complexes. Although the feasibility of this approach has been demonstrated in principle by recent results, many of the small-molecule modulators known to date have not been found by rational design approaches. In large part this is due to the challenges that one faces in dealing with protein binding epitopes compared to, e.g., enzyme binding pockets. Recent advances in the understanding of the energetics and dynamics of protein binding interfaces and methodological developments in the field of structure-based drug design methods may open up a way to apply rational design approaches also for finding protein-protein interaction modulators. These advances and developments include (I) computational approaches to dissect binding interfaces in terms of energetic contributions of single residues (to identify "hot spot" residues), (II) prediction of potential binding sites from unbound protein structures, (III) recognition of allosteric binding sites as alternatives to directly targeting interfaces, (IV) docking approaches that consider protein flexibility and improved descriptions of the solvent influence on electrostatic interactions, and (V) data-driven docking approaches. Here, we will summarize these developments with a particular emphasis on their applicability to screen for or design small-molecule modulators of protein-protein interactions.

Binding Sites↗

PubMed2002

Inhibition of trypsin by cowpea thionin: characterization, molecular modeling, and docking.

Higher plants produce several families of proteins with toxic properties, which act as defense compounds against pests and pathogens. The thionin family represents one family and comprises low molecular mass cysteine-rich proteins, usually basic and distributed in different plant tissues. Here, we report the purification and characterization of a new thionin from cowpea (Vigna unguiculata) with proteinase inhibitory activity. Cowpea thionin inhibits trypsin, but not chymotrypsin, binding with a stoichiometry of 1:1 as shown with the use of mass spectrometry. Previous annotations of thionins as proteinase inhibitors were based on their erroneous identification as homologues of Bowman-Birk family inhibitors. Molecular modeling experiments were used to propose a mode of docking of cowpea thionin with trypsin. Consideration of the dynamic properties of the cowpea thionin was essential to arrive at a model with favorable interface characteristics comparable with structures of trypsin-inhibitor complexes determined by X-ray crystallography. In the final model, Lys11 occupies the S1 specificity pocket of trypsin as part of a canonical style interaction.

Amino Acid Sequence↗

PubMed2004

Specificity of Amaranthus leucocarpus syn. hypocondriacus lectin for O-glycopeptides.

Amaranthus leucocarpus syn. hypochondriacus lectin (ALL) has been shown to be specific for N-acetyl-D-galactosamine (GalNAc). In this work, we determined a value of 1.0 x 10(-2) M for the association constant of ALL for GalNAc, calculated using fluorescence spectroscopy assays. Using neoglycopeptides obtained by in vitro O-glycosylation, we determined the main features of O-glycopeptides recognized by ALL using molecular dynamics simulations, capillary electrophoresis, and ELISA. Neo-glycopeptides were obtained by in vitro O-glycosylation reaction using microsomal preparations of murine thymocytes, human gastric fundus and colonic mucosa. ELISA assays were performed with peroxidase-labeled murine monoclonal IgG2, kappa light chain (5D4) antibodies against ALL. Among the in vitro neoglycopeptides, only those of TTSAPTTS containing GalNAc at Thr in #2 and #6 reacted with ALL. Neither the TTSAPTTS glycopeptide, containing a unique GalNAc residue at Thr in #2, nor others (with more than two GalNAc residues) interacted with the lectin. Computational docking assays of the lower energy conformers for interactions between glycopeptides and lectins confirmed that ALL recognized GalNAc residues when they are spaced out in glycan structures, whereas GalNAc residues arranged in clusters prevented interaction with the lectin, indicating that ALL is specific for a special GalNAc-containing motif found in different O-glycoproteins.

Acetylgalactosamine↗