PubMed Health⌕ Search

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 703 records · Page 39Linked to original sources
PubMed2006

Improved lead-finding for kinase targets using high-throughput docking.

Protein kinases represent a major class of drug targets for the pharmaceutical industry, and the identification of kinase inhibitors with novel and diverse chemotypes is therefore a high priority. Virtual screening methods are a primary source for the discovery of lead molecules for drug development, with high-throughput docking algorithms being among the most extensively used of these methods. A number of developments in virtual screening technologies have resulted in more effective computer-based compound screening of potential kinase inhibitors. These developments include the following: (i) conformational search methods for pose generation; (ii) improvements in the prediction of protein-ligand binding energy through the use of scoring functions; (iii) the use of interaction filters for identifying ligand poses with known binding determinants; and (iv) the impact of binding site flexibility upon high-throughput docking success. This review discusses the application of these methods in the context of the discovery of kinase adenosine triphosphate antagonists.

Animals↗
PubMed2004

Binding crevice for TT-232 in a homology model of type 1 somatostatin receptor.

Somatostatin receptor type 1 was modelled based on the atomic structure of bovine rhodopsin. Possible ways of binding interaction between somatostatin receptor type 1 and TT-232, a cycloheptapeptide analogue of somatostatin with broad therapeutic potential, were analysed by molecular docking. The twelve TT-232 conformations, obtained by NMR measurements in H(2)O-D(2)O mixture, were similar, disclosing a consensus backbone conformation. Several residues interacting with TT-232, such as Val133, Asp137 (helix 3), Arg197 (helix 4), Phe287, Gln291, Asn294 (helix 6), Ser305, and Tyr313 (helix 7), were found. In accordance, in vitro binding experiments indicated high-affinity binding of TT-232 to (125)I labelled somatostatin sites in brain membranes. The single binding crevice obtained by docking may allow the design and discovery of new peptidomimetics of TT-232 in the future.

Amino Acid Sequence↗
PubMed2006

DOCKGROUND resource for studying protein-protein interfaces.

MOTIVATION: Public resources for studying protein interfaces are necessary for better understanding of molecular recognition and developing intermolecular potentials, search procedures and scoring functions for the prediction of protein complexes. RESULTS: The first release of the DOCKGROUND resource implements a comprehensive database of co-crystallized (bound-bound) protein-protein complexes, providing foundation for the upcoming expansion to unbound (experimental and simulated) protein-protein complexes, modeled protein-protein complexes and systematic sets of docking decoys. The bound-bound part of DOCKGROUND is a relational database of annotated structures based on the Biological Unit file (Biounit) provided by the RCSB as a separated file containing probable biological molecule. DOCKGROUND is automatically updated to reflect the growth of PDB. It contains 67,220 pairwise complexes that rely on 14,913 Biounit entries from 34,778 PDB entries (January 30, 2006). The database includes a dynamic generation of non-redundant datasets of pairwise complexes based either on the structural similarity (SCOP classification) or on user-defined sequence identity. The growing DOCKGROUND resource is designed to become a comprehensive public environment for developing and validating new methodologies for modeling of protein interactions. AVAILABILITY: DOCKGROUND is available at http://dockground.bioinformatics.ku.edu. The current first release implements the bound-bound part.

Binding Sites↗
PubMed1993

3D-structures of the lipase from Rhizomucor miehei at different temperatures and computer modelling of a complex of the lipase with trilaurylglycerol.

The lipase of Rhizomucor miehei was measured at two different wavelengths (synchroton lambda = 1.06 A and at a rotating anode lambda = 1.542 A) and two temperatures (293 and 110 K). The structure could be solved using the C alpha-coordinates from the Brookhaven data base (code name 1TGL). Both structures were refined to R-values of 18.7% (synchrotron 293 K) and 20.0% (rotating anode 110 K) at a resolution of 8-2.3 and 8-2.5 A, respectively. Both structures are almost identical to the original 1TGL data set. The side chain positions of both crystal structures differ mainly in the parts of the molecule with relatively high temperature factors in the crystallographic refinement. The low-temperature structure forms more hydrogen bridges than the room temperature structure. From the results of Brzozowski et al. (1991), a preliminary model of the active RML was constructed. A docking of trilaurylglycerol led to a first model of the complex of RML/substrate.

Amino Acid Sequence↗
PubMed2006

Blind docking of drug-sized compounds to proteins with up to a thousand residues.

Blind docking was introduced for the detection of possible binding sites and modes of peptide ligands by scanning the entire surface of protein targets. In the present study, the method is tested on a group of drug-sized compounds and proteins with up to a thousand amino acid residues. Both proteins from complex structures and ligand-free proteins were used as targets. Robustness, limitations and future perspectives of the method are discussed. It is concluded that blind docking can be used for unbiased mapping of the binding patterns of drug candidates.

Binding Sites↗
PubMed2002

Docking multiple conformations of a flexible ligand into a protein binding site using NMR restraints.

A method is described for docking a large, flexible ligand using intra-ligand conformational restraints from exchange-transferred NOE (etNOE) data. Numerous conformations of the ligand are generated in isolation, and a subset of representative conformations is selected. A crude model of the protein-ligand complex is used as a template for overlaying the selected ligand structures, and each complex is conformationally relaxed by molecular mechanics to optimize the interaction. Finally, the complexes were assessed for structural quality. Alternative approaches are described for the three steps of the method: generation of the initial docking template; selection of a subset of ligand conformations; and conformational sampling of the complex. The template is generated either by manual docking using interactive graphics or by a computational grid-based search of the binding site. A subset of conformations from the total number of peptides calculated in isolation is selected based on either low energy and satisfaction of the etNOE restraints, or a cluster analysis of the full set. To optimize the interactions in the complex, either a restrained Monte Carlo-energy minimization (MCM) protocol or a restrained simulated annealing (SA) protocol were used. This work produced 53 initial complexes of which 8 were assessed in detail. With the etNOE conformational restraints, all of the approaches provide reasonable models. The grid-based approach to generate an initial docking template allows a large volume to be sampled, and as a result, two distinct binding modes were identified for a fifteen-residue peptide binding to an enzyme active site.

Anion Exchange Protein 1, Erythrocyte↗
PubMed2000

The third-dimensional structure of the complex between an Fv antibody fragment and an analogue of the main immunogenic region of the acetylcholine receptor: a combined two-dimensional NMR, homology, and molecular modeling approach.

Binding of autoantibodies to the acetylcholine receptor (AChR) plays a major role in the autoimmune disease Myasthenia gravis (MG). In this paper, we propose a structure model of a putative immunocomplex that gives rise to the reduction of functional AChR molecules during the course of MG. The model complex consists of the [G(70), Nle(76)] decapeptide analogue of the main immunogenic region (MIR), representing the major antigenic epitope of AChR, and the single chain Fv fragment of monoclonal antibody 198, a potent MG autoantibody. The structure of the complexed decapeptide antigen [G(70), Nle(76)]MIR was determined using two-dimensional nmr, whereas the antibody structure was derived by means of homology modeling. The final complex was constructed using calculational docking and molecular dynamics. We termed this approach "directed modeling," since the known peptide structure directs the prestructured antibody binding site to its final conformation. The independently derived structures of the peptide antigen and antibody binding site already showed a high degree of surface complementarity after the initial docking calculation, during which the peptide was conformationally restrained. The docking routine was a soft algorithm, applying a combination of Monte Carlo simulation and energy minimization. The observed shape complementarity in the docking process suggested that the structure assessments already led to anti-idiotypic conformations of peptide antigen and antibody fragment. Refinement of the complex by dynamic simulation yielded improved surface adaptation by small rearrangements within antibody and antigen. The complex presented herein was analyzed in terms of antibody-antigen interactions, properties of contacting surfaces, and segmental mobility. The structural requirements for AChR complexation by autoantibodies were explored and compared with experimental data from alanine scans of the MIR peptides. The analysis revealed that the N-terminal loop of the peptide structure, which is indispensable for antibody recognition, aligns three hydrophobic groups in a favorable arrangement leading to the burial of 40% of the peptide surface in the binding cleft upon complexation. These data should be valuable in the rational design of an Fv mutant with much improved affinity for the MIR and AChR to be used in therapeutic approaches in MG.

Amino Acid Sequence↗
PubMed2003

MD simulations of anthrax edema factor: calmodulin complexes with mutations in the edema factor "switch a" region and docking of 3'-deoxy ATP into the adenylyl cyclase active site of wild-type and mutant edema factor variants.

Bacillus anthracis, a spore-forming infectious bacterium, produces an exotoxin, called the edema factor (EF), that functions in part by disrupting internal signalling pathways. When complexed with human host cell calmodulin (CaM), EF becomes an active adenylyl cyclase, producing the internal signal substance cyclic-AMP in an uncontrolled fashion. Recently, the crystal structures for uncomplexed EF and EF:CaM complexes in the presence and absence of a substrate analog (3'-deoxy-ATP), were reported. EF mutational studies have implicated a number of residues important in CaM binding and/or in the generation of the adenylyl cyclase active site, formed by the movements of the EF switch A, B and C regions upon CaM binding. Here we report on the results of molecular dynamics (MD) simulations on two EF:CaM complexes, one containing wild-type EF and the other containing EF in which a cluster of residues in the switch A region (L523, K525, Q526 and V529) have been mutated to alanine. The switch A mutations cause a large increase in the flexibility of the switch C region, the rupture of a number of EF-CaM interactions, an expansion of the carboxyl-terminal domain of CaM, and a change in the Ca(2+) ion binding abilities of the CaM that is in complex with EF. The results indicate the importance of the mutated switch A residues in maintaining a compact EF:CaM complex that appears to be a prerequisite for the generation of a fully-functional adenylyl cyclase active site. The effects of mutating key residues (K346, K353, H577, E588, D590 and N639) in the active site region of EF (to alanine) on the ability of EF to bind the 3'-deoxy-ATP substrate analog were also examined. Active-site residue substitutions at positions 583 (N583A) and 577 (H577A) were found to be particularly disruptive for the placement of the adenine ring moiety into the position found in the x-ray crystal structure of the ligand-protein complex.

Adenylyl Cyclases↗
PubMedDate not supplied

VISCANA: visualized cluster analysis of protein-ligand interaction based on the ab initio fragment molecular orbital method for virtual ligand screening.

We have developed a visualized cluster analysis of protein-ligand interaction (VISCANA) that analyzes the pattern of the interaction of the receptor and ligand on the basis of quantum theory for virtual ligand screening. Kitaura et al. (Chem. Phys. Lett. 1999, 312, 319-324.) have proposed an ab initio fragment molecular orbital (FMO) method by which large molecules such as proteins can be easily treated with chemical accuracy. In the FMO method, a total energy of the molecule is evaluated by summation of fragment energies and interfragment interaction energies (IFIEs). In this paper, we have proposed a cluster analysis using the dissimilarity that is defined as the squared Euclidean distance between IFIEs of two ligands. Although the result of an ordered table by clustering is still a massive collection of numbers, we combine a clustering method with a graphical representation of the IFIEs by representing each data point with colors that quantitatively and qualitatively reflect the IFIEs. We applied VISCANA to a docking study of pharmacophores of the human estrogen receptor alpha ligand-binding domain (57 amino acid residues). By using VISCANA, we could classify even structurally different ligands into functionally similar clusters according to the interaction pattern of a ligand and amino acid residues of the receptor protein. In addition, VISCANA could estimate the correct docking conformation by analyzing patterns of the receptor-ligand interactions of some conformations through the docking calculation.

Binding Sites↗
PubMed2001

Accurate prediction of the bound conformation of galanthamine in the active site of Torpedo californica acetylcholinesterase using molecular docking.

The alkaloid (-)-galanthamine is known to produce significant improvement of cognitive performances in patients with the Alzheimer's disease. Its mechanism of action involves competitive and reversible inhibition of acetylcholinesterase (AChE). Herein, we correctly predict the orientation and conformation of the galanthamine molecule in the active site of AChE from Torpedo californica (TcAChE) using a combination of rigid docking and flexible geometry optimization with a molecular mechanics force field. The quality of the predicted model is remarkable, as indicated by the value of the RMS deviation of approximately 0.5A when compared with the crystal structure of the TcAChE-galanthamine complex. A molecular model of the complex between TcAChE and a galanthamine derivative, SPH1107, with a long chain substituent on the nitrogen has been generated as well. The side chain of this ligand is predicted to extend along the enzyme active site gorge from the anionic subsite, at the bottom, to the peripheral anionic site, at the top. The docking procedure described in this paper can be applied to produce models of ligand-receptor complexes for AChE and other macromolecular targets of drug design.

Acetylcholinesterase↗
PubMed1998

Structural features of the combining site region of Erythrina corallodendron lectin: role of tryptophan 135.

The role of Trp 135 and Tyr 108 in the combining site of Erythrina corallodendron lectin (ECorL) was investigated by physicochemical characterization of mutants obtained by site-directed mutagenesis, hemagglutination-inhibition studies, and molecular modeling, including dynamics simulations. The findings demonstrate that Trp 135 in ECorL: (1) is required for the tight binding of Ca2+ and Mn2+ to the lectin because mutation of this residue into alanine results in loss of these ions upon dialysis and concomitant reversible inactivation of the mutant; (2) contributes to the high affinity of methyl alpha-N-dansylgalactosaminide (MealphaGalNDns) to the lectin; and (3) is solely responsible for the fluorescence energy transfer between the aromatic residues of the lectin and the dansyl group in the ECorL-MealphaGalNDns complex. Docking of MealphaGalNDns into the combining site of the lectin reveals that the dansyl moiety is parallel with the indole of Trp 135, as required for efficient fluorescence energy transfer, in one of the two possible conformations that this ligand assumes in the bound state. In the W135A mutant, which still binds MealphaGalNDns strongly, the dansyl group may partially insert itself into the place formerly occupied by Trp 135, a process that from dynamics simulations does not appear to be energetically favored unless the loop containing this residue assumes an open conformation. However, a small fraction of the W135A molecules must be able to bind MealphaGalNDns in order to explain the relatively high affinity, as compared to galactose, still remaining for this ligand. A model for the molecular events leading to inactivation of the W135A mutant upon demetallization is also presented in which the cis-trans isomerization of the Ala 88-Asp 89 peptide bond, observed in high-temperature dynamics simulations, appears not to be a required step.

Acetylgalactosamine↗
PubMed2006

AClAP, Autonomous hierarchical agglomerative Cluster Analysis based protocol to partition conformational datasets.

MOTIVATION: Sampling the conformational space is a fundamental step for both ligand- and structure-based drug design. However, the rational organization of different molecular conformations still remains a challenge. In fact, for drug design applications, the sampling process provides a redundant conformation set whose thorough analysis can be intensive, or even prohibitive. We propose a statistical approach based on cluster analysis aimed at rationalizing the output of methods such as Monte Carlo, genetic, and reconstruction algorithms. Although some software already implements clustering procedures, at present, a universally accepted protocol is still missing. RESULTS: We integrated hierarchical agglomerative cluster analysis with a clusterability assessment method and a user independent cutting rule, to form a global protocol that we implemented in a MATLAB metalanguage program (AClAP). We tested it on the conformational space of a quite diverse set of drugs generated via Metropolis Monte Carlo simulation, and on the poses we obtained by reiterated docking runs performed by four widespread programs. In our tests, AClAP proved to remarkably reduce the dimensionality of the original datasets at a negligible computational cost. Moreover, when applied to the outcomes of many docking programs together, it was able to point to the crystallographic pose. AVAILABILITY: AClAP is available at the "AClAP" section of the website http://www.scfarm.unibo.it.

Algorithms↗
PubMed2005

Molecular docking of the highly hypolipidemic agent alpha-asarone with the catalytic portion of HMG-CoA reductase.

Docking experiments using a number of published crystal structures of HMG-CoA reductase with the potent hypocholesterolemic agent alpha-asarone are described. The results indicate that alpha-asarone binds in the enzyme's active site. The methoxy groups play a key role in the binding and probably also in its biological activity, as shown by extensive SAR studies reported for analogues of alpha-asarone. The docking results will be valuable for the structure-based design of novel hypolipidemic agents.

Allylbenzene Derivatives↗
PubMedDate not supplied

CoMFA 3D-QSAR analysis of HIV-1 RT nonnucleoside inhibitors, TIBO derivatives based on docking conformation and alignment.

HIV-1 RT is one of the key enzymes in the duplication of HIV-1. Inhibitors of HIV-1 RT are classified as nonnucleoside RT inhibitors (NNRTIs) and nucleoside analogues. NNRTIs bind in a region not associated with the active site of the enzyme. Within the NNRTI category, there is a set of inhibitors commonly referred to as TIBO inhibitors. Fifty TIBO inhibitors were used in the work to build 3-D QSAR models. The two known crystal structures of complexes are used to investigate and validate the docking protocol. The results show that the docking simulations reproduce the crystal complexes very well with RMSDs of approximately 1 A and approximately 0.6 A for 1REV and 1COU, respectively. The alignment of molecules and "active" conformation selection are the key to a successful 3D-QSAR model by CoMFA. The flexible docking (Autodock3) was used on determination of "active" conformation and molecular alignment, and CoMFA and CoMSIA were used to develop 3D-QSAR models of 50 TIBOs in the work. The 3D-QSAR models demonstrate a good ability to predict the activity of studied compounds (r2 = 0.972, 0.944, q2 = 0.704, 0.776). It is shown that the steric and electrostatic properties predicted by CoMFA contours can be related to the binding structure of the complex. The results demonstrate that the combination of ligand-based and receptor-based modeling is a powerful approach to build 3D-QSAR models.

Benzodiazepines↗
PubMed2006

Crystal structures of cytochrome c(L) and methanol dehydrogenase from Hyphomicrobium denitrificans: structural and mechanistic insights into interactions between the two proteins.

Methanol dehydrogenase (Hd-MDH) and its physiological electron acceptor, cytochrome c(L) (Hd-Cyt c(L)), isolated from a methylotrophic denitrifying bacterium, Hyphomicrobium denitrificans A3151, have been kinetically and structurally characterized; the X-ray structures of Hd-MDH and Hd-Cyt c(L) have been determined using molecular replacement at 2.5 and 2.0 A resolution, respectively. To explain the mechanism for electron transfer between these proteins, the dependence of MDH activity on the concentration of Hd-Cyt c(L) has been investigated at pH 4.5-7.0. The Michaelis constant for Hd-Cyt c(L) shows the smallest value (approximately 0.3 microM) at pH 5.5. The pseudo-first-order rate constant (k(obs)) of the reduction of Hd-Cyt c(L) exhibits a hyperbolic concentration dependence of Hd-MDH at pH 5.5, although k(obs) linearly increases at pH 6.5. These findings indicate formation of a transient complex between these proteins during an electron transfer event. Hd-MDH (148 kDa) is a large tetrameric protein with an alpha(2)beta(2) subunit composition, showing a high degree of structural similarity with other MDHs. Hd-Cyt c(L) (19 kDa) exhibiting the alpha-band at 550.7 nm has a unique C-terminal region involving a disulfide bond between Cys47 and Cys165. Moreover, there is a pair of Hd-Cyt c(L) monomers related with a pseudo-2-fold axis of symmetry in the asymmetric unit, and the two monomers tightly interact with each other through three hydrogen bonds. This configuration is the first example in the studies of cytochrome c as the physiological electron acceptor for MDH. The docking simulation between the coupled Hd-Cyt c(L) molecules and the heterotetrameric Hd-MDH molecule has been carried out.

Alcohol Oxidoreductases↗
PubMed2006

Molecular characterization of a disease associated conformational epitope on GAD65 recognised by a human monoclonal antibody b96.11.

Autoantibodies to the 65kDa isoform of glutamate decarboxylase (GAD65) are associated with type I diabetes and recognise highly conformational epitope(s) that remain to be defined. The human recombinant Fab from mAb b96.11 inhibits binding of most GAD65 antibody positive sera from patients and its epitope has previously been localized to the middle region of GAD65. Recent studies indicate that b96.11 antibody specificity predicts the risk of developing type 1 diabetes in prediabetic individuals. We describe the use homology modelling, protein-protein docking simulations and biopanning of random peptide phage displayed libraries with b96.11 to predict contact amino acids on the interface of GAD65/Fab b96.11 complex. Further analysis by in vitro mutagenesis of GAD65 followed by radioimmunoprecipitation refined the amino acids contributing to the b96.11 epitope. Our studies show an interface characterized by a protruding antibody-combining site centered on the long heavy chain CDR3 loop of Fab b96.11 establishing interactions with the critical residue Phe(344) in the core of the epitope on GAD65, surrounded by charged sites within (375)RK(376) and (305)DER(307). The epitope requires residues from both middle and the C-terminal domains, and is the first precise definition of an epitope on GAD65. The nature of the b96.11 epitope leads to considerations of potential structural variations for differences in antigenicity between the isoforms GAD65 and GAD67. The study shows the utility of using a combination of in silico techniques and experimental data for molecular characterization and localization of conformational epitopes for which crystal structures are lacking.

Amino Acid Sequence↗
PubMed2005

Approaching the CAPRI challenge with an efficient geometry-based docking.

The last 3 rounds (3-5) of CAPRI included a wide range of docking targets. Several targets were especially challenging, since they involved large-scale movements and symmetric rearrangement, while others were based on homology models. We have approached the targets with a variety of geometry-based docking algorithms that include rigid docking, symmetric docking, and flexible docking with symmetry constraints. For all but 1 docking target, we were able to submit at least 1 acceptable quality prediction. Here, we detail for each target the prediction methods used and the specific biological data employed, and supply a retrospective analysis of the results. We highlight the advantages of our techniques, which efficiently exploit the geometric shape complementarity properties of the interaction. These enable them to run only few minutes on a standard PC even for flexible docking, thus proving their scalability toward computational genomic scale experiments. We also outline the major required enhancements, such as the introduction of side-chain position refinement and the introduction of flexibility for both docking partners.

Algorithms↗