PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “molecular docking”

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 397 records · Page 22Linked to original sources

trans-Parinaric acid as a versatile spectroscopic label to study ligand binding properties of bovine beta-lactoglobulin.

Advantageous spectroscopic properties of the plant derived polyunsaturated trans-parinaric acid (tPnA) was demonstrated in obtaining valuable data on the ligand binding characteristics of the lipocalin member bovine beta-lactoglobulin A (BLG-A). Titration of the protein with tPnA resulted in the appearance of an intense negative induced circular dichroism (CD) band and bathochromic shift of the ultraviolet (UV) peak of the ligand. The extrinsic optical activity was interpreted by the chiral contribution of the allylic axial CH bonds of tPnA to the pi-pi(*) transition of the planar tetraene chromophore. Analysis of the series of induced CD curves obtained by CD titration experiment indicated the complexation of a single ligand molecule to a uniform protein binding site. Additionally, the dramatic increase of fluorescence intensity of the lactoglobulin bound ligand suggested the hydrophobic nature of the binding site. CD and fluorescence titration data were utilized to calculate the binding constant (K(a)) of which high value ( approximately 10(6)M(-1)) refers to strong protein association of tPnA. pH dependent reversible dis- and reappearance of the induced CD signal unambigously proved the inclusion of tPnA into the central hydrophobic cavity of the lactoglobulin governed by the protonation induced conformational movement of the EF loop at the opening of the calyx. This conclusion was supported and complemented by molecular docking calculations.

Animals↗

Docking of hydrophobic ligands with interaction-based matching algorithms.

MOTIVATION: Matching of chemical interacting groups is a common concept for docking and fragment placement algorithms in computer-aided drug design. These algorithms have been proven to be reliable and fast if at least a certain number of hydrogen bonds or salt bridges occur. However, the algorithms typically run into problems if hydrophobic fragments or ligands should be placed. In order to dock hydrophobic fragments without significant loss of computational efficiency, we have extended the interaction model and placement algorithms in our docking tool FlexX. The concept of multi-level interactions is introduced into the algorithms for automatic selection and placement of base fragments. RESULTS: With the multi-level interaction model and the corresponding algorithmic extensions, we were able to improve the overall performance of FlexX significantly. We tested the approach with a set of 200 protein-ligand complexes taken from the Brookhaven Protein Data Bank (PDB). The number of test cases which can be docked within 1.5 A RMSD from the crystal structure can be increased from 58 to 64%. The performance gain is paid for by an increase in computation time from 73 to 91 s on average per protein-ligand complex. AVAILABILITY: The FlexX molecular docking software is available for UNIX platforms IRIX, Solaris and Linux. See http://cartan.gmd.de/FlexX for additional information.

Algorithms↗

A computational comparison of the atomic models of the actomyosin interface.

Several atomic models of the actomyosin interface have been proposed based on the docking together of their component structures using electron microscopy and resonance energy-transfer measurements. Although these models are in approximate agreement in the location of the binding interfaces when myosin is tightly bound to actin, their relationships to molecular docking simulations based on computational free-energy calculations are investigated here. Both rigid-docking and flexible-docking conformational search strategies were used to identify free-energy minima at the interfaces between atomic models of myosin and actin. These results suggest that the docking model produced by resonance energy-transfer data is closer to a free-energy minimum at the interface than are the available atomic models based on electron microscopy. The conformational searches were performed using both scallop and chicken skeletal muscle myosins and identified similarly oriented actin-binding interfaces that serve to validate that these models are at the global minimum. These results indicate that the existing docking models are close to but not precisely at the lowest-energy initial contact site for strong binding between myosin and actin that should represent an initial contact between the two proteins; therefore, conformational changes are likely to be important during the transition to a strongly bound complex.

Actins↗

Combining NMR spectral and structural data to form models of polychlorinated dibenzodioxins, dibenzofurans, and biphenyls binding to the AhR.

A three-dimensional quantitative spectrometric data-activity relationship (3D-QSDAR) modeling technique which uses NMR spectral and structural information that is combined in a 3D-connectivity matrix has been developed. A 3D-connectivity matrix was built by displaying all possible assigned carbon NMR chemical shifts, carbon-to-carbon connections, and distances between the carbons. Two-dimensional 13C-13C COSY and 2D slices from the distance dimension of the 3D-connectivity matrix were used to produce a relationship among the 2D spectral patterns for polychlorinated dibenzofurans, dibenzodioxins, and biphenyls (PCDFs, PCDDs, and PCBs respectively) binding to the aryl hydrocarbon receptor (AhR). We refer to this technique as comparative structural connectivity spectral analysis (CoSCoSA) modeling. All CoSCoSA models were developed using forward multiple linear regression analysis of the predicted 13C NMR structure-connectivity spectral bins. A CoSCoSA model for 26 PCDFs had an explained variance (r2) of 0.93 and an average leave-four-out cross-validated variance (q(2)4) of 0.89. A CoSCoSA model for 14 PCDDs produced an r2 of 0.90 and an average leave-two-out cross-validated variance (q(2)2) of 0.79. One CoSCoSA model for 12 PCBs gave an r2 of 0.91 and an average q(2)2 of 0.80. Another CoSCoSA model for all 52 compounds had an r2 of 0.85 and an average q(2)2 of 0.52. Major benefits of CoSCoSA modeling include ease of development since the technique does not use molecular docking routines.

Algorithms↗

Identification of small molecule inhibitors of anthrax lethal factor.

The virulent spore-forming bacterium Bacillus anthracis secretes anthrax toxin composed of protective antigen (PA), lethal factor (LF) and edema factor (EF). LF is a Zn-dependent metalloprotease that inactivates key signaling molecules, such as mitogen-activated protein kinase kinases (MAPKK), to ultimately cause cell death. We report here the identification of small molecule (nonpeptidic) inhibitors of LF. Using a two-stage screening assay, we determined the LF inhibitory properties of 19 compounds. Here, we describe six inhibitors on the basis of a pharmacophoric relationship determined using X-ray crystallographic data, molecular docking studies and three-dimensional (3D) database mining from the US National Cancer Institute (NCI) chemical repository. Three of these compounds have K(i) values in the 0.5-5 microM range and show competitive inhibition. These molecular scaffolds may be used to develop therapeutically viable inhibitors of LF.

Animals↗

Mapping of possible binding sequences of two beta-sheet breaker peptides on beta amyloid peptide of Alzheimer's disease.

Aggregation of amyloid peptide (Abeta) has been identified as a major feature of the pathogenesis of Alzheimer's disease. Increased risk for disease is associated with increased formation of polymerized Abeta. Inhibition of formation of toxic (aggregated) form of Abeta is one of the therapeutic possibilities. Beta sheet breaker peptides (BSBs) fulfill the requirements of an effective inhibitor. After having attached to the Abeta molecules, BSBs can prevent aggregation of Abeta to polymeric forms (aggregates). In the present study, we performed molecular modelling of complex formation between Abeta and two BSB peptides. Our aim was to find proper binding sequences for the BSB peptides on Abeta and characterize them. A dimeric model of Abeta was also used to study the interaction of BSBs with the aggregated forms of Abeta and find the sequences responsible for the polymerization process. A fast and efficient computational method: molecular docking was used for the afore-mentioned purposes.

Amino Acid Sequence↗

Binding conformers searching method for ligands according to the structures of their receptors and its application to thrombin inhibitors.

AIM: To develop a method of finding binding conformers for ligands according to the three-dimensional structures of their receptors. METHODS: Combining the systematic search method of ligand with the molecular docking approach of ligand fitting into its receptor, we developed a binding conformer searching method for ligands. RESULTS: The binding conformers of phosphonopeptidyl thrombin inhibitors were recognized. The binding (interaction) energies between these inhibitors and thrombin were calculated with molecular mechanical method. CONCLUSION: Both of the total binding energies and steric binding energies have good correlations with the inhibitory activities of these thrombin inhibitors, demonstrating that our approach is reasonable. It can also be used to explain the inhibition mechanism of thrombin interacting with these inhibitors.

Binding Sites↗

Studying the roles of W86, E202, and Y337 in binding of acetylcholine to acetylcholinesterase using a combined molecular dynamics and multiple docking approach.

A combined molecular dynamics simulation and multiple ligand docking approach is applied to study the roles of the anionic subsite residues (W86, E202, Y337) in the binding of acetylcholine (ACh) to acetylcholinesterase (AChE). We find that E202 stabilizes docking of ACh via electrostatic interactions. However, we find no significant electrostatic contribution from the aromatic residues. Docking energies of ACh to mutant AChE show a more pronounced effect because of size/shape complementarity. Mutating to smaller residues results in poorer binding, both in terms of docking energy and statistical docking probability. Besides separating out electrostatics by turning off the partial charges from each residue and comparing it with the native, the mutations in this study are W86F, W86A, E202D, E202Q, E202A, Y337F, and Y337A. We also find that all perturbations result in a significant reduction in binding of extended ACh in the catalytically productive orientation. This effect is primarily caused by a small shift in preferred position of the quaternary tail.

Acetylcholine↗

Docking of 4-oxalocrotonate tautomerase substrates: implications for the catalytic mechanism.

The enzyme 4-oxalocrotonate tautomerase catalyzes the ketonization of dienols, which after further processing become intermediates in the Krebs cycle. The enzyme uses a general acid-base mechanism for proton transfer: the amino-terminal proline has been shown to function as the catalytic base and Arg39 has been implicated as the catalytic acid. We report the results of molecular docking simulations of 4-oxalocrotonate tautomerase with two substrates, 2-hydroxymuconate and 5-carboxymethyl-2-hydroxymuconate. pKa calculations are also performed for the free enzyme. The predicted binding mode of 2-hydroxymuconate is in agreement with experimental data. A model for the binding mode of 5-carboxymethyl-2-hydroxymuconate is proposed which explains the lower catalytic efficiency of the enzyme toward this substrate. The pKa predictions and docking simulations support residue Arg39 as the general acid for the enzyme catalysis.

Aldehydes↗

Distributed automated docking of flexible ligands to proteins: parallel applications of AutoDock 2.4.

AutoDock 2.4 predicts the bound conformations of a small, flexible ligand to a nonflexible macromolecular target of known structure. The technique combines simulated annealing for conformation searching with a rapid grid-based method of energy evaluation based on the AMBER force field. AutoDock has been optimized in performance without sacrificing accuracy; it incorporates many enhancements and additions, including an intuitive interface. We have developed a set of tools for launching and analyzing many independent docking jobs in parallel on a heterogeneous network of UNIX-based workstations. This paper describes the current release, and the results of a suite of diverse test systems. We also present the results of a systematic investigation into the effects of varying simulated-annealing parameters on molecular docking. We show that even for ligands with a large number of degrees of freedom, root-mean-square deviations of less than 1 A from the crystallographic conformation are obtained for the lowest-energy dockings, although fewer dockings find the crystallographic conformation when there are more degrees of freedom.

Azepines↗

Evaluating docked complexes with the HINT exponential function and empirical atomic hydrophobicities.

Methods that predict geometries of ligands binding to receptor molecules can facilitate ligand discovery and yield information on the factors governing complementarity. Here, the use of atomic hydrophobicities in evaluating binding modes has been examined with four ligand-receptor complexes of known structure. In each system, hundreds of hypothetical binding orientations were generated with DOCK and evaluated using the HINT (Hydropathic INTeractions) exponential function and atomic hydrophobic constants. In three of the four systems, the experimental binding mode received the best HINT score; in the fourth system, the experimental binding mode scored only slightly lower than a similar, apparently reasonable orientation. The HINT function may be generally useful as a scoring method in molecular docking.

Binding Sites↗

Flavonoid-3',5'-hydroxylase from Phalaenopsis: a novel member of cytochrome P450s, its cDNA cloning, endogenous expression and molecular modeling.

A novel cDNA for the flavonoid-3',5'-hydroxylase (F35H) gene was cloned from petals of Phalaenopsis, and designated as Phf35h (accession number DQ148458 in GenBank/EMBL/DDBJ). The genomic clone of Phf35h was isolated by a PCR-based strategy. Nucleotide sequence analysis revealed that its genomic clone contains one intron and an open reading frame encoding a polypeptide of 507 amino acid residues. Southern hybridization analysis indicated the presence of a single gene coding for Phf35h. RT-PCR analysis showed that the Phf35h mRNA is transcribed in late phase of petal development, which is concomitant with the appearance of anthocyanins in petal tissue. The transcript is abundant in the purple petals but not in leaves or roots. The three-dimensional model of PhF3'5'H protein is classified into an alpha-domain which contains most of the alpha-helices with three small beta-sheets, a beta-domain that contains the larger beta-sheets with three small alpha-helices by homology modeling. The substrate-binding site for dihydrokaempferol on PhF3'5'H protein was determined by molecular docking algorithm. A highly conserved HPPTPLSLPH sequence was predicted to contact the aromatic ring of dihydrokaempferol.

Amino Acid Sequence↗

Studying enzyme binding specificity in acetylcholinesterase using a combined molecular dynamics and multiple docking approach.

A combined molecular dynamics simulation and multiple ligand docking approach is applied to study the binding specificity of acetylcholinesterase (AChE) with its natural substrate acetylcholine (ACh), a family of substrate analogues, and choline. Calculated docking energies are well correlated to experimental k(cat)/K(M) values, as well as to experimental binding affinities of a related series of TMTFA inhibitors. The "esteratic" and "anionic" subsites are found to act together to achieve substrate binding specificity. We find that the presence of ACh in the active site of AChE not only stabilizes the setup of the catalytic triad but also tightens both subsites to achieve better binding. The docking energy gained from this induced fit is 0.7 kcal/mol for ACh. For the binding of the substrate tailgroup to the anionic subsite, both the size and the positive charge of the tailgroup are important. The removal of the positive charge leads to a weaker binding of 1 kcal/mol loss in docking energy. Substituting each tail methyl group with hydrogen results in both an incremental loss in docking energy and also a decrease in the percentage of structures docked in the active site correctly set up for catalysis.

Acetophenones↗

Structure-based discovery of inhibitors of thymidylate synthase.

A molecular docking computer program (DOCK) was used to screen the Fine Chemical Directory, a database of commercially available compounds, for molecules that are complementary to thymidylate synthase (TS), a chemotherapeutic target. Besides retrieving the substrate and several known inhibitors, DOCK proposed putative inhibitors previously unknown to bind to the enzyme. Three of these compounds inhibited Lactobacillus casei TS at submillimolar concentrations. One of these inhibitors, sulisobenzone, crystallized with TS in two configurations that differed from the DOCK-favored geometry: a counterion was bound in the substrate site, which resulted in a 6 to 9 angstrom displacement of the inhibitor. The structure of the complexes suggested another binding region in the active site that could be exploited. This region was probed with molecules sterically similar to sulisobenzone, which led to the identification of a family of phenolphthalein analogs that inhibit TS in the 1 to 30 micromolar range. These inhibitors do not resemble the substrates of the enzyme. A crystal structure of phenolphthalein with TS shows that it binds in the target site in a configuration that resembles the one suggested by DOCK.

Amino Acid Sequence↗

Parallel solid-phase synthesis and evaluation of inhibitors of Streptomyces coelicolor type II dehydroquinase.

A series of 1-substituted and 4-substituted benzyl analogues of the known inhibitor (1S,3R,4R)-1,3,4-trihydroxy-5-cyclohexene-1-carboxylic acid has been synthesized and tested as inhibitors of Streptomyces coelicolor type II dehydroquinase. The solid-phase syntheses of 18 new analogues are reported. The most potent inhibitor, 2-nitrobenzyloxy analogue 5i, has K(i) of 8 microM, more than 30 times lower than the K(M) of the substrate and approximately 4 times more potent than the original inhibitor. The binding modes of the synthesized analogues in the active site were studied by molecular docking with GOLD 2.0.

Benzyl Compounds↗

Molecular modeling as a powerful technique for understanding small-large molecules interactions.

In the present review we summarize recent work, aimed at a better understanding of the interactions in macromolecule ligand complexes, performed by means of computational tools such as pseudoreceptor generation, molecular docking, conformational search and energy minimization. While the first approach has been applied when the three-dimensional structural properties of the biological target were unknown, the remaining protocols exploited the knowledge of the overall structure of the involved macromolecules and their active sites. Molecular modeling techniques were used in the cases reported to study and propose macromolecular binding sites and to predict their interactions with bioactive conformers of the ligands.

Anti-HIV Agents↗

A novel computational analysis of ligand-induced conformational changes in the ATP binding sites of cyclin dependent kinases.

Protein kinases in general are known to be very flexible macromolecules. In this article, the conformational plasticity of the ATP binding site in cyclin dependent kinases is analyzed. Movement of the two lysine residues lining the ATP binding site are shown to play a major role in the conformational variability of the site. Linear models are developed to identify and quantify ligand properties that maximally influence the lysine side chain conformations. A few simple properties of the ligands are shown to account for more than 70% of the variation in the lysine conformations. The results are validated using test data and molecular simulation studies. Illustrative applications of the results of this analysis to finding the appropriate crystal structure for molecular docking and binding mode predictions of novel ligands are provided. This work provides a new approach to quantify ligand-induced conformational changes in the active sites of flexible proteins and to find the appropriate crystal structure for docking novel ligands.

Adenosine Triphosphate↗