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 559 records · Page 31Linked to original sources

PubMed1999

Exploring the conformational diversity of loops on conserved frameworks.

Loops are structurally variable regions, but the secondary structural elements bracing loops are often conserved. Motifs with similar secondary structures exist in the same and different protein families. In this study, we made an all-PDB-based analysis and produced 495 motif families accessible from the Internet. Every motif family contains some variable loops spanning a common framework (a pair of secondary structures). The diversity of loops and the convergence of frameworks were examined. In addition, we also identified 119 loops with conformational changes in different PDB files. These materials can give some directions for functional loop design and flexible docking.

Computer Simulation↗

PubMed2003

Lipid bilayer simulations of CXCR4 with inverse agonists and weak partial agonists.

CXCR4 is a G protein-coupled receptor (GPCR) that has multiple critical functions in normal and pathologic physiology that include regulation of the metastatic behavior of mammary carcinoma, and utilization as a coreceptor for infection by T-tropic strains of human immunodeficiency virus-1. Molecular dynamic simulations of the rhodopsin-based homology model of CXCR4 were performed in a solvated lipid bilayer to reproduce the microenvironment of this integral membrane protein. The amino acids in CXCR4 necessary for interaction with an inverse agonist, T140, and a weak partial agonist, AMD3100, identified by alanine scanning mutants, were spatially consistent when computationally docked. Whereas T140 binds residues in extracellular domains and regions of the hydrophobic core proximal to the cell surface, amino acids in the central hydrophobic core are critical to binding of AMD3100. The physical localization of T140 binding to CXCR4 by biochemical analyses corroborated the molecular and computational approaches. The structural basis for the interaction of T140 and AMD3100 with CXCR4 confirms that the mechanisms used by these agents are different. This complementary utilization of molecular, physical, and computation analysis provides a powerful approach to elucidate GPCR conformation.

Amino Acid Substitution↗

PubMed2000

Oligomerization reduces heparin affinity but enhances receptor binding of fibroblast growth factor 2.

The biological response of cells to fibroblast growth factors (FGFs) depends on heparan sulphate glycosaminoglycans sharing particular structural motifs. Heparin induced FGF dimerization has been suggested to mediate receptor dimerization and activation. Here we demonstrate that heparin-derived oligosaccharides that promote receptor binding and activation specifically induce the dimerization of basic FGF (FGF2). These heparin-induced dimers of FGF2 acquire high affinity for receptor binding and are biologically active. Using biotinylated FGF2 bound to immobilized streptavidin gradually saturated with biotin, enabled a quantitative analysis of heparin-dependent and heparin-independent FGF2 monomers and oligomers. Streptavidin induced FGF2 dimers bind and activate FGF receptors only in the presence of heparin. An excess of streptavidin, forcing biotin-FGF2 into monomers, reduces receptor binding and blocks FGF-dependent cell proliferation. All these suggest predominant receptor binding and activation by heparin associated FGF2 oligomers. Unexpectedly, heparin induced dimers and higher order oligomers lose most of their affinity towards heparin. Direct binding of soluble FGF receptors (FGFRs) to either monomers or dimers of FGF2, immobilized on heparin, confirm the preferred association of FGFRs with dimers of FGF2. Computerized molecular docking predicts a cis-oriented FGF2 dimer, stabilized by heparin, which complies with all the experimental data.

Animals↗

PubMed2006

Crystal structure of neoculin: insights into its sweetness and taste-modifying activity.

Although the majority of sweet compounds are of low molecular mass, several proteins are known to elicit sweet taste responses in humans. The fruit of Curculigo latifolia contains a heterodimeric protein, neoculin, which has both sweetness and a taste-modifying activity that converts sourness to sweetness. Here, we report the crystal structure of neoculin at 2.76A resolution. This is the first well-defined tertiary structure of a taste-modifying protein of this kind. The overall structure is quite similar to those of monocot mannose-binding lectins. However, crucial topological differences are observed in the C-terminal regions of both subunits. In both subunits of neoculin, the C-terminal tails turn up to form loops fixed by inter-subunit disulfide bonds that are not observed in the lectins. Indeed, the corresponding regions of the lectins stretch straight over the surface of another subunit. Such a C-terminal structural feature as is observed in neoculin results in a decrease in subunit-subunit interactions. Moreover, distribution of electrostatic potential on the surface of neoculin is unique and significantly different from those of the lectins, particularly in the basic subunit (NBS). We have found that there is a large cluster composed of six basic residues on the surface of NBS, and speculate that it might be involved in the elicitation of sweetness and/or taste-modifying activity of neoculin. Molecular dynamics simulation based on the crystallography results suggests that neoculin may adopt a widely "open" conformation at acidic pH, while unprotonated neoculin at neutral pH is in a "closed" conformation. Based on these simulations and the generation of a docking model between neoculin and the sweet-taste receptor, T1R2-T1R3, we propose the hypothesis that neoculin is in dynamic equilibrium between open and closed states, and that the addition of an acid shifts the equilibrium to the open state, allowing ligand-receptor interaction.

Amino Acid Sequence↗

PubMed1999

Molecular model of the interaction between nimesulide and human cyclooxygenase-2.

The cyclooxygenase-2 (COX-2) isoenzyme is a key target for COX-2-selective non-steroidal anti-inflammatory drugs (NSAIDs). An important difference in binding of nimesulide compared with non-selective NSAIDs appears to involve the amino acid at position 523 of the enzyme. Replacement of valine with isoleucine at this position provides access to a binding site that is larger in COX-2 than in COX-1. Nimesulide appears to exploit this enlarged binding site for establishing a number of favourable contacts with the enzyme that lead to selective inhibition of COX-2. We made these conclusions from a three-dimensional molecular model of the active site of human COX-2, constructed using the X-ray coordinates of COX-1 from sheep seminal vesicles and COX-2 from mouse fibroblasts as templates, with the aid of sequence alignment methods and molecular modelling techniques. The resulting model was refined, and the active site was probed for regions of steric and electrostatic complementarity for ligand binding. Docking studies were then undertaken with many different nimesulide conformers, a family of which could establish very favourable interactions with the NSAID binding site of human COX-2 by exploiting the extra space made available by the isoleucine/valine replacement. The stability of the resulting complexes was studied by simulating molecular dynamics.

Cyclooxygenase 1↗

PubMed2003

Molecular modeling of p-chlorophenoxyacetic acid binding to the CLC-0 channel.

Molecular simulation techniques were applied to predict the interaction of the CLC-0 Cl(-) channel and the channel-blocking molecule p-chlorophenoxyacetic acid (CPA). A three-dimensional model of the CLC-0 channel was constructed on the basis of the homology with the bacterial Cl(-) channel StCLC, the structure of which has been solved by X-ray crystallography. Docking of the CPA molecule was obtained by using a geometric recognition algorithm, yielding 5000 possible conformations. By restraining the simulation to those conformations in which CPA is near the intracellular mouth of the channel, the CPA-protein complex models were reduced to three sets of conformations, which are interconvertible within 2 ns when molecular dynamics is applied to the system. Point mutations of CLC-0 at three different positions predicted to interact with CPA in these configurations did, however, not greatly alter CPA inhibition, suggesting a deeper final binding location. In the model, binding of CPA to a more internal position in the ionic pathway was obtained by applying a constant force vector to CPA, pushing it toward the center of the channel. This technique allowed us to outline the possible intrachannel pathway of CPA and to describe qualitatively the binding sites and energy barriers of this pathway. The consistency of the obtained models and the experimental data indicates that the CLC-0-CPA complex model is reasonable and can be used in further biological studies, such as rational design of blocking agents of and mutagenesis of CLC Cl(-) channels.

2,4-Dichlorophenoxyacetic Acid↗

PubMed1996

Investigation of the conformational behaviour of permethylated cyclodextrins by molecular modelling.

Conformations of manually built native and permethylated alpha-, beta-, and gamma-cyclodextrins (CD) were investigated using various computer assisted molecular modelling methods. Calculations were carried by applying the MM+ and the Tripos force field. The influences of atomic charges on the macrocyclic conformations during the optimization procedure were analyzed. The permethylation of hydroxyl groups of cyclodextrins changes bond and torsion angles between the glucose monomers and of the primary substituents. A method to determine the diameters of the cyclodextrin cavity by a modelling approach is described. It is shown that due to permethylation the larger cavity opening is increased and the primary substituents are canted outwards. As a consequence, the torus shape of the molecule changes, which is an important feature for docking and fitting studies.

Algorithms↗

PubMed1999

Ligand solvation in molecular docking.

Solvation plays an important role in ligand-protein association and has a strong impact on comparisons of binding energies for dissimilar molecules. When databases of such molecules are screened for complementarity to receptors of known structure, as often occurs in structure-based inhibitor discovery, failure to consider ligand solvation often leads to putative ligands that are too highly charged or too large. To correct for the different charge states and sizes of the ligands, we calculated electrostatic and non-polar solvation free energies for molecules in a widely used molecular database, the Available Chemicals Directory (ACD). A modified Born equation treatment was used to calculate the electrostatic component of ligand solvation. The non-polar component of ligand solvation was calculated based on the surface area of the ligand and parameters derived from the hydration energies of apolar ligands. These solvation energies were subtracted from the ligand-receptor interaction energies. We tested the usefulness of these corrections by screening the ACD for molecules that complemented three proteins of known structure, using a molecular docking program. Correcting for ligand solvation improved the rankings of known ligands and discriminated against molecules with inappropriate charge states and sizes.

Aniline Compounds↗

PubMed2000

Deciphering common failures in molecular docking of ligand-protein complexes.

Common failures in predicting crystal structures of ligand-protein complexes are investigated for three ligand-protein systems by a combined thermodynamic and kinetic analysis of the binding energy landscapes. Misdocked predictions in ligand-protein docking are classified as 'soft' and 'hard' failures. While a soft failure arises when the search algorithm is unable to find the global energy minimum corresponding to the crystal structure, a hard failure results from a flaw of the energy function to qualify the crystal structure as the predicted lowest energy conformation in docking simulations. We find that neither the determination of a single structure with the lowest energy nor finding the most common binding mode is sufficient to predict crystal structures of the complexes, which belong to the category of hard failures. In a proposed hierarchical approach, structural similarity clustering of the conformations, generated from equilibrium simulations with the simplified energy function, is followed by energy refinement with the AMBER force field. This protocol, that involves a hierarchy of energy functions, resolves some common failures in ligand-protein docking and detects crystallographic binding modes that were not found during docking simulations.

Crystallography↗

PubMed2003

Biomimetic recognition and immunochemical assay of ligand binding to collagen.

PURPOSE: The objectives of this investigation were to analyze the adhesive interactions between type 1 collagen and two hydrophilic monomer primers (ligands) for dentin bonding through computer modeling and a novel immunochemical binding assay method. MATERIALS AND METHODS: The hydrophilic monomers studied included 2-hydroxyethyl methacrylate (HEMA) and 2-acryloyloxyethyl phosphate (PA). For computer modeling studies, a triple helical model structure of collagen fibril extracted from a public domain Protein Data Bank was used. The ligand conformations were modeled and optimized by Sybyl, and their interactions with the triple helical collagen structure in a solvent environment of water were simulated by AutoDock software. The effect of ligand binding to subsequent monoclonal antibody binding was also studied using an immunochemical binding assay method developed by us. The computer simulation results and binding assay results were analyzed for relationships pertinent to collagen-ligand interaction in bonding of dentin primers to demineralized dentin. RESULTS: Computer docking results indicate that ligand binding to collagen is favored to occur at cavity sites on the collagen molecular surfaces, where steric and electrostatic effects may play a critical role in mediating van der Waals and Coulombic interactions between ligand and receptor molecules. Prior ligand exposure of collagen reduces subsequent antibody binding during in vitro experiments. Differences in antibody binding were observed both as a function of ligand type and its concentration. CONCLUSION: Both steric complementarity and electrostatic complementarity conditions were observed in the docking site selection. Under conditions of ligand binding, antibody binding is diminished as a function of ligand structure and its concentration. The results suggest that approaches combining computer modeling and in vitro binding assay methods are powerful tools in evaluating dentinal adhesion at the atomic and/or molecular level.

Acrylates↗

PubMed2003

Molecular modelling probes: docking and scoring.

A general introduction to molecular modelling techniques in the area of protein-ligand interactions is given. Methods covered range from binding-site analysis to statistical treatment of sets of ligands. The main focus of this paper is on docking and scoring. After an outline of the main concepts, two specific application examples are given.

Algorithms↗

PubMed2003

Use of anti-neoepitope antibodies for the analysis of degradative events in cartilage and the molecular basis for neoepitope specificity.

Degradation of the cartilage proteoglycan, aggrecan, is an essential aspect of normal growth and development, and of joint pathology. The roles of different proteolytic enzymes in this process can be determined from the sites of cleavage in the aggrecan core protein, which generates novel termini (neoepitopes). Antibodies specific for the different neoepitopes generated by such cleavage events provide powerful tools with which to analyse these processes. The same approach can be used to differentiate the processed, active forms of proteases from their inactive pro-forms. Since the proteolytic processing of these enzymes requires the removal of the inhibitory pro-region, it also results in the generation of N-terminal neoepitopes. Using the newborn rat long bone as a model system, it was shown that the active form of ADAMTS-4 [ADAM (a disintegrin and metalloproteinase) with thrombospondin motifs-4], but not ADAMTS-5, co-localizes with the aggrecan cleavage neoepitopes known to be produced by this metalloproteinase. Thus, in long bone growth, aggrecan turnover seems to be dependent on ADAMTS-4 activity. To demonstrate the molecular basis of the specificity of anti-neoepitope antibodies, the Fv region of a monoclonal antibody specific for a neoepitope generated by the ADAMTS-4-mediated cleavage of aggrecan has been modelled and the binding of the peptide epitope simulated. In the docked structure, the N-terminus of the peptide antigen is clearly buried in the binding-site cavity. The absence of an open cleft makes it impossible for the intact substrate to pass through the binding site, providing a rationale for the specificity of this class of antibodies.

Amino Acid Sequence↗

PubMed1993

A computer model of the interleukin-4/receptor complex.

Interleukin-4 is a member of the cytokine family, a group of related messenger proteins which collectively help to moderate and control the immune response. It is believed that the folding topology of the beta-sheets of the interleukin-4 receptor (IL4R) is the same as that seen in the crystal structure of CD4. Although the sequence identity is low, homology modeling techniques have been used to model the IL4R structure from CD4. Refinement by molecular dynamics leads to a suggested structure which has been docked to interleukin-4 (IL4). Several residues of apparent importance for binding are identified.

Amino Acid Sequence↗

PubMed1999

Molecular dynamics simulations of P450 BM3--examination of substrate-induced conformational change.

Cytochrome P450 BM3, of bacterial origin, is one of only five isozymes of the ubiquitous family of over 400 metabolizing heme proteins with a known crystal structure and only one of two with both substrate-free and substrate-bound forms determined. P450 BM3 is of particular interest since it has a similar function and similar substrates as mammalian P450s particularly of the 4A subfamily. Thus, the extent to which the substrate-free form of P450 BM3 undergoes a conformational change upon binding of a typical fatty acid substrate, palmitoleic acid, has been the subject of recent active experimental effort. Surprisingly, direct examination of the substrate-free (pdb2hpd.ent and pdb2bmh.ent) and substrate-bound (pdb1fag.ent) forms do not provide a clear answer to this question. The main reason for this ambiguity is that the two substrate-free monomers reported in the crystal structures themselves have significantly different conformations from each other, one with a more open substrate-access channel than the other. Since there is no way to tell to which substrate-free form the substrate binds, the effect of substrate binding cannot be deduced directly from comparisons of the experimental substrate-bound and substrate-free forms. The computational studies reported here have been designed to more robustly establish the effect of substrate binding on this isozyme. Specifically, molecular dynamics simulations were performed for each of the two substrate-free forms found in the asymmetric unit of the X-ray structure and for the two corresponding substrate-bound forms, constructed by docking palmitloeic acid into each of them. Comparisons of the results showed that palmitoleic acid binding had little effect on the conformation of the more closed substrate-free form of P450 BM3. By contrast, in the more open substrate-free form, this same substrate induced a closing of the entrance to the substrate-binding channel. The MD averaged structure of these two complexes obtained from docking of pamitoleic acid into the two asymmetric units of the substrate-free form were also compared to that obtained starting with the X-ray structure of the substrate-bound form. These results taken together led to the conclusion that, if indeed the substrate induces conformational changes in P450 BM3, the mouth of the substrate-access channel first closes down in response to the presence of the substrate, followed by rotation of the F-G domain to further optimize the P450 BM3-substrate interaction that would occur at a later stage.

Bacterial Proteins↗

PubMed2006

Retrotransposon sequence variation in four asexual plant species.

Transposable elements (TEs) can be viewed as genetic parasites that persist in populations due to their capacity for increase in copy number and the inefficacy of selection against them. A corollary of this hypothesis is that TEs are more likely to spread within sexual populations and be eliminated or inactivated within asexual populations. While previous work with animals has shown that asexual taxa may contain less TE diversity than sexual taxa, comparable work with plants has been lacking. Here we report the results of a study of Ty1/copia, Ty3/gypsy, and LINE-like retroelement diversity in four asexual plant species. Retroelement-like sequences, with a high degree of conservation both within and between species, were isolated from all four species. The sequences correspond to several previously annotated retroelement subfamilies. They also exhibit a pattern of nucleotide substitution characterized by an excess of synonymous substitutions, suggestive of a history of purifying selection. These findings were compared with retroelement sequence evolution in sexual plant taxa. One likely explanation for the discovery of conserved TE sequences in the genomes of these asexual taxa is simply that asexuality within these taxa evolved relatively recently, such that the loss and breakdown of TEs is not yet detectable through analysis of sequence diversity. This explanation is examined by conducting stochastic simulation of TE evolution and by using published information to infer rough estimates of the ages of asexual taxa.

Amino Acid Sequence↗

PubMed2006

Binding mode analysis of the NADH cofactor in nitric oxide reductase: a theoretical study.

Nitric oxide reductase (Nor), a member of the cytochrome P450 superfamily, takes part in the denitrification process of fungi by reducing NO to N(2)O. Evidence indicates that Nor binds NADH, source of the reducing equivalents of the reaction, within its large hydrophilic ligand binding cavity on the distal side of heme and receives electrons directly from the cofactor. Here we present a binding mode analysis of the structure of the Nor-NO-NADH complex, performed in three steps. The NADH cofactor was first docked into the enzyme interior using the Monte Carlo multiple minimum algorithm, refined by low-mode conformational search and the final arrangement was obtained in a 5ns NPT molecular dynamics simulation. The NADH cofactor, in our results, is positioned--by Arg174, Lys291, Asp393 and several water molecules--within reactive distance of the NO binding spot suggesting a direct hydride shift mechanism between the two. The catalytically required water molecule is captured by NADH and the cofactor not only retains the suggested H-bonded proton transfer pathway between the active site and the solvent, but provides structural restraint for its members. We also found that direct interaction is formed between the cofactor and propionate A of the heme group, which flips from the proximal to the distal side of the heme plane in order to become an H-bonding partner of NADH. The role of Arg64 and Glu71 was suggested to be fixing the residues of the translocated helix B' to their new position.

Binding Sites↗

PubMed2006

Homology modeling and molecular dynamics study of chorismate synthase from Shigella flexneri.

Shigellosis is a major public health problem in many developing countries. Antibiotic therapy can reduce the severity of the dysentery and prevent potentially lethal complication. However, owing to the increased resistance to most of the widely used and inexpensive antibiotics, there is an urgent need for new antibacterial agents, particularly those that act on novel targets. Chorismate synthase (CS) is a key enzyme in the shikimic acid pathway, which is essential for the synthesis of aromatic amino acids in bacteria. As an anti-bacterial drug target, CS has been well validated. A homology model of Shigella-CS with the flavin mononucleotide (FMN) binding was constructed using the crystal structure of CS from other species. The substrate 5-enolpyruvylshikimate 3-phosphate (EPSP) was subsequently docked into the active site based on previous theoretical studies. Molecular dynamics (MD) was used to refine the starting ternary model. The model was well conserved during the 1.8 ns MD simulation with the equilibrium root mean square deviation (RMSD) value of 3.5 angstrom. The substrate binding energy was calculated and the electrostatic energy was found to be the most important term for binding. Decomposition of binding energies revealed that R129, R125, R327, R134 and R48 are important residues involved in substrate binding, which is useful for further site-directed mutagenesis experiments. In the absence of crystal structure, our study provides an early insight into the structure of CS from Shigella flexneri and its binding to the substrate and cofactor, thus facilitating the inhibitor design.

Amino Acid Sequence↗

PubMed1987

Molecular dynamics of a cytochrome c-cytochrome b5 electron transfer complex.

Cytochrome c and cytochrome b5 form an electrostatically associated electron transfer complex. Computer models of this and related complexes that were generated by docking the x-ray structures of the individual proteins have provided insight into the specificity and mechanism of electron transfer reactions. Previous static modeling studies were extended by molecular dynamics simulations of a cytochrome c-cytochrome b5 intermolecular complex. The simulations indicate that electrostatic interactions at the molecular interface results in a flexible association complex that samples alternative interheme geometries and molecular conformations. Many of these transient geometries appear to be more favorable for electron transfer than those formed in the initial model complex. Of particular interest is a conformational change that occurred in phenylalanine 82 of cytochrome c that allowed the phenyl side chain to bridge the two cytochrome heme groups.

Computer Graphics↗