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 1,027 records · Page 57Linked to original sources

PubMed1999

Three-dimensional models of alpha(2A)-adrenergic receptor complexes provide a structural explanation for ligand binding.

We have compared bacteriorhodopsin-based (alpha(2A)-AR(BR)) and rhodopsin-based (alpha(2A)-AR(R)) models of the human alpha(2A)-adrenengic receptor (alpha(2A)-AR) using both docking simulations and experimental receptor alkylation studies with chloroethylclonidine and 2-aminoethyl methanethiosulfonate hydrobromide. The results indicate that the alpha(2A)-AR(R) model provides a better explanation for ligand binding than does our alpha(2A)-AR(BR) model. Thus, we have made an extensive analysis of ligand binding to alpha(2A)-AR(R) and engineered mutant receptors using clonidine, para-aminoclonidine, oxymetazoline, 5-bromo-N-(4, 5-dihydro-1H-imidazol-2-yl)-6-quinoxalinamine (UK14,304), and norepinephrine as ligands. The representative docked ligand conformation was chosen using extensive docking simulations coupled with the identification of favorable interaction sites for chemical groups in the receptor. These ligand-protein complex studies provide a rational explanation at the atomic level for the experimentally observed binding affinities of each of these ligands to the alpha(2A)-adrenergic receptor.

Animals↗

PubMed1999

Theoretical study of the electrostatically driven step of receptor-G protein recognition.

This study proposes a theoretical model describing the electrostatically driven step of the alpha 1 b-adrenergic receptor (AR)-G protein recognition. The comparative analysis of the structural-dynamics features of functionally different receptor forms, i.e., the wild type (ground state) and its constitutively active mutants D142A and A293E, was instrumental to gain insight on the receptor-G protein electrostatic and steric complementarity. Rigid body docking simulations between the different forms of the alpha 1 b-AR and the heterotrimeric G alpha q, G alpha s, G alpha i1, and G alpha t suggest that the cytosolic crevice shared by the active receptor and including the second and the third intracellular loops as well as the cytosolic extension of helices 5 and 6, represents the receptor surface with docking complementarity with the G protein. On the other hand, the G protein solvent-exposed portions that recognize the intracellular loops of the activated receptors are the N-terminal portion of alpha 3, alpha G, the alpha G/alpha 4 loop, alpha 4, the alpha 4/beta 6 loop, alpha 5, and the C-terminus. Docking simulations suggest that the two constitutively active mutants D142A and A293E recognize different G proteins with similar selectivity orders, i.e., G alpha q approximately equal to G alpha s > G alpha i >> G alpha t. The theoretical models herein proposed might provide useful suggestions for new experiments aiming at exploring the receptor-G protein interface.

Amino Acid Sequence↗

PubMed2000

Identification of novel small molecule ligands that bind to tetanus toxin.

Tetanus toxin belongs to a family of clostridial protein neurotoxins for which there are no known antidotes. Another closely related member of this family, botulinum toxin, is being used with increasing frequency by physicians to treat severe muscle disorders. Botulinum toxin has also been produced in large quantities by terrorists for use as a biological weapon. To identify small molecule ligands that might bind to the targeting domain of tetanus and botulinum toxins and to facilitate the design of inhibitors and new reagents for their detection, molecular docking calculations were used to screen a large database of compounds for their potential to bind to the C fragment of tetanus toxin. Eleven of the predicted ligands were assayed by electrospray ionization mass spectrometry (ESI-MS) for binding to the tetanus toxin C fragment, and five ligands (45%) were found to bind to the protein. One of these compounds, doxorubicin, was observed to have strong hydrophobic interactions with the C fragment. To check the ligands for their ability to compete with ganglioside binding, each was also tested using a GT1b liposome assay. Doxorubicin was the only ligand found to competitively bind the tetanus toxin C fragment with an appreciable binding constant (9.4 microM).

Binding Sites↗

PubMed2000

The crystal structure of dihydrofolate reductase from Thermotoga maritima: molecular features of thermostability.

Two high-resolution structures have been obtained for dihydrofolate reductase from the hyperthermophilic bacterium Thermotoga maritima in its unliganded state, and in its ternary complex with the cofactor NADPH and the inhibitor, methotrexate. While the overall fold of the hyperthermophilic enzyme is closely similar to monomeric mesophilic dihydrofolate reductase molecules, its quaternary structure is exceptional, in that T. maritima dihydrofolate reductase forms a highly stable homodimer. Here, the molecular reasons for the high intrinsic stability of the enzyme are elaborated and put in context with the available data on the physical parameters governing the folding reaction. The molecule is extremely rigid, even with respect to structural changes during substrate binding and turnover. Subunit cooperativity can be excluded from structural and biochemical data. Major contributions to the high intrinsic stability of the enzyme result from the formation of the dimer. Within the monomer, only subtle stabilizing interactions are detectable, without clear evidence for any of the typical increments of thermal stabilization commonly reported for hyperthermophilic proteins. The docking of the subunits is optimized with respect to high packing density in the dimer interface, additional salt-bridges and beta-sheets. The enzyme does not show significant structural changes upon binding its coenzyme, NADPH, and the inhibitor, methotrexate. The active-site loop, which is known to play an important role in catalysis in mesophilic dihydrofolate reductase molecules, is rearranged, participating in the association of the subunits; it no longer participates in catalysis.

Amino Acid Sequence↗

PubMed2005

Modeling oligomers with Cn or Dn symmetry: application to CAPRI target 10.

The abundance of oligomeric proteins makes them a frequent target for structure prediction. However, homologous proteins sometimes adopt different oligomerization states, rendering the prediction of structures of whole oligomers beyond the scope of comparative modeling. This obstacle can be overcome by combining comparative modeling of the single subunit of an oligomer with docking techniques, designed for predicting subunit-subunit interfaces. We present here algorithms for predicting the structures of homo-oligomers with C(n) or D(n) (n > 2) symmetry. The prediction procedure includes a symmetry-restricted docking step followed by a C(n) or D(n) oligomer-forming step, in which the dimers from the docking step are assembled to oligomers. The procedure is applied to each of the crystallographically independent subunits in 8 C(n) and 3 D(n) oligomers, producing very accurate predictions. It is further applied to a single monomer of the tick-borne encephalitis virus coat protein E (Target 10 of the CAPRI experiment). The predicted trimer ranked 30, obtained via rigid-body geometric-hydrophobic docking followed by C(n) oligomer formation, is very similar to the experimentally observed trimer formed by domain II of this protein. Furthermore, the predicted trimer formed from the separated domain I is also close to the experimental structure.

Algorithms↗

PubMed1988

Comparison of model and nuclear magnetic resonance structures for the human inflammatory protein C5a.

The model structure previously proposed for human C5a, based upon the crystal structure of the homologous protein human C3a, is compared to the solution structure of human C5a recently determined by nuclear magnetic resonance (NMR) methods in our laboratory. The general folding and helix topography of the C5a protein were modeled very well. The N-terminus, which is disordered in the C3a crystal, was correctly predicted in the C5a model both as to its being a helix and as to its docking site on the rest of the molecule. On the other hand, the NMR data show that the biologically important C-terminal residues are disordered in solution, unlike the model and the C3a crystal structure where this region was helical.

Complement C5↗

PubMed2005

Carvedilol blocks cardiac KATP and KG but not IK1 channels by acting at the bundle-crossing regions.

We examined the effects of carvedilol on cardiac inwardly rectifying K(+) (Kir) channels, i.e., ATP-sensitive (K(ATP)), G-protein-activated (K(G)) and background (I(K1)) Kir channels. We found that carvedilol effectively inhibits K(ATP) and K(G), but not I(K1) channels. Carvedilol inhibits K(ATP) channels reconstituted in HEK293 cells with Kir6.2 lacking the C-terminal 26 amino acids (Kir6.2DeltaC26), suggesting that carvedilol acts in the channel pore. A sequence comparison of the three channels revealed that a cysteine residue, C166, in the inner helix of Kir6.2 is conserved in both Kir6.xs (K(ATP)) and Kir3.xs (K(G)), but not in Kir2.xs (I(K1)). The mutation of this residue (C166A) made Kir6.2DeltaC26 resistant to the drug. Homology modeling and docking simulation suggested that interaction between carvedilol and the pore could be located at the cytosolic portion of the inner helix (bundle-crossing region) containing C166. This study shows that carvedilol blocks specific groups of Kir channels by interacting with the bundle-crossing region.

Amino Acid Sequence↗

PubMed2004

Structural model of the UbcH5B/CNOT4 complex revealed by combining NMR, mutagenesis, and docking approaches.

The protein CNOT4 possesses an N-terminal RING finger domain that acts as an E3 ubiquitin ligase and specifically interacts with UbcH5B, a ubiquitin-conjugating enzyme. The structure of the CNOT4 RING domain has been solved and the amino acids important for the binding to UbcH5B have been mapped. Here, the residues of UbcH5B important for the binding to CNOT4 RING domain were identified by NMR chemical shift perturbation experiments, and these data were used to generate structural models of the complex with the program HADDOCK. Together with the NMR data, additional biochemical data were included in a second docking, and comparisons of the resulting model with the structure of the c-Cbl/UbcH7 complex reveal some significant differences, notably at specific residues, and give structural insights into the E2/E3 specificity.

Amino Acid Sequence↗

PubMed1997

Three-dimensional modelling of human cytochrome P450 1A2 and its interaction with caffeine and MeIQ.

The three-dimensional modelling of proteins is a useful tool to fill the gap between the number of sequenced proteins and the number of experimentally known 3D structures. However, when the degree of homology between the protein and the available 3D templates is low, model building becomes a difficult task and the reliability of the results depends critically on the correctness of the sequence alignment. For this reason, we have undertaken the modelling of human cytochrome P450 1A2 starting by a careful analysis of several sequence alignment strategies (multiple sequence alignments and the TOPITS threading technique). The best results were obtained using TOPITS followed by a manual refinement to avoid unlikely gaps. Because TOPITS uses secondary structure predictions, several methods that are available for this purpose (Levin, Gibrat, DPM, NnPredict, PHD, SOPM and NNSP) have also been evaluated on cytochromes P450 with known 3D structures. More reliable predictions on alpha-helices have been obtained with PHD, which is the method implemented in TOPITS. Thus, a 3D model for human cytochrome P450 1A2 has been built using the known crystal coordinates of P450 BM3 as the template. The model was refined using molecular mechanics computations. The model obtained shows a consistent location of the substrate recognition segments previously postulated for the CYP2 family members. The interaction of caffeine and a carcinogenic aromatic amine (MeIQ), which are characteristic P450 1A2 substrates, has been investigated. The substrates were solvated taking into account their molecular electrostatic potential distributions. The docking of the solvated substrates in the active site of the model was explored with the AUTODOCK programme, followed by molecular mechanics optimisation of the most interesting complexes. Stable complexes were obtained that could explain the oxidation of the considered substrates by cytochrome P450 1A2 and could offer an insight into the role played by water molecules.

Amino Acid Sequence↗

PubMed2000

Three-dimensional modeling of and ligand docking to vitamin D receptor ligand binding domain.

The ligand binding domain of the human vitamin D receptor (VDR) was modeled based on the crystal structure of the retinoic acid receptor. The ligand binding pocket of our VDR model is spacious at the helix 11 site and confined at the beta-turn site. The ligand 1alpha, 25-dihydroxyvitamin D(3) was assumed to be anchored in the ligand binding pocket with its side chain heading to helix 11 (site 2) and the A-ring toward the beta-turn (site 1). Three residues forming hydrogen bonds with the functionally important 1alpha- and 25-hydroxyl groups of 1alpha,25-dihydroxyvitamin D(3) were identified and confirmed by mutational analysis: the 1alpha-hydroxyl group is forming pincer-type hydrogen bonds with S237 and R274 and the 25-hydroxyl group is interacting with H397. Docking potential for various ligands to the VDR model was examined, and the results are in good agreement with our previous three-dimensional structure-function theory.

Amino Acid Sequence↗

PubMed2005

The structural comparison of the bacterial PepX and human DPP-IV reveals sites for the design of inhibitors of PepX activity.

X-prolyl dipeptidyl aminopeptidases (X-PDAP) are enzymes catalysing the release of dipeptides from the amino termini of polypeptides containing a proline or an alanine at the penultimate position. Involved in various mammalian regulation processes, as well as in chronic human diseases, they have been proposed to play a role in pathogenicity for Streptococci. We compared the structure of X-PDAP from Lactococcus lactis (PepX) with its human counterpart DPP-IV. Despite very different overall folds, the residues most implicated for X-PDAP activity are conserved in the same positions and orientations in both enzymes, thus defining a structural signature for the X-PDAP specificity that crosses the species frontiers of evolution. Starting from this observation, we tested some inhibitors of DPP-IV on PepX activity, for which no specific inhibitor is known. We thus found that PepX was highly sensitive to valine-pyrrolidide with a KI of 9.3 microm, close to that reported in DPP-IV inhibition. We finally used the structure of PepX from L. lactis as a template for computer-based homology modeling of PepX from the pathogenic Streptococcus gordonii. Docking simulations of valine-pyrrolidide into the active site of PepX led to the identification of key residues for a rational drug design against PepX from Streptococci. These results could have applications in human health giving new perspectives to the struggle against pathogens.

Adenosine Deaminase↗

PubMed2002

The ATP-Mg2+ binding site and cytoplasmic domain interactions of Na+,K+-ATPase investigated with Fe2+-catalyzed oxidative cleavage and molecular modeling.

This work utilizes Fe(2+)-catalyzed cleavages and molecular modeling to obtain insight into conformations of cytoplasmic domains and ATP-Mg(2+) binding sites of Na(+),K(+)-ATPase. In E(1) conformations the ATP-Fe(2+) complex mediates specific cleavages at 712VNDS (P domain) and near 440VAGDA (N domain). In E(2)(K), ATP-Fe(2+) mediates cleavages near 212TGES (A domain), near 440VAGDA, and between residues 460-490 (N domain). Cleavages at high ATP-Fe(2+) concentrations do not support suggestions for two ATP sites. A new reagent, fluorescein-DTPA, has been synthesized. The fluorescein-DTPA-Fe(2+) complex mediates cleavages similar to those mediated by ATP-Fe(2+). The data suggest the existence of N to P domain interactions in E(1)Na, with bound ATP-Fe(2+) or fluorescein-DPTA-Fe(2+), A-N, and A-P interactions in E(2)(K), and provide testable constraints for model building. Molecular models based on the Ca(2+)-ATPase structure are consistent with the predictions. Specifically, high-affinity ATP-Mg(2+) binding in E(1) is explained with the N domain tilted ca. 80 degrees toward the P domain, by comparison with well-separated N and P domains in the Ca-ATPase crystal structure. With ATP-Mg(2+) docked, bound Mg(2+) is close to both D710 (in 710DGVNDS) and D443 (in 440VAGDASE). D710 is known to be crucial for Mg(2+) binding. The cleavage and modeling data imply that D443 could also be a candidate for Mg(2+) binding. Comparison of E(1).ATP,Mg(2+) and E(2) models suggests an explanation of the high or low ATP affinities, respectively. We propose a scheme of ATP-Mg(2+) and Mg(2+) binding and N, P, and A domain interactions in the different conformations of the catalytic cycle.

Adenosine Triphosphate↗

PubMed2001

Structural analysis of BAG1 cochaperone and its interactions with Hsc70 heat shock protein.

BAG-family proteins share a conserved protein interaction region, called the 'BAG domain', which binds and regulates Hsp70/Hsc70 molecular chaperones. This family of cochaperones functionally regulates signal transducing proteins and transcription factors important for cell stress responses, apoptosis, proliferation, cell migration and hormone action. Aberrant overexpression of the founding member of this family, BAG1, occurs in human cancers. In this study, a structure-based approach was used to identify interacting residues in a BAG1--Hsc70 complex. An Hsc70-binding fragment of BAG1 was shown by multidimensional NMR methods to consist of an antiparallel three-helix bundle. NMR chemical shift experiments marked surface residues on the second (alpha 2) and third (alpha 3) helices in the BAG domain that are involved in chaperone binding. Structural predictions were confirmed by site-directed mutagenesis of these residues, resulting in loss of binding of BAG1 to Hsc70 in vitro and in cells. Molecular docking of BAG1 to Hsc70 and mutagenesis of Hsc70 marked the molecular surface of the ATPase domain necessary for interaction with BAG1. The results provide a structural basis for understanding the mechanism by which BAG proteins link molecular chaperones and cell signaling pathways.

Adenosine Triphosphatases↗

PubMed2005

Synechocystis ferredoxin/ferredoxin-NADP(+)-reductase/NADP+ complex: Structural model obtained by NMR-restrained docking.

Ferredoxin (Fd) and ferredoxin-NADP(+)-reductase (FNR) are two terminal physiological partners of the photosynthetic electron transport chain. Based on a nuclear magnetic resonance (NMR)-restrained-docking approach, two alternative structural models of the Fd-FNR complex in the presence of NADP+ are proposed. The protein docking simulations were performed with the software BiGGER. NMR titration revealed a 1:1 stoichiometry for the complex and allowed the mapping of the interacting residues at the surface of Fd. The NMR chemical shifts were encoded into distance constraints and used with theoretically calculated electronic coupling between the redox cofactors to propose experimentally validated docked complexes.

Amino Acid Sequence↗

PubMed2004

A pharmacophore docking algorithm and its application to the cross-docking of 18 HIV-NNRTI's in their binding pockets.

The docking of small molecules into the binding site of a target protein is an important but difficult step in structure-based drug design. The performance of a docking algorithm is usually evaluated by re-docking ligands into their native binding sites. We have explored the cross-docking of 18 HIV-NNRTIs (non-nucleoside inhibitors of HIV reverse transcriptase) of which the ligand-protein structure has been determined: each of the 18 ligands was docked into each of the 18 binding sites. The docking algorithms studied are an energy-based simulated annealing algorithm and a novel pharmacophore docking algorithm. It turns out that the energy-based docking of the ligands into non-native pockets is far less successful than the docking into their native pockets. The results can be improved by using explicit pharmacophore information, and by docking a ligand into a panel of protein structures and selecting the ligand-protein combination with the lowest interaction energy as the final result.

Algorithms↗

PubMed1994

The conformation of the sialyl Lewis X ligand changes upon binding to E-selectin.

High-field NMR spectroscopy has been used to study the complex formed by the tetrasaccharide sialyl Lewis X and its receptor, E-selectin. Transferred NOEs demonstrate a specific interaction between the protein and ligand and enable measurement of the dissociation constant for the complex to be between approximately 1.1 and 2.0 mM. Differences between Overhauser spectra for free and bound sialyl Lewis X highlight a conformational change upon binding. This can be pinpointed to a change in the torsion angle of the glycosidic link between the sialyl and galactosyl residues and used to select a likely "bound" conformation from four low-energy species. Docking the bound form of sialyl Lewis X onto a model of the lectin domain of E-selectin suggests that the conformational change upon binding results primarily from steric interactions.

Carbohydrate Sequence↗

PubMed2000

Inhibitors of kinesin activity from structure-based computer screening.

Kinesin motor proteins use ATP hydrolysis for transport along microtubules in the cell. We sought to identify small organic ligands to inhibit kinesin's activity. Candidate molecules were identified by computational docking of commercially available compounds using the computer program DOCK. Compounds were docked at either of two sites, and a selection was tested for inhibition of microtubule-stimulated ATPase activity. Twenty-two submillimolar inhibitors were identified. Several inhibitors appeared to be competitive for microtubule binding and not for ATP binding, and three compounds showed 50% inhibition down to single-digit micromolar levels. Most inhibitors grouped into four distinct classes (fluoresceins, phenolphthaleins, anthraquinones, and naphthylene sulfonates). We measured the binding of one inhibitor, rose bengal lactone (RBL), to kinesin (dissociation constant 2.5 microM) by its increase in steady-state fluorescence anisotropy. The RBL binding site on kinesin was localized by fluorescent resonance energy transfer (FRET) using a donor fluorophore (coumarin) covalently attached at unique, surface-exposed cysteine residues engineered at positions 28, 149, 103, 220, or 330. RBL was found to bind in its original docked site: the pocket cradled by loop 8 and beta-strand 5 in kinesin's three-dimensional structure. These results confirm this region's role in microtubule binding and identify this pocket as a novel binding site for kinesin inhibition.

Adenosine Triphosphatases↗

PubMedDate not supplied

Docking of combinatorial peptide libraries into a broadly cross-reactive human IgM.

A monoclonal IgM cryoglobulin with diverse binding behavior was isolated from a patient (Mez) with Waldenström's macroglobulinemia. It gave very high titers in the binding of combinatorially synthesized libraries of peptides ranging in size from two to eight residues. The crystal structure of Mez Fv revealed that the binding site was divided into two cavities of unequal volumes with dimensions and chemical properties that were compatible with the binding of peptides. Access to this unique combination of structural information and peptide binding data led us to carry out Mez-peptide docking simulations to gain insight into the Mez binding propensities. In the present article, the results for docking of five peptide libraries are combined with discussions of the methods and approximations involved in the docking process. We analyze the origins of peptide binding affinity for Mez IgM in terms of its cross-reactivity and its structural preferences.

Antibodies, Monoclonal↗