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

PubMed2004

Mapping protein matrix cavities in human cytoglobin through Xe atom binding.

Cytoglobin is the fourth recognized globin type, almost ubiquitously distributed in human tissues; its function is still poorly understood. Cytoglobin displays a core region of about 150 residues, structurally related to hemoglobin and myoglobin, and two extra segments, about 20 residues each, at the N- and C-termini. The core region hosts a large apolar cavity, held to provide a ligand diffusion pathway to/from the heme, and/or ligand temporary docking sites. Here we report the crystal structure (2.4A resolution, R-factor 19.1%) of a human cytoglobin mutant bearing the CysB2(38) --> Ser and CysE9(83) --> Ser substitutions (CYGB*), treated under pressurized xenon. Three Xe atoms bind to the heme distal site region of CYGB* mapping the protein matrix apolar cavity. Despite the conserved globin fold, the cavity found in CYGB* is structured differently from those recognized to play a functional role in myoglobin, neuroglobin, truncated hemoglobins, and Cerebratulus lacteus mini-hemoglobin.

Binding Sites↗

PubMed2005

Active site binding modes of curcumin in HIV-1 protease and integrase.

Structure models for the interaction of curcumin with HIV-1 integrase (IN) and protease (PR) were investigated using computational docking. Curcumin was found to bind preferentially in similar ways to the active sites of both IN and PR. For IN, the binding site is formed by residues Asp64, His67, Thr66, Glu92, Thr93, Asp116, Ser119, Asn120, and Lys159. Docked curcumin contacts the catalytic residues adjacent to Asp116 and Asp64, and near the divalent metal (Mg2+). In the PR docking, the curcumin structure fitted well to the active site, interacting with residues Asp25, Asp29, Asp30, Gly27', Asp29', and Asp30'. The results suggest that o-hydroxyl and/or keto-enol structures are important for both IN and PR inhibitory actions. The symmetrical structure of curcumin seems to play an important role for binding to the PR protein, whereas the keto-enol and only one side of the terminal o-hydroxyl showed tight binding to the IN active site.

Binding Sites↗

PubMed1998

Clarification of the binding mode of teleocidin and benzolactams to the Cys2 domain of protein kinase Cdelta by synthesis of hydrophobically modified, teleocidin-mimicking benzolactams and computational docking simulation.

Phorbol esters (12-O-tetradecanoylphorbol 13-acetate; TPA) and teleocidins are known to be potent tumor promoters and to activate protein kinase C (PKC) by binding competitively to the enzyme. The relationship between the chemical structures and the activities of these compounds has attracted much attention because of the marked structural dissimilarities. The benzolactam 5, with an eight-membered lactam ring and benzene ring instead of the nine-membered lactam ring and indole ring of teleocidins, reproduces the active ring conformation and biological activities of teleocidins. Herein we describe the synthesis of benzolactams with hydrophobic substituents at various positions. Structure-activity data indicate that the existence of a hydrophobic region between C-2 and C-9 and the steric factor at C-8 play critical roles in the appearance of biological activities. We also computationally simulated the docking of teleocidin and the modified benzolactam molecules to the Cys2 domain structure observed in the crystalline complex of PKCdelta with phorbol 13-acetate. Teleocidin and benzolactams fitted well into the same cavity as phorbol 13-acetate. Of the three functional groups hydrogen-bonding to the protein, two hydrogen-bonded with protein atoms in common with phorbol 13-acetate, but the third one hydrogen-bonded with a different protein atom from that in the case of phorbol 13-acetate. The model explains well the remarkable difference in activity between 5 and its analogue having a bulky substituent at C-8.

Binding Sites↗

PubMed1983

Experimental and computer graphics simulation analyses of the DNA interaction of 1,8-bis-(2-diethylaminoethylamino)-anthracene-9,10-dione, a compound modelled on doxorubicin.

The crystal structure of the anthraquinone derivative 1,8-bis-(2-diethylaminoethylamino)-anthracene-9,10-dione has been established. This compound was prepared as a potential DNA-intercalating agent based on the proven intercalators doxorubicin and mitoxantrone. Its DNA-binding properties have been examined experimentally by spectroscopic, thermal denaturation and ccc-DNA unwinding techniques: the results are consistent with an intercalative mode of binding to DNA. Computer graphics stimulation of the intercalative docking of this compound into the self-complementary dimer of d(CpG) has provided a minimum energy geometrical arrangement for the bound drug in the intercalation site comparable to that for proflavine when intercalated into the same d(CpG) model system. Entry of the compound into the site can only occur via the major groove.

Animals↗

PubMed2005

AMBER force field implementation of the boronate function to simulate the inhibition of beta-lactamases by alkyl and aryl boronic acids.

Boronic acids are a very appealing class of serine proteases inhibitors whose rational design suffers, in spite of their therapeutic potential, from the lack of boron-related parameters in force fields commonly used for proteins. We introduced bonded, non-bonded and point charges in the MacroModel/Amber force field, as well as GB/SA solvation parameters, to model boronic acids as tetrahedral adducts formed after protease's serine Ogamma coordination. With the aim to check the implemented force field, flexible docking studies were performed on three crystallographic complexes of beta-lactamases with boronic acids that output the crystallographic conformation of the complexes as the global minimum energy structure. Although the used approach was basic, nevertheless the resultant force field seems to be efficient and suitable for the structure-based design of new boronic inhibitors of beta-lactamases.

Boronic Acids↗

PubMed2003

Markovian negentropies in bioinformatics. 1. A picture of footprints after the interaction of the HIV-1 Psi-RNA packaging region with drugs.

MOTIVATION: Many experts worldwide have highlighted the potential of RNA molecules as drug targets for the chemotherapeutic treatment of a range of diseases. In particular, the molecular pockets of RNA in the HIV-1 packaging region have been postulated as promising sites for antiviral action. The discovery of simpler methods to accurately represent drug-RNA interactions could therefore become an interesting and rapid way to generate models that are complementary to docking-based systems. RESULTS: The entropies of a vibrational Markov chain have been introduced here as physically meaningful descriptors for the local drug-nucleic acid complexes. A study of the interaction of the antibiotic Paromomycin with the packaging region of the RNA present in type-1 HIV has been carried out as an illustrative example of this approach. A linear discriminant function gave rise to excellent discrimination among 80.13% of interacting/non-interacting sites. More specifically, the model classified 36/45 nucleotides (80.0%) that interacted with paromomycin and, in addition, 85/106 (80.2%) footprinted (non-interacting) sites from the RNA viral sequence were recognized. The model showed a high Matthews' regression coefficient (C = 0.64). The Jackknife method was also used to assess the stability and predictability of the model by leaving out adenines, C, G, or U. Matthews' coefficients and overall accuracies for these approaches were between 0.55 and 0.68 and 75.8 and 82.7, respectively. On the other hand, a linear regression model predicted the local binding affinity constants between a specific nucleotide and the aforementioned antibiotic (R2 = 0.83,Q2 = 0.825). These kinds of models may play an important role either in the discovery of new anti-HIV compounds or in the elucidation of their mode of action. AVAILABILITY: On request from the corresponding author (humbertogd@cbq.uclv.edu.cu or humbertogd@navegalia.com).

Binding Sites↗

PubMed1998

Fragment-based modeling of NAD binding to the catalytic subunits of diphtheria and pertussis toxins.

We describe a novel application of a fragment-based ligand docking technique; similar methods are commonly applied to the de novo design of ligands for target protein binding sites. We have used several new flexible docking and superposition tools, as well as a more conventional rigid-body (fragment) docking method, to examine NAD binding to the catalytic subunits of diphtheria (DT) and pertussis (PT) toxins, and to propose a model of the NAD-PT complex. Docking simulations with the rigid NAD fragments adenine and nicotinamide revealed that the low-energy dockings clustered in three distinct sites on the two proteins. Two of the sites were common to both fragments and were related to the structure of NAD bound to DT in an obvious way; however, the adenine subsite of PT was shifted relative to that of DT. We chose adenine/nicotinamide pairs of PT dockings from these clusters and flexibly superimposed NAD onto these pairs. A Monte Carlo-based flexible docking procedure and energy minimization were used to refine the modeled NAD-PT complexes. The modeled complex accounts for the sequence and structural similarities between PT and DT and is consistent with many results that suggest the catalytic importance of certain residues. A possible functional role for the structural difference between the two complexes is discussed.

Amino Acid Sequence↗

PubMed2002

Inhibition of cathepsin G by 2-amino-3,1-benzoxazin-4-ones: kinetic investigations and docking studies.

A series of benzoxazinones was used to investigate the interaction of human cathepsin G with acyl-enzyme inhibitors. With respect to the primary specificity of cathepsin G, inhibitors with hydrophobic or basic residues at position 2 were included in the study. Parameters of the enzyme acylation and deacylation were determined by slow-binding kinetics in the presence of a chromogenic substrate. For selected inhibitors, the time course of the enzyme-catalyzed conversion of the inhibitors was followed. This approach was suitable to elucidate a rate-determining deacylation step. Docking simulations of the noncovalent enzyme-inhibitor complexes were performed and several clusters were analyzed for each inhibitor. The amino acids of the active site that participate in the binding of the inhibitors were determined. The arrangements in several clusters of an inhibitor were not uniform with respect to the orientation by which the inhibitor was bound in the S(1) pocket. Docking of the basic piperazino derivatives 6 and 10 indicated an interaction with Glu 226 at the bottom of the S(1) specificity pocket. The (N-methyl)benzylamino derivative 1 showed the strongest acylation rate (k(on)=1200 M(-1) s(-1)), which was attributed to a high extent of pseudo-productive orientations of the noncovalent preassociation complex.

Benzoxazoles↗

PubMed2006

The betaI/betaIII-tubulin isoforms and their complexes with antimitotic agents. Docking and molecular dynamics studies.

Both microtubule destabilizer and stabilizer agents are important molecules in anticancer therapy. In particular, paclitaxel has been demonstrated to be effective for the treatment of ovarian, breast, and nonsmall cell lung carcinomas. It has been shown that emergence of resistance against this agent correlates with an increase in the relative abundance of tubulin isoform betaIII and that the more recently discovered IDN5390 can be effectively used once resistance has emerged. In this paper, we analyze the binding modes of these antimitotic agents to type I and III isoforms of beta-tubulin by computational methods. Our results are able to provide a molecular explanation of the experimental data. Using the same protocol, we could also show that no preference for any of the two isoforms can be detected for epothilone A, a potentially very interesting drug for which no data about the emergence of resistance is currently available. Our analysis provides structural insights about the recognition mode and the stabilization mechanism of these antimitotic agents and provides useful suggestions for the design of more potent and selective antimitotic agents.

Algorithms↗

PubMed1997

Simulation of alternative binding modes in a structure-based QSAR study of HIV-1 protease inhibitors.

We have used a published set of inhibitors of HIV-1 protease to build a COMBINE-type structure-based QSAR model with good predictive ability (r2 = 0.90, q2 = 0.69). Since the compounds in the training series exhibit most of their structural variability on one-half of the pseudosymmetrical binding cavity and only one binding orientation was explored for each molecule, the model describes mainly the effect of the structural changes on interactions involving only one-half of the binding cavity (pockets S1' and S2'). Thus, the model cannot be expected to give accurate predictions for new compounds exhibiting structural variation in both halves. The model does in fact show a tendency to underpredict slightly the biological activity of the molecules in the external test set. In an attempt to improve the quality of the model, both possible orientations of the ligands are now considered so that structural variation takes place in all binding pockets. One possibility would have been to build an additional set of complexes with the inhibitors docked in a reversed orientation. The alternative we have explored, however, consists of manipulating the data matrix describing the interaction energies so that each row is duplicated and the order of the variables in the duplicated rows is swapped between subunits. This simple approach has produced a new model that is similar in quality to the original model (r2 = 0.89, q2 = 0.64) but lacks the tendency to underpredict the activity of the compounds in the external set. Moreover, since equivalent residues are assigned equivalent weights, the model is insensitive to ligand orientation and is easier to interpret.

Artificial Intelligence↗

PubMed2001

Antifreeze protein from shorthorn sculpin: identification of the ice-binding surface.

Shorthorn sculpins, Myoxocephalus scorpius, are protected from freezing in icy seawater by alanine-rich, alpha-helical antifreeze proteins (AFPs). The major serum isoform (SS-8) has been reisolated and analyzed to establish its correct sequence. Over most of its length, this 42 amino acid protein is predicted to be an amphipathic alpha-helix with one face entirely composed of Ala residues. The other side of the helix, which is more heterogeneous and hydrophilic, contains several Lys. Computer simulations had suggested previously that these Lys residues were involved in binding of the peptide to the [11-20] plane of ice in the <-1102> direction. To test this hypothesis, a series of SS-8 variants were generated with single Ala to Lys substitutions at various points around the helix. All of the peptides retained significant alpha-helicity and remained as monomers in solution. Substitutions on the hydrophilic helix face at position 16, 19, or 22 had no obvious effect, but those on the adjacent Ala-rich surface at positions 17, 21, and 25 abolished antifreeze activity. These results, with support from our own modeling and docking studies, show that the helix interacts with the ice surface via the conserved alanine face, and lend support to the emerging idea that the interaction of fish AFPs with ice involves appreciable hydrophobic interactions. Furthermore, our modeling suggests a new N terminus cap structure, which helps to stabilize the helix, whereas the role of the lysines on the hydrophilic face may be to enhance solubility of the protein.

Alanine↗

PubMed2006

Modeling subtype-selective agonists binding with alpha4beta2 and alpha7 nicotinic acetylcholine receptors: effects of local binding and long-range electrostatic interactions.

The subtype-selective binding of 14 representative agonists with alpha4beta2 and alpha7 nicotinic acetylcholine receptors (nAChRs) has been studied by performing homology modeling, molecular docking, geometry optimizations, and microscopic and phenomenological binding free energy calculations. All of the computational results demonstrate that the subtype selectivity of the agonists binding with alpha4beta2 and alpha7 7 nAChRs is affected by both local binding and long-range electrostatic interactions between the receptors and the protonated structures of the agonists. The effects of the long-range electrostatic interactions are mainly due to the distinct difference in the net charge of the ligand-binding domain between the two nAChR subtypes. For the alpha4beta2-selective agonists examined, the microscopic binding modes with the alpha4beta2 nAChR are very similar to the corresponding modes with the alpha7 nAChR, and therefore, the subtype selectivity of these agonists binding with alpha4beta2 and alpha7 nAChRs is dominated by the long-range electrostatic interactions. For the alpha7-selective agonists, their microscopic binding modes with the alpha7 nAChR are remarkably different from those with the alpha4beta2 nAChR so that the local binding (including the hydrogen bonding and cation-pi interactions) with the alpha7 nAChR is much stronger than that with the alpha4beta2 nAChR. The calculated phenomenological binding free energies are in good agreement with available experimental data for the relative binding free energies concerning the subtype selectivity of agonists binding with the two different nAChR subtypes. The fundamental insights obtained in the present study should be valuable for future rational design of potential therapeutic agents targeted to specific nAChR subtypes.

Amino Acid Sequence↗

PubMed1994

Rigid-body docking with mutant constraints of influenza hemagglutinin with antibody HC19.

An automatic docking algorithm has been applied to the modeling of the complex between hemagglutinin from influenza virus and the Fab fragment of a monoclonal antibody raised against this antigen. We have introduced here the use of biochemical information provided by mutants of hemagglutinin. The docking procedure finds a small number of candidate solutions where three sites of escape mutations are buried and form hydrogen bonds in the interface. The localization of the epitope is improved by additional biochemical data about mutants that do not affect antibody binding. Five candidate solutions with low energy, reasonably well-packed interfaces, and six to ten hydrogen bonds are compatible with mutant information. One of the five stands out as generally better than the others from these points of views.

Algorithms↗

PubMed2003

A further investigation of the cytochrome b5-cytochrome c complex.

The interaction of reduced rabbit cytochrome b(5) with reduced yeast iso-1 cytochrome c has been studied through the analysis of (1)H-(15)N HSQC spectra, of (15)N longitudinal ( R(1)) and transverse ( R(2)) relaxation rates, and of the solvent exchange rates of protein backbone amides. For the first time, the adduct has been investigated also from the cytochrome c side. The analysis of the NMR data was integrated with docking calculations. The result is that cytochrome b(5) has two negative patches capable of interacting with a single positive surface area of cytochrome c. At low protein concentrations and in equimolar mixture, two different 1:1 adducts are formed. At high concentration and/or with excess cytochrome c, a 2:1 adduct is formed. All the species are in fast exchange on the scale of differences in chemical shift. By comparison with literature data, it appears that the structure of one 1:1 adduct changes with the origin or primary sequence of cytochrome b(5).

Animals↗

PubMed2005

Docking prediction using biological information, ZDOCK sampling technique, and clustering guided by the DFIRE statistical energy function.

We entered the CAPRI experiment during the middle of Round 4 and have submitted predictions for all 6 targets released since then. We used the following procedures for docking prediction: (1) the identification of possible binding region(s) of a target based on known biological information, (2) rigid-body sampling around the binding region(s) by using the docking program ZDOCK, (3) ranking of the sampled complex conformations by employing the DFIRE-based statistical energy function, (4) clustering based on pairwise root-mean-square distance and the DFIRE energy, and (5) manual inspection and relaxation of the side-chain conformations of the top-ranked structures by geometric constraint. Reasonable predictions were made for 4 of the 6 targets. The best fraction of native contacts within the top 10 models are 89.1% for Target 12, 54.3% for Target 13, 29.3% for Target 14, and 94.1% for Target 18. The origin of successes and failures is discussed. .

Algorithms↗

PubMed1997

Modeling of the three-dimensional structure of the human melanocortin 1 receptor, using an automated method and docking of a rigid cyclic melanocyte-stimulating hormone core peptide.

A model is presented of the melanocortin 1 receptor (MC1R), constructed by use of an unbiased, objective method. The model is created directly from data derived from multiple sequence analysis, a low-resolution EM-projection map of rhodopsin, and the approximate membrane thickness. The model agrees well with available data concerning natural mutations of MC1Rs occurring in different species. A model is also presented of the most rigid ligand for this receptor, the cyclic pentapeptide cHFRWG, shown docked in the receptor model. The receptor-ligand complex model agrees well with available experimental data. The ligand is located between transmembrane region 1 (TM1), TM2, TM3, TM6, and TM7 of the receptor. Multiple interactions occur between ligand and receptor, including interactions with Leu-48 (TM1), Ser-52 (TM1), Glu-55 (TM1), Asn-91 (TM2), Glu-94 (TM2), Thr-95 (TM2) Ile-98 (TM2), Asp-121 (TM3), Thr-124 (TM3), Phe-257 (TM6), Phe-283 (TM7), Asn-290 (TM7), and Asp-294 (TM7) of the receptor.

Amino Acid Sequence↗

PubMed2005

A structural model of an amyloid protofilament of transthyretin.

A docking-and-alignment protocol was devised in order to build amyloid protofilaments of Transthyretin (TTR), starting from partially disrupted TTR monomeric subunits and based on experimentally available information. The docking approach is driven by a combination of shape complementarity and energetic criteria, and uses constraints derived from experimental data obtained for the fibrillar state. The dimeric structures obtained were then subjected to an alignment scheme followed by clustering analysis, producing a collection of protofilaments with distinct geometric properties. The selected protofilament model presented here does agree with known experimental data and general amyloid properties; it is formed by two extended continuous beta-sheets with the beta-strands perpendicular to the main axis of the protofilament and a helical twist with a period of approximately 48 beta-strands. This TTR proto-filament model may be an important step in the understanding of the molecular mechanisms of TTR aggregation, as well as, a valuable instrument in drug design strategies against amyloid diseases.

Amyloid↗

PubMed1995

Identification of agonist binding sites of vasopressin and oxytocin receptors.

The present study aims at delineating residues in the vasopressin/oxytocin receptor family responsible for the high affinity binding of the hormone. Therefore, we have constructed a computer-generated 3 dimensional model of the rat V1a vasopressin receptor subtype which allowed us to propose residues likely to be involved in agonist binding. Among these residues, several are highly conserved in the receptor family. They were selected for site-directed mutagenesis on the basis of putative direct interaction with bound ligands. The present model and experimental results led us to conclude that the hormone is docked in a pocket completely buried in the transmembrane core of the receptor. Large polar residues, such as glutamine and lysine, located in transmembrane regions 2,3,4 and 6 are involved in the binding of the neurohypophysial hormone. Since all the mutated residues are highly conserved in AVP and OT receptors, we propose that the agonist binding site is similar in all members of the receptor family; only minor changes were found in antagonist potencies, suggesting that agonist and antagonist binding sites do not completely overlap.

Amino Acid Sequence↗