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

PubMed2005

Approaching the CAPRI challenge with an efficient geometry-based docking.

The last 3 rounds (3-5) of CAPRI included a wide range of docking targets. Several targets were especially challenging, since they involved large-scale movements and symmetric rearrangement, while others were based on homology models. We have approached the targets with a variety of geometry-based docking algorithms that include rigid docking, symmetric docking, and flexible docking with symmetry constraints. For all but 1 docking target, we were able to submit at least 1 acceptable quality prediction. Here, we detail for each target the prediction methods used and the specific biological data employed, and supply a retrospective analysis of the results. We highlight the advantages of our techniques, which efficiently exploit the geometric shape complementarity properties of the interaction. These enable them to run only few minutes on a standard PC even for flexible docking, thus proving their scalability toward computational genomic scale experiments. We also outline the major required enhancements, such as the introduction of side-chain position refinement and the introduction of flexibility for both docking partners.

Algorithms↗

PubMed2005

Scoring docking models with evolutionary information.

We have developed methods for the extraction of evolutionary information from multiple sequence alignments for use in the study of the evolution of protein interaction networks and in the prediction of protein interaction. For Rounds 3, 4, and 5 of the CAPRI experiment, we used scores derived from the analysis of multiple sequence alignments to submit predictions for 7 of the 12 targets. Our docking models were generated with Hex and GRAMM, but all our predictions were selected using methods based on multiple sequence alignments and on the available experimental evidence. With this approach, we were able to predict acceptable level models for 4 of the targets, and for a fifth target, we located the residues involved in the binding surface. Here we detail our successes and highlight several of the limitations and problems that we faced while dealing with particular docking cases.

Algorithms↗

PubMed2006

Design, synthesis, and evaluation of new type of L-amino acids containing pyridine moiety as nitric oxide synthase inhibitor.

New amino acids 7-12 were designed and synthesized as candidate inhibitors of human nitric oxide synthase (NOS). The 2-aminopyridine-containing l-amino acids 8 had potent inhibitory activity toward all of the human NOS isozymes. However, the regioisomers 9 and 10, and 2-methylpyridine-containing compound 11 had much lower inhibitory activity. Human NOS isozymes were also inhibited by 7, which lacks an amino group on the pyridine moiety. A computational docking study was carried out to investigate the mechanism of the inhibitory effect.

Amino Acids↗

PubMed2005

Peptide-cleaving catalyst selective for peptide deformylase.

A peptide-cleaving catalyst selective for peptide deformylase (PDF) was obtained from a library containing about 15 000 catalyst candidates. The catalyst cleaved the polypeptide backbone of PDF at Gln(152)-Arg(153). Docking simulations suggested multiple modes of interactions in the complex formed between the catalyst and PDF.

Amides↗

PubMed2006

Docking of photosystem I subunit C using a constrained geometric simulation.

The elucidation of assembly pathways of multi-subunit protein complexes is a problem of great interest in structural biology and biomolecular modeling. In this study, we use a new computer algorithm for the simulation of large-scale motion in proteins to dock the subunit PsaC onto Photosystem I. We find that a complicated docking pathway involving multiple conformational changes can be quickly simulated by actively targeting only a few residues at a time to their target positions. Simulations for two possible docking scenarios are explored, and experimental approaches to distinguish between them are discussed.

Algorithms↗

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↗

PubMed1997

Molecular modeling of the neurophysin I/oxytocin complex.

Neurophysins I and II (NPI and NPII) act in the neurosecretory granules as carrier proteins for the neurophyseal hormones oxytocin (OT) and vasopressin (VP), respectively. The NPI/OT functional unit, believed to be an (NPI/OT)2 heterotetramer, was modeled using low-resolution structure information, viz. the C alpha carbon atom coordinates of the homologous NPII/dipeptide complex (file 1BN2 in the Brookhaven Protein Databank) as a template. Its all-atom representation was obtained using standard modeling tools available within the INSIGHT/Biopolymer modules supplied by Biosym Technologies Inc. A conformation of the NPI-bound OT, similar to that recently proposed in a transfer NOE experiment, was docked into the ligand-binding site by a superposition of its Cys1-Tyr2 fragment onto the equivalent portion of the dipeptide in the template. The starting complex for the initial refinements was prepared by two alternative strategies, termed Model I and Model II, each ending with a approximately 100 ps molecular dynamics (MD) simulation in water using the AMBER 4.1 force field. The free homodimer NPI2 was obtained by removal of the two OT subunits from their sites, followed by a similar structure refinement. The use of Model I, consisting of a constrained simulated annealing, resulted in a structure remarkably similar to both the NPII/dipeptide complex and a recently published solid-state structure of the NPII/OT complex. Thus, Model I is recommended as the method of choice for the preparation of the starting all-atom data for MD. The MD simulations indicate that, both in the homodimer and in the heterotetramer, the 3(10)-helices demonstrate an increased mobility relative to the remaining body of the protein. Also, the C-terminal domains in the NPI2 homodimer are more mobile than the N-terminal ones. Finally, a distinct intermonomer interaction is identified, concentrated around its most prominent, although not unique, contribution provided by an H-bond from Ser25 O gamma in one NPI unit to Glu81 O epsilon in the other unit. This interaction is present in the heterotetramer (NPI/OT)2 and absent or weak in the NPI2 homodimer. We speculate that this interaction, along with the increased mobility of the 3(10)-helices and the carboxy domains, may contribute to the allosteric communication between ligand binding and NPI dimerization.

Allosteric Regulation↗

PubMed2005

From structure to function: methods and applications.

The rapid increase in experimental data along with recent progress in computational methods has brought modern biology a step closer toward solving one of the most challenging problems: prediction of protein function. Comprehension of protein function at its most basic level requires understanding of molecular interactions. Currently, it is becoming universally accepted that the scale of the accumulated data for analysis and for prediction necessitate highly efficient computational tools with appropriate application capabilities. The review presents the up-to-date advances in computational methods for structural pattern discovery and for prediction of molecular associations. We focus on their applications toward a range of biological problems and highlight the advantages of the combination of these methods and their integration with biological experiments. We provide examples, synergistically merging structural modeling, rigid and flexible structural alignment and detection of conserved structural patterns and docking (rigid and flexible with hinge-bending movements). We hope the review will lead to a broader utilization of computational methods, and their cross-fertilization with experiment.

Algorithms↗

PubMed2005

Different environments for a realistic simulation of GPCRs-application to the M2 muscarinic receptor.

A model of the human M(2) muscarinic receptor was taken as an example for a class A G-protein coupled receptor to explore the influence of different environments in a molecular dynamics simulation (MDS) on the protein structure. The most commonly used environment is the vacuum, although it is very unnatural for a transmembrane protein. As an alternative a membrane-like system, consisting of a lipophilic central layer and two aqueous flanking layers, was tested. The most realistic system that can be applied is a phospholipid bilayer with a surrounding physiological sodium chloride solution. From all systems good protein structures were received, nevertheless clear differences between the systems were detected in the structural comparison of the models. Subsequently it was analyzed whether the observed structural differences influence ligand binding. For this purpose the antagonist (S)-scopolamine was docked into the binding cavity, which is well known by many reported single and multiple point mutations. As expected from the observed structural variations triggered by the type of environment employed in MDS, also differences in the binding mode of (S)-scopolamine were detected, all contacts, however, which are known to be important were found.

Chemistry, Pharmaceutical↗

PubMed2005

An effective simulation of aqueous micellar aggregates by computational models.

We have computationally studied the interaction modes, localization and orientation of a benzene (Bz) molecule on the surface of micelles formed by cetyltrimethylammonium salts CTAX. Experimental 1H-NMR data on complexation shifts induced by Bz on the polar head hydrogens and on the adjacent methylene hydrogens of CTAX have been interpreted using a computational approach that combines an automatic molecular docking procedure with a calculation module that accounts for NMR complexation shifts due to ring current diamagnetic anisotropy. Three different models were used to reduce the complexity of the micellar system. Computational results, in good agreement with available experimental data, point to a preferential localization of the Bz molecule along the CTAX alkyl tail, about 3.9 angstroms away from the charged nitrogen. The Bz molecular plane is predicted perpendicular to the C-H bonds of the alkyl tail. The good results obtained with the simplest model suggest that it could be used to study more complex systems involving surfactants endowed with molecular recognition or catalytic abilities.

Catalysis↗

PubMed2001

FlexE: efficient molecular docking considering protein structure variations.

Side-chain or even backbone adjustments upon docking of different ligands to the same protein structure, a phenomenon known as induced fit, are frequently observed. Sometimes point mutations within the active site influence the ligand binding of proteins. Furthermore, for homology derived protein structures there are often ambiguities in side-chain placement and uncertainties in loop modeling which may be critical for docking applications. Nevertheless, only very few molecular docking approaches have taken into account such variations in protein structures. We present the new software tool FlexE which addresses the problem of protein structure variations during docking calculations. FlexE can dock flexible ligands into an ensemble of protein structures which represents the flexibility, point mutations, or alternative models of a protein. The FlexE approach is based on a united protein description generated from the superimposed structures of the ensemble. For varying parts of the protein, discrete alternative conformations are explicitly taken into account, which can be combinatorially joined to create new valid protein structures.FlexE was evaluated using ten protein structure ensembles containing 105 crystal structures from the PDB and one modeled structure with 60 ligands in total. For 50 ligands (83 %) FlexE finds a placement with an RMSD to the crystal structure below 2.0 A. In all cases our results are of similar quality to the best solution obtained by sequentially docking the ligands into all protein structures (cross docking). In most cases the computing time is significantly lower than the accumulated run times for the single structures. FlexE takes about five and a half minutes on average for placing one ligand into the united protein description on a common workstation. The example of the aldose reductase demonstrates the necessity of considering protein structure variations for docking calculations. We docked three potent inhibitors into four protein structures with substantial conformational changes within the active site. Using only one rigid protein structure for screening would have missed potential inhibitors whereas all inhibitors can be docked taking all protein structures into account.

Aldehyde Reductase↗

PubMed2001

Molecular simulation of the interaction of kappa-conotoxin-PVIIA with the Shaker potassium channel pore.

Molecular simulation techniques were appplied to predict the interaction of the voltage-dependent Shaker potassium channel with the channel-blocking toxin kappa-conotoxin-PVIIA (PVIIA). A structural thee-dimensional model of the extracellular vestibule of the potassium channel was constructed based on structural homologies with the bacterial potassium channel Kcsa, whose structure has been solved by X-ray crystallography. The docking of the PVIIA molecule was obtained by a geometric recognition algorithm, yielding 100 possible conformations. A series of residue-residue distance restraints, predicted from mutation-cycle experiments, were used to select a small set of a plausible channel-toxin complex models among the resulting possible conformations. The four final conformations, with similar characteristics, can explain most of the single-point mutation experiments done with this system. The models of the Shaker-PVIIA interaction predict two clusters of amino acids, critical for the binding of the toxin to the channel. The first cluster is the amino acids R2, I3, Q6 and K7 that form the plug of the toxin that interacts with the entrance to the selectivity filter of the channel. The second cluster of residues, R22, F23, N24 and K25, interacts with a channel region near to the external entrance of the pore vestibule. The consistency of the obtained models and the experimental data indicate that the Shaker-PVIIA complex model is reasonable and can be used in further biological studies such as the rational design of blocking agents of potassium channels and the mutagenesis of both toxins and potassium channels.

Algorithms↗

PubMed1992

Molecular modeling of beta-cyclodextrin complexes with nootropic drugs.

The geometry and structural features of the inclusion complexes of beta-cyclodextrin (beta-CD) with the chiral antiamnesic drugs (+/-)-1-benzyl-4-hydroxymethylpyrrolidin-2-one (WEB-1868). (+/-)-1-benzenesulfonyl-5-ethoxypyrrolidin-2-one (RU-35929), and (+/-)-1-(3-pyridinlysulfonyl)-5-ethoxypyrrolidin-2-one (RU-47010) were studied by the molecular modeling method (MacroModel interactive computer program). Docking procedures yielded the most stable complexes, which showed the aromatic ring of the guests inside the cavity and the pyrrolidinone ring out from the side of the beta-CD secondary hydroxyl groups. The binding energies were essentially due to hydrogen-bonded structures involving the C=O group of the guests. Selective interactions allowed chiral discrimination, and accordingly, separate beta-CD complexes of the R and S enantiomers of each guest compound were studied. The almost round beta-CD structure, in all the cases, assumed an elliptic shape on passing from the isolated molecule to the docked complex. The optimized structures and conformations of beta-CD and its inclusion compounds showed acceptable general agreement with information from proton nuclear magnetic resonance studies.

Computer Simulation↗

PubMed2006

Molecular evolution of adrenoceptors and dopamine receptors: implications for the binding of catecholamines.

We derived homology models for all human catecholamine-binding GPCRs (CABRs; the alpha-1, alpha-2, and beta-adrenoceptors and the D1-type and D2-type dopamine receptor) using the bovine rhodopsin-11-cis-retinal X-ray structure. Interactions were predicted from the endogenous ligands norepinephrine or dopamine and from the binding site and were used to optimize receptor-ligand interactions. Similar binding modes in the complexes agree with a large "binding core" conserved across the CABRs, that is, D3.32, V(I)3.33, T3.37, S5.42, S(A/C)5.43, S5.46, F6.51, F6.52, and W6.48. Model structures and docking simulations suggest that extracellular loop 2 could provide a common attachment point for the ligands' beta-hydroxyl via a hydrogen bond donated by the main-chain NH group of residue xl2.52. The modeled CABRs and docking modes are in good agreement with published experimental studies. Complementarity between the ligand and the binding site suggests that the bovine rhodopsin structure is a suitable template for modeling agonist-bound CABRs.

Adrenergic Agonists↗

PubMed2005

ATTRACT: protein-protein docking in CAPRI using a reduced protein model.

Protein-protein complex structures have been predicted for CAPRI Rounds 3 and 5 using a reduced protein model. Proteins are represented by up to 3 pseudoatoms per amino acid. The docking approach termed ATTRACT is based on energy minimization in translational and rotational degrees of freedom of one protein with respect to another protein. The reduced protein model allows one to perform systematic docking minimization of many thousand start structures in reasonable computer time. Flexibility of critical surface side-chains can be accounted for by a multiple conformational copy approach. The multicopy approach allows simultaneous adjustment of side-chain conformations and optimization of translational and rotational degrees of freedom of one protein with respect to the partner during docking. For 3 (Targets 8, 14, and 19) out of 5 CAPRI targets, the approach resulted in predictions in close agreement with experiment [root-mean-square deviation (RMSD) of backbone atoms within 10 A of the protein-protein interface < 1.8 A]. The comparison of predicted and experimental structures of the CAPRI targets indicates that besides local conformational changes (e.g., changes in side-chain conformations), global conformational changes of the protein backbone can be critical for complex formation. These conformational changes not accounted for during docking are a likely reason for the unrealistic predictions in 2 cases (Targets 9 and 18).

Algorithms↗

PubMed2000

Modeling the interactions of a peptide-major histocompatibility class I ligand with its receptors. I. Recognition by two alpha beta T cell receptors.

A three-dimensional model of the complex between an Influenza Hemagglutinin peptide, Ha255-262, and its restricting element, the mouse major histocompatibility complex (MHC) class I molecule, Kk, was built by homology modeling and subsequently refined by simulated annealing and restrained molecular dynamics. Next, three-dimensional models of two different T cell receptors (TCRs) both specific for the Ha255-262/Kk complex were generated based on previously published TCR X-ray structures. Finally, guided by the recently published X-ray structures of ternary TCR/peptide/MHC-I complexes, the TCR models were successfully docked into the Ha255-262/Kk model. We have previously used a systematic and exhaustive panel of 144 single amino acid substituted analogs to analyze both MHC binding and T cell recognition of the parental viral peptide. This large body of experimental data was used to evaluate the models. They were found to account well for the experimentally obtained data, lending considerable support to the proposed models and suggesting a universal docking mode for alpha beta TCRs to MHC-peptide complexes. Such models may also be useful in guiding future rational experimentation.

Amino Acid Sequence↗

PubMed2006

Solvated docking: introducing water into the modelling of biomolecular complexes.

MOTIVATION: Interfacial water, which plays an important role in mediating biomolecular interactions, has been neglected in the modelling of biomolecular complexes. METHODS: We present a solvated docking approach that explicitly accounts for the presence of water in protein-protein complexes. Our solvated docking protocol is based on the concept of the first encounter complex in which a water layer is present in-between the molecules. It mimics the pathway from this initial complex towards the final assembly in which most waters have been expelled from the interface. Docking is performed from solvated biomolecules and waters are removed in a biased Monte Carlo procedure based on water-mediated contact propensities obtained from an analysis of high-resolution crystal structures. RESULTS: We demonstrate the feasibility of this approach for protein-protein complexes representing both 'wet' and 'dry' interfaces. Solvated docking leads to improvements both in quality and scoring. Water molecules are recovered that closely match the ones in the crystal structures. AVAILABILTY: Solvated docking will be made available in the future release of HADDOCK version 2.0 (http://www.nmr.chem.uu.nl/haddock).

Binding Sites↗

PubMed2004

Coupled motions between pore and voltage-sensor domains: a model for Shaker B, a voltage-gated potassium channel.

A high-resolution crystal structure of KvAP, an archeabacterial voltage-gated potassium (Kv) channel, complexed with a monoclonal Fab fragment has been recently determined. Based on this structure, a mechanism for the activation (opening) of Kv channels has been put forward. This mechanism has since been criticized, suggesting that the resolved structure is not representative of the family of voltage-gated potassium channels. Here, we propose a model of the transmembrane domain of Shaker B, a well-characterized Kv channel, built by homology modeling and docking calculations. In this model, the positively charged S4 helices are oriented perpendicular to the membrane and localized in the groove between segments S5 and S6 of adjacent subunits. The structure and the dynamics of the full atomistic model embedded in a hydrated lipid bilayer were investigated by means of two large-scale molecular dynamics simulations under transmembrane-voltage conditions known to induce, respectively, the resting state (closed) and the activation (opening) of voltage-gated channels. Upon activation, the model undergoes conformational changes that lead to an increase of the hydration of the charged S4 helices, correlated with an upward translation and a tilting of the latter, concurrently with movements of the S5 helices and the activation gate. Although small, these conformational changes ultimately result in an alteration of the ion-conduction pathway. Our findings support the transporter model devised by Bezanilla and collaborators, and further underline the crucial role played by internal hydration in the activation of the channel.

Computer Simulation↗