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

PubMed1995

First-principles calculation of the folding free energy of a three-helix bundle protein.

The folding and unfolding of a three-helix bundle protein were explored with molecular-dynamics simulations, cluster analysis, and weighted-histogram techniques. The folding-unfolding process occurs by means of a "folding funnel," in which a uniform and broad distribution of conformational states is accessible outside of the native manifold. This distribution narrows near a transition region and becomes compact within the native manifold. Key thermodynamic steps in folding include initial interactions around the amino-terminal helix-turn-helix motif, interactions between helices I and II, and, finally, the docking of helix III onto the helix I-II subdomain. A metastable minimum in the calculated free-energy surface is observed at approximately 1.5 times the native volume. Folding-unfolding thermodynamics are dominated by the opposing influences of protein-solvent energy, which favors unfolding, and the overall entropy, which favors folding by means of the hydrophobic effect.

Amino Acid Sequence↗

PubMed1997

Molecular docking to ensembles of protein structures.

Until recently, applications of molecular docking assumed that the macromolecular receptor exists in a single, rigid conformation. However, structural studies involving different ligands bound to the same target biomolecule frequently reveal modest but significant conformational changes in the target. In this paper, two related methods for molecular docking are described that utilize information on conformational variability from ensembles of experimental receptor structures. One method combines the information into an "energy-weighted average" of the interaction energy between a ligand and each receptor structure. The other method performs the averaging on a structural level, producing a "geometry-weighted average" of the inter-molecular force field score used in DOCK 3.5. Both methods have been applied in docking small molecules to ensembles of crystal and solution structures, and we show that experimentally determined binding orientations and computed energies of known ligands can be reproduced accurately. The use of composite grids, when conformationally different protein structures are available, yields an improvement in computational speed for database searches in proportion to the number of structures.

Computer Simulation↗

PubMed2006

The dynamics of the MgATP-driven closure of MalK, the energy-transducing subunit of the maltose ABC transporter.

The nucleotide binding domains (NBDs) are the energy supplying subunits of ATP-binding cassette (ABC) proteins. They power transport by binding and hydrolyzing ATP. Tracing the pathway between different conformational states of the NBDs during ATP binding, hydrolysis, and release has, however, proven difficult. We have used molecular dynamics simulations to study the ATP-driven association of the NBDs of the maltose ABC transporter, MalK, based on the crystal structures of its open and semiopen dimers. When MgATP was introduced into the binding pockets, the semiopen dimer transitioned to a closed conformation, whereas the open dimer evolved to a semiopen state. In the absence of docked MgATP, however, the twin NBDs of both the open and semiopen starting configurations drifted further apart. Both the presence of MgATP and direct cross-interface protein-protein hydrogen bonds, primarily involving the D-loop, quite likely play a key role in initiating closure. The simulations of the MgATP-docked semiopen form indicate that completion of closure is driven mainly by cross-interface contacts between the gamma-phosphate of ATP and residues in the signature motif. Our simulations also give insight into possible interactions of MalK with the regulatory proteins MalT and enzyme IIA(glc).

ATP-Binding Cassette Transporters↗

PubMed2001

Molecular model of the Escherichia coli Na1/H1 antiporter NhaA.

A three-dimensional electron density projection map of the ion-coupled membrane protein Escherichia coli Na+/H+ antiporter (NhaA) was recently published. Based on this projection map, and previous biophysical studies determining the assignment of the 12 transmembrane alpha-helices (TMHs), a three-dimensional molecular model of the NhaA was constructed, using interactive molecular graphics and energy calculations. The diuretic drug, amiloride, was docked into the model and putative interacting amino acids were identified. The model suggests that the pH dependent activity of NhaA may be explained by charge changes in the intracellular loop between TMH8 and TMH9 which alter the positions of TMHs 4, 5 and 11 relative to each other, such that a pore area of the transporter protein is opened.

Bacterial Proteins↗

PubMed2005

Calculation of the Maxwell stress tensor and the Poisson-Boltzmann force on a solvated molecular surface using hypersingular boundary integrals.

The electrostatic interaction among molecules solvated in ionic solution is governed by the Poisson-Boltzmann equation (PBE). Here the hypersingular integral technique is used in a boundary element method (BEM) for the three-dimensional (3D) linear PBE to calculate the Maxwell stress tensor on the solvated molecular surface, and then the PB forces and torques can be obtained from the stress tensor. Compared with the variational method (also in a BEM frame) that we proposed recently, this method provides an even more efficient way to calculate the full intermolecular electrostatic interaction force, especially for macromolecular systems. Thus, it may be more suitable for the application of Brownian dynamics methods to study the dynamics of protein/protein docking as well as the assembly of large 3D architectures involving many diffusing subunits. The method has been tested on two simple cases to demonstrate its reliability and efficiency, and also compared with our previous variational method used in BEM.

Algorithms↗

PubMed2005

An NMR-based docking model for the physiological transient complex between cytochrome f and cytochrome c6.

The physiological transient complex between cytochrome f (Cf) and cytochrome c(6) (Cc(6)) from the cyanobacterium Nostoc sp. PCC 7119 has been analysed by NMR spectroscopy. The binding constant at low ionic strength is 8 +/- 2 mM(-1), and the binding site of Cc(6) for Cf is localized around its exposed haem edge. On the basis of the experimental data, the resulting docking simulations suggest that Cc(6) binds to Cf in a fashion that is analogous to that of plastocyanin but differs between prokaryotes and eukaryotes.

Cytochromes c6↗

PubMed2005

Unbinding pathways of an agonist and an antagonist from the 5-HT3 receptor.

The binding sites of 5-HT3 and other Cys-loop receptors have been extensively studied, but there are no data on the entry and exit routes of ligands for these sites. Here we have used molecular dynamics simulations to predict the pathway for agonists and antagonists exiting from the 5-HT3 receptor binding site. The data suggest that the unbinding pathway follows a tunnel at the interface of two subunits, which is approximately 8 A long and terminates approximately 20 A above the membrane. The exit routes for an agonist (5-HT) and an antagonist (granisetron) were similar, with trajectories toward the membrane and outward from the ligand binding site. 5-HT appears to form many hydrogen bonds with residues in the unbinding pathway, and experiments show that mutating these residues significantly affects function. The location of the pathway is also supported by docking studies of granisetron, which show a potential binding site for granisetron on the unbinding route. We propose that leaving the binding pocket along this tunnel places the ligands close to the membrane and prevents their immediate reentry into the binding pocket. We anticipate similar exit pathways for other members of the Cys-loop receptor family.

Binding Sites↗

PubMed1997

Computer simulation study on the conversion of spirolactone to enones over H-Y zeolite.

We report here the results of computer modeling studies and quantum chemical calculations on the spirolactone-to-enone conversion reaction over the zeolite catalysts, especially H-Y zeolite. We studied the adsorption mode of the molecules inside the supercage of H-Y and the mechanism of electron transfer between organic molecules and the framework of zeolite H-Y by density functional theory (DFT) calculations. Because the organic molecules considered in the present study are less symmetrical we docked the molecule inside the supercage of H-Y and energy minimization was then applied to these docked structures to yield representative low-energy binding sites for the molecules within the host structure. The interaction energy results show that the major interaction is between the methylene hydrogen of the molecule and the oxygen of the framework. The molecular electrostatic potential maps show that the ketonic oxygen of the reactant molecules abstract proton from Brönsted acid site. Thus the mechanism proposed by DFT calculation matches well with the experimental observations.

Computer Simulation↗

PubMed2005

A Brownian dynamics study of the effects of cytochrome f structure and deletion of its small domain in interactions with cytochrome c6 and plastocyanin in Chlamydomonas reinhardtii.

The availability of seven different structures of cytochrome f (cyt f) from Chlamydomonas reinhardtii allowed us, using Brownian dynamics simulations, to model interactions between these molecules and their redox partners, plastocyanin (PC) and cytochrome c6 (cyt c6) in the same species to study the effect of cyt f structure on its function. Our results showed that different cyt f structures, which are very similar, produced different reaction rates in interactions with PC and cyt c6. We were able to attribute this to structural differences among these molecules, particularly to a small flexible loop between A-184 and G-191 (which has some of the highest crystallographic temperature factors in all of the cyt f structures) on the cyt f small domain. We also showed that deletion of the cyt f small domain affected cyt c6 more than PC, due to their different binding positions on cyt f. One function of the small domain in cyt f may be to guide PC or cyt c6 to a uniform dock with cyt f, especially due to electrostatic interactions with K-188 and K-189 on this domain. Our results could serve as a good guide for future experimental work on these proteins to understand better the electron transfer process between them. Also, these results demonstrated the sensitivity and the power of the Brownian dynamics simulations in the study of molecular interactions.

Algorithms↗

PubMed1999

Comparison of two implementations of the incremental construction algorithm in flexible docking of thrombin inhibitors.

A set of 32 known thrombin inhibitors representing different chemical classes has been used to evaluate the performance of two implementations of incremental construction algorithms for flexible molecular docking: DOCK 4.0 and FlexX 1.5. Both docking tools are able to dock 10-35% of our test set within 2 A of their known, bound conformations using default sampling and scoring parameters. Although flexible docking with DOCK or FlexX is not able to reconstruct all native complexes, it does offer a significant improvement over rigid body docking of single, rule-based conformations, which is still often used for docking of large databases. Docking of sets of multiple conformers of each inhibitor, obtained with a novel protocol for diverse conformer generation and selection, yielded results comparable to those obtained by flexible docking. Chemical scoring, which is an empirically modified force field scoring method implemented in DOCK 4.0, outperforms both interaction energy scoring by DOCK and the Böhm scoring function used by FlexX in rigid and flexible docking of thrombin inhibitors. Our results indicate that for reliable docking of flexible ligands the selection of anchor fragments, conformational sampling and currently available scoring methods still require improvement.

Algorithms↗

PubMed2006

The dipeptide H-Trp-Glu-OH shows highly antagonistic activity against PPARgamma: bioassay with molecular modeling simulation.

The peroxisome proliferator-activated receptor gamma (PPARgamma) is an important therapeutic drug target for several conditions, including diabetes, inflammation, dyslipidemia, hypertension, and cancer. It is shown that an antagonist or partial agonist of PPARgamma has attractive potential applications in the discovery of novel antidiabetic agents that may retain efficacious insulin-sensitizing properties and minimize potential side effects. In this work, the dipeptide H-Trp-Glu-OH (G3335) was discovered to be a novel PPARgamma antagonist. Biacore 3000 results based on the surface plasmon resonance (SPR) technique showed that G3335 exhibits a highly specific binding affinity against PPARgamma (K(D) = 8.34 microM) and is able to block rosiglitazone, a potent PPARgamma agonist, in the stimulation of the interaction between the PPARgamma ligand-binding domain (LBD) and RXRalpha-LBD. Yeast two-hybrid assays demonstrated that G3335 exhibits strong antagonistic activity (IC50 = 8.67 microM) in perturbing rosiglitazone in the promotion of the PPARgamma-LBD-CBP interaction. Moreover, in transactivation assays, G3335 was further confirmed as an antagonist of PPARgamma in that G3335 could competitively bind to PPARgamma against 0.1 microM rosiglitazone to repress reporter-gene expression with an IC50 value of 31.9 muM. In addition, homology modeling and molecular-docking analyses were performed to investigate the binding mode of PPARgamma-LBD with G3335 at the atomic level. The results suggested that residues Cys285, Arg288, Ser289, and His449 in PPARgamma play vital roles in PPARgamma-LBD-G3335 binding. The significance of Cys285 for PPARgamma-LBD-G3335 interaction was further demonstrated by PPARgamma point mutation (PPARgamma-LBD-Cys285Ala). It is hoped our current work will provide a powerful approach for the discovery of PPARgamma antagonists, and that G3335 might be developed as a possible lead compound in diabetes research.

Animals↗

PubMed2006

Dimerization of the lutropin receptor: insights from computational modeling.

A computational approach based upon rigid-body docking, ad hoc filtering, and cluster analysis has been carried out to predict likely interfaces in LHR homodimers. Quaternary structure predictions emphasize the role of helices 4, 5 and 6, with prominence to helix 4, in mediating inter-monomer interactions. Intermolecular interactions essentially involve the transmembrane domains rather than the hydrophilic loops and do not implicate disulfide bridges.Collectively, molecular dynamics simulations on the isolated receptor and computational modeling of LHR homodimerization suggest that mutation-induced LHR activation favors H4-H4 contacts involving the highly conserved W491 from both the receptors monomers.

Amino Acid Sequence↗

PubMed2001

Collagen-ligand interaction in dentinal adhesion: computer visualization and analysis.

The objective of this study was to characterize the interactions of selected ligand molecules with collagen structure through computer visualization of the reacting molecules and the resulting complexes. Five ligand molecules were studied. They were 2-Hydroxyethyl methacrylate, Glutaraldehyde-HEMA adduct, Glyceryl dimethacrylate, Methacryloyloxyethyl maleate and Acryloyloxyethyl citraconate. These ligands were selected with oxygen as a common heteroatom for a reactive or functional site. Energy minimized 3-D structures of the molecules were generated by Sybyl molecular modeling software. The structures were subjected to a systematic conformational search, yielding conformations of the molecules with a common recognition site with both steric and electrostatic complementarity to appropriate receptor sites in a type I collagen molecular structure. The ligands were also docked to collagen receptor by autodock procedures and the receptor sites where docking occurred were evaluated. The energy of the molecules and their complexes with collagen was evaluated and compared. The computer visualization results reveal that steric complementarity between receptor sites in collagen and optimally configured ligands may be the basis of micromechanical bonding between collagen and the ligands. Typically, ligands docked on the cavities of collagen molecular surface and wrapped around the cavities which follow the helical turns of the collagen macromolecule. In addition, analysis of electrostatic potential features revealed electrostatic complementarity as an additional source of interaction. Hydrogen bonds between ligands and collagen molecule were detected in the complexes of several of the conformations of all the ligands. Thus computer simulation studies show that steric and electrostatic complementarity and consequent interactions form the potential basis of binding between dentin adhesive ligands and type I collagen.

Adhesiveness↗

PubMedDate not supplied

Docking ligands to vasopressin and oxytocin receptors via genetic algorithm.

The aim of the study was to computer-dock selected ligands to neurophyseal receptors in order to identify amino acid residues responsible for ligand-receptor interactions. To this aim, reliable oxytocin receptor (OTR) and arginine-vasopressin receptor (V1aR/V2R) models were built. The OTR-selective agonist [Thr4,Gly7]OT, the OTR-selective cyclohexapeptide antagonist L-366,948 and OT itself were docked via genetic algorithm to OTR, V1aR, and V2R and relaxed using a constrained simulated annealing protocol. For the analysis of receptor/ligand interactions a subset of initial conformations was chosen using energetic and steric criteria. All three ligands seem to prefer similar modes of binding to the receptors, manifested by repetitive residues of the receptors which directly interact with the ligands. Taking into account that many aspects of mechanisms of G protein-coupled receptor (GPCR) action are still unsolved, the results obtained with the docking simulations may propose future experimental research, especially in site-directed mutagenesis analysis and searching for key amino acid residues responsible for drug activities.

Algorithms↗

PubMed1999

Inhibition of tumor necrosis factor-alpha (TNF-alpha)/TNF-alpha receptor binding by structural analogues of suramin.

Suramin, a symmetrical polysulfonated urea derivative, promotes the dissociation of trimeric human tumor necrosis factor-alpha (TNF-alpha) into biologically inactive subunits and prevents the interaction of TNF-alpha with its cellular receptors. The aim of this work was to identify compounds structurally related to suramin which inhibit the binding of TNF-alpha to its receptor. Molecular modeling studies were performed on suramin and TNF-alpha molecules and likely interaction sites were identified in the docked complex. On this basis, Evans blue, trypan blue, sulfonazo III, beryllon II, and 1,3,6-naphthalenetrisulfonic acid trisodium salt were identified as polysulfonated compounds endowed, to various extents, with the structural characteristics responsible for interaction with TNF-alpha. N,N-bis(3,5-di-tert-butylphenyl)-3,4,9,10-perylenedicarboximide was used as an unrelated structure. The capacity of these molecules to inhibit the binding of TNF-alpha with its receptor p55 was tested in vitro by means of a specific immunoenzymatic assay using suramin as reference compound. Evans blue and trypan blue inhibited TNF-alpha/p55 binding with an IC50 of 0.75 and 1.00 mM, respectively (suramin IC50: 0.65 mM); no effect was observed with the other molecules. Molecular modeling analyses on Evans blue and trypan blue docked into the TNF-alpha molecule support these experimental results by demonstrating that these compounds share with suramin a similar binding mode to TNF-alpha. The results of this work provide a new insight into and useful hints for the design of new chemical entities endowed with a potent and selective activity on TNF-alpha.

Antigens, CD↗

PubMed2002

Protein--protein docking of electron transfer complexes: cytochrome c oxidase and cytochrome c.

Electron transferring protein complexes form only transiently and the crystal structures of electron transfer protein--protein complexes involving cytochrome c could so far be determined only for the pairs of yeast cytochrome c peroxidase (CcP) with iso-1-cytochrome c (iso-1-cyt c) and with horse heart cytochrome c (cyt c). This article presents models from computational docking for complexes of cytochrome c oxidase (COX) from Paracoccus denitrificans with horse heart cytochrome c, and with its physiological counterpart cytochrome c552 (c552). Initial docking is performed with the FTDOCK program, which permits an exhaustive search of translational and rotational space. A filtering procedure is then applied to reduce the number of complexes to a manageable number. In a final step of structural and energetic refinement, the complexes are optimized by rigid-body energy minimization with the molecular mechanics package CHARMM. This methodology was first tested on the CcP:iso-1-cyt c complex, in which the complex with the lowest CHARMM energy has an RMSD from the crystal structure of only 1.8 A (C(alpha) carbon atoms). Notably, the crystal conformation has an even lower energy. The same procedure was then applied to COX:cyt c and COX:c552. The lowest-energy COX:cyt c complex is very similar to a docking model previously described for the complex of bovine cytochrome c oxidase with horse heart cytochrome c. For the COX:c552 complex, cytochrome c552 is found in two different orientations, depending on whether it is docked against COX from a two-subunit or from a four-subunit crystal structure, respectively. Both conformations are discussed critically in the light of the available experimental data.

Algorithms↗

PubMed2004

Altered inactivation pathway of factor Va by activated protein C in the presence of heparin.

Inactivation of factor Va (FVa) by activated protein C (APC) is a predominant mechanism in the down-regulation of thrombin generation. In normal FVa, APC-mediated inactivation occurs after cleavage at Arg306 (with corresponding rate constant k'306) or after cleavage at Arg506 (k506) and subsequent cleavage at Arg306 (k306). We have studied the influence of heparin on APC-catalyzed FVa inactivation by kinetic analysis of the time courses of inactivation. Peptide bond cleavage was identified by Western blotting using FV-specific antibodies. In normal FVa, unfractionated heparin (UFH) was found to inhibit cleavage at Arg506 in a dose-dependent manner. Maximal inhibition of k506 by UFH was 12-fold, with the secondary cleavage at Arg306 (k306) being virtually unaffected. In contrast, UFH stimulated the initial cleavage at Arg306 (k'306) two- to threefold. Low molecular weight heparin (Fragmin) had the same effects on the rate constants of FVa inactivation as UFH, but pentasaccharide did not inhibit FVa inactivation. Analysis of these data in the context of the 3D structures of APC and FVa and of simulated APC-heparin and FVa-APC complexes suggests that the heparin-binding loops 37 and 70 in APC complement electronegative areas surrounding the Arg506 site, with additional contributions from APC loop 148. Fewer contacts are observed between APC and the region around the Arg306 site in FVa. The modeling and experimental data suggest that heparin, when bound to APC, prevents optimal docking of APC at Arg506 and promotes association between FVa and APC at position Arg306.

Blotting, Western↗

PubMed2005

Theoretically predicted structures of plasma membrane Ca(2+)-ATPase and their susceptibilities to oxidation.

Oxidative damage to the plasma membrane Ca(2+)-ATPase (PMCA) appears to contribute to the decreased clearance of intracellular Ca(2+) in the neurons of aged brain, possibly contributing to its vulnerability to numerous age-related diseases such as Alzheimer's disease. The precise sites of oxidative susceptibility have not been identified. However, it is known that calmodulin (CaM) protects the purified PMCA against oxidative inactivation, perhaps via conformational restructuring of the protein through dissociation of a 20 residue domain (C20W) in the C-terminal region that function as a CaM-binding site. In order to postulate likely oxidation sites and the mechanism underlying the protection offered by CaM, we have generated a three-dimensional model of PMCA via a combination of homology/comparative modeling, threading, protein-protein docking, and guidance from prior biochemical and analytical studies. The resulting model was validated based on surface polarity/hydrophobicity profiling, standard ProCheck, WhatIF, and PROVE checks, as well as comparison with empirical structure-function observations. This model was then used to identify likely oxidation sites by comparing time-averaged solvent accessibility of potentially oxidizable surface residues as measured from molecular dynamics simulations of intact PMCA and the PMCA sequence from which C20W has been deleted. The resulting model complex has permitted us to identify three amino acids whose solvent accessibility is greatly reduced by the C20W dissociation: Tyr 589, Met 622, and Met 831.

Adenosine Triphosphate↗