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 505 records · Page 28Linked to original sources

PubMed1996

Scoring noncovalent protein-ligand interactions: a continuous differentiable function tuned to compute binding affinities.

Exploitation of protein structures for potential drug leads by molecular docking is critically dependent on methods for scoring putative protein-ligand interactions. An ideal function for scoring must exhibit predictive accuracy and high computational speed, and must be tolerant of variations in the relative protein-ligand molecular alignment and conformation. This paper describes the development of an empirically derived scoring function, based on the binding affinities of protein-ligand complexes coupled with their crystallographically determined structures. The function's primary terms involve hydrophobic and polar complementarity, with additional terms for entropic and solvation effects. The issue of alignment/conformation dependence was solved by constructing a continuous differentiable nonlinear function with the requirement that maxima in ligand conformation/alignment space corresponded closely to crystallographically determined structures. The expected error in the predicted affinity based on cross-validation was 1.0 log unit. The function is sufficiently fast and accurate to serve as the objective function of a molecular-docking search engine. The function is particularly well suited to the docking problem, since it has spatially narrow maxima that are broadly accessible via gradient descent.

Bacterial Proteins↗

PubMed2003

Molecular modeling of the neurohypophyseal receptor/atosiban complexes.

The neurohypophyseal nonapeptide hormone oxytocin (OT) is the strongest uterotonic substance known and is responsible for the initiation of labor. Conversely, oxytocin antagonists blocking uterine OT receptor can suppress uterus contraction. In this paper we describe a computer simulated docking pertinent to affinity of an oxytocin antagonist atosiban towards OT receptor, versus vasopressin V1a and V2 receptors.

Amino Acids↗

PubMed2006

Docking protein domains in contact space.

BACKGROUND: Many biological processes involve the physical interaction between protein domains. Understanding these functional associations requires knowledge of the molecular structure. Experimental investigations though present considerable difficulties and there is therefore a need for accurate and reliable computational methods. In this paper we present a novel method that seeks to dock protein domains using a contact map representation. Rather than providing a full three dimensional model of the complex, the method predicts contacting residues across the interface. We use a scoring function that combines structural, physicochemical and evolutionary information, where each potential residue contact is assigned a value according to the scoring function and the hypothesis is that the real configuration of contacts is the one that maximizes the score. The search is performed with a simulated annealing algorithm directly in contact space. RESULTS: We have tested the method on interacting domain pairs that are part of the same protein (intra-molecular domains). We show that it correctly predicts some contacts and that predicted residues tend to be significantly closer to each other than other pairs of residues in the same domains. Moreover we find that predicted contacts can often discriminate the best model (or the native structure, if present) among a set of optimal solutions generated by a standard docking procedure. CONCLUSION: Contact docking appears feasible and able to complement other computational methods for the prediction of protein-protein interactions. With respect to more standard docking algorithms it might be more suitable to handle protein conformational changes and to predict complexes starting from protein models.

Algorithms↗

PubMed2006

Insights into the o-acetylation reaction of hydroxylated heterocyclic amines by human arylamine N-acetyltransferases: a computational study.

A computational study was performed to better understand the differences between human arylamine N-acetyltransferase (NAT) 1 and 2. Homology models were constructed from available crystal structures, and comparisons of the active site residues 125, 127, and 129 for these two enzymes provide insight into observed substrate differences. The NAT2 model provided a basis for understanding how some of the common polymorphisms may affect the structure of this protein. Molecular dynamics simulations of the human NAT models and the template structure (NAT from Mycobacterium smegmatis) were performed and showed the models to be stable and reasonable. Docking studies of hydroxylated heterocyclic amines in the models of NAT1 and NAT2 probed the differences exhibited by these two proteins with mutagenic agents. The hydroxylated heterocyclic amines were only able to fit into the NAT2 active site, and an alternative binding site by the phosphate-binding loop was found using our models and will be discussed. Quantum mechanical calculations on the O-acetylation reaction of the hydroxylated heterocyclic amines N-OH MeIQx and N-OH PhIP show that the reaction coordinates differ for these two compounds, but the activation barrier separating the reactant from the product are both low. The results of this study suggest that common polymorphisms in human NAT2 are distant from the active site and are more likely to destabilize the enzyme than affect catalysis. Additionally, the quantum mechanical calculations show that the observed differences in mutagenic activity between N-OH MeIQx and N-OH PhIP are not related to their acetylation reaction with NAT.

Acetylation↗

PubMed2003

Modeling and docking of the three-dimensional structure of the human melanocortin 4 receptor.

A three-dimensional structure of the human melanocortin 4 receptor (hMC4R) is constructed in this study using a computer-aided molecular modeling approach. Human melanocortin 4 receptor is a G Protein-Coupled Receptor (GPCR). We structurally aligned transmembrane helices with bovine rhodopsin transmembrane domains, simulated both intracellular and extracellular loop domains on homologous loop regions in other proteins of known 3D structure and modeled the C terminus on the corresponding part of bovine rhodopsin. Then tandem minimization and dynamics calculations were run to refine the crude structure. The simulative model was tested by docking with a triplet peptide (RFF) ligand. It was found that the ligand is located among transmembrane regions TM3, TM4, TM5, and TM6 of hMC4R. In consistence with mutational and biochemical data, binding site is mainly formed as a hydrophobic and negatively charged pocket. The model constructed here might provide a structural framework for making rational predictions in relevant fields.

Amino Acid Sequence↗

PubMed2005

Synthesis, biological evaluation, and modeling studies of inhibitors aimed at the malarial proteases plasmepsins I and II.

The increasing resistance of the malarial parasite to antimalarial drugs is a major contributor to the reemergence of the disease and increases the need for new drug targets. The two aspartic proteases, plasmepsins I and II, from Plasmodium falciparum have recently emerged as potential targets. In an effort to inhibit these hemoglobinases, a series of inhibitors encompassing a basic hydroxyethylamine transition state isostere as a central fragment were prepared. The synthesized compounds were varied in the P1' position and exhibited biological activities in the range of 31 to >2000 nM. To try to rationalize the results, molecular docking and 3D-QSAR analysis were used.

Animals↗

PubMed2007

Interaction of malaria parasite-inhibitory antibodies with the merozoite surface protein MSP1(19) by computational docking.

Merozoite surface protein 1 (MSP1) of the malaria parasite Plasmodium falciparum is an important vaccine candidate antigen. Antibodies specific for the C-terminal maturation product, MSP1(19), have been shown to inhibit erythrocyte invasion and parasite growth. Specific monoclonal antibodies react with conformational epitopes contained within the two EGF-like domains that constitute the antigen MSP1(19). To gain greater insight into the inhibitory process, the authors selected two strongly inhibitory antibodies (designated 12.8 and 12.10) and modeled their structures by homology. Computational docking was used to generate antigen-antibody complexes and a selection filter based on NMR data was applied to obtain plausible models. Molecular Dynamics simulations of the selected complexes were performed to evaluate the role of specific side chains in the binding. Favorable complexes were obtained that complement the NMR data in defining specific binding sites. These models can provide valuable guidelines for future experimental work that is devoted to the understanding of the action mechanism of invasion-inhibitory antibodies.

Amino Acid Sequence↗

PubMed2006

A new Bacillus pasteurii urease inhibitor from Euphorbia decipiens.

Inhibition of Bacillus pasteurii urease enzyme by 3,7,15-tri-O-acetyl-5-O-nicotinoyl-13,14-dihydroxymyrsinol (1), a diterpene ester with a myrsinol-type skeleton, isolated from Euphorbia decipiens Boiss. and Buhse, was un-competitive consistent with the molecular docking results. The Ki value was 117.40 +/- 0.7 microM.

Bacillus↗

PubMedDate not supplied

Molecular docking towards drug discovery.

Fueled by advances in molecular structure determination, tools for structure-based drug design are proliferating rapidly. Lead discovery through searching of ligand databases with molecular docking techniques represents an attractive alternative to high-throughout random screening. The size of commercial databases imposes severe computational constraints on molecular docking, compromising the level of calculational detail permitted for each putative ligand. We describe alternative philosophies for docking which effectively address this challenge. With respect to the dynamic aspects of molecular recognition, these strategies lie along a spectrum of models bounded by the Lock-and-Key and Induced-Fit theories for ligand binding. We explore the potential of a rigid model in exploiting species specificity and of a tolerant model in predicting absolute ligand binding affinity. Current molecular docking methods are limited primarily by their ability to rank docked complexes; we therefore place particular emphasis on this aspect of the problem throughout our validation of docking strategies.

Animals↗

PubMed2006

Synthesis and evaluation of N-acetyl-L-tyrosine based compounds as PPARalpha selective activators.

The development of type 2 diabetes in obese individuals is linked to lipid accumulation in non-adipose tissues. A series of N-acetyl-L-tyrosine derivatives were synthesized and evaluated for PPAR transactivation. Compounds 4d and 4f were found to show better PPARalpha transactivation as compared to PPARgamma. Molecular docking analysis was carried out to study their important interactions with the active site of PPARalpha.

Binding Sites↗

PubMed2006

Molten globule formation in apomyoglobin monitored by the fluorescent probe Nile Red.

The interaction of nile red (NR) with apomyoglobin (ApoMb) in the native (pH 7) and molten globule (pH 4) states was investigated using experimental and computational methods. NR binds to hydrophobic locations in ApoMb with higher affinity (K(d) = 25 +/- 5 microM) in the native state than in the molten globule state (K(d) = 52 +/- 5 microM). In the molten globule state, NR is located in a more hydrophobic environment. The dye does not bind to the holoprotein, suggesting that the binding site is located at the heme pocket. In addition to monitoring steady-state properties, the fluorescence emission of NR is capable of tracking submillisecond, time-resolved structural rearrangements of the protein, induced by a nanosecond pH jump. Molecular dynamics simulations were run on ApoMb at neutral pH and at pH 4. The structure obtained for the molten globule state is consistent with the experimentally available structural data. The docking of NR with the crystal structure shows that the ligand binds into the binding pocket of the heme group, with an orientation bringing the planar ring system of NR to overlap with the position of two of the heme porphyrin rings in Mb. The docking of NR with the ApoMb structure at pH 4 shows that the dye binds to the heme pocket with a slightly less favorable binding energy, in keeping with the experimental K(d) value. Under these conditions, NR is positioned in a different orientation, reaching a more hydrophobic environment in agreement with the spectroscopic data.

Animals↗

PubMed2003

Molecular modeling and metabolic studies of the interaction of catechol-O-methyltransferase and a new nitrocatechol inhibitor.

Catechol-O-methyltransferase (COMT, EC 2.1.1.6) plays a central role in the metabolic inactivation of neurotransmitters and neuroactive xenobiotics possessing a catechol motif. 1-(3,4-Dihydroxy-5-nitrophenyl)-2-phenyl-ethanone (BIA 3-202) is a novel nitrocatechol-type inhibitor of COMT, the potential clinical benefit of which is currently being evaluated in the treatment of Parkinson's disease. In the present work we characterize the molecular interactions of BIA 3-202 within the active site of COMT and discuss their implication on the regioselectivity of metabolic O-methylation. Unrestrained flexible-docking simulations suggest that the solution structure of this complex is better described as an ensemble of alternative binding modes, in contrast to the well defined bound configuration revealed by the X-ray structures of related nitrocatechol inhibitors, co-crystallized with COMT. The docking results wherein presented are well supported by experimental evidence, where the pattern of in vitro enzymatic O-methylation and O-demethylation reactions are analyzed. We propose a plausible explanation for the paradoxical in vivo regioselectivity of O-methylation of BIA 3-202, as well as of its related COMT inhibitor tolcapone. Both compounds undergo in vivo O-methylation by COMT at either meta or para catechol hydroxyl groups. However, results herein presented suggest that, in a subsequent step, the p-O-methyl derivatives are selectively demethylated by a microsomal enzyme system. The overall balance is the accumulation of the m-O-methylated metabolites over the para-regioisomers. The implications for the general recognition of nitrocatechol-type inhibitors by COMT and the regioselectivity of their metabolic O-methylation are discussed.

Acetophenones↗

PubMed2001

Rapid atomic density methods for molecular shape characterization.

Two methods for rapid characterization of molecular shape are presented. Both techniques are based on the density of atoms near the molecular surface. The Fast Atomic Density Evaluation (FADE) algorithm uses fast Fourier transforms to quickly estimate densities. The Pairwise Atomic Density Reverse Engineering (PADRE) method derives modified density measures from the relationship between atomic density and total potentials. While many shape-characterization techniques define shape relative to a surface, the descriptors returned by FADE and PADRE can measure local geometry from points within the three-dimensional space surrounding a molecule. The methods can be used to find crevices and protrusions near the surface of a molecule and to test shape complementarity at the interface between docking molecules.

Algorithms↗

PubMed2005

Classification of protein complexes based on docking difficulty.

Based on the results of several groups using different docking methods, the key properties that determine the expected success rate in protein-protein docking calculations are measures of conformational change, interface area, and hydrophobicity. A classification of protein complexes in terms of these measures provides a prediction of docking difficulty. This classification is used to study the targets of the CAPRI docking experiment. Results show that targets with a moderate expected difficulty were indeed predicted well by a number of groups, whereas the use of additional a priori information was necessary to obtain good results for some very difficult targets. The analysis indicates that CAPRI and other relatively large-scale docking studies represent very important steps toward understanding the capabilities and limitations of current protein-protein docking methods.

Algorithms↗

PubMed2006

A model structure for the heterodimer apoA-IMilano-apoA-II supports its peculiar susceptibility to proteolysis.

In this study, we propose a structure for the heterodimer between apolipoprotein A-I(Milano) and apolipoprotein A-II (apoA-I(M)-apoA-II) in a synthetic high-density lipoprotein (HDL) containing L-alpha-palmitoyloleoyl phosphatidylcholine. We applied bioinformatics/computational tools and procedures, such as molecular docking, molecular and essential dynamics, starting from published crystal structures for apolipoprotein A-I and apolipoprotein A-II. Structural and energetic analyses onto the simulated system showed that the molecular dynamics produced a stabilized synthetic HDL. The essential dynamic analysis showed a deviation from the starting belt structure. Our structural results were validated by limited proteolysis experiments on HDL from apoA-I(M) carriers in comparison with control HDL. The high sensitivity of apoA-I(M)-apoA-II to proteases was in agreement with the high root mean-square fluctuation values and the reduction in secondary structure content from molecular dynamics data. Circular dichroism on synthetic HDL containing apoA-I(M)-apoA-II was consistent with the alpha-helix content computed on the proposed model.

Apolipoprotein A-I↗

PubMed2006

A new methodology for the simulation of flexible protein-ligand interactions.

A methodology has been developed for the simulation of induced fit between a ligand and its target protein. It utilizes constrained molecular dynamics where atoms determined to be immobile from difference distance matrix studies are fixed. Application of this methodology to HIV-1 reverse transcriptase (RT) as the example target protein has demonstrated its robustness. Short simulation times are sufficient to achieve good refinement of docking poses resulting from exchange of structurally dissimilar inhibitors between crystal structures.

Alkynes↗

PubMed1990

Does anthramycin bind to Z-DNA as well as to B-DNA? A molecular mechanics study.

Molecular mechanics simulations have been presented on the covalent complexes between anthramycin and the decanucleotides d(GCGCGCGCGC)2 and d(GCGCGTGCGC).d(GCGCACGCGC) in the novel Z form, with the drug docked in the minor groove. The simulations predict that anthramycin binds to Z-DNA almost as well as it does to B-DNA, but with a few key differences in the structural aspects of binding. In particular, the 5' orientation of the drug is preferred in the Z form, contrasting with the preference for 3' orientation in the B form and the drug has a left-handed twist. While no experimental studies have been published on the binding of pyrrolo(1,4)benzodiazepines to left-handed forms of DNA, the energy refined stereochemically satisfactory models provide valuable information which will hopefully simulate high resolution 2-D NMR/NOE studies in solution.

Anthramycin↗

PubMed2002

Modeling the structure of agitoxin in complex with the Shaker K+ channel: a computational approach based on experimental distance restraints extracted from thermodynamic mutant cycles.

Computational methods are used to determine the three-dimensional structure of the Agitoxin (AgTx2)-Shaker complex. In a first stage, a large number of models of the complex are generated using high temperature molecular dynamics, accounting for side chain flexibility with distance restraints deduced from thermodynamic analysis of double mutant cycles. Four plausible binding mode candidates are found using this procedure. In a second stage, the quality and validity of the resulting complexes is assessed by examining the stability of the binding modes during molecular dynamics simulations with explicit water molecules and by calculating the binding free energies of mutant proteins using a continuum solvent representation and comparing with experimental data. The docking protocol and the continuum solvent model are validated using the Barstar-Barnase and the lysozyme-antibody D1.2 complexes, for which there are high-resolution structures as well as double mutant data. This combination of computational methods permits the identification of two possible structural models of AgTx2 in complex with the Shaker K+ channel, additional structural analysis providing further evidence in favor of a single model. In this final complex, the toxin is bound to the extracellular entrance of the channel along the pore axis via a combination of hydrophobic, hydrogen bonding, and electrostatic interactions. The magnitude of the buried solvent accessible area corresponding to the protein-protein contact is on the order of 1000 A with roughly similar contributions from each of the four subunits. Some side chains of the toxin adopt different conformation than in the experimental solution structure, indicating the importance of an induced-fit upon the formation of the complex. In particular, the side chain of Lys-27, a residue highly conserved among scorpion toxins, points deep into the pore with its positively charge amino group positioned at the outer binding site for K+. Specific site-directed mutagenesis experiments are suggested to verify and confirm the structure of the toxin-channel complex.

Animals↗