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 829 records · Page 46Linked to original sources

PubMed2002

Inhibition of trypsin by cowpea thionin: characterization, molecular modeling, and docking.

Higher plants produce several families of proteins with toxic properties, which act as defense compounds against pests and pathogens. The thionin family represents one family and comprises low molecular mass cysteine-rich proteins, usually basic and distributed in different plant tissues. Here, we report the purification and characterization of a new thionin from cowpea (Vigna unguiculata) with proteinase inhibitory activity. Cowpea thionin inhibits trypsin, but not chymotrypsin, binding with a stoichiometry of 1:1 as shown with the use of mass spectrometry. Previous annotations of thionins as proteinase inhibitors were based on their erroneous identification as homologues of Bowman-Birk family inhibitors. Molecular modeling experiments were used to propose a mode of docking of cowpea thionin with trypsin. Consideration of the dynamic properties of the cowpea thionin was essential to arrive at a model with favorable interface characteristics comparable with structures of trypsin-inhibitor complexes determined by X-ray crystallography. In the final model, Lys11 occupies the S1 specificity pocket of trypsin as part of a canonical style interaction.

Amino Acid Sequence↗

PubMed2004

Specificity of Amaranthus leucocarpus syn. hypocondriacus lectin for O-glycopeptides.

Amaranthus leucocarpus syn. hypochondriacus lectin (ALL) has been shown to be specific for N-acetyl-D-galactosamine (GalNAc). In this work, we determined a value of 1.0 x 10(-2) M for the association constant of ALL for GalNAc, calculated using fluorescence spectroscopy assays. Using neoglycopeptides obtained by in vitro O-glycosylation, we determined the main features of O-glycopeptides recognized by ALL using molecular dynamics simulations, capillary electrophoresis, and ELISA. Neo-glycopeptides were obtained by in vitro O-glycosylation reaction using microsomal preparations of murine thymocytes, human gastric fundus and colonic mucosa. ELISA assays were performed with peroxidase-labeled murine monoclonal IgG2, kappa light chain (5D4) antibodies against ALL. Among the in vitro neoglycopeptides, only those of TTSAPTTS containing GalNAc at Thr in #2 and #6 reacted with ALL. Neither the TTSAPTTS glycopeptide, containing a unique GalNAc residue at Thr in #2, nor others (with more than two GalNAc residues) interacted with the lectin. Computational docking assays of the lower energy conformers for interactions between glycopeptides and lectins confirmed that ALL recognized GalNAc residues when they are spaced out in glycan structures, whereas GalNAc residues arranged in clusters prevented interaction with the lectin, indicating that ALL is specific for a special GalNAc-containing motif found in different O-glycoproteins.

Acetylgalactosamine↗

PubMed2004

The three-dimensional structure of Arabidopsis thaliana O-methyltransferase predicted by homology-based modelling.

O-methylation of flavonoid compounds is an important enzymatic reaction since it not only reduces the chemical reactivity of their phenolic hydroxyl groups but also increases their lipophilicity and, hence, their intracellular compartmentation. Several genes encoding flavonoid O-methyltransferases (OMTs) have been isolated and characterized both at the molecular and biochemical levels. In contrast with mammalian enzymes, plant OMTs exhibit narrow substrate specificities as well as position-specific activities, so that the homology comparison, derived using programs such as BLAST can not provide sufficient information on the enzyme function or its substrate preference. In order to study these characteristics, therefore, another approach, homology-based modelling is being carried out. We report here the determination of the 3-D structure of Arabidopsis thaliana O-methyltransferase, AtOMT1 as well as its dynamics when complexed with its substrate. The predicted structure obtained by homology-based modelling is conserved during molecular dynamics simulations. AtOMT1 exhibits a structure similar to that of caffeic acid O-methyltransferase, COMT when the latter was used as a template. Whereas COMT includes 20 alpha-helices and nine beta-sheets, AtOMT1 has 16 and 9, respectively. Although the homology between both proteins is higher than 77% and all amino acids surrounding the active sites, except one residue, are similar in their primary sequences, the two proteins exhibit different substrate preferences. The differences in substrate specificity may be explained on the basis of the predicted structures of the protein and its complex with the substrate. In addition, docking the substrate into the active site of the protein allowed the study of the structural change of the active site on the dihedral angle distribution of the residues surrounding the active site.

Amino Acid Sequence↗

PubMed2000

Research on anti-HIV-1 agents. Investigation on the CD4-Suradista binding mode through docking experiments.

Sulfonated distamycin (Suradista) derivatives exhibit anti-HIV-1 activity by inhibiting the binding of the viral envelope glycoprotein gp120 to its receptor (CD4). With the aim to propose a possible binding mode between Suradistas and the CD4 macromolecule, molecular docking experiments, followed by energy minimization of the complexes thus obtained, were performed. Computational results show that ligand binding at the CD4 surface involves two or three positively charged regions of the macromolecule, in agreement with the results of X-ray crystallographic analysis of a ternary complex (CD4/gp120/neutralizing antibody) recently reported in the literature. Our findings account well for the structure-activity relationship found for Suradista compounds.

Anti-HIV Agents↗

PubMed1998

Aminoglycoside binding to the hammerhead ribozyme: a general model for the interaction of cationic antibiotics with RNA.

A variety of drugs inhibit biological key processes by binding to a specific RNA component. We focus here on the well-analysed hammer-head ribozyme RNA that is inhibited by aminoglycoside antibiotics, a process considered as a paradigm for studying drug/RNA interactions. With insight gained from molecular dynamics simulations of the ribozyme in the presence of Mg2+ identified by crystallography and of aminoglycosides in solution, a general model for aminoglycoside binding to RNA is proposed. A striking structurally based complementarity between the charged ammonium groups of the aminoglycosides and the metal binding sites in the hammerhead was uncovered. Despite dynamical flexibility of the aminoglycosides, several of the intramolecular distances between the charged ammonium groups of the drugs were found to be rather constant. Intramolecular ammonium distances of the aminoglycosides span ranges similar to the interionic distances between Mg2+ in the hammerhead. Successful docking of aminoglycosides to the hammerhead ribozyme could be achieved by positioning the ammonium groups at the sites occupied by Mg2+. The covalently linked ammonium groups of the aminoglycosides are thus able to complement in space the negative electrostatic potential created by a three-dimensional RNA fold. Consequently, it is suggested that aminoglycoside-derived sugars could constitute a basic set of yardstick synthons ideal for rational and combinatorial synthesis of drugs targeted at biologically relevant RNA folds.

Anti-Bacterial Agents↗

PubMed1998

Flexible ligand docking using conformational ensembles.

Molecular docking algorithms suggest possible structures for molecular complexes. They are used to model biological function and to discover potential ligands. A present challenge for docking algorithms is the treatment of molecular flexibility. Here, the rigid body program, DOCK, is modified to allow it to rapidly fit multiple conformations of ligands. Conformations of a given molecule are pre-calculated in the same frame of reference, so that each conformer shares a common rigid fragment with all other conformations. The ligand conformers are then docked together, as an ensemble, into a receptor binding site. This takes advantage of the redundancy present in differing conformers of the same molecule. The algorithm was tested using three organic ligand protein systems and two protein-protein systems. Both the bound and unbound conformations of the receptors were used. The ligand ensemble method found conformations that resembled those determined in X-ray crystal structures (RMS values typically less than 1.5 A). To test the method's usefulness for inhibitor discovery, multi-compound and multi-conformer databases were screened for compounds known to bind to dihydrofolate reductase and compounds known to bind to thymidylate synthase. In both cases, known inhibitors and substrates were identified in conformations resembling those observed experimentally. The ligand ensemble method was 100-fold faster than docking a single conformation at a time and was able to screen a database of over 34 million conformations from 117,000 molecules in one to four CPU days on a workstation.

Algorithms↗

PubMed2003

Discovery of novel low molecular weight inhibitors of IMPDH via virtual needle screening.

Novel, low molecular weight inhibitors of IMPDH have been discovered through the application of a validated virtual screening protocol. A series of 21 IMPDH inhibitors were used to validate the docking procedure. Application of this procedure to the selection of compounds for screening from an in-house database resulted in a 50-fold reduction in the size of the screening set (3425 to 74 compounds) and gave a hit-rate of 10% on biological evaluation.

Binding Sites↗

PubMed2006

3D-QSAR and molecular docking studies on pyrazolopyrimidine derivatives as glycogen synthase kinase-3beta inhibitors.

Glycogen synthase kinase-3 (GSK-3), a serine/threonine kinase, is a fascinating enzyme with diverse biological actions in intracellular signaling systems, making it an emerging target for diseases such as diabetes mellitus, cancer, chronic inflammation, bipolar disorders and Alzheimer's disease. It is important to inhibit GSK-3 selectively and the net effect of the GSK-3 inhibitors thus should be target specific, over other phylogenetically related kinases such as CDK-2. In the present work, we have carried out three-dimensional quantitative structure activity relationship (3D-QSAR) studies on novel class of pyrazolopyrimidine derivatives as GSK-3 inhibitors reported to have improved cellular activity. Docked conformation of the most active molecule in the series, which shows desirable interactions in the receptor, was taken as template for alignment of the molecules. Statistically significant CoMFA and CoMSIA models were generated using 49 molecules in training set. By applying leave-one-out (LOO) cross-validation study, r(cv)2 values of 0.53 and 0.48 for CoMFA and CoMSIA, respectively and non-cross-validated (r(ncv)2) values of 0.98 and 0.92 were obtained for CoMFA and CoMSIA models, respectively. The predictive ability of CoMFA and CoMSIA models was determined using a test set of 12 molecules which gave predictive correlation coefficients (r(pred)2) of 0.47 and 0.48, respectively, indicating good predictive power. Based upon the information derived from CoMFA and CoMSIA contour maps, we have identified some key features that explain the observed variance in the activity and have been used to design new pyrazolopyrimidine derivatives. The designed molecules showed better binding affinity in terms of estimated docking scores with respect to the already reported systems; hence suggesting that newly designed molecules can be more potent and selective towards GSK-3beta inhibition.

Binding, Competitive↗

PubMed2003

Role of the histidine triad-like motif in nucleotide hydrolysis by the rotavirus RNA-packaging protein NSP2.

Octamers formed by the nonstructural protein NSP2 of rotavirus are proposed to function as molecular motors in the packaging of the segmented double-stranded RNA genome. The octamers have RNA binding, helix unwinding, and Mg(2+)-dependent NTPase activities and play a crucial role in assembly of viral replication factories (viroplasms). Comparison of x-ray structures has revealed significant structural homology between NSP2 and the histidine triad (HIT) family of nucleotidyl hydrolases, which in turn has suggested the location of the active site for NTP hydrolysis in NSP2. Consistent with the structural predictions, we show here using site-specific mutagenesis and ATP docking simulations that the active site for NTP hydrolysis is localized to residues within a 25-A-deep cleft between the C- and N-terminal domains of the NSP2 monomer. Although lacking the precise signature HIT motif (HØHØHØØ where Ø is a hydrophobic residue), our analyses demonstrate that histidines (His(221) and His(225)) represent critical residues of the active site. Similar to events occurring during nucleotide hydrolysis by HIT proteins, NTP hydrolysis by NSP2 was found to produce a short lived phosphorylated intermediate. Evaluation of the biological importance of the NTPase activity of NSP2 by transient expression in mammalian cells showed that such activity has no impact on the ability of NSP2 to induce the hyperphosphorylation of NSP5 or to interact with NSP5 to form viroplasm-like structures. Hence the NTPase activity of NSP2 probably has a role subsequent to the formation of viroplasms, consistent with its suspected involvement in RNA packaging and/or replication.

Adaptor Proteins, Signal Transducing↗

PubMed2004

Testing a flexible-receptor docking algorithm in a model binding site.

Sampling receptor flexibility is challenging for database docking. We consider a method that treats multiple flexible regions of the binding site independently, recombining them to generate different discrete conformations. This algorithm scales linearly rather than exponentially with the receptor's degrees of freedom. The method was first evaluated for its ability to identify known ligands of a hydrophobic cavity mutant of T4 lysozyme (L99A). Some 200000 molecules of the Available Chemical Directory (ACD) were docked against an ensemble of cavity conformations. Surprisingly, the enrichment of known ligands from among a much larger number of decoys in the ACD was worse than simply docking to the apo conformation alone. Large decoys, accommodated in the larger cavity conformations sampled in the ensemble, were ranked better than known small ligands. The calculation was redone with an energy correction term that considered the cost of forming the larger cavity conformations. Enrichment improved, as did the balance between high-ranking large and small ligands. In a second retrospective test, the ACD was docked against a conformational ensemble of thymidylate synthase. Compared to docking against individual enzyme conformations, the flexible receptor docking approach improved enrichment of known ligands. Including a receptor conformational energy weighting term improved enrichment further. To test the method prospectively, the ACD database was docked against another cavity mutant of lysozyme (L99A/M102Q). A total of 18 new compounds predicted to bind this polar cavity and to change its conformation were tested experimentally; 14 were found to bind. The bound structures for seven ligands were determined by X-ray crystallography. The predicted geometries of these ligands all corresponded to the observed geometries to within 0.7A RMSD or better. Significant conformational changes of the cavity were observed in all seven complexes. In five structures, part of the observed accommodations were correctly predicted; in two structures, the receptor conformational changes were unanticipated and thus never sampled. These results suggest that although sampling receptor flexibility can lead to novel ligands that would have been missed when docking a rigid structure, it is also important to consider receptor conformational energy.

Algorithms↗

PubMed2002

On the role of the crystal environment in determining protein side-chain conformations.

The role of crystal packing in determining the observed conformations of amino acid side-chains in protein crystals is investigated by (1) analysis of a database of proteins that have been crystallized in different unit cells (space group or unit cell dimensions) and (2) theoretical predictions of side-chain conformations with the crystal environment explicitly represented. Both of these approaches indicate that the crystal environment plays an important role in determining the conformations of polar side-chains on the surfaces of proteins. Inclusion of the crystal environment permits a more sensitive measurement of the achievable accuracy of side-chain prediction programs, when validating against structures obtained by X-ray crystallography. Our side-chain prediction program uses an all-atom force field and a Generalized Born model of solvation and is thus capable of modeling simple packing effects (i.e. van der Waals interactions), electrostatic effects, and desolvation, which are all important mechanisms by which the crystal environment impacts observed side-chain conformations. Our results are also relevant to the understanding of changes in side-chain conformation that may result from ligand docking and protein-protein association, insofar as the results reveal how side-chain conformations change in response to their local environment.

Computer Simulation↗

PubMed2001

Homology modelling and molecular dynamics studies of human placental tissue protein 13 (galectin-13).

The primary structure of the newly sequence analysed placental tissue protein 13 (PP13) was highly homologous to several members of the beta-galactoside-binding S-type lectin (galectin) family. By homology modelling, the three-dimensional structure of PP13 was built based on high-resolution crystal structures of homologues and also their characteristic 'jellyroll' fold was found in the case of PP13. Our model has been deposited in the Brookhaven Protein Data Bank. By multiple sequence alignment and structure-based secondary structure prediction, we underlined the structural similarity of PP13 with its homologues. The secondary structure of PP13 was identical with 'proto-type' galectins consisting of a five- and a six-stranded beta-sheet, joined by two alpha-helices, and galectins' highly conserved carbohydrate-recognition domain (CRD) was also present in PP13. Of the eight consensus residues in the CRD, four identical and three conservatively substituted were shared by PP13. By docking simulations PP13 possessed sugar-binding activity with highest affinity to N-acetyllactosamine and lactose typical of most galectins. All ligands were docked into the putative CRD of PP13. Based on several lines of evidence discussed in this paper demonstrating that PP13 is a novel galectin, PP13 was also designated galectin-13. These computational results provide some new insights into the possible role and importance of PP13 in various processes of the human body and can be of help in the initial steps of further functional research.

Amino Acid Sequence↗

PubMed2004

Electron paramagnetic resonance spectroscopy and molecular modelling of the interaction of myelin basic protein (MBP) with calmodulin (CaM)-diversity and conformational adaptability of MBP CaM-targets.

The classic 18.5 kDa isoform of murine myelin basic protein (mMBP) has been shown to bind calmodulin (CaM) strongly and specifically in vitro. Here, we have used site-directed spin labelling (SDSL) and electron paramagnetic resonance (EPR) spectroscopy to map more precisely the sites of interaction of recombinant mMBP (rmMBP) with CaM. On the basis of these and previous experimental data, and the predictions of CaM-binding motifs using the Calmodulin Target Database (), three main segments of MBP were suggested for the interaction. The first site is located at the C-terminus; the second one lies in the central portion of the protein and forms an amphipathic alpha-helix in reconstituted myelin-mimetic systems; the third is quite close to the N-terminus. The murine Golli-MBP isoform J37 has also been shown to bind CaM in vitro, and an interaction site was predicted in the N-terminal Golli-specific portion of the protein. From these four segments, we selected peptide fragments of 12-14 residues in length, chosen on the bases of their amphipathicity and CaM-target characteristics. We modelled each of these peptides as alpha-helices, and performed docking simulations to investigate their interactions with the CaM peptide-binding tunnel. Different yet almost equally favourable CaM-binding modes were found for each of them. The experimental SDSL/EPR and theoretical modelling results were in good agreement, and supported the conjecture that there are several plausible CaM-binding sites in MBP, that could be induced into an alpha-helical conformation by their interaction with CaM and account for strong immobilisation of spin-labeled residues in all three segments. Phosphorylation and deimination were also emulated and simulated for known sites of MBP post-translational modification. The results obtained confirmed the appropriate utilisation of simple residue substitutions to mimic the natural modifications, and demonstrated molecular mechanisms by which MBP-CaM interactions could be modulated in vivo.

Amino Acid Sequence↗

PubMed2005

MEDock: a web server for efficient prediction of ligand binding sites based on a novel optimization algorithm.

The prediction of ligand binding sites is an essential part of the drug discovery process. Knowing the location of binding sites greatly facilitates the search for hits, the lead optimization process, the design of site-directed mutagenesis experiments and the hunt for structural features that influence the selectivity of binding in order to minimize the drug's adverse effects. However, docking is still the rate-limiting step for such predictions; consequently, much more efficient algorithms are required. In this article, the design of the MEDock web server is described. The goal of this sever is to provide an efficient utility for predicting ligand binding sites. The MEDock web server incorporates a global search strategy that exploits the maximum entropy property of the Gaussian probability distribution in the context of information theory. As a result of the global search strategy, the optimization algorithm incorporated in MEDock is significantly superior when dealing with very rugged energy landscapes, which usually have insurmountable barriers. This article describes four different benchmark cases that span a diverse set of different types of ligand binding interactions. These benchmarks were compared with the use of the Lamarckian genetic algorithm (LGA), which is the major workhorse of the well-known AutoDock program. These results demonstrate that MEDock consistently converged to the correct binding modes with significantly smaller numbers of energy evaluations than the LGA required. When judged by a threshold of the number of energy evaluations consumed in the docking simulation, MEDock also greatly elevates the rate of accurate predictions for all benchmark cases. MEDock is available at http://medock.csie.ntu.edu.tw/ and http://bioinfo.mc.ntu.edu.tw/medock/.

Algorithms↗

PubMed2004

The parsley plastocyanin-turnip cytochrome f complex: a structurally distorted but kinetically functional acidic patch.

In general, inter-protein electron transfer proceeds via the formation of transient complexes. The initial stage of the interaction between plastocyanin (PCu) and cytochrome f (cyt f ) from plants is mediated by complementary electrostatics. Given the diffuse nature of its acidic patch, parsley PCu is an atypical example of a plant PCu. The interaction of this PCu with turnip cyt f was investigated by stopped-flow kinetics, NMR spectroscopy and protein-docking simulations. We show that, despite the altered acidic patch, parsley PCu is as efficient as spinach PCu in accepting electrons from cyt f, over the physiological range of ionic strength. At high ionic strength, the rate constant for the reaction of cyt f with parsley PCu is twice that of the spinach protein. This difference in reactivity is attributed to variations in the hydrophobic patch of parsley PCu. The results of NMR studies and protein-docking simulations indicate that parsley PCu and its spinach analogue adopt different orientations in their complexes with cyt f.

Amino Acid Sequence↗

PubMed2006

Design, synthesis, biological evaluation, and molecular modeling studies of TIBO-like cyclic sulfones as non-nucleoside HIV-1 reverse transcriptase inhibitors.

TIBO- and TBO-like sulfone derivatives 1 and 2 were designed, synthesized, and tested for their ability to block the replication cycle of HIV-1 in infected cells. The anti-HIV-1 activities of sulfones 3, which were intermediates in the syntheses of 1 and 2, were also evaluated. Surprisingly, the sulfone analogues of TIBO R82913 (compounds 1) were inactive, whereas interesting results were obtained for truncated derivatives 2. Compound 2 w was the most potent among this series in cell-based assays (EC50=0.07 microM, CC50>200 microM, SI>2857). It was twofold less potent than R82913, but more selective. An X-ray crystallographic analysis was carried out to establish the absolute configuration of 2 w and its enantiomer 2 x, which were obtained by semipreparative HPLC of 2 v, one of the most potent racemates. Compounds 1-3 were proven to target HIV-1 RT. In fact, representative derivatives inhibited recombinant HIV-1 RT in vitro at concentrations similar to those active in cell-based assays. 3D QSAR studies and docking simulations were developed on TIBO- and TBO-like sulfone derivatives to rationalize their anti-HIV-1 potencies and to predict the activity of novel untested sulfone derivatives. Predictive 3D QSAR models were obtained with a receptor-based alignment by docking of TIBO- and TBO-like derivatives into the NNBS of RT.

Cell Line↗

PubMed2005

Beta2 subunit contribution to 4/7 alpha-conotoxin binding to the nicotinic acetylcholine receptor.

The structures of acetylcholine-binding protein (AChBP) and nicotinic acetylcholine receptor (nAChR) homology models have been used to interpret data from mutagenesis experiments at the nAChR. However, little is known about AChBP-derived structures as predictive tools. Molecular surface analysis of nAChR models has revealed a conserved cleft as the likely binding site for the 4/7 alpha-conotoxins. Here, we used an alpha3beta2 model to identify beta2 subunit residues in this cleft and investigated their influence on the binding of alpha-conotoxins MII, PnIA, and GID to the alpha3beta2 nAChR by two-electrode voltage clamp analysis. Although a beta2-L119Q mutation strongly reduced the affinity of all three alpha-conotoxins, beta2-F117A, beta2-V109A, and beta2-V109G mutations selectively enhanced the binding of MII and GID. An increased activity of alpha-conotoxins GID and MII was also observed when the beta2-F117A mutant was combined with the alpha4 instead of the alpha3 subunit. Investigation of A10L-PnIA indicated that high affinity binding to beta2-F117A, beta2-V109A, and beta2-V109G mutants was conferred by amino acids with a long side chain in position 10 (PnIA numbering). Docking simulations of 4/7 alpha-conotoxin binding to the alpha3beta2 model supported a direct interaction between mutated nAChR residues and alpha-conotoxin residues 6, 7, and 10. Taken together, these data provide evidence that the beta subunit contributes to alpha-conotoxin binding and selectivity and demonstrate that a small cleft leading to the agonist binding site is targeted by alpha-conotoxins to block the nAChR.

Amino Acid Sequence↗

PubMedDate not supplied

[Optical biosensor study and molecular modeling of interactions between cytochrome P450cam and cytochrome b5].

The formation of complexes of cytochrome P450cam (P450cam) with full-length cytochrome b5 (d-b5) and its tryptic water-soluble fragment (t-b5) was analyzed using a two-channel IAsys+ optical biosensor. It was found that t-b5 can form complexes with P450cam, while d-b5 does not interact with P450cam. The involvement of amine groups of P450cam in the complex formation was demonstrated. The temperature dependence of t-b5(im)/P450cam complex formation was measured. The association rate constant (k(on)) increased with temperature, while the dissociation rate constant (k(off)) practically remained unchanged. It was concluded that hydrophobic interactions play a key role in the complex formation, while electrostatic interactions are significant for complex stabilization. Based on temperature dependence the activation energy, enthalpy and entropy of complex formation were calculated. It was shown that the entropy component plays a key role in t-b5(im)/P450cam interaction. Computer modeling of P450cam/t-b5 and P450cam/d-b5 interactions was carried out. Using the method of molecular docking some hypotheses of protein-protein complexes were advanced and the best ones were selected based on geometric complementarity, calculated binding energy and probability of electron tunneling between proteins. The computer modeling has shown that only P450cam and t-b5 can form the stable complex. These results are in good agreement with the experimental data obtained with the optical biosensor.

Amino Acid Sequence↗