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,009 records · Page 56Linked to original sources

PubMed2004

Evaluation of comparative cytochrome P450 2B4 model by photoaffinity labeling.

A homology model of rabbit CYP 2B4 was constructed on the basis of the crystallographic structure of truncated mammalian CYP 2C5/3 and bacterial soluble CYPs. To validate the CYP 2B4 homology model photoaffinity labeling was employed. Three probes (I-III) containing a photo-labile azido-group and an amino-group on opposite ends of the molecule were designed for photoaffinity labeling of the CYP 2B4 in increasing distance from the heme iron. Spectroscopic data proved probes I (the shortest) and II (a middle sized) to be coordinated with the heme iron via their amino-groups in the enzyme active center while the probe III (the longest) was not bound in this way. This binding orientation of probes I and II is in accordance with the model predicting ion-pairing of the negatively charged side chain of CYP 2B4 Asp 105 and a positively charged nitrogen located in an appropriate position in structures of probes I and II, only. The lack of heme binding of the probe III is clear from its docking into the CYP 2B4 model since no Asp 105 ion-pairing is possible. The target of photoactivated probe II, Arg 197, in a distance of about 16.5 A from the heme iron, exactly matches the position of that amino acid residue, predicted from the CYP 2B4 homology model. Moreover, using this technique, a substrate access channel has been identified. To assess the predicted substrate-binding pocket, an interaction of a specific CYP 2B4 substrate, diamantane, was examined. In "silico" docking revealed strong binding of diamantane in an orientation allowing experimentally observed C4-hydroxylation. Our homology model of CYP 2B4 is thus consistent with experimental metabolic and photoaffinity labeling data.

Amino Acid Sequence↗

PubMed1998

Covalent sequestration of the nitrogen mustard mechlorethamine by metallothionein.

The research reported here demonstrates covalent binding to the metal-binding protein metallothionein (MT) by the therapeutic nitrogen mustard mechlorethamine. The most surprising aspect of this interaction is the selectivity of the alkylating agent for specific residues of MT. A combination of MS and proteolytic and enzymatic methods was used to deduce specific locations of mechlorethamine alkylation. These experiments indicated that alkylation occurs predominantly in the carboxyl domain of MT, with one molecule of mechlorethamine covalently cross-linking two cysteine residues. Electrospray MS revealed the retention of all seven metal ions in the cross-linked MT/mechlorethamine adducts, highlighting the uniqueness of this protein. Computerized docking experiments supported the hypothesis that selective binding precedes selective alkylation, and the structure of the drug indicates the minimal structural requirements for this binding. These results support the idea that MT overexpressed in tumor cells contributes to the inactivation of anticancer drugs.

Alkylating Agents↗

PubMed1999

Molecular docking studies on interaction of diverse retinol structures with human alcohol dehydrogenases predict a broad role in retinoid ligand synthesis.

Some members of the human alcohol dehydrogenase (ADH) family possess retinol dehydrogenase activity and may thus function in production of the active nuclear receptor ligand retinoic acid. Many diverse natural forms of retinol exist including all-trans-retinol (vitamin A(1)), 9-cis-retinol, 3,4-didehydroretinol (vitamin A(2)), 4-oxo-retinol, and 4-hydroxy-retinol as well as their respective carboxylic acid derivatives which are active ligands for retinoid receptors. This raises the question of whether ADHs can accommodate all these different retinols and thus participate in the activation of several retinoid ligands. The crystal structures of human ADH1B and ADH4 provide the opportunity to examine their active sites for potential binding to many diverse retinol structures using molecular docking algorithms. The criteria used to score successful docking included achievement of distances of 1.9-2.4 A between the catalytic zinc and the hydroxyl oxygen of retinol and 3.2-3.6 A between C-4 of the coenzyme NAD and C-15 of retinol. These distances are sufficient to enable hydride transfer during the oxidation of an alcohol to an aldehyde. By these criteria, all-trans-retinol, 4-oxo-retinol, and 4-hydroxy-retinol were successfully docked to both ADH1B and ADH4. However, 9-cis-retinol and 3,4-didehydroretinol, which have more restrictive conformations, were successfully docked to only ADH4 which possesses a wider active site than ADH1B and more easily accommodates the C-19 methyl group. Furthermore, docking of all retinols was more favorable in the active site of ADH4 rather than ADH1B as measured by force field and contact scores. These findings suggest that ADH1B has a limited capacity to metabolize retinols, but that ADH4 is well suited to function in the metabolism of many diverse retinols and is predicted to participate in the synthesis of the active ligands all-trans-retinoic acid, 9-cis-retinoic acid, 3, 4-didehydroretinoic acid, 4-oxo-retinoic acid, and 4-hydroxy-retinoic acid.

Alcohol Dehydrogenase↗

PubMed2005

Terephthalamide derivatives as mimetics of helical peptides: disruption of the Bcl-x(L)/Bak interaction.

A series of Bcl-x(L)/Bak antagonists, based on a terephthalamide scaffold, was designed to mimic the alpha-helical region of the Bak peptide. These molecules showed favorable in vitro activities in disrupting the Bcl-x(L)/Bak BH3 domain complex (terephthalamides 9 and 26, K(i) = 0.78 +/- 0.07 and 1.85 +/- 0.32 microM, respectively). Extensive structure-affinity studies demonstrated a correlation between the ability of terephthalamide derivatives to disrupt Bcl-x(L)/Bak complex formation and the size of variable side chains on these molecules. Treatment of human HEK293 cells with the terephthalamide derivative 26 resulted in disruption of the Bcl-x(L)/Bax interaction in whole cells with an IC(50) of 35.0 microM. Computational docking simulations and NMR experiments suggested that the binding cleft for the BH3 domain of the Bak peptide on the surface of Bcl-x(L) is the target area for these synthetic inhibitors.

Biomimetic Materials↗

PubMed2005

Solution structure of a common substrate mimetic of cyclooxygenase-downstream synthases bound to an engineered thromboxane A2 synthase using a high-resolution NMR technique.

Understanding the docking mechanism of the common substrate, prostaglandin H(2) (PGH(2)), into the active sites of different cyclooxygenase(COX)-downstream synthases is a key step toward uncovering the molecular basis of the isomerization of PGH(2) to different prostanoids. A high-resolution NMR spectroscopy was used to determine the conformational changes and solution 3D structure of U44069, a PGH(2) analogue, bound to one of the COX-downstream synthases-an engineered thromboxane A(2) synthase (TXAS). The dynamic binding was clearly observed by (1)D NMR titration. The detailed conformational change and 3D structure of U44069 bound to the TXAS were demonstrated by 2D (1)H NMR experiments using transferred NOEs. Through the assignments for the 2D (1)H NMR spectra, TOCSY, DQF-COSY, NOESY, and the structural calculations based on the NOE constraints, they demonstrated that the widely open conformation with a triangle shape of the free U44069 changed to a compact structure with an oval shape when bound to the TXAS. The putative substrate-binding pocket of the TXAS model fits the conformation of the TXAS-bound U44069 appropriately, but could not fit the free form of U44069. It was the first to provide structural information for the dynamic docking of the PGH(2) mimic of the TXAS in solution, and to imply that PGH(2) undergoes conformational changes when bound to different COX-downstream synthases, which may play important roles in the isomerization of PGH(2) to different prostanoids. The NMR technique can be used as a powerful tool to determine the conformations of PGH(2) bound to other COX-downstream synthases.

Binding Sites↗

PubMed2003

Homology modeling of the central catalytic domain of insertion sequence ISLC3 isolated from Lactobacillus casei ATCC 393.

The tertiary structure of the central catalytic domain of insertion sequence ISLC3 isolated from Lactobacillus casei ATCC 393 was predicted using the homology modeling approach. The novel insertion sequence was isolated by us from the template bacteriophage phiA3 of L.casei ATCC 393. The number of amino acid residues of the ISLC3 central catalytic domain was 116 and was treated as the query sequence. There were five Web-available threading methods used to find some primary structure templates for the query sequence. These primary templates were further screened using the SWISS-MODEL Protein Modeling Server and the default parameter settings therein to give six final structure templates. All of these final structure templates were the integrase (IN) protein of retroviruses. Multiple sequence alignment using these IN sequences against the query one revealed the signature DDE motif. Based on the structures of these final templates, the structure of the query sequence was constructed using the InsightII/Discover/Homology programs. A metal ion, Mg(2+), was inserted into the center of the putative catalytic pocket formed by the DDE residues of the predicted structure in the final rounds of refinement by molecular dynamics (MD) simulations. The structure with a metal ion included was designated with Mg and that without a metal ion was designated free Mg. The average exposed surface area of some hydrophobic residues of both the predicted free Mg and with Mg structures were computed and compared with those computed for the six structure templates. Whereas the predicted with Mg structure was slightly more exposed than the predicted free Mg structure, the former appeared to be more stable than the latter, as revealed by the lower conformation energy recorded for the former during the structure refinement by MD simulations. To verify further the predicted structures, the coordinates of both predicted structures were fed into the ERRAT Protein Verification Server. It was found that the quality of the predicted with Mg structure was much better than that of the free Mg structure. The validation results also indicated that regions of the predicted with Mg structure that can be rejected at the 95% confidence level were approximately 20% whereas those which can be rejected at the same level for the six structure templates were approximately 10%. The predicted with Mg structure was also docked into a short oligonucleotide representing the substrate of the ISLC3 transposase using the DOCK_4.0.2 program. It was found that both Glu140 and Asp68 residues of the DDE motif of the predicted with Mg structure were able to form hydrogen bonds with the DNA substrate, which was similar to what was observed in a docking study using the retrovirus IN 1asu and its DNA substrate.

Amino Acid Sequence↗

PubMed1996

Knowledge-based modeling of a legume lectin and docking of the carbohydrate ligand: the Ulex europaeus lectin I and its interaction with fucose.

Ulex europaeus isolectin I is specific for fucose-containing oligosaccharide such as H type 2 trisaccharide alpha-L-Fuc (1-->2) beta-D-Gal (1-->4) beta-D-GlcNAc. Several legume lectins have been crystallized and modeled, but no structural data are available concerning such fucose-binding lectin. The three-dimensional structure of Ulex europaeus isolectin I has been constructed using seven legume lectins for which high-resolution crystal structures were available. Some conserved water molecules, as well as the structural cations, were taken into account for building the model. In the predicted binding site, the most probable locations of the secondary hydroxyl groups were determined using the GRID method. Several possible orientations could be determined for a fucose residue. All of the four possible conformations compatible with energy calculations display several hydrogen bonds with Asp-87 and Ser-132 and a stacking interaction with Tyr-220 and Phe-136. In two orientations, the O-3 and O-4 hydroxyl groups of fucose are the most buried ones, whereas two other, the O-2 and O-3 hydroxyl groups are at the bottom of the site. Possible docking modes are also studied by analysis of the hydrophobic and hydrophilic surfaces for both the ligand and the protein. The SCORE method allows for a quantitative evaluation of the complementarity of these surfaces, on the basis of molecular lipophilicity calculations. The predictions presented here are compared with known biochemical data.

Artificial Intelligence↗

PubMed2000

Experimental design for analysis of complex kinetics using surface plasmon resonance.

Using BIAcore surface plasmon resonance technology, we found that the real-time association kinetics of Fabs specific for hen egg-white lysozyme did not conform to a 1:1 Langmuir association model. Heterogeneity of the components is not the source of the complex kinetics. Informed by independent structural data suggesting conformational flexibility differences among these antibodies, we chose global mathematical analysis based on a two-phase model, consistent with the encounter-docking view of protein-protein associations. Experimental association times (T(a)) from 2 to 250 min revealed that initial dissociation rates decreased with increasing T(a), confirming a multiphasic association. The relationship between observed dissociation rate and T(a) is characteristic of each antibody-antigen complex. We define a new parameter, T(50), the time at which the encounter and final complexes are of equimolar concentration. The observed T(50) is a function of analyte concentration and the encounter and docking rate constants. Simulations showed that when the ligand is saturated at high analyte concentrations, T(50) reaches a minimum value, T(50)(MIN), which can be used to compare antigen-antibody complexes. For high-affinity complexes with rapid rearrangement to a stable complex, T(50)(MIN) approaches T(1/2) of the rearrangement forward rate constant. We conclude that experiments with a range of T(a) are essential to assess the nature of the kinetics, regardless of whether a two-state or 1:1 model is applicable. We suggest this strategy because each T(a) potentially reveals a different distribution of molecular states; for two-step analysis, a range of T(a) that brackets T(50) is optimal.

Animals↗

PubMed2003

Algorithms for computational solvent mapping of proteins.

Computational mapping methods place molecular probes (small molecules or functional groups) on a protein surface to identify the most favorable binding positions by calculating an interaction potential. We have developed a novel computational mapping program called CS-Map (computational solvent mapping of proteins), which differs from earlier mapping methods in three respects: (i) it initially moves the ligands on the protein surface toward regions with favorable electrostatics and desolvation, (ii) the final scoring potential accounts for desolvation, and (iii) the docked ligand positions are clustered, and the clusters are ranked on the basis of their average free energies. To understand the relative importance of these factors, we developed alternative algorithms that use the DOCK and GRAMM programs for the initial search. Because of the availability of experimental solvent mapping data, lysozyme and thermolysin are considered as test proteins. Both DOCK and GRAMM speed up the initial search, and the combined algorithms yield acceptable mapping results. However, the DOCK-based approaches place the consensus site farther from its experimentally determined position than CS-Map, primarily because of the lack of a solvation term in the initial search. The GRAMM-based program also finds the correct consensus site for thermolysin. We conclude that good sampling is the most important requirement for successful mapping, but accounting for desolvation and clustering of ligand positions also help to reduce the number of false positives.

2-Propanol↗

PubMed2000

The kinetic mechanism of the dual phosphorylation of the ATF2 transcription factor by p38 mitogen-activated protein (MAP) kinase alpha. Implications for signal/response profiles of MAP kinase pathways.

The mitogen-activated protein kinases (MAPKs) are a family of enzymes conserved among eukaryotes that regulate cellular activities in response to numerous external signals. They are the terminal component of a three-kinase cascade that is evolutionarily conserved and whose arrangement appears to offer considerable flexibility in encompassing the diverse biological situations for which they are employed. Although multistep protein phosphorylation within mitogen-activated protein kinase (MAPK) cascades can dramatically influence the sensitivity of signal propagation, an investigation of the mechanism of multisite phosphorylation by a MAPK has not been reported. Here we report a kinetic examination of the phosphorylation of Thr-69 and Thr-71 of the glutathione S-transferase fusion protein of the trans-activation domain of activating transcription factor-2 (GST-ATF2-(1-115)) by p38 MAPKalpha (p38alpha) as a model system for the phosphorylation of ATF2 by p38alpha. Our experiments demonstrated that GST-ATF2-(1-115) is phosphorylated in a two-step distributive mechanism, where p38alpha dissociates from GST-ATF2-(1-115) after the initial phosphorylation of either Thr-69 or Thr-71. Whereas p38alpha showed similar specificity for Thr-71 and Thr-69 in the unphosphorylated protein, it displayed a marked difference in specificity toward the mono-phosphoisomers. Phosphorylation of Thr-71 had no significant effect on the rate of Thr-69 phosphorylation, but Thr-69 phosphorylation reduced the specificity, k(cat)/K(M), of p38alpha for Thr-71 by approximately 40-fold. Computer simulation of the mechanism suggests that the activation of ATF2 by p38alpha in vivo is essentially Michaelian and provides insight into how the kinetics of a two-step distributive mechanism can be adapted to modulate effectively the sensitivity of a signal transduction pathway. This work also suggests that whereas MAPKs utilize docking interactions to bind substrates, they can be weak and transient in nature, providing just enough binding energy to promote the phosphorylation of a specific substrate.

Activating Transcription Factor 2↗

PubMed1999

Probing the binding of indolactam-V to protein kinase C through site-directed mutagenesis and computational docking simulations.

Protein kinase C (PKC) comprises a family of ubiquitous enzymes transducing signals by the lipophilic second messenger sn-1, 2-diacylglycerol (DAG). Teleocidin and its structurally simpler congener indolactam-V (ILV) bind to PKC with high affinity. In this paper, we report our computational docking studies on ILV binding to PKC using an automatic docking computer program, MCDOCK. In addition, we used site-directed mutagenesis to assess the quantitative contribution of crucial residues around the binding site of PKC to the binding affinity of ILV to PKC. On the basis of the docking studies, ILV binds to PKC in its cis-twist conformation and forms a number of optimal hydrogen bond interactions. In addition, the hydrophobic groups in ILV form "specific" hydrophobic interactions with side chains of a number of conserved hydrophobic residues in PKC. The predicted binding mode for ILV is entirely consistent with known structure-activity relationships and with our mutational analysis. Our mutational analysis establishes the quantitative contributions of a number of conserved residues to the binding of PKC to ILV. Taken together, our computational docking simulations and analysis by site-directed mutagenesis provide a clear understanding of the interaction between ILV and PKC and the structural basis for design of novel, high-affinity, and isozyme-selective PKC ligands.

Animals↗

PubMed2005

Q-SiteFinder: an energy-based method for the prediction of protein-ligand binding sites.

MOTIVATION: Identifying the location of ligand binding sites on a protein is of fundamental importance for a range of applications including molecular docking, de novo drug design and structural identification and comparison of functional sites. Here, we describe a new method of ligand binding site prediction called Q-SiteFinder. It uses the interaction energy between the protein and a simple van der Waals probe to locate energetically favourable binding sites. Energetically favourable probe sites are clustered according to their spatial proximity and clusters are then ranked according to the sum of interaction energies for sites within each cluster. RESULTS: There is at least one successful prediction in the top three predicted sites in 90% of proteins tested when using Q-SiteFinder. This success rate is higher than that of a commonly used pocket detection algorithm (Pocket-Finder) which uses geometric criteria. Additionally, Q-SiteFinder is twice as effective as Pocket-Finder in generating predicted sites that map accurately onto ligand coordinates. It also generates predicted sites with the lowest average volumes of the methods examined in this study. Unlike pocket detection, the volumes of the predicted sites appear to show relatively low dependence on protein volume and are similar in volume to the ligands they contain. Restricting the size of the pocket is important for reducing the search space required for docking and de novo drug design or site comparison. The method can be applied in structural genomics studies where protein binding sites remain uncharacterized since the 86% success rate for unbound proteins appears to be only slightly lower than that of ligand-bound proteins. AVAILABILITY: Both Q-SiteFinder and Pocket-Finder have been made available online at http://www.bioinformatics.leeds.ac.uk/qsitefinder and http://www.bioinformatics.leeds.ac.uk/pocketfinder

Algorithms↗

PubMed2003

Ligand-induced conformational changes: improved predictions of ligand binding conformations and affinities.

A computational docking strategy using multiple conformations of the target protein is discussed and evaluated. A series of low molecular weight, competitive, nonpeptide protein tyrosine phosphatase inhibitors are considered for which the x-ray crystallographic structures in complex with protein tyrosine phosphatase 1B (PTP1B) are known. To obtain a quantitative measure of the impact of conformational changes induced by the inhibitors, these were docked to the active site region of various structures of PTP1B using the docking program FlexX. Firstly, the inhibitors were docked to a PTP1B crystal structure cocrystallized with a hexapeptide. The estimated binding energies for various docking modes as well as the RMS differences between the docked compounds and the crystallographic structure were calculated. In this scenario the estimated binding energies were not predictive inasmuch as docking modes with low estimated binding energies corresponded to relatively large RMS differences when aligned with the corresponding crystal structure. Secondly, the inhibitors were docked to their parent protein structures in which they were cocrystallized. In this case, there was a good correlation between low predicted binding energy and a correct docking mode. Thirdly, to improve the predictability of the docking procedure in the general case, where only a single target protein structure is known, we evaluate an approach which takes possible protein side-chain conformational changes into account. Here, side chains exposed to the active site were considered in their allowed rotamer conformations and protein models containing all possible combinations of side-chain rotamers were generated. To evaluate which of these modeled active sites is the most likely binding site conformation for a certain inhibitor, the inhibitors were docked against all active site models. The receptor rotamer model corresponding to the lowest estimated binding energy is taken as the top candidate. Using this protocol, correct inhibitor binding modes could successfully be discriminated from proposed incorrect binding modes. Moreover, the ranking of the estimated ligand binding energies was in good agreement with experimentally observed binding affinities.

Binding Sites↗

PubMed2001

3-Aryl[1,2,4]triazino[4,3-a]benzimidazol-4(10H)-ones: a new class of selective A1 adenosine receptor antagonists.

Radioligand binding assays using bovine cortical membrane preparations and biochemical in vitro studies revealed that various 3-aryl[1,2,4]triazino[4,3-a]benzimidazol-4(10H)-one (ATBI) derivatives, previously reported by us as ligands of the central benzodiazepine receptor (BzR) (Primofiore, G.; et al. J. Med. Chem. 2000, 43, 96-102), behaved as antagonists at the A1 adenosine receptor (A1AR). Alkylation of the nitrogen at position 10 of the triazinobenzimidazole nucleus conferred selectivity for the A1AR vs the BzR. The most potent ligand of the ATBI series (10-methyl-3-phenyl[1,2,4]triazino[4,3-a]benzimidazol-4(10H)-one 12) displayed a Ki value of 63 nM at the A1AR without binding appreciably to the adenosine A2A and A3 nor to the benzodiazepine receptor. Pharmacophore-based modeling studies in which 12 was compared against a set of well-established A1AR antagonists suggested that three hydrogen bonding sites (HB1 acceptor, HB2 and HB3 donors) and three lipophilic pockets (L1, L2, and L3) might be available to antagonists within the A1AR binding cleft. According to the proposed pharmacophore scheme, the lead compound 12 engages interactions with the HB2 site (via the N2 nitrogen) as well as with the L2 and L3 sites (through the pendant and the fused benzene rings). The results of these studies prompted the replacement of the methyl with more lipophilic groups at the 10-position (to fill the putative L1 lipophilic pocket) as a strategy to improve A1AR affinity. Among the new compounds synthesized and tested, the 3,10-diphenyl[1,2,4]triazino[4,3-a]benzimidazol-4(10H)-one (23) was characterized by a Ki value of 18 nM which represents a 3.5-fold gain of A1AR affinity compared with the lead 12. A rhodopsin-based model of the bovine adenosine A1AR was built to highlight the binding mode of 23 and two well-known A1AR antagonists (III and VII) and to guide future lead optimization projects. In our docking simulations, 23 receives a hydrogen bond (via the N1 nitrogen) from the side chain of Asn247 (corresponding to the HB1 and HB2 sites) and fills the L1, L2, and L3 lipophilic pockets with the 10-phenyl, 3-phenyl, and fused benzene rings, respectively.

Amino Acid Sequence↗

PubMed2005

Design and synthesis of inositolphosphoglycan putative insulin mediators.

The binding modes of a series of molecules, containing the glucosamine (1-->6) myo-inositol structural motif, into the ATP binding site of the catalytic subunit of cAMP-dependent protein kinase (PKA) have been analysed using molecular docking. These calculations predict that the presence of a phosphate group at the non-reducing end in pseudodisaccharide and pseudotrisaccharide structures properly orientate the molecule into the binding site and that pseudotrisaccharide structures present the best shape complementarity. Therefore, pseudodisaccharides and pseudotrisaccharides have been synthesised from common intermediates using effective synthetic strategies. On the basis of this synthetic chemistry, the feasibility of constructing small pseudotrisaccharide libraries on solid-phase using the same intermediates has been explored. The results from the biological evaluation of these molecules provide additional support to an insulin-mediated signalling system which involves the intermediacy of inositolphosphoglycans as putative insulin mediators.

Adenosine Triphosphate↗

PubMed2004

Direct photoaffinity labeling by dolastatin 10 of the amino-terminal peptide of beta-tubulin containing cysteine 12.

Tubulin with bound [5-3H]dolastatin 10 was exposed to ultraviolet light, and 8-10% of the bound drug cross-linked to the protein, most of it specifically. The primary cross-link was to the peptide spanning amino acid residues 2-31 of beta-tubulin, but the specific amino acid could not be identified. Indirect studies indicated that cross-link formation occurred between cysteine 12 and the thiazole moiety of dolastatin 10. An equipotent analog of dolastatin 10, lacking the thiazole ring, did not form an ultraviolet light-induced cross-link to beta-tubulin. Preillumination of tubulin with ultraviolet light, known to induce cross-link formation between cysteine 12 and exchangeable site nucleotide, inhibited the binding of [5-3H]dolastatin 10 and cross-link formation more potently than it inhibited the binding of colchicine or vinblastine to tubulin. Conversely, binding of dolastatin 10 to tubulin inhibited formation of the cross-link between cysteine 12 and the exchangeable site nucleotide. Dithiothreitol inhibited formation of the beta-tubulin/dolastatin 10 cross-link but not the beta-tubulin/exchangeable site nucleotide cross-link. Modeling studies revealed a highly favored binding site for dolastatin 10 at the + end of beta-tubulin in proximity to the exchangeable site GDP. Computational docking of an energy-minimized dolastatin 10 conformation at this site placed the thiazole ring of dolastatin 10 8-9 A from the sulfur atom of cysteine 12. Dolastatin 15 and cryptophycin 1 could also be docked into positions that overlapped more extensively with the docked dolastatin 10 than with each other. This result was consistent with the observed binding properties of these peptides.

Animals↗

PubMed2006

Aggregation of cyclodextrins as an important factor to determine their complexation behavior.

Here, we report a study on the complexation behavior of carotenoids with cyclodextrins (CDs) using solubility experiments and molecular-modelling methods. Carotenoids are an important group of naturally occurring dyes found in vegetables and fruits. Their antioxidant property has initiated investigations on their possible use as drugs. However, carotenoids are lipophilic molecules with very little inherent aqueous solubility. Cyclodextrin complexation has been widely used in order to increase the potential applications of hydrophobic compounds. Thus, the aim of our investigation was to design carotenoids with enhanced water solubility by cyclodextrin complexation. Molecular modelling of carotenoid-cyclodextrin complexes with a 1 : 1 stoichiometry successfully explained the experimentally observed capability of beta-cyclodextrins (beta-CDs) to form complexes with carotenoids as opposed to alpha-cyclodextrins (alpha-CDs) and gamma-cyclodextrins (gamma-CDs). Furthermore, molecular-dynamics calculations revealed that the aggregation properties of CD derivatives significantly influence their complexation behavior. Our docking calculations showed that RAMEB (random methylated beta-CD) is the beta-CD derivative that possesses the lowest tendency to aggregate. Solubility experiments yielded the same results, namely, RAMEB complexes possess the best water solubility. Our results showed that complexation of a ligand not buried inside of the CD cavity is dependent on two factors: i) the geometry of the inclusion part of the complex; ii) the self-aggregation property of the CD itself. The lower affinity the CDs possess for self-aggregation, the more likely are they involved in interactions with carotenoids. These results suggest that self-aggregation of CDs should be considered as an important parameter determining complexation in general.

Antioxidants↗

PubMed2006

Computer-aided analysis of the interactions of glutamine synthetase with its inhibitors.

Mechanism of inhibition of glutamine synthetase (EC 6.3.1.2; GS) by phosphinothricin and its analogues was studied in some detail using molecular modeling methods. Among three possible conformations of phosphinothricin in the active site of GS, this compatible with binding mode of methionine sulfoximine, determined recently by crystallography, was found to be energetically favored. Basing on these results eleven inhibitors of GS were docked into its active site. Taking into consideration that phosphinothricin acts as suicide inhibitor, which is due to phosphorylation by the enzyme, seven of studied analogues were additionally analyzed in their phosphorylated forms. All the inhibitor-enzyme complexes were evaluated quantitatively by using eight scoring functions implemented in Insight and Sybyl program packages and significant correlation between the obtained scores and experimental pK(i) values was achieved. Computed surface charge distribution for five selected inhibitors in both free and phosphorylated forms and their comparison with electronic structure of enzymatic reaction transition state allowed us to determine important electronic features required to construct potent inhibitors of glutamine synthetase.

Aminobutyrates↗