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 685 records · Page 38Linked to original sources

PubMed2002

Molecular simulation of interaction between estrogen receptor and selective estrogen receptor modulators.

AIM: To study the mechanism of interaction between a series of potent racemic selective estrogen receptor modulators (SERM) and estrogen receptors (ER). METHODS: Active conformations of these conformationally restricted raloxifene analogues in binding pocket were determined by molecular mechanics. The interactive energies between ligand and receptor were calculated by docking program. RESULTS: Both R and S configurations of these SERM were accommodated by the binding pocket of ER. The hydroxy group of compounds forms hydrogen bonds with amino acid residues of ER and the phenolic group mimics the A ring of estradiol. The most potential compounds were those with two hydroxy groups and accommodated by binding pocket in S configuration with phenolic group at C(16) imitating A ring of estradiol. CONCLUSION: Chiral center conferred little effect on the binding affinity of these conformationally restricted raloxifene analogues. The hydroxy group(s) play(s) a critical role to the orientation of compounds in active pocket of ER and the binding between ligand and receptor.

Binding Sites↗

PubMed1999

The role of waters in docking strategies with incremental flexibility for carbohydrate derivatives: heat-labile enterotoxin, a multivalent test case.

Molecular docking studies of carbohydrate derivatives in protein binding sites are often challenging because of water-mediated interactions and the inherent flexibility of the many terminal hydroxyl groups. Using the recognition process between heat-labile enterotoxin from Escherichia coli and ganglioside GM1 as a paradigm, we developed a modeling protocol that includes incremental conformational flexibility of the ligand and predicted water interactions. The strategy employs a modified version of the Monte Carlo docking program AUTODOCK and water affinity potentials calculated with GRID. After calibration of the protocol on the basis of the known binding modes of galactose and lactose to the toxin, blind predictions were made for the binding modes of four galactose derivatives: lactulose, melibionic acid, thiodigalactoside, and m-nitrophenyl-alpha-galactoside. Subsequent crystal structure determinations have demonstrated that our docking strategy can predict the correct binding modes of carbohydrate derivatives within 1.0 A from experiment. In addition, it is shown that repeating the docking simulations in each of the seemingly identical binding sites of the multivalent toxin increases the chance of finding the correct binding mode.

Binding Sites↗

PubMed2005

The role and application of in silico docking in chemical genomics research.

In silico docking techniques are being used to investigate the complementarity at the molecular level of a ligand and a protein target. As such, docking studies can be used to identify the structural features that are important for binding and for in silico screening efforts in which suitable binding partners can be identified. Here we describe a practical approach for setting up docking simulations using different docking programs. We also cover the analysis and rescoring of the obtained docking poses. Possible pitfalls in the docking studies are discussed and hints are provided to resolve commonly occurring problems.

Binding Sites↗

PubMed1984

Solvent accessible surface area and excluded volume in proteins. Analytical equations for overlapping spheres and implications for the hydrophobic effect.

An analytical formula has been derived for the calculation of the solvent accessible surface area of a protein molecule or equivalently the surface area exterior to an arbitrary number of overlapping spheres. The directional derivative of this function with respect to atomic co-ordinates is provided to facilitate minimization procedures used with molecular docking algorithms and energy calculations. An analytical formula for the calculation of the volume enclosed within the accessible surface, the excluded volume, is also derived. Although the area function is not specific to the structures of proteins, the derivation was motivated by the need for a computationally feasible simulation of the hydrophobic effect in proteins. A computer program using the equations for area has been tested and has had limited application to the docking of protein alpha-helices. Possible relationships of the solvent excluded volume to hydrophobic interaction free energy and transfer free energy of solute molecules are derived from the statistical mechanics of solution.

Mathematics↗

PubMed2006

ADP_EM: fast exhaustive multi-resolution docking for high-throughput coverage.

MOTIVATION: Efficient fitting tools are needed to take advantage of a fast growth of atomic models of protein domains from crystallography or comparative modeling, and low-resolution density maps of larger molecular assemblies. Here, we report a novel fitting algorithm for the exhaustive and fast overlay of partial high-resolution models into a low-resolution density map. The method incorporates a fast rotational search based on spherical harmonics (SH) combined with a simple translational scanning. RESULTS: This novel combination makes it possible to accurately dock atomic structures into low-resolution electron-density maps in times ranging from seconds to a few minutes. The high-efficiency achieved with simulated and experimental test cases preserves the exhaustiveness needed in these heterogeneous-resolution merging tools. The results demonstrate its efficiency, robustness and high-throughput coverage. AVAILABILITY: http://sbg.cib.csic.es/Software/ADP_EM. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Algorithms↗

PubMed2006

Empirical rules facilitate the search for binding sites on protein surfaces.

Computational surface screening of 3D protein structures is a valuable means of finding possible docking sites for substrates, effectors and similar molecules. It can be improved by considering properties of molecules which are known to bind to protein surfaces, and thus reflect the required properties of binding sites. In-depth studies are available on drugs and lead compounds as binding partners with statistically assured properties. Here we present a simple strategy for finding binding sites, which is based on the empirical rule-of-five by Lipinski et al. for oral drugs and the rule-of-three by Congreve et al. for leads. The fast automated search with the new C-code TRIDOCK yields a preliminary set of sites, thus facilitating further investigation by visual, comparative and quantitative work. Possible binding sites are tagged by pseudo-atoms added to the structure file for molecular graphical evaluation. Usually, the strategy yields not just a few single sites, but an accumulation of several sites in known substrate binding pockets. Clusters are also found at known or putative protein-protein docking interfaces. A comparison of the activated and inactivated form of the GTPase Ras reveals clear differences and identifies a niche, which is possibly a suitable new target for compounds that bind specifically to activated Ras.

Bacteriorhodopsins↗

PubMed1987

Comparison of computer modelling and X-ray results of the binding of a pyrazole derivative to liver alcohol dehydrogenase.

The binding to liver alcohol dehydrogenase of the inhibitor 2,4-(4-pyrazolyl)-butylisothiourea has been studied both by modelling experiments using computer graphics with interactive energy minimization and by X-ray crystallographic structure determination. For the modelling experiments, we used the program system TOM, which was developed in our laboratory as an extension of the program FRODO. Different strategies for using computer graphics with interactive energy minimization were tested. Two essentially different binding modes were found. One of these was favoured from energy minimizations using a potential energy function which was the sum of a Coulomb interaction term and two different van der Waals' interaction terms for non-bonded and torsional interactions. This binding mode was close to the crystallographic observed structure. The results show that flexibility of both ligand and receptor side-chains as well as main-chain conformations are important for docking to the active site of liver alcohol dehydrogenase.

Alcohol Dehydrogenase↗

PubMed2002

Computer modeling of three-dimensional structure of DNA-packaging RNA (pRNA) monomer, dimer, and hexamer of Phi29 DNA packaging motor.

A striking common feature in the maturation of all linear double-stranded DNA viruses is that their lengthy genome is translocated with remarkable velocity into the limited space within a preformed protein shell and packaged into near crystalline density. A DNA-translocating motor, powered by ATP hydrolysis, accomplishes this task, which would otherwise be energetically unfavorable. DNA-packaging RNA, pRNA, forms a hexameric complex to serve as a vital component of the DNA translocating motor of bacterial virus Phi29. The sequential action of six pRNA ensures continual function in the DNA translocation process. The Phi29 motor has been assembled with purified components synthesized by chemical or biotechnological approaches and is able to pump the viral DNA into the protein shell in vitro. pRNA dimers are the building blocks of the hexamer. The computer models of the three-dimensional structure of the motor was constructed based on experimental data derived from photoaffinity cross-linking by psoralen, phenphi (cis-Rh(1,10-phenanthroline)(9,10-phenan-threnequinone diimine)Cl(2)(+)), and azidophenacyl; chemical modification and chemical modification interference with dimethyl sulfate, 1-cyclohexyl-3-(2-morpholinoethyl)carbodiimide metho-p-toluene sulfonate, and kethoxal; complementary modification; and nuclease probing by single- and double-stranded specific RNases. The shapes of these computer models are very similar to the published pRNA images of cryo-atomic force microscopy. pRNA hexamer docking with the connector crystal structure reveals a very impressive match with the available biochemical, genetic, and physical data.

Bacteriophages↗

PubMed2006

Structural basis for tetraspanin functions as revealed by the cryo-EM structure of uroplakin complexes at 6-A resolution.

Tetraspanin uroplakins (UPs) Ia and Ib, together with their single-spanning transmembrane protein partners UP II and IIIa, form a unique crystalline 2D array of 16-nm particles covering almost the entire urothelial surface. A 6 A-resolution cryo-EM structure of the UP particle revealed that the UP tetraspanins have a rod-shaped structure consisting of four closely packed transmembrane helices that extend into the extracellular loops, capped by a disulfide-stabilized head domain. The UP tetraspanins form the primary complexes with their partners through tight interactions of the transmembrane domains as well as the extracellular domains, so that the head domains of their tall partners can bridge each other at the top of the heterotetramer. The secondary interactions between the primary complexes and the tertiary interaction between the 16-nm particles contribute to the formation of the UP tetraspanin network. The rod-shaped tetraspanin structure allows it to serve as stable pilings in the lipid sea, ideal for docking partner proteins to form structural/signaling networks.

Animals↗

PubMed2006

CIRSE: a solvation energy estimator compatible with flexible protein docking and design applications.

We present the Coordinate Internal Representation of Solvation Energy (CIRSE) for computing the solvation energy of protein configurations in terms of pairwise interactions between their atoms with analytic derivatives. Currently, CIRSE is trained to a Poisson/surface-area benchmark, but CIRSE is not meant to fit this benchmark exclusively. CIRSE predicts the overall solvation energy of protein structures from 331 NMR ensembles with 0.951+/-0.047 correlation and predicts relative solvation energy changes between members of individual ensembles with an accuracy of 15.8+/-9.6 kcal/mol. The energy of individual atoms in any of CIRSE's 17 types is predicted with at least 0.98 correlation. We apply the model in energy minimization, rotamer optimization, protein design, and protein docking applications. The CIRSE model shows some propensity to accumulate errors in energy minimization as well as rotamer optimization, but these errors are consistent enough that CIRSE correctly identifies the relative solvation energies of designed sequences as well as putative docked complexes. We analyze the errors accumulated by the CIRSE model during each type of simulation and suggest means of improving the model to be generally useful for all-atom simulations.

Models, Molecular↗

PubMed2005

Refinement of protein cores and protein-peptide interfaces using a potential scaling approach.

Refinement of side chain conformations in protein model structures and at the interface of predicted protein-protein or protein-peptide complexes is an important step during protein structural modelling and docking. A common approach for side chain prediction is to assume a rigid protein main chain for both docking partners and search for an optimal set of side chain rotamers to optimize the steric fit. However, depending on the target-template similarity in the case of comparative protein modelling and on the accuracy of an initially docked complex, the main chain template structure is only an approximation of a realistic target main chain. An inaccurate rigid main chain conformation can in turn interfere with the prediction of side chain conformations. In the present study, a potential scaling approach (PS-MD) during a molecular dynamics (MD) simulation that also allows the inclusion of explicit solvent has been used to predict side chain conformations on semi-flexible protein main chains. The PS-MD method converges much faster to realistic protein-peptide interface structures or protein core structures than standard MD simulations. Depending on the accuracy of the protein main chain, it also gives significantly better results compared with the standard rotamer search method.

Amino Acid Sequence↗

PubMed2005

Novel, highly potent adenosine deaminase inhibitors containing the pyrazolo[3,4-d]pyrimidine ring system. Synthesis, structure-activity relationships, and molecular modeling studies.

This study reports the synthesis of a number of 1- and 2-alkyl derivatives of the 4-aminopyrazolo[3,4-d]pyrimidine (APP) nucleus and their evaluation as inhibitors of ADA from bovine spleen. The 2-substituted aminopyrazolopyrimidines proved to be potent inhibitors, most of them exhibiting K(i) values in the nanomolar/subnanomolar range. In this series the inhibitory activity is enhanced with the increase in length of the alkyl chain, reaching a maximum with the n-decyl substituent. Insertion of a 2'-hydroxy group in the n-decyl chain gave 3k, whose (R)-isomer displayed the highest inhibitory potency of the series (K(i) 0.053 nM), showing an activity 2 orders of magnitude higher than that of (+)-EHNA (K(i) 1.14 nM), which was taken as the reference standard. Docking simulations of aminopyrazolopyrimidines into the ADA binding site were also performed, to rationalize the structure-activity relationships of this class of inhibitors.

Adenosine Deaminase Inhibitors↗

PubMedDate not supplied

A thermostable sugar-binding protein from the Archaeon Pyrococcus horikoshii as a probe for the development of a stable fluorescence biosensor for diabetic patients.

In this work is presented the first attempt to develop an innovative ultrastable protein-based biosensor for blood glucose detections. The gene of a putative thermostable sugar-binding protein has been cloned and expressed in E. coli. The recombinant protein has been purified to homogeneity by thermoprecipitation and affinity chromatography steps. The recombinant protein is a monomer with an apparent molecular weight of 55,000 as judged by gel filtration and sodium dodecyl sulfate polyacrylamide gel eletrophoresis. Circular dichroism experiments showed that the protein possesses a secondary structure content rich in alpha-helices and beta-structures and that the protein is highly stable as investigated in the range of temperature between 20 and 95 degrees C. Fluorescence spectroscopy experiments demonstrated that the recombinant protein binds glucose with a dissociation constant of about 10 mM, a concentration of sugar very close to the concentration of glucose present in the human blood. A docking simulation on the modeled structure of the protein confirms its ability to bind glucose and proposes possible modifications to improve the affinity for glucose and/or its detection. The obtained results suggest the use of the protein as a probe for a stable glucose biosensor.

Archaea↗

PubMed2000

Modeling of kappa-opioid receptor/agonists interactions using pharmacophore-based and docking simulations.

The interaction of the kappa-opioid receptor with arylacetamide and benzomorphan derivatives acting as agonists was modeled through pharmacophore-based and docking calculations. Potentially bioactive conformations of representative ligands (U-50,488 and its benzo-fused analogues 4 and 6 for arylacetamides and MPCB for benzomorphans) were identified by systematic conformational analysis and docked into a 3D model of the kappa-receptor. The obtained complexes, refined by energy-minimization and molecular dynamics, were evaluated for their consistency with structure-activity relationships and site-directed mutagenesis data. The following interactions are hypothesized to govern the ligand-receptor recognition process: (i) a salt bridge between the Asp138 carboxylate and the protonated nitrogen of the bound agonist; (ii) a hydrogen bond donated by the Tyr312 hydroxyl to the carbonyl oxygen of arylacetamides and MPCB; (iii) hydrophobic interactions established by the dichlorophenyl moiety of arylacetamides and the pendant phenyl ring of MPCB with the surrounding side chains of Tyr312, Leu224, Leu295, and Ala298; (iv) a pi-stacking contact between the Tyr312 side chain and the phenyl ring of arylacetamides; (v) a hydrogen bond linking the His291 imidazole ring to the phenolic hydroxy group featured by typical benzomorphans and the arylacetamides 4 and 6.

3,4-Dichloro-N-methyl-N-(2-(1-pyrrolidinyl)-cycloh↗

PubMed2005

Improving CAPRI predictions: optimized desolvation for rigid-body docking.

The ICM Docking and Interface Side-Chain Optimization (ICM-DISCO) showed promising predictive results during the first CAPRI experiment by successfully finding medium- or high-accuracy models in 3 of the 7 targets. A key factor was the ability to recognize near-native rigid-body geometries in a relatively low number of alternative docking poses, together with the successful refinement of the rigid-body docking interfaces. Since then, we have focused on improving the scoring function to optimally discriminate the near-native rigid-body conformations. For that, we have defined a new desolvation descriptor for rigid-body docking, based on atomic solvation parameters (ASPs) derived from octanol-water transfer experiments. This and other new approaches have been gradually incorporated into our docking procedure during our participation on the second CAPRI experiment. Overall, we produced reasonable models for 8 of the 9 official targets. Especially encouraging were those cases in which a homology model of 1 of the subunits had to be used during the docking simulations. And not less gratifying has been the successful prediction of antibody-antigen targets in a completely automatic, unrestrained fashion. In summary, our success rate (89%) shows a consistent improvement over the previous CAPRI rounds, and suggests that a correct desolvation description is key for improved protein-protein docking predictions.

Algorithms↗

PubMed1989

A computational approach to the mechanism of self-cleavage of hammerhead RNA.

Extensive minimization and dynamics computational studies of the hammerhead structural domain of the virusoid of lucerne transient streak virus have been carried out. The following observations at the self-cleavage position are derived from the resulting three-dimensional structure: (i) the cytosine base is at the surface and does not interact with another base (it is free to move), and (ii) the ribose-phosphate backbone is forced to take an abrupt turn since it bridges stems I and III, and this turn points the pro-S and pro-R oxygens of the phosphate to the inward side of the hammerhead. These structural features are independent of the hammerhead being open or closed and allow an unencumbered 3'- to 2'-endo conformational change of the ribose with the resultant creation of an unusual stereochemistry that allows a direct in-line self-cleavage reaction. In the closed hammerhead structure, interactions of stems I and II create a vacancy into which the catalytic hydrated Mg(II) may be docked on labile phosphate. This opening is not present if stems I and II are shortened.

Base Sequence↗

PubMed1994

Monte Carlo docking of protein-DNA complexes: incorporation of DNA flexibility and experimental data.

A Monte Carlo simulation program (MONTY) has been developed to dock proteins onto DNA. Protein and DNA interact via square-well potentials for hydrogen bond and van der Waals interactions. The effect of the inclusion of DNA flexibility and experimentally derived restraints has been tested on members of the helix-turn-helix family of DNA binding proteins. Unwinding and bending the DNA double helix improves the number of correctly retrieved hydrogen bonds in simulations starting from the 434 cro protein monomer complexed with a standard B-DNA OR1 half-site. Agreement with phosphate ethylation interference and mutagenesis data is rewarded with energy bonuses. This protocol was tested on protein-DNA complexes of 434 cro, lac headpiece and a mutant lac headpiece resembling the gal repressor headpiece with the recognition helices in correct and reversed orientations in the DNA major groove. The inclusion of experimental data gives an improved convergence of the correctly oriented structures and allows for an easier discrimination between correctly and incorrectly docked complexes.

Algorithms↗

PubMed2004

A Brownian dynamics study of the interaction of Phormidium laminosum plastocyanin with Phormidium laminosum cytochrome f.

The interaction of Phormidium laminosum plastocyanin (PC) with P. laminosum cytochrome f (cyt f) was studied using Brownian dynamics (BD) simulations. Few complexes and a low rate of electron transfer were observed for wild-type PC. Increasing the positive electrostatic field on PC by the addition of a Zn(2+) ion in the neighborhood of D44 and D45 on PC (as found in crystal structure of plastocyanin) increased the number of complexes formed and the calculated rates of electron transfer as did PC mutations D44A, D45A, E54A, and E57A. Mutations of charged residues on Phormidium PC and Phormidium cyt f were used to map binding sites on both proteins. In both the presence and absence of the Zn(2+) ion, the following residues on PC interact with cyt f: D44, D45, K6, D79, R93, and K100 that lie in a patch just below H92 and Y88 and D10, E17, and E70 located on the upper portion of the PC molecule. In the absence of the Zn(2+) ion, K6 and K35 on the top of the PC molecule also interact with cyt f. Cyt f residues involved in binding PC, in the absence of the Zn(2+) ion, include E165, D187, and D188 that are located on the small domain of cyt f. The orientation of PC in the complexes was quite random in accordance with NMR results. In the presence of the Zn(2+) ion, K53 and E54 in the lower patch of the PC molecule also interact with cyt f and PC interacts with E86, E95, and E123 on the large domain of cyt f. Also, the orientation of PC in the complexes was much more uniform than in the absence of the Zn(2+) ion. The difference may be due to both the larger electrostatic field and the greater asymmetry of the charge distribution on PC observed in the presence of the Zn(2+) ion. Hydrophobic interactions were also observed suggesting a model of cyt f-PC interactions in which electrostatic forces bring the two molecules together but hydrophobic interactions participate in stabilizing the final electron-transfer-active dock.

Binding Sites↗