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 1,045 records · Page 58Linked to original sources

PubMed1994

Molecular modeling of adenosine receptors. The ligand binding site on the rat adenosine A2A receptor.

The amino acid sequence of the rat adenosine A2A receptor and the atomic coordinates of bacteriorhodopsin were combined to generate a three-dimensional model for the adenosine A2A receptor. This model consists of seven amphipathic alpha-helices, forming a pore that is rather hydrophilic compared to the hydrophobic outside of the protein. Subsequently, a highly potent and selective ligand for this receptor, 2-(cyclohexylmethylidinehydrazino)adenosine (SHA 174), was docked into this cavity. A binding site is proposed that takes into account the conformational characteristics of the ligand. Moreover, it involves two histidine residues that were shown to be important for ligand coordination from chemical modification studies. Subsequently, the deduced binding site was used to model other potent ligands, including 8-(3-chlorostyryl)caffeine, a new A2-selective antagonist, that could all be accommodated consistent with earlier biochemical and pharmacological findings. Finally, some thoughts on how adenosine receptor activation might proceed are put forward, based on structural analogies with the enzyme family of serine proteases.

Adenosine↗

PubMed2005

Combining protein modeling and 6D-QSAR. Simulating the binding of structurally diverse ligands to the estrogen receptor.

We present a concept for the in silico simulation of adverse effects triggered by drugs and chemicals. The underlying philosophy combines flexible docking (software Yeti) for the identification of the binding mode(s) and 6D-QSAR (software Quasar) for their quantification. The results obtained for 106 diverse molecules binding to the estrogen receptor (q2 = 0.903; p2 = 0.885) suggest that our approach is suitable for the identification of an endocrine-disrupting potential associated with drugs and chemicals.

Binding Sites↗

PubMedDate not supplied

A combination of molecular dynamics and docking calculations to explore the binding mode of ADS-J1, a polyanionic compound endowed with anti-HIV-1 activity.

The HIV-1 entry process is an important target for the design of new pharmaceuticals for the multidrug therapy of AIDS. A lot of polyanionic compounds, such as polysulfonated and polysulfated, are reported in the literature for their ability to block early stages of HIV-1 replication. Several studies have been performed to elucidate the mechanism of the anti-HIV-1 activity of sulfated polysaccharides and polyanions in general, including binding to cell surface CD4 and interfering with the gp120-coreceptor interaction. Here, we show molecular modeling investigations on ADS-J1, a polyanionic compound with anti-HIV activity that is able to interfere with gp120-coreceptor interactions. Agreeing with experimental data, computer simulations suggested that the V3 loop of gp120 was the preferential binding site for ADS-J1 onto HIV-1. Moreover, mutations induced by the inhibitor significantly changed the stereoelectronic properties of the gp120 surface, justifying a marked drop in the affinity of ADS-J1 toward an ADS-J1-resistant HIV-1 strain.

Anti-HIV Agents↗

PubMed2000

Automated docking of alpha-(1-->4)- and alpha-(1-->6)-linked glucosyl trisaccharides and maltopentaose into the soybean beta-amylase active site.

The Lamarckian genetic algorithm of AutoDock 3.0 was used to dock alpha-maltotriose, methyl alpha-panoside, methyl alpha-isopanoside, methyl alpha-isomaltotrioside, methyl alpha-(6(1)-alpha-glucopyranosyl)-maltoside, and alpha-maltopentaose into the closed and, except for alpha-maltopentaose, into the open conformation of the soybean beta-amylase active site. In the closed conformation, the hinged flap at the mouth of the active site closes over the substrate. The nonreducing end of alpha-maltotriose docks preferentially to subsites -2 or +1, the latter yielding nonproductive binding. Some ligands dock into less optimal conformations with the nonreducing end at subsite -1. The reducing-end glucosyl residue of nonproductively-bound alpha-maltotriose is close to residue Gln194, which likely contributes to binding to subsite +3. In the open conformation, the substrate hydrogen-bonds with several residues of the open flap. When the flap closes, the substrate productively docks if the nonreducing end is near subsites -2 or -1. Trisaccharides with alpha-(1-->6) bonds do not successfully dock except for methyl alpha-isopanoside, whose first and second glucosyl rings dock exceptionally well into subsites -2 and -1. The alpha-(1-->6) bond between the second and third glucosyl units causes the latter to be improperly positioned into subsite +1; the fact that isopanose is not a substrate of beta-amylase indicates that binding to this subsite is critical for hydrolysis.

Algorithms↗

PubMed2006

The potential of P1 site alterations in peptidomimetic protease inhibitors as suggested by virtual screening and explored by the use of C-C-coupling reagents.

A synthetic concept is presented that allows the construction of peptide isostere libraries through polymer-supported C-acylation reactions. A phosphorane linker reagent is used as a carbanion equivalent; by employing MSNT as a coupling reagent, the C-acylation can be conducted without racemization. Diastereoselective reduction was effected with L-selectride. The reagent linker allows the preparation of a norstatine library with full variation of the isosteric positions including the P1 side chain that addresses the protease S1 pocket. Therefore, the concept was employed to investigate the P1 site specificity of peptide isostere inhibitors systematically. The S1 pocket of several aspartic proteases including plasmepsin II and cathepsin D was modeled and docked with approximately 500 amino acid side chains. Inspired by this virtual screen, a P1 site mutation library was designed, synthesized, and screened against three aspartic proteases (plasmepsin II, HIV protease, and cathepsin D). The potency of norstatine inhibitors was found to depend strongly on the P1 substituent. Large, hydrophobic residues such as biphenyl, 4-bromophenyl, and 4-nitrophenyl enhanced the inhibitory activity (IC50) by up to 70-fold against plasmepsin II. In addition, P1 variation introduced significant selectivity, as up to 9-fold greater activity was found against plasmepsin II relative to human cathepsin D. The active P1 site residues did not fit into the crystal structure; however, molecular dynamics simulation suggested a possible alternative binding mode.

Aminocaproates↗

PubMed2005

Biologically enhanced sampling geometric docking and backbone flexibility treatment with multiconformational superposition.

An efficient biologically enhanced sampling geometric docking method is presented based on the FTDock algorithm to predict the protein-protein binding modes. The active site data from different sources, such as biochemical and biophysical experiments or theoretical analyses of sequence data, can be incorporated in the rotation-translation scan. When discretizing a protein onto a 3-dimensional (3D) grid, a zero value is given to grid points outside a sphere centered on the geometric center of specified residues. In this way, docking solutions are biased toward modes where the interface region is inside the sphere. We also adopt a multiconformational superposition scheme to represent backbone flexibility in the proteins. When these procedures were applied to the targets of CAPRI, a larger number of hits and smaller ligand root-mean-square deviations (RMSDs) were obtained at the conformational search stage in all cases, and especially Target 19. With Target 18, only 1 near-native structure was retained by the biologically enhanced sampling geometric docking method, but this number increased to 53 and the least ligand RMSD decreased from 8.1 A to 2.9 A after performing multiconformational superposition. These results were obtained after the CAPRI prediction deadlines.

Algorithms↗

PubMed2003

Analysis and optimization of structure-based virtual screening protocols (1): exploration of ligand conformational sampling techniques.

Ligand conformational flexibility has long been recognized as an important issue in virtual screening (VS). To this end, a number of different methodologies have been adapted to tackle the problem. Many of said techniques were originally designed for ligand derived pharmacophore screens, but have subsequently been fashioned for application within structure-based virtual screening (SVS). A popular adaptation is the pre-calculation of diverse ligand conformations for subsequent docking in target active sites. In this paper, we study a number of the software programs currently being used in conformer generation, analyzing their ability to regenerate known ligand binding conformations. The implications of these studies are discussed, from the perspective of VS in general and SVS in particular.

Computational Biology↗

PubMed1995

Solution structure of a band 3 peptide inhibitor bound to aldolase: a proposed mechanism for regulating binding by tyrosine phosphorylation.

Human erythrocyte band 3 inhibits glycolytic enzymes, including aldolase, by binding these cytoplasmic enzymes at its N-terminus. Phosphorylation of Y8 disrupts inhibition, and there is evidence that in vivo glycolysis levels in erythrocytes are regulated in part by a phosphorylation/dephosphorylation signaling pathway. The structural basis for control by phosphorylation has been investigated by NMR studies on a complex between aldolase and a synthetic peptide corresponding to the first 15 residues of band 3 (MEELQDDYEDMMEEN-NH2). The structure of this band 3 peptide (B3P) when it is bound to rabbit muscle aldolase was determined using the exchange-transferred nuclear Overhauser effect (ETNOE). Two hundred NMR structures for B3P were generated by simulated annealing molecular dynamics with NMR-derived distance restraints and excluding electrostatic terms. Twenty structures were further refined against a force field including full partial charges. The important conformational feature of B3P in the bound state is a folded loop structure involving residues 4-9 and M12 that surrounds Y8 and is stabilized by a hydrophobic cluster with the ring of Y8 sandwiched between the methyl groups of L4 and M12. Differential line broadening indicates that this loop structure binds aldolase in a relatively specific manner, while terminal regions are structurally heterogeneous. To better understand B3P inhibition of aldolase and the mechanism of phosphorylation control, a complex was modeled by docking B3P into the active site of aldolase and optimizing the fit using restrained molecular dynamics and energy minimization. The B3P loop is complementary in conformation to the beta-barrel central core containing the aldolase active site residues. Binding is electrostatic in nature with numerous ionic and hydrogen-bonding interactions involving several conserved lysine and arginine residues of aldolase. How phosphorylation of band 3 could disrupt inhibition was considered by modeling a phosphoryl moiety onto Y8 of B3P. An energetic analysis with respect to rigid phosphate rotation suggests that aldolase inhibition is reversed primarily because of electrostatic repulsion between B3P residues that destabilizes the B3P loop formed in the complex. This proposed intramolecular mechanism for blocking protein--protein association by electrostatic repulsion with the phosphoryl group may be applicable to other protein--protein signaling complexes.

Amino Acid Sequence↗

PubMed2001

A comparative docking study and the design of potentially selective MMP inhibitors.

As part of a program aimed at the design and synthesis of constrained MMP inhibitors, a survey of the reported X-ray and NMR structures of MMP/inhibitor complexes was performed, revealing mutations of key amino acids at different subsites between MMPs. A comparative study of fully automated docking programs AutoDock and DOCK in closely approximating the X-ray crystal structures of ten selected MMP inhibitors was performed. AutoDock proved to be highly reliable, efficient and predictive for a set of inhibitors with less than six atom types.

Binding Sites↗

PubMed2000

Docking of peptide-T onto the D1 domain of the CD4 receptor.

Peptide T (pepT) is a segment of the human immunodeficiency virus (HIV) envelope protein gp120. The peptide competitively binds to the CD4 receptor of a subset of peripheral T lymphocytes and inhibits binding of gp120. Previous studies of this laboratory allowed the assessment of a bioactive form of the peptide and a pharmacophore for the peptide-receptor interaction. In the present study the proposed bioactive form of pepT and its (4-8) segment, the smallest pepT fragment shown to retain full activity, were docked onto the D1 domain of the CD4 receptor. The bioactive conformation of the peptides complements well a cleft on the surface of the CD4 receptor, shown to be the attachment site of gp120 from site directed mutagenesis experiments. These studies provide an improved description of the ligand-receptor pharmacophore.

Amygdalin↗

PubMed2004

Identification of binding epitope of a monoclonal antibody (Z12) against human TNF-alpha using computer modeling and deletion mutant technique.

The genes of the heavy and light chain variable region (VH, VL) of Z12 antibody against hTNF-alpha were cloned, and according to the translated sequence of amino acids, the spatial structures of VH and VL domains were modeled by using homology-based modeling method, followed by constructing the whole three-dimensional structure of Fv fragment. The complex model of Fv interacting with hTNF-alpha was gained with computer-guided molecular docking method, based on which, it was predicted that the epitope recognized by Z12 was from 141 to 146 of hTNF-alpha. hTNF-alpha molecule was divided into two fragments of N-terminal region from 1 to 91 and C-terminal region from 92 to 157 with prokaryotic expression. The measured results suggested that the antigenic epitope recognized by Z12 antibody was located in the C-terminal region 92-157 of hTNF-alpha, proving the predicted result reliable preliminarily. Further experimental results showed that after hTNF-a 141-146 residues were deleted, Z12 antibody almost lost the ability to recognize the mutant, suggesting that the amino acid residues from 141 to 146 of hTNF-alpha were specially recognized by Z12 antibody.

Animals↗

PubMed2001

The hairpin ribozyme substrate binding-domain: a highly constrained D-shaped conformation.

The two domains of the hairpin ribozyme-substrate complex, usually depicted as straight structural elements, must interact with one another in order to form an active conformation. Little is known about the internal geometry of the individual domains in an active docked complex. Using various crosslinking and structural approaches in conjunction with molecular modeling (constraint-satisfaction program MC-SYM), we have investigated the conformation of the substrate-binding domain in the context of the active docked ribozyme-substrate complex. The model generated by MC-SYM showed that the domain is not straight but adopts a bent conformation (D-shaped) in the docked state of the ribozyme, indicating that the two helices bounding the internal loop are closer than was previously assumed. This arrangement rationalizes the observed ability of hairpin ribozymes with a circularized substrate-binding strand to cleave a circular substrate, and provides essential information concerning the organization of the substrate in the active conformation. The internal geometry of the substrate-binding strand places G8 of the substrate-binding strand near the cleavage site, which has allowed us to predict the crucial role played by this nucleotide in the reaction chemistry.

Base Sequence↗

PubMed2005

Biological and molecular properties of a new alpha(v)beta3/alpha(v)beta5 integrin antagonist.

The aim of the present study was to identify specific alpha(v)beta3/alpha(v)beta5 integrin antagonists active on tumor-induced angiogenesis. To this purpose, in vitro integrin-binding assays were used to screen a library of conformationally constrained bicyclic lactam Arg-Gly-Asp-containing pseudopeptides. The results identified ST1646 as a high-affinity specific ligand for alpha(v)beta3 and alpha(v)beta5 integrins with negligible interacting with alpha5beta1 integrin. In all the assays, ST1646 was equipotent to or more potent than the well-characterized integrin antagonists c(RGDfV) and cyclo(Arg-Gly-Asp-d-Phe-[NMe]Val) (EMD121974). In the chorioallantoic membrane assay, topical administration of ST1646 was able to prevent the angiogenic responses elicited by recombinant fibroblast growth factor-2 or vascular endothelial growth factor. In addition, systemic administration of ST1646 in mice exerted a significant antiangiogenic activity on neovascularization triggered by mammary carcinoma MDA-MB435 cells implanted s.c. in a dorsal air sac via a (Millipore Filter Corporation, Bedford, MA) chamber. Moreover, ST1646 delivery via an osmotic pump inhibited the growth and vascularization of tumor xenografts originating from the injection of alpha(v)beta3/alpha(v)beta5-expressing human ovarian carcinoma cells in nude mice. In agreement with the biochemical and pharmacologic studies, Monte Carlo/Stochastic Dynamics simulation showed that the bicyclic scaffold in ST1646 forced the compound to assume a preferred conformation superimposable to the X-ray conformation of alpha(v)beta3-bound EMD121974. Accordingly, computer-docking studies indicated that the ST1646-alpha(v)beta3 integrin complex maintains the ligand-receptor distances and interactions observed in the crystalline EMD121974-alpha(v)beta3 integrin complex. Taken together, these observations indicate that ST1646 represents a dual alpha(v)beta3/alpha(v)beta5 integrin antagonist with interesting biochemical and biological features to be tested in cancer therapy.

Animals↗

PubMed1996

CORCEMA evaluation of the potential role of intermolecular transferred NOESY in the characterization of ligand-receptor complexes.

We report a theoretical characterization of the intermolecular transferred NOESY (inter-TrNOESY) between ligands and receptor macromolecules that bind reversibly, using a COmplete Relaxation and Conformational Exchange MAtrix (CORCEMA) theory developed in our laboratory. We examine the dependence of inter-TrNOESY on the dissociation constant, off-rate, ligand-to-receptor ratio, and distance variations between protons of interacting species within the complex. These factors are analyzed from simulations on two model systems: (i) neuraminidase complexed to a transition-state analogue; and (ii) thermolysin complexed to a leucine-based inhibitor. The latter case utilizes a three-state model of interaction to simulate the effect of hinge-bending motions on the inter-TrNOESY. Our calculations suggest a potential role for inter-TrNOESY (when observable) and CORCEMA analysis in properly docking the ligand within the active site, and in refining the conformation of the ligand-receptor (active-site) complex. These findings have implications on the structure-based design of ligands (e.g., inhibitors) reversibly binding to receptors (e.g., enzymes).

Binding Sites↗

PubMed2006

Biochemical profiling in silico--predicting substrate specificities of large enzyme families.

A general high-throughput method for in silico biochemical profiling of enzyme families has been developed based on covalent docking of potential substrates into the binding sites of target enzymes. The method has been tested by systematically docking transition state--analogous intermediates of 12 substrates into the binding sites of 20 alpha/beta hydrolases from 15 homologous families. To evaluate the effect of side chain orientations to the docking results, 137 crystal structures were included in the analysis. A good substrate must fulfil two criteria: it must bind in a productive geometry with four hydrogen bonds between the substrate and the catalytic histidine and the oxyanion hole, and a high affinity of the enzyme-substrate complex as predicted by a high docking score. The modelling results in general reproduce experimental data on substrate specificity and stereoselectivity: the differences in substrate specificity of cholinesterases toward acetyl- and butyrylcholine, the changes of activity of lipases and esterases upon the size of the acid moieties, activity of lipases and esterases toward tertiary alcohols, and the stereopreference of lipases and esterases toward chiral secondary alcohols. Rigidity of the docking procedure was the major reason for false positive and false negative predictions, as the geometry of the complex and docking score may sensitively depend on the orientation of individual side chains. Therefore, appropriate structures have to be identified. In silico biochemical profiling provides a time efficient and cost saving protocol for virtual screening to identify the potential substrates of the members of large enzyme family from a library of molecules.

Binding Sites↗

PubMed2002

Docking studies reveal a selective binding of D-penicillamine to the transactivator protein of human immunodeficiency virus type 1.

DOCK and Affinity studies were carried out to study the binding of D- and L-penicillamine to the transactivator protein (tat) of human immunodeficiency virus type 1 (HIV-1). These studies reveal a selective binding of D-penicillamine to the cysteine-rich region covering amino acid residues 20-38 of the tat protein. A careful analysis of the components of the binding energy of the D- and L-isomers reveals that the D-isomer has a more favorable van der Waals interaction resulting from an optimal placement of the dimethylthiomethyl side chain in the binding site. This observation matches the experimental data that D-penicillamine is a more potent inhibitor of tat-mediated transactivation than the L-isomer. The docking and experimental data offer an interesting approach to design structural molecules with potential application to block signal functions of the tat protein in HIV-1 pathogenesis.

Binding Sites↗

PubMed2004

Benzamides and benzamidines as specific inhibitors of epidermal growth factor receptor and v-Src protein tyrosine kinases.

The benzamides 1 and the benzamidines 2 as well as the cyclic benzamidines 3 were designed and synthesized as the mimics of 4-anilinoquinazolines for an inhibitor of EGFR tyrosine kinase. The specific inhibitions of EGFR tyrosine kinase were observed in the benzamides 1c and 1d, and the benzamidine 2a, whereas the specific inhibitions of v-Src kinase were observed in the benzamide 1j and the benzamidine 2d at a 10microg/mL concentration of compounds. The cyclic benzamidines 3a and 3b showed potent kinase inhibition of EGFR at a 1.0microg/mL concentration. According to the docking simulation using the X-ray structure of EGFR kinase domain in complex with erlotinib, the LigScore2 scoring function value of erlotinib was calculated as 5.61, whereas that of the benzamide 1c was 5.05. In a similar manner, the LigScore2 value of the cyclic benzamidine 3a was calculated as 5.10.

Benzamides↗

PubMed2006

Pheromone discrimination by the pheromone-binding protein of Bombyx mori.

Pheromone-binding proteins are postulated to contribute to the exquisite specificity of the insect's olfactory system, acting as a filter by preferentially binding only one of the components of the natural pheromone. Here, we investigated the possible discrimination of the two very similar components of the natural pheromone gland from the silk moth, Bombyx mori, bombykol and bombykal, by the only pheromone-binding protein (BmorPBP) known to be expressed in the pheromone-detecting sensilla. Free-energy calculations and virtual docking indicate that both bombykol and bombykal bind to BmorPBP with similar affinity. In addition, in vitro competitive binding assays showed that both bombykol and bombykal were bound by BmorPBP with nearly the same high affinity. While BmorPBP might filter out other physiologically irrelevant compounds hitting the sensillar lymph, discrimination between the natural pheromone compounds must be achieved by molecular interactions with their cognate receptors.

Alkadienes↗