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 955 records · Page 53Linked to original sources

PubMed2003

Prediction of 5-HT3 receptor agonist-binding residues using homology modeling.

5-HT(3) receptors demonstrate significant structural and functional homology to other members of the Cys-loop ligand-gated ion channel superfamily. The extracellular domains of these receptors share similar sequence homology (approximately 20%) with Limnaea acetylcholine binding protein, for which an x-ray crystal structure is available. We used this structure as a template for computer-based homology modeling of the 5-HT(3) receptor extracellular domain. AutoDock software was used to dock 5-HT into the putative 5-HT(3) receptor ligand-binding site, resulting in seven alternative energetically favorable models. Residues located no more than 5 A from the docked 5-HT were identified for each model; of these, 12 were found to be common to all seven models with five others present in only certain models. Some docking models reflected the cation-pi interaction previously demonstrated for W183, and data from these and other studies were used to define our preferred models.

Amino Acid Sequence↗

PubMed2005

Crystal structure of a hyperthermophilic archaeal acylphosphatase from Pyrococcus horikoshii--structural insights into enzymatic catalysis, thermostability, and dimerization.

Acylphosphatases catalyze the hydrolysis of the carboxyl-phosphate bond in acyl phosphates. Although acylphosphatase-like sequences are found in all three domains of life, no structure of acylphosphatase has been reported for bacteria and archaea so far. Here, we report the characterization of enzymatic activities and crystal structure of an archaeal acylphosphatase. A putative acylphosphatase gene (PhAcP) was cloned from the genomic DNA of Pyrococcus horikoshii and was expressed in Escherichia coli. Enzymatic parameters of the recombinant PhAcP were measured using benzoyl phosphate as the substrate. Our data suggest that, while PhAcP is less efficient than other mammalian homologues at 25 degrees C, the thermophilic enzyme is fully active at the optimal growth temperature (98 degrees C) of P. horikoshii. PhAcP is extremely stable; its apparent melting temperature was 111.5 degrees C and free energy of unfolding at 25 degrees C was 54 kJ mol(-)(1). The 1.5 A crystal structure of PhAcP adopts an alpha/beta sandwich fold that is common to other acylphosphatases. PhAcP forms a dimer in the crystal structure via antiparallel association of strand 4. Structural comparison to mesophilic acylphosphatases reveals significant differences in the conformation of the L5 loop connecting strands 4 and 5. The extreme thermostability of PhAcP can be attributed to an extensive ion-pair network consisting of 13 charge residues on the beta sheet of the protein. The reduced catalytic efficiency of PhAcP at 25 degrees C may be due to a less flexible active-site residue, Arg20, which forms a salt bridge to the C-terminal carboxyl group. New insights into catalysis were gained by docking acetyl phosphate to the active site of PhAcP.

Acid Anhydride Hydrolases↗

PubMedDate not supplied

Fully automated flexible docking of ligands into flexible synthetic receptors using forward and inverse docking strategies.

The prediction of the structure of host-guest complexes is one of the most challenging problems in supramolecular chemistry. Usual procedures for docking of ligands into receptors do not take full conformational freedom of the host molecule into account. We describe and apply a new docking approach which performs a conformational sampling of the host and then sequentially docks the ligand into all receptor conformers using the incremental construction technique of the FlexX software platform. The applicability of this approach is validated on a set of host-guest complexes with known crystal structure. Moreover, we demonstrate that due to the interchangeability of the roles of host and guest, the docking process can be inverted. In this inverse docking mode, the receptor molecule is docked around its ligand. For all investigated test cases, the predicted structures are in good agreement with the experiment for both normal (forward) and inverse docking. Since the ligand is often smaller than the receptor and, thus, its conformational space is more restricted, the inverse docking approach leads in most cases to considerable speed-up. By having the choice between two alternative docking directions, the application range of the method is significantly extended. Finally, an important result of this study is the suitability of the simple energy function used here for structure prediction of complexes in organic media.

Algorithms↗

PubMed1991

Modeling and conformation analysis of beta-cyclodextrin complexes.

A series of beta-cyclodextrin complexes containing various guest molecules was studied using computer-aided molecular modeling and conformation analysis techniques. The geometry of each complex was studied using crystallographic data. The positions of the glycosidic O4 atoms indicate that the beta-cyclodextrin molecules are elliptically distorted. This distortion can be related to the van der Waals volume of the guest molecules. This correlation is different for aromatic and non-aromatic guest compounds. Rigid body docking experiments demonstrated that in crystal structures the guest molecule occupies a position in the cavity of nearly minimum interaction energy when there are no other molecules having interactions with the guest molecule. From the crystallographic data several rules could be deduced which seem to determine the conformation of beta-cyclodextrin molecules in complexes. A procedure was developed to construct beta-cyclodextrin molecules that are able to encompass guest molecules having a given van der Waals volume.

Carbohydrate Conformation↗

PubMed1997

The cannabinoid receptor: computer-aided molecular modeling and docking of ligand.

A three-dimensional model of human cannabinoid receptor is constructed using computer-aided molecular modeling techniques. The helices of bacteriorhodopsin were used as the initial template to construct the transmembrane helices. The extracellular and intracellular loops were added using the SYBYL molecular modeling package. The extracellular N terminus was modeled on the basis of its similarity to rat oncomodulin. Similarly, the C terminus was constructed on the basis of similarity to bovine prothrombin fragment 1. The final structure was refined by several runs of minimization and dynamics calculation using the CHARMm package. delta 9-Tetra hydrocannabinol was docked into the internal cavity using the AUTODOCK program. Our study snows that there may be a calcium-binding site in the extracellular N terminus of this receptor. The ligand binds mainly to a hydrophobic site, which consists of residues Met-240, Trp-241 (TMH-4), Trp-356, Leu-359, Leu-360 (TMH-6), and Ala-283 (TMH-5). Its phenolic hydroxyl group forms a hydrogen bond with the carboxy group of Ala-198 (TMH-3). The results of modeling agree well with experimental QSAR studies.

Amino Acid Sequence↗

PubMed2005

Finding biologically relevant protein domain interactions: conserved binding mode analysis.

Proteins evolved through the shuffling of functional domains, and therefore, the same domain can be found in different proteins and species. Interactions between such conserved domains often involve specific, well-determined binding surfaces reflecting their important biological role in a cell. To find biologically relevant interactions we developed a method of systematically comparing and classifying protein domain interactions from the structural data. As a result, a set of conserved binding modes (CBMs) was created using the atomic detail of structure alignment data and the protein domain classification of the Conserved Domain Database. A conserved binding mode is inferred when different members of interacting domain families dock in the same way, such that their structural complexes superimpose well. Such domain interactions with recurring structural themes have greater significance to be biologically relevant, unlike spurious crystal packing interactions. Consequently, this study gives lower and upper bounds on the number of different types of interacting domain pairs in the structure database on the order of 1000-2000. We use CBMs to create domain interaction networks, which highlight functionally significant connections by avoiding many infrequent links between highly connected nodes. The CBMs also constitute a library of docking templates that may be used in molecular modeling to infer the characteristics of an unknown binding surface, just as conserved domains may be used to infer the structure of an unknown protein. The method's ability to sort through and classify large numbers of putative interacting domain pairs is demonstrated on the oligomeric interactions of globins.

Algorithms↗

PubMed2005

Tetrazolyl isoxazole amino acids as ionotropic glutamate receptor antagonists: synthesis, modelling and molecular pharmacology.

Two 3-(5-tetrazolylmethoxy) analogues, 1a and 1b, of (RS)-2-amino-3-(3-hydroxy-5-methyl-4-isoxazolyl)propionic acid (AMPA), a selective AMPA receptor agonist, and (RS)-2-amino-3-(5-tert-butyl-3-hydroxy-4-isoxazolyl)propionic acid (ATPA), a GluR5-preferring agonist, were synthesized. Compounds 1a and 1b were pharmacologically characterized in receptor binding assays, and electrophysiologically on homomeric AMPA receptors (GluR1-4), homomeric (GluR5 and GluR6) and heteromeric (GluR6/KA2) kainic acid receptors, using two-electrode voltage-clamped Xenopus laevis oocytes expressing these receptors. Both analogues proved to be antagonists at all AMPA receptor subtypes, showing potencies (Kb=38-161 microM) similar to that of the AMPA receptor antagonist (RS)-2-amino-3-[3-(carboxymethoxy)-5-methyl-4-isoxazolyl]propionic acid (AMOA) (Kb=43-76 microM). Furthermore, the AMOA analogue, 1a, blocked two kainic acid receptor subtypes (GluR5 and GluR6/KA2), showing sevenfold preference for GluR6/KA2 (Kb=19 microM). Unlike the iGluR antagonist (S)-2-amino-3-[5-tert-butyl-3-(phosphonomethoxy)-4-isoxazolyl]propionic acid [(S)-ATPO], the corresponding tetrazolyl analogue, 1b, lacks kainic acid receptor effects. On the basis of docking to a crystal structure of the isolated extracellular ligand-binding core of the AMPA receptor subunit GluR2 and a homology model of the kainic acid receptor subunit GluR5, we were able to rationalize the observed structure-activity relationships.

Animals↗

PubMed2004

Accurate detection of protein:ligand binding sites using molecular dynamics simulations.

Accurate prediction of location of cavities and surface grooves in proteins is important, as these are potential sites for ligand binding. Several currently available programs for cavity detection are unable to detect cavities near the surface or surface grooves. In the present study, an optimized molecular dynamics based procedure is described for detection and quantification of interior cavities as well as surface pockets. This is based on the observation that the mobility of water in such pockets is significantly lower than that of bulk water. The algorithm efficiently detects surface grooves that are sites of protein-ligand and protein-protein interaction. The algorithm was also used to substantially improve the performance of an automated docking procedure for docking monomers of nonobligate protein-protein complexes. In addition, it was applied to predict key residues involved in the binding of the E. coli toxin CcdB with its inhibitor. Predictions were subsequently validated by mutagenesis experiments.

Algorithms↗

PubMed1998

The surface-charge asymmetry and dimerisation of cytochrome c550 from Paracoccus denitrificans--implications for the interaction with cytochrome c peroxidase.

The implications of the dimeric state of cytochrome c550 for its binding to Paracoccus cytochrome c peroxidase and its delivery of the two electrons required to restore the active enzyme during catalysis have been investigated. The amino acid sequence of cytochrome c550 of Paracoccus denitrificans strain LMD 52.44 was determined and showed 21 differences from that of strain LMD 22.21. Based on the X-ray structure of the latter, a structure for the cytochrome c550 monomer from strain 52.44 is proposed and a dipole moment of 945 debye was calculated with an orientation close to the exposed haem edge. The behaviour of the cytochrome on molecular-exclusion chromatography is indicative of an ionic strength-dependent monomer (15 kDa)/dimer (30 kDa) equilibrium that can also be detected by 1H-NMR spectroscopy. The apparent mass of 50 kDa observed at very low ionic strength was consistent with the presence of a strongly asymmetric dimer. This was confirmed by cross-linking studies, which showed that a cross-linked species of mass 30 kDa on SDS behaved with an apparent mass of 50 kDa on molecular-exclusion chromatography. A programme which carried out and evaluated molecular docking of two monomers to give a dimer generated a most probable dimer in which the monomer dipoles lay almost antiparallel to each other. The resultant dipole moment of the dimer is therefore small. Although this finding calls into question the possibility of preorientation of a strongly asymmetrically charged cytochrome as it collides with a redox partner, the stoichiometry of complex formation with cytochrome c peroxidase as studied by 1H-NMR spectroscopy shows that it is the monomer that binds.

Amino Acid Sequence↗

PubMedDate not supplied

Virtual screening of novel noncovalent inhibitors for SARS-CoV 3C-like proteinase.

The SARS coronavirus 3C-like proteinase is considered as a potential drug design target for the treatment of severe acute respiratory syndrome (SARS). Owing to the lack of available drugs for the treatment of SARS, the discovery of inhibitors for SARS coronavirus 3C-like proteinase that can potentially be optimized as drugs appears to be highly desirable. We have built a "flexible" three-dimensional model for SARS 3C-like proteinase by homology modeling and multicanonical molecular dynamics method and used the model for virtual screening of chemical databases. After Dock procedures, strategies including pharmocophore model, consensus scoring, and "drug-like" filters were applied in order to accelerate the process and improve the success rate of virtual docking screening hit lists. Forty compounds were purchased and tested by HPLC and colorimetric assay against SARS 3C-like proteinase. Three of them including calmidazolium, a well-known antagonist of calmodulin, were found to inhibit the enzyme with an apparent K(i) from 61 to 178 microM. These active compounds and their binding modes provide useful information for understanding the binding sites and for further selective drug design against SARS and other coronavirus.

Binding Sites↗

PubMed2002

Progress toward virtual screening for drug side effects.

The development and application of a computational protocol for conducting virtual screens of drug side interactions is described. A conventional drug-docking algorithm (AutoDock) is used to conduct two separate studies. First, a series of docking simulations is performed by using guanosine diphosphate and adenosine diphosphate as prototype drugs with the goal of successfully differentiating known receptors from a large number of bait receptors. Using the electrostatic energy of the purine ring as a basis for discrimination allows the correct identification of receptors in blind studies with 100% specificity and 94% sensitivity. In a second study, similar methodology is used to investigate the binding of clinically relevant inhibitors (Gleevec, purvalanol A, and hymenialdisine) to a variety of protein kinase targets. Overall, excellent agreement between experimental and predicted preferences for kinase targets is obtained. An important conclusion from the latter study is that homology-modeled structures of putative receptors may reasonably be used as targets for docking when directly solved crystal structures are not available. The prospects for routine application of the methodology as a means of identifying potential side interactions of candidate drugs are discussed.

Adenosine Diphosphate↗

PubMed1996

Validation of the use of intermolecular NOE constraints for obtaining docked structures of protein-ligand complexes.

The use of intermolecular NOEs for docking a small ligand molecule into its target protein has been investigated with the aim of determining the effectiveness and methodology of this type of NOE docking calculation. A high-resolution X-ray structure of a protein-ligand complex has been used to simulate loose distance constraints of varying degrees of quality, typical of those estimated from experimental NOE intensities. These simulated data were used to examine the effect of the number, distribution and representation of the experimental constraints on the precision and accuracy of the calculated structures. A standard simulated annealing protocol was used, as well as a more novel method based on rigid-body dynamics. The results showed some analogies with those from similar studies on complete protein NMR structure determinations, but it was found that more constraints per torsion angle are required to define docked structures of similar quality. The effectiveness of different NOE-constraint averaging methods was explored and the benefits of using 'R-6 averaging' rather than 'centre averaging' with small sets of NOE constraints were shown. The starting protein structure used in docking calculations was obtained from previous X-ray or NMR structure studies on a related complex. The effects on the calculated conformations of introducing structural differences into the binding site of the initial protein structure were also considered.

Binding Sites↗

PubMed2001

Protein ligand docking based on empirical method for binding affinity estimation.

An empirical protein-ligand binding affinity estimation method, SCORE, was incorporated into a popular docking program, DOCK4. The combined program, ScoreDock, was used to reconstruct the 200 protein-ligand complex structures and found to give good results for the complexes with high binding affinities. A quality assessment method for docking results from ScoreDock was developed based on the whole test set and tested by additionally selected complexes. The method significantly improves the docking accuracy and was shown to be reliable in docking quality assessment. As a docking tool in structural based drug design, ScoreDock can screen out final hits directly based on the predicted negative logarithms of dissociation equilibrium constants of protein-ligand complexes, and can explicitly deal with structure water molecules, as well as metal atoms.

Computer Simulation↗

PubMed2002

Protein-protein docking with multiple residue conformations and residue substitutions.

The protein docking problem has two major aspects: sampling conformations and orientations, and scoring them for fit. To investigate the extent to which the protein docking problem may be attributed to the sampling of ligand side-chain conformations, multiple conformations of multiple residues were calculated for the uncomplexed (unbound) structures of protein ligands. These ligand conformations were docked into both the complexed (bound) and unbound conformations of the cognate receptors, and their energies were evaluated using an atomistic potential function. The following questions were considered: (1) does the ensemble of precalculated ligand conformations contain a structure similar to the bound form of the ligand? (2) Can the large number of conformations that are calculated be efficiently docked into the receptors? (3) Can near-native complexes be distinguished from non-native complexes? Results from seven test systems suggest that the precalculated ensembles do include side-chain conformations similar to those adopted in the experimental complexes. By assuming additivity among the side chains, the ensemble can be docked in less than 12 h on a desktop computer. These multiconformer dockings produce near-native complexes and also non-native complexes. When docked against the bound conformations of the receptors, the near-native complexes of the unbound ligand were always distinguishable from the non-native complexes. When docked against the unbound conformations of the receptors, the near-native dockings could usually, but not always, be distinguished from the non-native complexes. In every case, docking the unbound ligands with flexible side chains led to better energies and a better distinction between near-native and non-native fits. An extension of this algorithm allowed for docking multiple residue substitutions (mutants) in addition to multiple conformations. The rankings of the docked mutant proteins correlated with experimental binding affinities. These results suggest that sampling multiple residue conformations and residue substitutions of the unbound ligand contributes to, but does not fully provide, a solution to the protein docking problem. Conformational sampling allows a classical atomistic scoring function to be used; such a function may contribute to better selectivity between near-native and non-native complexes. Allowing for receptor flexibility may further extend these results.

Amino Acid Substitution↗

PubMedDate not supplied

Interactions of some PGHS-2 selective inhibitors with the PGHS-1: an automated docking study by BioDock.

The automated stochastic docking procedure BioDock has been applied to a series of inhibitors of PGH synthase, the key enzyme in the synthesis of eicosanoids from arachidonic acid. Some PGHS-2 selective inhibitors have been docked to the structure of the ovine PGHS-1 enzyme, as recently obtained by means of X-ray crystallographic analysis, in order to highlight possible structural bases for selectivity.

Computer Simulation↗

PubMedDate not supplied

A technique to study molecular recognition in drug design: preliminary application of free energy derivatives to inhibition of a malarial cysteine protease.

We present molecular dynamics studies on model complexes of inhibitors of a malarial cysteine protease. The initial model for such complexes came from the model building of the protein using its homology with other cysteine proteases and calculations using DOCK to generate new lead compounds. Some of the initial model-built structures were quite stable for 100 psec of dynamics; others moved significantly from their model-built orientation. We also calculated the free energy derivatives at each atom in the inhibitor, both in water and in the binding site. The results of these calculations suggest directions for the design of new, more potent enzyme inhibitors and agree qualitatively with some of the experimental findings. Nonetheless, we stress that we have only used this methodology in an interpretive rather than a predictive manner.

Binding Sites↗

PubMed1998

Mechanism of action of chromogranin A on catecholamine release: molecular modeling of the catestatin region reveals a beta-strand/loop/beta-strand structure secured by hydrophobic interactions and predictive of activity.

A novel fragment of chromogranin A, known as 'catestatin' (bovine chromogranin A344-364), inhibits catecholamine release from chromaffin cells and noradrenergic neurons by acting as a non-competitive nicotinic cholinergic antagonist, and may therefore constitute an endogenous autocrine feedback regulator of sympathoadrenal activity. To characterize how this activity depends on the peptide's structure, we searched for common 3-dimensional motifs for this primary structure or its homologs. Catestatin's primary structure bore significant (29-35.5% identity, general alignment score 44-57) sequence homology to fragment sequences within three homologs of known 3-dimensional structures, based on solved X-ray crystals: 8FAB, IPKM, and 2IG2. Each of these sequences exists in nature as a beta-strand/loop/beta-strand structure, stabilized by hydrophobic interactions between the beta-strands. The catestatin structure was stable during molecular dynamics simulations. The catestatin loop contains three Arg residues, whose electropositive side chains form the terminus of the structure, and give rise to substantial uncompensated charge asymmetry in the molecule. A hydrophobic moment plot revealed that catestatin is the only segment of chromogranin A predicted to contain amphiphilic beta-strand. Circular dichroism in the far ultraviolet showed substantial (63%) beta-sheet structure, especially in a hydrophobic environment. Alanine-substitution mutants of catestatin established a crucial role for the three central arginine residues in the loop (Arg351, Arg353, and Arg358), though not for two arginine residues in the strand region toward the amino-terminus. [125I]Catestatin bound to Torpedo membranes at a site other than the nicotinic agonist binding site. When the catestatin structure was 'docked' with the extracellular domain of the Torpedo nicotinic cholinergic receptor, it interacted principally with the beta and delta subunits, in a relatively hydrophobic region of the cation pore extracellular orifice, and the complex of ligand and receptor largely occluded the cation pore, providing a structural basis for the non-competitive nicotinic cholinergic antagonist properties of the peptide. We conclude that a homology model of catestatin correctly predicts actual features of the peptide, both physical and biological. The model suggests particular spatial and charge features of the peptide which may serve as starting points in the development of non-peptide mimetics of this endogenous nicotinic cholinergic antagonist.

Animals↗

PubMedDate not supplied

Miglustat (NB-DNJ) works as a chaperone for mutated acid beta-glucosidase in cells transfected with several Gaucher disease mutations.

Gaucher disease (GD) is a disorder of glycosphinglipid metabolism caused by deficiency of lysosomal acid beta-glucosidase (GC), resulting in progressive deposition of glucosylceramide in macrophages. The glucose analogue, N-butyl-deoxynojirimycin (NB-DNJ, Miglustat), is an inhibitor of the ceramide-specific glucosyltransferase (CSG) which catalyzes the first step of glycosphingolipids biosynthesis and is currently approved for the oral treatment of type 1 GD. Using site-directed mutagenesis, we constructed plasmids containing wild-type and several mutations in glucocerebrosidase (GBA) gene. The plasmids were transfected into COS-7 cells and stable transfected cell lines were obtained by geneticin (G418) selection. Cells were cultured during 6 days with medium with or without 10 microM NB-DNJ. The addition of NB-DNJ to COS-7 cell medium leads to 1.3-, 2.1-, 2.3-, 3.6-, and 9.9-fold increase in the activity of S364R, wild-type, N370S, V15M, and M123T GC, respectively. However, no significant changes were observed in the activity of the L444P, L336P, and S465del mutated proteins, but a small decrease in the rare P266L variant was observed. These results suggest that NB-DNJ, in addition to the inhibitory effect on CSG, also works as a "chemical chaperone", increasing the activity of acid beta-glucosidase of wild-type and several GC mutated proteins, including the most frequent N370S mutation. The specific location of the Miglustat binding site in GC is unknown. Potential binding sites in the enzyme have been searched for using computational molecular docking. The searching strategy identified three potential GC binding sites for Miglustat, one being the substrate-binding site of the enzyme, which was the best-ranked site by AutoDock program. Therefore, it is possible that Miglustat exerts its chaperoning activity on acid beta-glucosidase by acting as an inhibitor bound at the active site. This increase on the activity of the acid beta-glucosidase would imply that Miglustat is not only a substrate reducer but also an inhibitor of the GC degradation, with very promising clinical implications for the treatment of GD patients.

1-Deoxynojirimycin↗