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 991 records · Page 55Linked to original sources

PubMed2004

Inhibition of Golgi mannosidase II with mannostatin A analogues: synthesis, biological evaluation, and structure-activity relationship studies.

Mannostatin and aminocyclopentitetrol analogues with various substitutions at the amino function were synthesized. These compounds were tested as inhibitors of human Golgi and lysosomal alpha-mannosidases. Modification of the amine of mannostatin had only marginal effects, whereas similar modifications of aminocyclopentitetrol led to significantly improved inhibitors. Ab initio calculations and molecular docking studies were employed to rationalize the results. It was found that mannostatin and aminocyclopentitretrol could bind to Golgi alpha-mannosidase II in a similar mode to that of the known inhibitor swainsonine. However, due to the flexibility of the five-membered rings of these compounds, additional low-energy binding modes could be adopted. These binding modes may be relevant for the improved activities of the benzyl-substituted compounds. The thiomethyl moiety of mannostatin was predicted to make favorable hydrophobic interactions with Arg228 and Tyr727 that would possibly account for its greater inhibitory activity.

Amines↗

PubMed2001

[1,2,4]Triazino[4,3-a]benzimidazole acetic acid derivatives: a new class of selective aldose reductase inhibitors.

Acetic acid derivatives of [1,2,4]triazino[4,3-a]benzimidazole (TBI) were synthesized and tested in vitro and in vivo as a novel class of aldose reductase (ALR2) inhibitors. Compound 3, (10-benzyl[1,2,4]triazino[4,3-a]benzimidazol-3,4(10H)-dion-2-yl)acetic acid, displayed the highest inhibitory activity (IC(50) = 0.36 microM) and was found to be effective in preventing cataract development in severely galactosemic rats when administered as an eyedrop solution. All the compounds investigated were selective for ALR2, since none of them inhibited appreciably aldehyde reductase, sorbitol dehydrogenase, or glutathione reductase. The activity of 3 was lowered by inserting various substituents on the pendant phenyl ring, by shifting the acetic acid moiety from the 2 to the 3 position of the TBI nucleus, or by cleaving the TBI system to yield benzimidazolylidenehydrazines as open-chain analogues. A three-dimensional model of human ALR2 was built, taking into account the conformational changes induced by the binding of inhibitors such as zopolrestat, to simulate the docking of 3 into the enzyme active site. The theoretical binding mode of 3 was fully consistent with the structure-activity relationships in the TBI series and will guide the design of novel ALR2 inhibitors.

Acetates↗

PubMed1993

Molecular defect in factor IX Tokyo: substitution of valine-182 by alanine at position P2' in the second cleavage site by factor XIa resulting in impaired activation.

Utilizing polymerase chain reaction and directly sequencing the amplified exon 6 of the factor IX gene derived from a mild hemophilia Bm patient, we have identified a T to C mutation at nucleotide 20,525. This point mutation predicted a Val182 to Ala substitution in the abnormal factor IX molecule, designated as factor IX Tokyo. The patient manifested a low factor IX activity and a moderately prolonged ox-brain prothrombin time but a normal factor IX antigen level in plasma. Immunopurified factor IX derived from the patient was found to have a normal molecular weight but a reduced specific activity (23% of normal). Limited proteolysis by activated factor XI or by a snake venom-derived factor X-activating enzyme was considerably delayed, indicating the presence of structural alteration(s) most probably at or near the second enzyme-cleavage site. Once activated, however, factor IXa Tokyo was able to activate factor X normally and was inactivated by antithrombin III also in a normal fashion. The structural model of factor IXa and a docking model of factor IX and activated factor VII (factor VIIa) suggested that the Val182 to Ala substitution would not affect the local conformation of the catalytic domain. This mutation would rather loosen the fitness of the molecule into the substrate-binding pocket of factor VIIa due to a shorter side chain of the Ala substitution at the P2' position of the second cleavage site.

Adolescent↗

PubMed2006

Homology modeling and molecular docking study of translationally controlled tumor protein and artemisinin.

Translationally controlled tumor protein (TCTP), also known as histamine releasing factor (HRF), is found abundantly in different eukaryotic cell types. The sequence homology of TCTP between different species is very high, belonging to the MSS4/DSS4 superfamily of proteins. TCTP is involved in both cell growth and human late allergy reaction, as well as having a calcium binding property; however, its primary biological functions remain to be clearly elucidated. In regard to many possible functions, the TCTP of Plasmodium falciparum (Pf) is known to bind with an antimalarial agent, artemisinin, which is activated by heme. It is assumed that the endoperoxide-bridge of artemisinin is opened up by heme to form a free radical, which then eventually alkylates, probably to the Cys14 of PfTCTP. Study of the docking of artemisinin with heme, and subsequently with PfTCTP, was carried out to verify the above hypothesis on the basis of structural interactions. The three dimensional (3D) structure of PfTCTP was built by homology modeling, using the NMR structure of the TCTP of Schizosaccharomyces pombe as a template. The quality of the model was examined based on its secondary structure and biological function, as well as with the use of structure evaluating programs. The interactions between artemisinin, heme and PfTCTP were then studied using the docking program, FlexiDock. The center of the peroxide bond of artemisinin and the Fe of heme were docked within a short distance of 2.6A, implying the strong possibility of an interaction between the two molecules, as proposed. When the activated form of artemisinin was docked on the PfTCTP, the C4-radical of the drug faced towards the sulfur of Cys14 within a distance of 2.48A, again suggesting the possibility of alkylation having occurred. These results confirm the proposed mechanism of the antimalarial effect of artemisinin, which will provide a reliable method for establishing the mechanism of its biological activity using a molecular modeling study.

Antimalarials↗

PubMed2001

Improved scoring of ligand-protein interactions using OWFEG free energy grids.

A new approach to rapidly score protein-ligand interactions is tested on several protein-ligand systems. Results using this approach - the OWFEG free energy grid - are quite promising and are generally in better agreement with experiment (in some cases much better) than those obtained employing scoring techniques currently in wide use. The OWFEG free energy grid is generated from a one-window free energy perturbation MD simulation (Pearlman, D. A. J. Med. Chem. 1999, 42, 4313-4324). The OWFEG approach is applied to three protein systems: IMPDH, MAP kinase p38, and HIV-1 aspartyl protease. OWFEG scores are compared to experimental K(i) and IC50 data in each case. Empirical scoring functions applied to the same systems for comparison include ChemScore, Piecewise Linear Potential (PLP), and Dock energy score.

Databases, Factual↗

PubMed2001

A molecular modeling and 3D QSAR study of a large series of indole inhibitors of human non-pancreatic secretory phospholipase A2.

Automated docking allowing protein-based alignment was performed for a series of 188 indole inhibitors of the human non-pancreatic secretory phospholipase A2 (hnps-PLA2). All the substituted indoles were docked to the crystal structure of hnps-PLA2 and a three-dimensional QSAR model was then established using the CoMFA method. The set of 188 compounds was divided into two subsets, the first one constituting the training set (126 compounds), while the second constituted the test set (62 compounds). The established CoMFA model derived from the training set was then applied to the test set. A good correlation between predicted and experimental activity data allows to validate the 3D QSAR model. A second and global 3D QSAR including all the compounds was established, allowing the creation of the hnps-PLA2 pharmacophore.

Anti-Inflammatory Agents↗

PubMed1991

The three-dimensional structure of a DNA hairpin in solution two-dimensional NMR studies and structural analysis of d(ATCCTATTTATAGGAT).

The hairpin formed by d(ATCCTATTTATAGGAT) was studied by means of two-dimensional NMR spectroscopy and conformational analysis. Almost all 1H resonances of the stem region could be assigned, while the 1H and 31P spectra of the loop region were interpreted completely; this includes the stereospecific assignment of the H5' and H5" resonances. The derivation of the detailed loop structure was carried out in a stepwise fashion including some improved and new methods for structure determination from NMR data. In the first step, the mononucleotide structures were examined. The conformational space available to the mononucleotide was scanned systematically by varying the glycosidic torsion angle and pseudorotational parameters. Each generated conformer was tested against the experimental J coupling constants and NOE parameters. In the following stage, the structures of dinucleotides and longer fragments were derived. Inter-residue distances between protons were calculated by means of a procedure in which the simulated NOEs, obtained via a relaxation-matrix approach, were fitted to the experimental NOEs without the introduction of a molecular model. In addition, the backbone torsion angles beta, gamma and epsilon were deduced from homocoupling and heterocoupling constants. These data served as constraints in the next step, in which the loop sequence was subjected to a multi-conformer generation procedure. The resulting structures were tested against the mentioned constraints and disregarded if these constraints were violated. This yielded a family of structures for the loop region, confined to a relatively narrow conformational space. A representative conformation was subsequently docked on a B-type stem which fulfilled the structural constraints (derived from the NMR experiments for the stem region) to yield the hairpin structure. Results obtained from subsequent restrained-molecular-mechanics as well as free-molecular-mechanics calculations are in accordance with those obtained by means of the analysis described above. The structure of the hairpin loop is a compactly folded conformation and the first base of the central TTTA region forms a Hoogsteen T-A pair with the fourth base. This Hoogsteen base pair is stacked upon the sixth base pair of the B-type double-helical stem. The second base of the loop is folded into the minor groove, whereas the third base of the loop is partly stacked on the first and fourth bases. The phosphate backbone exhibits a sharp turn between the third and fourth nucleotides of the loop. The peculiar structure of this hairpin loop is discussed in relation to loop folding in DNA and RNA hairpins and in relation to a general model for loop folding.

Base Sequence↗

PubMed2005

Structure of the terminal oxygenase component of angular dioxygenase, carbazole 1,9a-dioxygenase.

Carbazole 1,9a-dioxygenase (CARDO) catalyzes the dihydroxylation of carbazole by angular position (C9a) carbon bonding to the imino nitrogen and its adjacent C1 carbon. This reaction is an initial degradation reaction of the carbazole degradation pathway by various bacterial strains. Only a limited number of Rieske non-heme iron oxygenase systems (ROSs) can catalyze this novel reaction, termed angular dioxygenation. Angular dioxygenation is also involved in the degradation pathways of carbazole-related compounds, dioxin, and CARDO can catalyze the angular dioxygenation for dioxin. CARDO consists of a terminal oxygenase component (CARDO-O), and the electron transport components, ferredoxin (CARDO-F) and ferredoxin reductase (CARDO-R). CARDO-O has a homotrimeric structure, and governs the substrate specificity of CARDO. Here, we have determined the crystal structure of CARDO-O of Janthinobacterium sp. strain J3 at a resolution of 1.95A. The alpha3 trimeric overall structure of the CARDO-O molecule roughly corresponds to the alpha3 partial structures of other terminal oxygenase components of ROSs that have the alpha3beta3 configuration. The CARDO-O structure is a first example of the terminal oxygenase components of ROSs that have the alpha3 configuration, and revealed the presence of the specific loops that interact with a neighboring subunit, which is proposed to be indispensable for stable alpha3 interactions without structural beta subunits. The shape of the substrate-binding pocket of CARDO-O is markedly different from those of other oxygenase components involved in naphthalene and biphenyl degradation pathways. Docking simulations suggested that carbazole binds to the substrate-binding pocket in a manner suitable for catalysis of angular dioxygenation.

Bacterial Proteins↗

PubMed2005

Electrical docking of microtubules for kinesin-driven motility in nanostructures.

We demonstrate localized electrical control of the docking of microtubules onto engineered kinesin-coated structures. After applying a voltage to a gold electrode, we observe an enhanced transport of microtubules from solution toward the surface and a subsequent increase of the amount of moving microtubule shuttles. Switching off the voltage leads to a partial detachment of microtubules from the surface. The surface coverage of microtubules, during both the docking and undocking events, follows an exponential time dependence. We provide a simple kinetic model, incorporating the equilibrium between free and surface-bound microtubules, that explains these data.

Binding Sites↗

PubMed2006

Prediction of binding affinities between the human amphiphysin-1 SH3 domain and its peptide ligands using homology modeling, molecular dynamics and molecular field analysis.

The SH3 domain of the human protein amphiphysin-1, which plays important roles in clathrin-mediated endocytosis, actin function and signaling transduction, can recognize peptide motif PXRPXR (X is any amino acid) with high affinity and specificity. We have constructed a complex structure of the amphiphysin-1 SH3 domain and a high-affinity peptide ligand PLPRRPPRA using homology modeling and molecular docking, which was optimized by molecular dynamics (MD). Three-dimensional quantitative structure-affinity relationship (3D-QSAR) analyses on the 200 peptides with known binding affinities to the amphiphysin-1 SH3 domain was then performed using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). The best CoMSIA model showed promising predictive power, giving good predictions for about 95% of the peptides in the test set (absolute prediction errors less than 1.0). It was used to validate peptide-SH3 binding structure and provide insight into the structural requirements for binding of peptides to SH3 domains. Finally, MD simulations were performed to analyze the interaction between the SH3 domain and another peptide GFPRRPPPRG that contains with the PXRPXsR (s represents residues with small side chains) motif. MD simulations demonstrated that the binding conformation of GFPRRPPPRG is quite different from that of PLPRRPPRAA especially the four residues at the C terminal, which may explain why the CoMSIA model cannot give good predictions on the peptides of the PXRPXsR motif. Because of its efficiency and predictive power, the 3D-QSAR model can be used as a scoring filter for predicting peptide sequences bound to SH3 domains.

Amino Acid Sequence↗

PubMed2004

Single-molecule enzymology of RNA: essential functional groups impact catalysis from a distance.

The hairpin ribozyme is a minimalist paradigm for studying RNA folding and function. In this enzyme, two domains dock by induced fit to form a catalytic core that mediates a specific backbone cleavage reaction. Here, we have fully dissected its reversible reaction pathway, which comprises two structural transitions (docking/undocking) and a chemistry step (cleavage/ligation), by applying a combination of single-molecule fluorescence resonance energy transfer (FRET) assays, ensemble cleavage assays, and kinetic simulations. This has allowed us to quantify the effects that modifications of essential functional groups remote from the site of catalysis have on the individual rate constants. We find that all ribozyme variants show similar fractionations into effectively noninterchanging molecule subpopulations of distinct undocking rate constants. This leads to heterogeneous cleavage activity as commonly observed for RNA enzymes. A modification at the domain junction additionally leads to heterogeneous docking. Surprisingly, most modifications not only affect docking/undocking but also significantly impact the internal chemistry rate constants over a substantial distance from the site of catalysis. We propose that a network of coupled molecular motions connects distant parts of the RNA with its reaction site, which suggests a previously undescribed analogy between RNA and protein enzymes. Our findings also have broad implications for applications such as the action of drugs and ligands distal to the active site or the engineering of allostery into RNA.

Catalysis↗

PubMed1997

Molecular modeling of the interaction of diagnostic radiopharmaceuticals with receptor proteins: m2 antagonist binding to the muscarinic m2 subtype receptor.

Models of the m2 muscarinic receptor have been built and acetylcholine and an antagonist of the quinuclidinyl benzilate family docked to the putative active site. We have incorporated aspects of homology, site-directed mutagenesis studies and structure-activity studies of specific lead compounds in the construction of our receptor models with a primary focus on the structure of the binding sites. We have observed a deep pocket binding of 5-BrQNT, suggesting a plausible explanation for the observation that agonists and antagonists do not bind competitively. The results of these computational studies are interpreted within the context of the observed in vitro results. Our goal is to assist in the development of subtype receptor selective radiopharmaceuticals for use in PET and SPECT.

Acetylcholine↗

PubMed2004

Binding studies and computer-aided modelling of macromolecule/odorant interactions.

Odorant-to-biopolymer (proteins and polysaccharides) binding properties influence the partition coefficients of odorant/matrix mixtures and depend on the molecular structure of the guest and the host. While air/solvent partition coefficients of odorants influence the flavor intensity of compounds in the headspace above the solvent, more complex mechanisms are proposed for the binding of a molecule to a food matrix. Odorant air/solvent partition coefficients are dependent on both the physico-chemical properties of odorant and solvent. Binding affinities for flavor compounds on various biopolymers can be estimated by calculation of physico-chemical descriptors for the odorants. Binding affinities of gamma- and delta-lactones (C7-C11) to bovine serum albumin (BSA) and beta-lactoglobulin (BLG) were investigated by ultracentrifugation and equilibrium-dialysis techniques. Quantitative structure-activity relationships (QSAR) of lactone binding on proteins (BLG and BSA) were performed by the measurement of lipophilicity and H-bond strength. Large differences in observed protein-binding properties for the various compounds clearly demonstrated that structure-activity relationship was significantly influenced by the lipophilicity of the odorant. If the structure of the receptor molecule is known, computational ligand-macromolecule docking experiments can be used to predict binding affinities for unknown compounds with the receptor molecule. A BLG-lactone binding position, not previously reported in the literature, has been identified and confirmed by competitive binding studies. A model has been developed to estimate the free energy of binding of odorants to biopolymers. Estimated free energies of binding of lactones with BLG from computational methods were in very good agreement with the experimentally obtained results.

Animals↗

PubMed2002

A fast empirical approach to binding free energy calculations based on protein interface information.

Three useful variables from the interfaces of 20 protein-protein complexes were investigated. These variables are the side-chain accessible number (N(b)), the number of hydrophilic pairs (N(pair)) and buried a polar solvent accessible surface areas (DeltaDeltaASA(apol)). An empirical model based on the three variables was developed to describe the free energy of protein associations. As the results show, the side-chain accessible numbers characterize the loss of side-chain conformational entropy of protein interactions and the effective empirical function presented here has great capability for estimating the binding free energy. It was found that the variables of interface information capture most of the significant features of protein-protein association. Also, we applied the model based on the variables as a rescoring function to docking simulations and found that it has the potential to distinguish the 'true' binding mode. It is clear that the simple and empirical scale developed here is an attractive target function for calculating binding free energy for various biological processes to rational protein design.

Algorithms↗

PubMed2002

Modeling of binding modes and inhibition mechanism of some natural ligands of farnesyl transferase using molecular docking.

Several natural inhibitors of farnesyl transferase have been reported in the literature: some compounds are competitive with farnesyl pyrophosphate (FPP), whereas other ones are competitive with Ras proteins, even though it is usually hard to highlight their inhibition mechanism, which is still unknown for several natural compounds. The aim of this work is to show that the molecular docking analysis can be successfully used to underline the inhibition mechanism of these natural compounds. First, the selected compounds were subjected to a detailed docking analysis, by means of BioDock, a program able to reveal the most likely binding mode for each ligand. By comparing these results with the binding sites for the natural substrates, earlier determined, it was possible to highlight the site specificity and the inhibition mechanism of the selected compounds. In addition, it is possible to relate the binding mode of these molecules with their lipole values, which is appreciably less for peptidomimetics than for FPP mimetic and reveals a straightforward method to predict and to understand the inhibition mechanism of these natural derivatives.

Alkyl and Aryl Transferases↗

PubMedDate not supplied

Creating artificial binding pocket boundaries to improve the efficiency of flexible ligand docking.

Traditionally, algorithms for binding site characterization or identification focus on the problem of identifying atoms within a macromolecule that might be responsible for ligand binding. In this manuscript, we focus on the binding pocket problem from a different perspective as a challenge of calculating an artificial binding pocket boundary that is sufficient to isolate binding pocket volume. The approach involves the calculation of a macromolecule encapsulating surface (MES) that separates binding pocket volume from outside space. We show that the MES can be used to increase the efficiency of flexible docking as implemented in AutoDock 3.0. The most significant improvement in docking efficiency is seen when the entire protein is searched and results show additional support for the use of AutoDock, in and of itself, as a feasible tool for binding-site identification for cases in which a protein ligand is known.

Algorithms↗

PubMed2004

On the rational design of substrate mimetics: The function of docking approaches for the prediction of protease specificities.

The behaviour of substrate mimetics in mediating the acceptance of nonspecific acyl moieties by proteases has been investigated as a direct function of their site-specific ester leaving groups. In this contribution we report on a computational approach to rationalise this interplay and to predict the power of a potential ester moiety to act as a suitable substrate mimetic for a given enzyme by means of an automated docking procedure. Investigations with seven distinct substrate mimetics and two proteases, subtilisin and chymotrypsin, show a clear correlation between the theoretically calculated binding energies DeltaE and the specificity constants k(cat)KM(-1) obtained from parallel hydrolysis kinetic studies. These results prove the general function of the docking approach as a rational model not only in predicting the general acceptance of a substrate mimetic in a qualitative manner, but also to provide reliable information on its individual specificity towards proteases.

Alanine↗

PubMed1999

A structure-function study of the C2 domain of cytosolic phospholipase A2. Identification of essential calcium ligands and hydrophobic membrane binding residues.

The C2 domain of cytosolic phospholipase A2 (cPLA2) is involved in the Ca2+-dependent membrane binding of this protein. To identify protein residues in the C2 domain of cPLA2 essential for its Ca2+ and membrane binding, we selectively mutated Ca2+ ligands and putative membrane-binding residues of cPLA2 and measured the effects of mutations on its enzyme activity, membrane binding affinity, and monolayer penetration. The mutations of five Ca2+ ligands (D40N, D43N, N65A, D93N, N95A) show differential effects on the membrane binding and activation of cPLA2, indicating that two calcium ions bound to the C2 domain have differential roles. The mutations of hydrophobic residues (F35A, M38A, L39A, Y96A, Y97A, M98A) in the calcium binding loops show that the membrane binding of cPLA2 is largely driven by hydrophobic interactions resulting from the penetration of these residues into the hydrophobic core of the membrane. Leu39 and Val97 are fully inserted into the membrane, whereas Phe35 and Tyr96 are partially inserted. Finally, the mutations of four cationic residues in a beta-strand (R57E/K58E/R59E/R61E) have modest and negligible effects on the binding of cPLA2 to zwitterionic and anionic membranes, respectively, indicating that they are not directly involved in membrane binding. In conjunction with our previous study on the C2 domain of protein kinase C-alpha (Medkova, M., and Cho, W. (1998) J. Biol. Chem. 273, 17544-17552), these results demonstrate that C2 domains are not only a membrane docking unit but also a module that triggers membrane penetration of protein and that individual Ca2+ ions bound to the calcium binding loops play differential roles in the membrane binding and activation of their parent proteins.

Amino Acid Sequence↗