SEARCH · PubMed Health

Results for “molecular docking”

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

PubMed2006

Computer-aided molecular design of highly potent HIV-1 RT inhibitors: 3D QSAR and molecular docking studies of efavirenz derivatives.

Ligand- and structure-based design approaches have been applied to an extended series of 74 efavirenz compounds effectively inhibiting wild type (WT) and mutant type (K103N) HIV-1 reverse transcriptase (RT). For ligand-based approach, three dimensional quantitative structure-activity relationship (3D-QSAR) methods, comparative molecular field analysis (CoMFA) and comparative similarity indices analysis (CoMSIA), were performed. The starting geometry of efavirenz was obtained from X-ray crystallographic data. The efavirenz derivatives were constructed and fully optimized by ab-initio molecular orbital method at HF/3-21G level. Reliable QSAR models for high predictive abilities were developed. Regarding WT and K103N inhibitions, CoMFA models with r2/cv = 0.651 and 0.678 and CoMSIA models with r2/cv = 0.662 and 0.743 were derived, respectively. The interpretation obtained from the models highlights different structural requirements for inhibition of WT and K103N HIV-1 RT. To elucidate potential binding modes of efavirenz derivatives in the binding pocket of WT and K103N HIV-1 RT, structure-based approach based on computational docking studies of selected efavirenz compounds were performed by using GOLD and FlexX programs. The results derived from docking analysis give additional information and further probe the inhibitor-enzyme interactions. The correlation of the results obtained from 3D QSAR and docking models validate each other and lead to better understanding of the structural requirements for the activity. Therefore, these integrated results are informative to provide key features and a helpful guideline for novel compound design active against HIV-1 RT.

Alkynes↗

PubMed2001

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↗

PubMed2005

Importance of accurate charges in molecular docking: quantum mechanical/molecular mechanical (QM/MM) approach.

The extent to which accuracy of electric charges plays a role in protein-ligand docking is investigated through development of a docking algorithm, which incorporates quantum mechanical/molecular mechanical (QM/MM) calculations. In this algorithm, fixed charges of ligands obtained from force field parameterization are replaced by QM/MM calculations in the protein environment, treating only the ligands as the quantum region. The algorithm is tested on a set of 40 cocrystallized structures taken from the Protein Data Bank (PDB) and provides strong evidence that use of nonfixed charges is important. An algorithm, dubbed "Survival of the Fittest" (SOF) algorithm, is implemented to incorporate QM/MM charge calculations without any prior knowledge of native structures of the complexes. Using an iterative protocol, this algorithm is able in many cases to converge to a nativelike structure in systems where redocking of the ligand using a standard fixed charge force field exhibits nontrivial errors. The results demonstrate that polarization effects can play a significant role in determining the structures of protein-ligand complexes, and provide a promising start towards the development of more accurate docking methods for lead optimization applications.

Algorithms↗

PubMed2006

Structure of the virus capsid protein VP1 of enterovirus 71 predicted by some homology modeling and molecular docking studies.

The homology modeling technique has been used to construct the structure of enterovirus 71 (EV 71) capsid protein VP1. The protein is consisted of 297 amino acid residues and treated as the target. The amino acid sequence identity between the target protein and sequences of template proteins 1EAH, 1PIV, and 1D4M searched from NCBI protein BLAST and WorkBench protein tools were 38, 37, and 36%, respectively. Based on these template structures, the protein model was constructed by using the InsightII/Homology program. The protein model was briefly refined by energy minimization and molecular dynamics (MD) simulation steps. The protein model was validated using some web available servers such as ERRAT, PROCHECK, PROVE, and PROSA2003. However, an inconsistency between the docking scores and the measured activity was observed for a series of EV 71 VP1 inhibitors synthesized by Shia et al. (J Med Chem 2002, 45, 1644) and docked into the binding pocket of the protein model using the DOCK 4.0.2 program. The protein model with an EV 71 VP1 inhibitor docked and engulfed was then refined further by some MD simulation steps in the presence of water molecules. The docking scores obtained for these inhibitors after such a MD refinement were well correlated with the activities. The structure-activity relationships for the ligand-protein model system was also analyzed using the GRID-VOLSURF programs and the corresponding noncrossvalidated and crossvalidated (by leave-one-out) r2 and q2 were 0.99 and 0.61, respectively. The hydrophobic nature of the binding pocket of the protein model was also examined using the GRID21 program. The possibility of improving the potency of the current series of EV 71 VP1 inhibitors was discussed based on all the studies presented.

Amino Acid Sequence↗

PubMed2003

A novel scoring function for molecular docking.

We present a novel scoring function for docking of small molecules to protein binding sites. The scoring function is based on a combination of two main approaches used in the field, the empirical and knowledge-based approaches. To calibrate the scoring function we used an iterative procedure in which a ligand's position and its score were determined self-consistently at each iteration. The scoring function demonstrated superiority in prediction of ligand positions in docking tests against the commonly used Dock, FlexX and Gold docking programs. It also demonstrated good accuracy of binding affinity prediction for the docked ligands.

Binding Sites↗

PubMed2004

Molecular docking of four beta-amyloid1-42 fragments on the alpha7 nicotinic receptor: delineating the binding site of the Abeta peptides.

Three-dimensional structures of the complexes between the Abeta(1-42) fragments Abeta(1-11), Abeta(10-20), Abeta(12-28), and Abeta(22-35) and the alpha7 nicotinic receptor were obtained with the aid of the ESCHER program. Furthermore, short high-temperature molecular dynamics simulations in vacuo were employed to relax the complexes and allow the peptides to accommodate in the binding site. The final models have shown that Abeta peptides do bind on the same site, which is delineated by loop C of one subunit and the loops 62-74 and G of the adjacent subunit on the receptor. This finding is supported by previous experimental and theoretical data, and should help one to obtain a better and more detailed structural information about the activity of the Abeta peptides and their repercussion in the disorders at molecular level, which are characteristic of the Alzheimer's disease.

Amyloid beta-Peptides↗

PubMed2003

Molecular docking studies of natural cholinesterase-inhibiting steroidal alkaloids from Sarcococca saligna.

Alkaloids isolated from Sarcococca saligna significantly inhibit acetyl- and butyrylcholinesterase enzyme, suggesting discovery of inhibitors for nervous-system disorders. Studying interactions with the active site of the AChE enzyme from Torpedo californica, we have identified hydrophobic interactions inside the aromatic gorge area as the major stabilizing factor in enzyme-inhibitor complexes of these alkaloids. Molecular Dynamics simulation of a predicted complex indicates that ligand binding does not extensively alter enzyme structure, but reduces flexibility at the gorge.

Acetylcholinesterase↗

PubMed1994

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↗

PubMed2006

Efficient molecular docking of NMR structures: application to HIV-1 protease.

Docking ligands into an ensemble of NMR conformers is essential to structure-based drug discovery if only NMR structures are available for the target. However, sequentially docking ligands into each NMR conformer through standard single-receptor-structure docking, referred to as sequential docking, is computationally expensive for large-scale database screening because of the large number of NMR conformers involved. Recently, we developed an efficient ensemble docking algorithm to consider protein structural variations in ligand binding. The algorithm simultaneously docks ligands into an ensemble of protein structures and achieves comparable performance to sequential docking without significant increase in computational time over single-structure docking. Here, we applied this algorithm to docking with NMR structures. The HIV-1 protease was used for validation in terms of docking accuracy and virtual screening. Ensemble docking of the NMR structures identified 91% of the known inhibitors under the criterion of RMSD < 2.0 A for the best-scored conformation, higher than the average success rate of single docking of individual crystal structures (66%). In the virtual screening test, on average, ensemble docking of the NMR structures obtained higher enrichments than single-structure docking of the crystal structures. In contrast, docking of either the NMR minimized average structure or a single NMR conformer performed less satisfactorily on both binding mode prediction and virtual screening, indicating that a single NMR structure may not be suitable for docking calculations. The success of ensemble docking of the NMR structures suggests an efficient alternative method for standard single docking of crystal structures and for considering protein flexibility.

Algorithms↗

PubMed2004

Lessons in molecular recognition: the effects of ligand and protein flexibility on molecular docking accuracy.

The key to success for computational tools used in structure-based drug design is the ability to accurately place or "dock" a ligand in the binding pocket of the target of interest. In this report we examine the effect of several factors on docking accuracy, including ligand and protein flexibility. To examine ligand flexibility in an unbiased fashion, a test set of 41 ligand-protein cocomplex X-ray structures were assembled that represent a diversity of size, flexibility, and polarity with respect to the ligands. Four docking algorithms, DOCK, FlexX, GOLD, and CDOCKER, were applied to the test set, and the results were examined in terms of the ability to reproduce X-ray ligand positions within 2.0A heavy atom root-mean-square deviation. Overall, each method performed well (>50% accuracy) but for all methods it was found that docking accuracy decreased substantially for ligands with eight or more rotatable bonds. Only CDOCKER was able to accurately dock most of those ligands with eight or more rotatable bonds (71% accuracy rate). A second test set of structures was gathered to examine how protein flexibility influences docking accuracy. CDOCKER was applied to X-ray structures of trypsin, thrombin, and HIV-1-protease, using protein structures bound to several ligands and also the unbound (apo) form. Docking experiments of each ligand to one "average" structure and to the apo form were carried out, and the results were compared to docking each ligand back to its originating structure. The results show that docking accuracy falls off dramatically if one uses an average or apo structure. In fact, it is shown that the drop in docking accuracy mirrors the degree to which the protein moves upon ligand binding.

Algorithms↗

PubMed2000

Deciphering common failures in molecular docking of ligand-protein complexes.

Common failures in predicting crystal structures of ligand-protein complexes are investigated for three ligand-protein systems by a combined thermodynamic and kinetic analysis of the binding energy landscapes. Misdocked predictions in ligand-protein docking are classified as 'soft' and 'hard' failures. While a soft failure arises when the search algorithm is unable to find the global energy minimum corresponding to the crystal structure, a hard failure results from a flaw of the energy function to qualify the crystal structure as the predicted lowest energy conformation in docking simulations. We find that neither the determination of a single structure with the lowest energy nor finding the most common binding mode is sufficient to predict crystal structures of the complexes, which belong to the category of hard failures. In a proposed hierarchical approach, structural similarity clustering of the conformations, generated from equilibrium simulations with the simplified energy function, is followed by energy refinement with the AMBER force field. This protocol, that involves a hierarchy of energy functions, resolves some common failures in ligand-protein docking and detects crystallographic binding modes that were not found during docking simulations.

Crystallography↗

PubMed2005

Granuphilin molecularly docks insulin granules to the fusion machinery.

The Rab27a effector granuphilin is specifically localized on insulin granules and is involved in their exocytosis. Here we show that the number of insulin granules morphologically docked to the plasma membrane is markedly reduced in granuphilin-deficient beta cells. Surprisingly, despite the docking defect, the exocytosis of insulin granules in response to a physiological glucose stimulus is significantly augmented, which results in increased glucose tolerance in granuphilin-null mice. The enhanced secretion in mutant beta cells is correlated with a decrease in the formation of the fusion-incompetent syntaxin-1a-Munc18-1 complex, with which granuphilin normally interacts. Furthermore, in contrast to wild-type granuphilin, its mutant that is defective in binding to syntaxin-1a fails to restore granule docking or the protein level of syntaxin-1a in granuphilin-null beta cells. Thus, granuphilin not only is essential for the docking of insulin granules but simultaneously imposes a fusion constraint on them through an interaction with the syntaxin-1a fusion machinery. These findings provide a novel paradigm for the docking machinery in regulated exocytosis.

Adenoviridae↗

PubMed1990

Distance-constrained molecular docking by simulated annealing.

An optimized method based on the principle of simulated annealing is presented for determining the relative position and orientation of interacting molecules. The spatial relationships of these molecules are described by intermolecular distance constraints between specific pairs of atoms, such as found in hydrogen bonds or from experimentally determined data. The method makes use of a random walk through six rotational and translational degrees of freedom where the constituent molecules are treated as rigid bodies. Van der Waals repulsions are used only to define a lower bound on distances between constrained atom pairs within the docking procedure. A cost function comprised of purely geometric constraints is optimized via simulated annealing, in order to search for the best orientation and position of the two molecules. Our docking procedure is applied to eight serine proteinase complexes from the Brookhaven Protein Data Bank. For each simulation 100 computations were performed. A typical docking computation requires only a few seconds of CPU time on a VAXserver 3500. The influence of the number of constraints on the final docked positions was studied. The sensitivity of the docking procedure to a ligand structure which is not well defined is also addressed. Possible applications of this method include using approximate distances incorporating complete energy functions.

Aspartic Acid Endopeptidases↗

PubMed2006

Homology modeling and molecular docking study of translationally controlled tumor protein and artemisinin.

Translationally controlled tumor protein (TCTP), also known as histamine releasing factor (HRF), is found abundantly in different eukaryotic cell types. The sequence homology of TCTP between different species is very high, belonging to the MSS4/DSS4 superfamily of proteins. TCTP is involved in both cell growth and human late allergy reaction, as well as having a calcium binding property; however, its primary biological functions remain to be clearly elucidated. In regard to many possible functions, the TCTP of Plasmodium falciparum (Pf) is known to bind with an antimalarial agent, artemisinin, which is activated by heme. It is assumed that the endoperoxide-bridge of artemisinin is opened up by heme to form a free radical, which then eventually alkylates, probably to the Cys14 of PfTCTP. Study of the docking of artemisinin with heme, and subsequently with PfTCTP, was carried out to verify the above hypothesis on the basis of structural interactions. The three dimensional (3D) structure of PfTCTP was built by homology modeling, using the NMR structure of the TCTP of Schizosaccharomyces pombe as a template. The quality of the model was examined based on its secondary structure and biological function, as well as with the use of structure evaluating programs. The interactions between artemisinin, heme and PfTCTP were then studied using the docking program, FlexiDock. The center of the peroxide bond of artemisinin and the Fe of heme were docked within a short distance of 2.6A, implying the strong possibility of an interaction between the two molecules, as proposed. When the activated form of artemisinin was docked on the PfTCTP, the C4-radical of the drug faced towards the sulfur of Cys14 within a distance of 2.48A, again suggesting the possibility of alkylation having occurred. These results confirm the proposed mechanism of the antimalarial effect of artemisinin, which will provide a reliable method for establishing the mechanism of its biological activity using a molecular modeling study.

Antimalarials↗

PubMed1997

CASP2 molecular docking predictions with the LIGIN software.

Seven docking predictions were made with the LIGIN program. In six cases the location of the binding pocket was identified correctly by systematically docking everywhere within the protein structure. In two cases the ligand was docked to within 1.8 A RMSD of the experimentally determined structure. LIGIN has not been optimized to deal with highly flexible ligands that dock at the surface of proteins. Consequently, in three cases the exposed part of the ligand was docked poorly, although the buried parts were docked well, and made similar atomic contacts with the protein as in the experimentally determined structure.

Amiloride↗

PubMed2006

Molecular docking studies on 4-thiazolidinones as HIV-1 RT inhibitors.

Flexible docking simulations were performed on two series of 4-thiazolidinones as HIV-1 reverse transcriptase (HIV-1 RT) inhibitors. This was done by analyzing the interaction of these compounds with the allosteric site of the HIV-1 reverse transcriptase enzyme. The binding scores for these compounds were also congruent with their anti-HIV activity. A good correlation between the predicted binding free energies and the experimentally observed inhibitory activities (EC(50)) suggest that the identified binding conformations of these inhibitors are reliable. The results of docking studies provide an insight into the pharmacophoric structural requirements for the HIV-1 RT inhibitory activity of this class of molecules.

Allosteric Site↗

PubMed2004

Performance of 3D-database molecular docking studies into homology models.

The performance of docking studies into protein active sites constructed by homology model building was investigated using CDK2 and factor VIIa screening data sets. When the sequence identity between model and template near the binding site area is greater than approximately 50%, roughly 5 times more active compounds are identified than would be found randomly. This performance is comparable to docking to crystal structures.

Binding Sites↗

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↗