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 541 records · Page 30Linked to original sources

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↗

PubMed2004

Modeling the E. coli 4-hydroxybenzoic acid oligoprenyltransferase ( ubiA transferase) and characterization of potential active sites.

4-hydroxybenzoate oligoprenyltransferase of E. coli, encoded in the gene ubiA, is an important key enzyme in the biosynthetic pathway to ubiquinone. It catalyzes the prenylation of 4-hydroxybenzoic acid in position 3 using an oligoprenyl diphosphate as a second substrate. Up to now, no X-ray structure of this oligoprenyltransferase or any structurally related enzyme is known. Knowledge of the tertiary structure and possible active sites is, however, essential for understanding the catalysis mechanism and the substrate specificity. With homology modeling techniques, secondary structure prediction tools, molecular dynamics simulations, and energy optimizations, a model with two putative active sites could be created and refined. One active site selected to be the most likely one for the docking of oligoprenyl diphosphate and 4-hydroxybenzoic acid is located near the N-terminus of the enzyme. It is widely accepted that residues forming an active site are usually evolutionary conserved within a family of enzymes. Multiple alignments of a multitude of related proteins clearly showed 100% conservation of the amino acid residues that form the first putative active site and therefore strongly support this hypothesis. However, an additional highly conserved region in the amino acid sequence of the ubiA enzyme could be detected, which also can be considered a putative (or rudimentary) active site. This site is characterized by a high sequence similarity to the aforementioned site and may give some hints regarding the evolutionary origin of the ubiA enzyme. Semiempirical quantum mechanical PM3 calculations have been performed to investigate the thermodynamics and kinetics of the catalysis mechanism. These results suggest a near S(N)1 mechanism for the cleavage of the diphosphate ion from the isoprenyl unit. The 4-hydroxybenzoic acid interestingly appears not to be activated as benzoate anion but rather as phenolate anion to allow attack of the isoprenyl cation to the phenolate, which appeared to be the rate limiting step of the whole process according to our quantum chemical calculations. Our models are a basis for developing inhibitors of this enzyme, which is crucial for bacterial aerobic metabolism. [figure: see text]. Structure of the model of ubiA oligoprenyltransferase derived from the photosynthetic reaction center (1PRC). Putative active amino acid residues and substrates are shown as capped sticks to describe their location and geometry in the putative active sites. The violet spheres identify Mg2+.

Amino Acid Sequence↗

PubMed2005

Coupling between catalytic site and collective dynamics: a requirement for mechanochemical activity of enzymes.

Growing evidence supports the view that enzymatic activity results from a subtle interplay between chemical kinetics and molecular motions. A systematic analysis is performed here to delineate the type and level of coupling between catalysis and conformational mechanics. The dynamics of a set of 98 enzymes representative of different EC classes are analyzed with the Gaussian network model (GNM) and compared with experimental data. In more than 70% of the examined enzymes, the global hinge centers predicted by the GNM are found to be colocalized with the catalytic sites experimentally identified. Low translational mobility (< 7%) is observed for the catalytic residues, consistent with the fine-tuned design of enzymes to achieve precise mechanochemical activities. Ligand binding sites, while closely neighboring catalytic sites, enjoy a moderate flexibility to accommodate the ligand binding. These findings could serve as additional criteria for assessing drug binding residues and could lessen the computational burden of substrate docking searches.

Binding Sites↗

PubMed2005

The relationship between the flexibility of proteins and their conformational states on forming protein-protein complexes with an application to protein-protein docking.

We investigate the extent to which the conformational fluctuations of proteins in solution reflect the conformational changes that they undergo when they form binary protein-protein complexes. To do this, we study a set of 41 proteins that form such complexes and whose three-dimensional structures are known, both bound in the complex and unbound. We carry out molecular dynamics simulations of each protein, starting from the unbound structure, and analyze the resulting conformational fluctuations in trajectories of 5 ns in length, comparing with the structure in the complex. It is found that fluctuations take some parts of the molecules into regions of conformational space close to the bound state (or give information about it), but at no point in the simulation does each protein as whole sample the complete bound state. Subsequent use of conformations from a clustered MD ensemble in rigid-body docking is nevertheless partially successful when compared to docking the unbound conformations, as long as the unbound conformations are themselves included with the MD conformations and the whole globally rescored. For one key example where sub-domain motion is present, a ribonuclease inhibitor, principal components analysis of the MD was applied and was also able to produce conformations for docking that gave enhanced results compared to the unbound. The most significant finding is that core interface residues show a tendency to be less mobile (by size of fluctuation or entropy) than the rest of the surface even when the other binding partner is absent, and conversely the peripheral interface residues are more mobile. This surprising result, consistent across up to 40 of the 41 proteins, suggests different roles for these regions in protein recognition and binding, and suggests ways that docking algorithms could be improved by treating these regions differently in the docking process.

Cluster Analysis↗

PubMed2005

Toward the supramolecular structure of collagen: a molecular dynamics approach.

The structure, stability, and conformational dynamics of an assembly of two pentameric bundles made of collagen-like triple helical segments are explored using 1.2 ns molecular dynamics simulations in three environments: 8.0% (v/v) formaldehyde/water solution, 1.4% (v/v) gallic acid/water solution, and pure water. Stable supramolecular arrangements, where the two collagen units are very close to each other at interacting distances, are identified via docking and energy minimization procedures. Analysis of the interaction with formaldehyde and gallic acid suggests that they perturb the protein in a similar way depending on hydrogen-bonding capability, hydrophobic association properties, and the size and concentration of the compound.

Amino Acids↗

PubMed2004

From conversion to aggregation: protofibril formation of the prion protein.

The ability to diagnose and treat prion diseases is limited by our current understanding of the conversion process of the protein from healthy to harmful isoform. Whereas the monomeric, benign species is well characterized, the misfolded conformations responsible for infectivity and neurodegeneration remain elusive. There is mounting evidence that fibrillization intermediates, or protofibrils, but not mature fibrils or plaques, are the pathogenic species in amyloid diseases. Here, we use molecular dynamics to simulate the conversion of the prion protein. Molecular dynamics simulation produces a scrapie prion protein-like conformation enriched in beta-structure that is in good agreement with available experimental data. The converted conformation was then used to model a protofibril by means of the docking of hydrophobic patches of the template structure to form hydrogen-bonded sheets spanning adjacent subunits. The resulting protofibril model provides a non-branching aggregate with a 3(1) axis of symmetry that is in good agreement with a wide variety of experimental data; importantly, it was derived from realistic simulation of the conversion process.

Amino Acid Sequence↗

PubMed2002

Lead discovery using molecular docking.

As the structures of more and more proteins and nucleic acids become available, molecular docking is increasingly considered for lead discovery. Recent studies consider the hit-rate enhancement of docking screens and the accuracy of docking structure predictions. As more structures are determined experimentally, docking against homology-modeled targets also becomes possible for more proteins. With more docking studies being undertaken, the 'drug-likeness' and specificity of docking hits is also being examined.

Animals↗

PubMed2005

Enhancing the accuracy of virtual screening: molecular dynamics with quantum-refined force fields.

A methodology aimed at improving the accuracy of current docking-scoring procedures is proposed, and validated through detailed tests of its performance in predicting the activity of HIV-1 protease inhibitors. This methodology is based on molecular dynamics simulations using a force field whose effective charges are refined by means of a novel procedure that relies on quantum-mechanical calculations and preserves the internal consistency of the parameterization scheme.

Drug Design↗

PubMed2005

Multiple ligand-binding modes in bacterial R67 dihydrofolate reductase.

R67 dihydrofolate reductase (DHFR), a bacterial plasmid-encoded enzyme associated with resistance to the drug trimethoprim, shows neither sequence nor structural homology with the chromosomal DHFR. It presents a highly symmetrical toroidal structure, where four identical monomers contribute to the unique central active-site pore. Two reactants (dihydrofolate, DHF), two cofactors (NADPH) or one of each (R67*DHF*NADPH) can be found simultaneously within the active site, the last one being the reactive ternary complex. As the positioning of the ligands has proven elusive to empirical determination, we addressed the problem from a theoretical perspective. Several potential structures of the ternary complex were generated using the docking programs AutoDock and FlexX. The variability among the final poses, many of which conformed to experimental data, prompted us to perform a comparative scoring analysis and molecular dynamics simulations to assess the stability of the complexes. Analysis of ligand-ligand and ligand-protein interactions along the 4 ns trajectories of eight different structures allowed us to identify important inter-ligand contacts and key protein residues. Our results, combined with published empirical data, clearly suggest that multipe binding modes of the ligands are possible within R67 DHFR. While the pterin ring of DHF and the nicotinamide ring of NADPH assume a stacked endo-conformation at the centre of the pore, probably assisted by V66, Q67 and I68, the tails of the molecules extend towards opposite ends of the cavity, adopting multiple configurations in a solvent rich-environment where hydrogen-bond interactions with K32 and Y69 may play important roles.

Binding Sites↗

PubMed2006

Flap opening dynamics in HIV-1 protease explored with a coarse-grained model.

We present a one-bead coarse-grained model that enables dynamical simulations of proteins on the time scale of tens of microseconds. The parameterization of the force field includes accurate conformational terms that allow for fast and reliable exploration of the configurational space. The model is applied to the dynamics of flap opening in HIV-1 protease. The experimental structure of the recently crystallized semi-open conformation of HIV-1 protease is well reproduced in the simulation, which supports the accuracy of our model. Thanks to very long simulations and extensive sampling of opening and closing events, we also investigate the thermodynamics and kinetics of the opening process. We have shown that the effect of the solvent slows down the dynamics to the experimentally observed time scales. The model is found to be reliable for application to substrate docking simulations, which are currently in progress.

HIV Protease↗

PubMed2000

Monte Carlo simulations of the chiral recognition of fenoprofen enantiomers by cyclomaltoheptaose (beta-cyclodextrin).

Differential complexation of fenoprofen enantiomers by cyclomaltoheptaose (beta-cyclodextrin) was investigated by Monte Carlo docking simulations. The chiral discrimination of (R)- and (S)-fenoprofen by beta-cyclodextrin was discussed in terms of the difference in the interaction energies and the patterns of molecular interactions. The interaction energies between each enantiomer of fenoprofen and beta-cyclodextrin were consistent with the reported experimental results that showed that the S isomer interacted preferentially with beta-cyclodextrin and was retained longer in a separation process than the R isomer. The thermodynamic preference of inclusion complex formation of (S)-fenoprofen could be explained by the orientation of the phenyl group attached to the chiral carbon, which provided closer contact and thus more favorable intermolecular interactions between the host and guest molecule. The results presented here would be very useful for the prediction of chiral recognition ability of beta-cyclodextrin.

Carbohydrate Conformation↗

PubMed2005

Molecular modeling of wild-type and D816V c-Kit inhibition based on ATP-competitive binding of ellipticine derivatives to tyrosine kinases.

The D816V activating mutation of the c-Kit kinase domain often causes human mastocytosis. Although inhibitors of wild-type c-Kit are known (e.g. STI-571), they are at least 10 times less active against the c-Kit mutant. Several derivatives of ellipticine (5,11-dimethyl-6H-pyrido[3,4-b]carbazole), substituted at positions 1, 2, 9, and 11, were found to inhibit purified D816V and wild-type c-Kit kinase domains with comparable potencies by competing with ATP binding. We investigated the difference between these inhibitors by modeling the D816V mutation in crystal structures of inactive and active c-Kit. Molecular dynamics simulations strongly suggested that the D816V point mutation shifts the conformational equilibrium of c-Kit kinase domain toward the active conformation. All ellipticine compounds were subsequently docked to the D816V mutant c-Kit model. The model provides possible explanations for the structure-activity relationships observed among ellipticine compounds, resulting in new insights into D816V c-Kit mutant inhibition.

Adenosine Triphosphate↗

PubMed2006

Finding a needle in a haystack: development of a combinatorial virtual screening approach for identifying high specificity heparin/heparan sulfate sequence(s).

We describe a combinatorial virtual screening approach for predicting high specificity heparin/heparan sulfate sequences using the well-studied antithrombin-heparin interaction as a test case. Heparan sulfate hexasaccharides were simulated in the 'average backbone' conformation, wherein the inter-glycosidic bond angles were held constant at the mean of the known solution values, irrespective of their sequence. Molecular docking utilized GOLD with restrained inter-glycosidic torsions and intra-ring conformations, but flexible substituents at the 2-, 3-, and 6-positions and explicit incorporation of conformational variability of the iduronate residues. The approach reproduces the binding geometry of the sequence-specific heparin pentasaccharide to within 2.5 A. Screening of a combinatorial virtual library of 6,859 heparin hexasaccharides using a dual filter strategy, in which predicted antithrombin affinity was the first filter and self-consistency of docking was the second, resulted in only 10 sequences. Of these, nine were found to bind antithrombin in a manner identical to the natural pentasaccharide, while a novel hexasaccharide bound the inhibitor in a unique but dramatically different geometry and orientation. This work presents the first approach on combinatorial library screening for heparin/heparan sulfate GAGs to determine high specificity sequences and opens up huge opportunities to investigate numerous other physiologically relevant GAG-protein interactions.

Antithrombins↗

PubMed2002

Mapping of stereoselective recognition sites on human serum transferrin by capillary electrophoresis and molecular modelling.

Stereoselective recognition of chiral compounds can be used for mapping of surface interaction sites on proteins. Iron-free human serum transferrin is a suitable chiral selector in capillary electrophoresis used in native form in solution. Separation of optical isomers of tryptophan-methylester, tryptophan-ethylester and tryptophan-butylester and various drugs were studied in capillary zone electrophoresis applying a distinct transferrin zone prior to sample injection. Changes in the electrophoretic patterns (i.e., in the migration properties) of the molecules reflected the possible interactions with the protein. The tryptophan derivatives and eight drugs possessed stereoselective interactions, seven drugs showed interactions without appreciable chiral separation, and the others did not present any direct complexation with the protein molecules. Molecular modelling was performed to characterize the binding areas at the iron binding site of iron-free transferrin. The docking of tryptophan derivatives on transferrin showed that the R-enantiomers possess a stronger complexation with transferrin, whereas the S-enantiomers are bound by weaker interactions, which is in excellent agreement with the capillary electrophoresis results, where the R-enantiomers were always retarded stronger by transferrin. A ranking of drugs by the lipo score parameter of the docking shows an accordance with the stereoselective interactions by the protein.

Binding Sites↗

PubMed1995

Molecular recognition of the inhibitor AG-1343 by HIV-1 protease: conformationally flexible docking by evolutionary programming.

BACKGROUND: An important prerequisite for computational structure-based drug design is prediction of the structures of ligand-protein complexes that have not yet been experimentally determined by X-ray crystallography or NMR. For this task, docking of rigid ligands is inadequate because it assumes knowledge of the conformation of the bound ligand. Docking of flexible ligands would be desirable, but requires one to search an enormous conformational space. We set out to develop a strategy for flexible docking by combining a simple model of ligand-protein interactions for molecular recognition with an evolutionary programming search technique. RESULTS: We have developed an intermolecular energy function that incorporates steric and hydrogen-bonding terms. The parameters in this function were obtained by docking in three different protein systems. The effectiveness of this method was demonstrated by conformationally flexible docking of the inhibitor AG-1343, a potential new drug against AIDS, into HIV-1 protease. For this molecule, which has nine rotatable bonds, the crystal structure was reproduced within 1.5 A root-mean-square deviation 34 times in 100 simulations, each requiring eight minutes on a Silicon Graphics R4400 workstation. The energy function correctly evaluates the crystal structure as the global energy minimum. CONCLUSIONS: We believe that a solution of the docking problem may be achieved by matching a simple model of molecular recognition with an efficient search procedure. The necessary ingredients of a molecular recognition model include only steric and hydrogen-bond interaction terms. Although these terms are not necessarily sufficient to predict binding affinity, they describe ligand-protein interactions faithfully enough to enable a docking program to predict the structure of the bound ligand. This docking strategy thus provides an important tool for the interdisciplinary field of rational drug design.

Biological Evolution↗

PubMed1995

Monte Carlo docking with ubiquitin.

The development of general strategies for the performance of docking simulations is prerequisite to the exploitation of this powerful computational method. Comprehensive strategies can only be derived from docking experiences with a diverse array of biological systems, and we have chosen the ubiquitin/diubiquitin system as a learning tool for this process. Using our multiple-start Monte Carlo docking method, we have reconstructed the known structure of diubiquitin from its two halves as well as from two copies of the uncomplexed monomer. For both of these cases, our relatively simple potential function ranked the correct solution among the lowest energy configurations. In the experiments involving the ubiquitin monomer, various structural modifications were made to compensate for the lack of flexibility and for the lack of a covalent bond in the modeled interaction. Potentially flexible regions could be identified using available biochemical and structural information. A systematic conformational search ruled out the possibility that the required covalent bond could be formed in one family of low-energy configurations, which was distant from the observed dimer configuration. A variety of analyses was performed on the low-energy dockings obtained in the experiment involving structurally modified ubiquitin. Characterization of the size and chemical nature of the interface surfaces was a powerful adjunct to our potential function, enabling us to distinguish more accurately between correct and incorrect dockings. Calculations with the structure of tetraubiquitin indicated that the dimer configuration in this molecule is much less favorable than that observed in the diubiquitin structure, for a simple monomer-monomer pair. Based on the analysis of our results, we draw conclusions regarding some of the approximations involved in our simulations, the use of diverse chemical and biochemical information in experimental design and the analysis of docking results, as well as possible modifications to our docking protocol.

Algorithms↗

PubMed1999

Exploring the conformational diversity of loops on conserved frameworks.

Loops are structurally variable regions, but the secondary structural elements bracing loops are often conserved. Motifs with similar secondary structures exist in the same and different protein families. In this study, we made an all-PDB-based analysis and produced 495 motif families accessible from the Internet. Every motif family contains some variable loops spanning a common framework (a pair of secondary structures). The diversity of loops and the convergence of frameworks were examined. In addition, we also identified 119 loops with conformational changes in different PDB files. These materials can give some directions for functional loop design and flexible docking.

Computer Simulation↗

PubMed2003

Lipid bilayer simulations of CXCR4 with inverse agonists and weak partial agonists.

CXCR4 is a G protein-coupled receptor (GPCR) that has multiple critical functions in normal and pathologic physiology that include regulation of the metastatic behavior of mammary carcinoma, and utilization as a coreceptor for infection by T-tropic strains of human immunodeficiency virus-1. Molecular dynamic simulations of the rhodopsin-based homology model of CXCR4 were performed in a solvated lipid bilayer to reproduce the microenvironment of this integral membrane protein. The amino acids in CXCR4 necessary for interaction with an inverse agonist, T140, and a weak partial agonist, AMD3100, identified by alanine scanning mutants, were spatially consistent when computationally docked. Whereas T140 binds residues in extracellular domains and regions of the hydrophobic core proximal to the cell surface, amino acids in the central hydrophobic core are critical to binding of AMD3100. The physical localization of T140 binding to CXCR4 by biochemical analyses corroborated the molecular and computational approaches. The structural basis for the interaction of T140 and AMD3100 with CXCR4 confirms that the mechanisms used by these agents are different. This complementary utilization of molecular, physical, and computation analysis provides a powerful approach to elucidate GPCR conformation.

Amino Acid Substitution↗