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 523 records · Page 29Linked to original sources

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↗

PubMed2000

Model peptide studies of sequence repeats derived from the intracrystalline biomineralization protein, SM50. II. Pro,Asn-rich tandem repeats.

In the biomineralization process, a number of Pro-rich proteins participate in the formation of three-dimensional supramolecular structures. One such protein superfamily, the Pro,Gly-rich sea urchin intracrystalline spicule matrix proteins, form protein-protein supramolecular assemblies that modify the microstructure of the inorganic mineral phase (calcite) within embryonic sea urchin spicules and adult sea urchin spines. These proteins represent a useful model for understanding Pro sequence usage and the resulting generation of extended or "open" structures for protein-protein and/or protein-crystal recognition. In the sea urchin spicule matrix protein, SM50 (Strongylocentrotus purpuratus), there exists an unusual 20-residue Pro,Asn-containing repeat, &bond;PNNPNNPNPNNPNNPNNPNPbond which links the upstream 15-residue C-terminal domain and the downstream 211-residue beta-spiral repeat domain. To define the structural preferences of this 20-residue repeat, we created a 20-residue N- and C-terminal "capped" peptidomimetic of this sequence. Using far-uv CD dichroism, CH(alpha) and alpha-(15)N conformational shifts, (3)J(NH-CHalpha) coupling constants, sequential d(NN(i, i + 1)) rotating frame nuclear Overhauser effect connectivities, d(alphaN(i, i + 1))/d(NN(i, i + 1)) intensity ratios, amide temperature shift coefficients, amide solvent exchange, and simulated annealing refinement protocols, we have determined that this 20-residue repeat motif adopts an extended "twist" structure consisting of turn- and coil-like regions. These findings are consistent with previous studies, which have shown that Pro-rich tandem repeats adopt extended, flexible structures in solution. We hypothesize that this 20-residue repeat may fulfill the role of a mineral-binding domain, a protein-protein docking domain, or as an internal "molecular spacer" for the SM50 protein during spicule biocomposite formation.

Amino Acid Motifs↗

PubMed2003

Strategy of computer-aided drug design.

Modern strategies of computer-aided drug design (CADD) are reviewed. The task of CADD in the pipeline of drug discovery is accelerating of finding the new lead compounds and their structure optimization for the following pharmacological tests. The main directions in CADD are based on the availability of the experimentally determined three-dimensional structure of the target macromolecule. If spatial structure is known the methods of structure-based drug design are used. In the opposite case the indirect methods of CADD based on the structures of known ligands (ligand-based drug design) are used. The interrelationship between the main directions of CADD is reviewed. The main CADD approaches of molecule de novo design and database mining are described. They include methods of molecular docking, de novo design, design of pharmacophore and quantity structure-activity relationship models. New ways and perspectives of CADD are discussed.

Computational Biology↗

PubMed1997

Construction of a 3D model of cytochrome P450 2B4.

A three-dimensional structural model of rabbit phenobarbital-inducible cytochrome P450 2B4 (LM2) was constructed by homology modeling techniques previously developed for building and evaluating a 3D model of the cytochrome P450choP isozyme. Four templates with known crystal structures including cytochrome P450cam, terp, BM-3 and eryF were used in multiple sequence alignments and construction of the cytochrome P450 2B4 coordinates. The model was evaluated for its overall quality using available protein analysis programs and found to be satisfactory. The model structure was stable at room temperature during a 140 ps unconstrained full protein molecular dynamics simulation. A putative substrate access channel and binding site were identified. Two different substrates, benzphetamine and androstenedione, that are metabolized by cytochrome P450 2B4 with pronounced product specificity were docked into the putative binding site. Two orientations were found for each substrate that could lead to the observed preferred products. Using a geometric fit method three regions on the surface of the model cytochrome P450 structure were identified as possible sites for interaction with cytochrome b5, a redox partner of P450 2B4. Residues that may interact with the substrates and with cytochrome b5 have been identified and mutagenesis studies are currently in progress.

Amino Acid Sequence↗

PubMed2000

Enzyme selectivity of new cyclooxygenase-2/5 lipoxygenase inhibitors using molecular modeling approach.

We have studied the conformational flexibility of three 5-keto-substituted 7-tert-butyl-2,3-dihydro-3,3-dimethylbenzofurans (DHDMBFs) which show dual cyclooxygenase (COX) and 5-lipoxygenase (LOX) inhibition and are potential candidates as antiinflammatory agents and analgesics. The conformations were studied by systematic search, molecular mechanics (MM) and simulated annealing molecular dynamics (SAMD) techniques. We also studied several structure based parameters and distribution of molecular electrostatic potential (MEP) around these molecules. All the three compounds were docked in the active cavity of cyclooxygenase-2 (COX-2) using graphical and energy grid search techniques. The complex geometries were optimized by MM. The results on conformational flexibility, inter-atomic distances and angles, MEP distribution and points of contacts with peptide side chains in active cavity have been used to understand the mechanistic cause of differential action of these molecules.

Arachidonate 5-Lipoxygenase↗

PubMed2003

Molecular determinants of KCNQ1 channel block by a benzodiazepine.

KCNQ1 channels underlie the slow delayed rectifier K+ current, mediate repolarization of cardiac action potentials, and are a potential therapeutic target for treatment of arrhythmia. (E)-(+)-N-[(3R)-2,3-dihydro-1-methyl-2-oxo-5-phenyl-1H-1,4-benzodiazepin-3-yl]-3-(2,4-dichlorophenyl)-2-propenamide [L-735821 (L-7)] is a potent blocker of KCNQ1 channels. Here we describe the structural determinants of KCNQ1 that are critical for high-affinity block by L-7 using site-directed mutagenesis to alter specific residues and voltage clamp to record channel currents in Xenopus laevis oocytes. Chimeric channels were constructed by combination of regions from L-7-sensitive KCNQ1 and L-7-insensitive KCNQ2 channel subunits. This approach localized the drug interaction site to the pore and S6 domains of KCNQ1. Substitution of single amino acids identified Thr-312 of the pore domain and Ile-337, Phe-339, Phe-340, and Ala-344 of the S6 domain as the most important molecular determinants of channel block. Some mutations also altered the inactivation properties of KCNQ1, but there was no correlation between extent of inactivation and sensitivity to block by L-7. Modeling was used to simulate the docking of L-7 to the KCNQ1 channel pore. The docking was consistent with our experimental data and predicts that L-7 blocks K+ conductance by physically precluding the occupancy of a K+ ion to a pore helix-coordinated site within the central hydrated cavity, a crucial step in ion permeation.

Amino Acid Sequence↗

PubMed2006

A combination of docking/dynamics simulations and pharmacophoric modeling to discover new dual c-Src/Abl kinase inhibitors.

A computational protocol was applied to identify molecular scaffolds untested toward the c-Src tyrosine kinase. A combination of docking and dynamics calculations allowed us to build three-dimensional models of the complexes between Src and several of its known inhibitors. Interactions most contributing to activity of the inhibitors, in terms of hydrogen bonds and hydrophobic contacts, were codified into pharmacophoric models that were in turn applied to perform a search of commercially available compounds within the Asinex database. As a result, we identified 1,3,4-thiadiazoles and pyrazolydine-3,5-diones showing inhibitory activity in the submicromolar range in a cell-free assay toward Src. Moreover, since several of the compounds used to generate pharmacophores were also known as Abl inhibitors, we tested the identified hits toward Abl tyrosine kinase, finding activity in the submicromolar range. Such biological data suggested that the computational protocol is an efficient tool for identifying new hits toward both Src and Abl.

CSK Tyrosine-Protein Kinase↗

PubMed2001

High-resolution probing of local conformational changes in proteins by the use of multiple labeling: unfolding and self-assembly of human carbonic anhydrase II monitored by spin, fluorescent, and chemical reactivity probes.

Two different spin labels, N-(1-oxyl-2,2,5,5-tetramethyl-3-pyrrolidinyl)iodoacetamide (IPSL) and (1-oxyl-2,2,5,5-tetramethylpyrroline-3-methyl) methanethiosulfonate (MTSSL), and two different fluorescent labels 5-((((2-iodoacetyl)amino)-ethyl)amino)naphtalene-1-sulfonic acid (IAEDANS) and 6-bromoacetyl-2-dimetylaminonaphtalene (BADAN), were attached to the introduced C79 in human carbonic anhydrase (HCA II) to probe local structural changes upon unfolding and aggregation. HCA II unfolds in a multi-step manner with an intermediate state populated between the native and unfolded states. The spin label IPSL and the fluorescent label IAEDANS reported on a substantial change in mobility and polarity at both unfolding transitions at a distance of 7.4-11.2 A from the backbone of position 79. The shorter and less flexible labels BADAN and MTSSL revealed less pronounced spectroscopic changes in the native-to-intermediate transition, 6.6-9.0 A from the backbone. At intermediate guanidine (Gu)-HCl concentrations the occurrence of soluble but irreversibly aggregated oligomeric protein was identified from refolding experiments. At approximately 1 M Gu-HCl the aggregation was found to be essentially complete. The size and structure of the aggregates could be varied by changing the protein concentration. EPR measurements and line-shape simulations together with fluorescence lifetime and anisotropy measurements provided a picture of the self-assembled protein as a disordered protein structure with a representation of both compact as well as dynamic and polar environments at the site of the molecular labels. This suggests that a partially folded intermediate of HCA II self-assembles by both local unfolding and intermolecular docking of the intermediates vicinal to position 79. The aggregates were determined to be 40-90 A in diameter depending on the experimental conditions and spectroscopic technique used.

Amino Acid Substitution↗

PubMed1995

Dramatic differences in the motions of the mouth of open and closed cytochrome P450BM-3 by molecular dynamics simulations.

Molecular dynamics trajectories were calculated separately for each of the two molecules in the asymmetric unit of the crystal structure of the hemoprotein domain of cytochrome P450BM-3. Each simulation was 200 ps in length and included a 10 A layer of explicit solvent. The simulated time-average structure of each P450BM-3 molecule is closer to its crystal structure than the two molecular dynamics time-averaged structures are to each other. In the crystal structure, molecule 2 has a more accessible substrate binding pocket than molecule 1, and this difference is maintained throughout the simulations presented here. In particular, the substrate docking regions of molecule 1 and molecule 2 diverge in the solution state simulations. The mouth of the substrate binding pocket is significantly more mobile in the simulation of molecule 2 than in the simulation of molecule 1. For molecule 1, the width of the mouth is only slightly larger than its X-ray value of 8.7 A and undergoes fluctuations of about 1 A. However, in molecule 2, the mouth of the substrate binding pocket is dramatically more open in the time-average molecular dynamics structure (14.7 A) than in the X-ray structure (10.9 A). Furthermore, this region of the protein undergoes large amplitude motions during the trajectory that are not seen in the trajectory of molecule 1, repeatedly opening and closing up to 7 A. Presumably, the binding of different substrates will induce the mouth region to adopt different conformations from within the wide range of structures that are accessible.

Binding Sites↗

PubMed1995

Principles of protein-protein recognition from structure to thermodynamics.

Specific recognition is illustrated by X-ray structures of protease-inhibitor, antigen-antibody and other high affinity complexes including five electron transfer complexes. We attempt to give a physical definition to affinity and specificity on the basis of these data. In a protein-protein complex, specific recognition results from the assembly of complementary surfaces into well-packed interfaces that cover about 1500 A2 and contain about ten hydrogen bonds. These interfaces are larger than between molecules in protein crystals, and smaller than between subunits in oligomeric proteins. We relate the size and chemical nature of interfaces in complexes to the thermodynamical parameters that characterize affinity: the heat capacity and free enthalpy (Gibbs energy) of dissociation at equilibrium, the activation free enthalpy for the dissociation reaction. The same structural and thermodynamical parameters are inadequate for representing the specificity of recognition. We propose instead to describe specificity with the help of statistical physics, and we illustrate the application of the random energy model to antigen-antibody recognition by analyzing results of computer simulations by docking.

Chemical Phenomena↗

PubMed2003

A structural model of agonist binding to the alpha3beta4 neuronal nicotinic receptor.

(alpha3)2(beta4)3 is the most abundant type of neuronal nicotinic ACh receptor (nAChR) mediating cholinergic actions on the autonomic nervous system. Studies to refine or devise drugs selectively acting on (alpha3)2(beta4)3 receptors would benefit from a detailed description of the hitherto unclear agonist-binding domain. The present study reports a three-dimensional model for the ligand-binding domain (LBD) of this receptor either in its unoccupied or agonist-bound conformation. The receptor model was based on the structure of the acetylcholine-binding protein (AChBP), and was obtained using molecular modelling techniques. ACh, nicotine and cytisine (full agonists), muscarine (a selective agonist of muscarinic ACh receptors) and the allosteric modulator eserine were docked into the binding pockets of the receptor model. Simulated agonist-receptor complexes were compared with the agonist-binding complex of the AChBP, as well as of the (alpha4)2(beta2)3 type of nAChR, which is the commonest in the brain. Agonist docking identified discrete amino-acid residues of the beta subunits important for pharmacological selectivity of nAChRs. The predicted affinity of muscarine for the nAChR was low, suggesting the present model to be suitable for effective discrimination of nicotinic agonist binding versus nonselective cholinergic binding. Furthermore, the current model outlined a potential binding site for the allosteric modulator eserine, the site of action of which has remained elusive. The present LBD model of the receptor in its free state provides a novel structural framework to interpret experimental observations and a useful template for future investigations to develop (alpha3)2(beta4)3-selective ligands.

Amino Acid Sequence↗

PubMed2001

Elucidation of the mode of substrate binding to hydroxynitrile lyase from Hevea brasiliensis.

The hydroxynitrile lyase from Hevea brasiliensis (Hb-HNL) is used as a catalyst in enantiospecific syntheses of alpha-hydroxynitriles from aldehydes and methyl-ketones. The catalyzed reaction represents one of the few industrially relevant examples of enzyme mediated C-C coupling reactions. In this work, we modeled Hb-HNL substrate complexes that have as yet proven inaccessible to experimental structure determination and were able to identify two binding modes for the natural substrate acetone cyanohydrin in docking simulations. Discrimination of the two alternatives was achieved by modeling complexes with two different chiral cyanohydrins followed by an analysis of the respective relative binding energies from molecular mechanics and thermodynamic integration. Only for one of the alternative binding modes the experimentally established S-selectivity of the enzyme was correctly predicted. Our results yielded further support for an enzymatic mechanism involving the catalytic triad Ser80, His235, and Asp207 as a general acid/base. A pivotal role was ascribed to Lys236, which seems to be crucial for enzymatic activity at low pH values. In addition, the modeling calculations provided possible explanations for the observed substrate and enantioselectivity of the enzyme that rationalize available mutational data and will be the basis for future protein engineering efforts.

Aldehyde-Lyases↗

PubMed1997

Flexible protein-ligand docking by global energy optimization in internal coordinates.

Eight protein-ligand complexes were simulated by using global optimization of a complex energy function, including solvation, surface tension, and side-chain entropy in the internal coordinate space of the flexible ligand and the receptor side chains [Abagyan, R.A., Totrov, M.M. J. Mol. Biol. 235: 983-1002, 1994]. The procedure uses two types of efficient random moves, a pseudobrownian positional move [Abagyan, R.A., Totrov, M.M., Kuznetsov, D.A. J. Comp. Chem. 15:488-506, 1994] and a Biased-Probability multitorsion move [Abagyan, R.A., Totrov, M.M. J. Mol. Biol. 235: 983-1002, 1994], each accompanied by full local energy minimization. The best docking solutions were further ranked according to the interaction energy, which included intramolecular deformation energies of both receptor and ligand, the interaction energy, surface tension, side-chain entropic contribution, and an electrostatic term evaluated as a boundary element solution of the Poisson equation with the molecular surface as a dielectric boundary. The geometrical accuracy of the docking solutions ranged from 30% to 70% according to the relative displacement error measure at a 1.5 A scale. Similar results were obtained when the explicit receptor atoms were replaced with a grid potential.

Computer Simulation↗

PubMed2003

Molecular simulation of a single-chain antibody against AChE to explore molecular basis of inhibitory effect of 3F3 McAb on enzyme activity.

AIM: To explore the molecular basis of the inhibitory effect of 3F3, a monoclonal antibody against acetylcholinesterase (AChE), by computer-aided molecular simulation. METHODS: The single-chain 3F3 antibody (Sc3F3) was designed by joining VH and VL via a flexible linker (Gly4Ser)3. The amino acid sequence of the recombinant Sc3F3 was then subjected to computer-aided molecular modeling, and docking with the antigen molecule AChE to mimic the immunoactive interaction in a three-dimensional fashion. RESULTS: The modeled structure of Sc3F3 manifested the common features of a classical antibody. Both VH and VL were composed of two ?-sheets and connecting loops. The docking profile of the action between Sc3F3 with AChE demonstrated the formation of a stable structure. The van der Waals force played an important role suggesting that the complex was formed mainly via hydrophobic interactions between Sc3F3 and AChE molecules. CONCLUSION: The spatial structure of the complex of Sc3F3 and AChE showed that Sc3F3 overlaid the entrance of the active center gorge of AChE blocking the access of substrate.

Acetylcholinesterase↗