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

X-ray structure of antistasin at 1.9 A resolution and its modelled complex with blood coagulation factor Xa.

The three-dimensional structure of antistasin, a potent inhibitor of blood coagulation factor Xa, from the Mexican leech Haementeria officinalis was determined at 1.9 A resolution by X-ray crystallography. The structure reveals a novel protein fold composed of two homologous domains, each resembling the structure of hirustasin, a related 55-residue protease inhibitor. However, hirustasin has a different overall shape than the individual antistasin domains, it contains four rather than two beta-strands, and does not inhibit factor Xa. The two antistasin domains can be subdivided into two similarly sized subdomains with different relative orientations. Consequently, the domain shapes are different, the N-terminal domain being wedge-shaped and the C-terminal domain flat. Docking studies suggest that differences in domain shape enable the N-terminal, but not C-terminal, domain of antistasin to bind and inhibit factor Xa, even though both have a very similar reactive site. Furthermore, a putative exosite binding region could be defined in the N-terminal domain of antistasin, comprising residues 15-17, which is likely to interact with a cluster of positively charged residues on the factor Xa surface (Arg222/Lys223/Lys224). This exosite binding region explains the specificity and inhibitory potency of antistasin towards factor Xa. In the C-terminal domain of antistasin, these exosite interactions are prevented due to the different overall shape of this domain.

Amino Acid Sequence↗

Synthesis, biological evaluation, and quantitative receptor docking simulations of 2-[(acylamino)ethyl]-1,4-benzodiazepines as novel tifluadom-like ligands with high affinity and selectivity for kappa-opioid receptors.

The synthesis and biological evaluation of a series of 2-substituted 5-phenyl-1,4-benzodiazepines, structurally related to tifluadom (5), the only benzodiazepine that acts simultaneously as a kappa-opioid agonist and a cholecystokinin-A (CCK-A) antagonist, are reported. The radioligand binding models used in these studies were [(125)I](BH)-CCK-8 in rat pancreas (CCK-A), [(3)H]-(MENLE(28,31))-cck-8 in guinea pig cerebral cortex (CCK-B), and [(3)H]U-69593 (kappa(1)), [(3)H]DAMGO (mu), and [(3)H]DADLE (delta) in guinea pig brain. All the title compounds were devoid of significant affinity for both CCK-A and CCK-B receptors, while some of them bound with nanomolar affinity and high selectivity for kappa-opioid receptors. In particular, the 2-thienyl derivative 7A(X = H) with a K(i) = 0.50 nM represents a clear improvement with respect to tifluadom, showing a comparable potency but higher selectivity. The application of computational simulations and linear regression analysis techniques to the complexes between guinea pig kappa (kappa(1))-receptor and the title compounds allowed the identification of the structural determinants for recognition and quantitative elucidation of the structure-affinity relationships in this class of receptors.

Amino Acid Sequence↗

A model for short alpha-neurotoxin bound to nicotinic acetylcholine receptor from Torpedo californica: comparison with long-chain alpha-neurotoxins and alpha-conotoxins.

Short-chain alpha-neurotoxins from snakes are highly selective antagonists of the muscle-type nicotinic acetylcholine receptors (nAChR). Although their spatial structures are known and abundant information on topology of binding to nAChR is obtained by labeling and mutagenesis studies, the accurate structure of the complex is not yet known. Here, we present a model for a short alpha-neurotoxin, neurotoxin II from Naja oxiana (NTII), bound to Torpedo californica nAChR. It was built by comparative modeling, docking and molecular dynamics using 1H NMR structure of NTII, cross-linking and mutagenesis data, cryoelectron microscopy structure of Torpedo marmorata nAChR [Unwin, N., 2005. Refined structure of the nicotinic acetylcholine receptor at 4A resolution. J. Mol. Biol. 346, 967-989] and X-ray structures of acetylcholine-binding protein (AChBP) with agonists [Celie, P.H., van Rossum-Fikkert, S.E., van Dijk, W.J., Brejc, K., Smit, A.B., Sixma, T.K., 2004. Nicotine and carbamylcholine binding to nicotinic acetylcholine receptors as studied in AChBP crystal structures. Neuron 41 (6), 907-914] and antagonists: alpha-cobratoxin, a long-chain alpha-neurotoxin [Bourne, Y., Talley, T.T., Hansen, S.B., Taylor, P., Marchot, P., 2005. Crystal structure of Cbtx-AChBP complex reveals essential interactions between snake alpha-neurotoxins and nicotinic receptors. EMBO J. 24 (8), 1512-1522] and alpha-conotoxin [Celie, P.H., Kasheverov, I.E., Mordvintsev, D.Y., Hogg, R.C., van Nierop, P., van Elk, R., van Rossum-Fikkert, S.E., Zhmak, M.N., Bertrand, D., Tsetlin, V., Sixma, T.K., Smit, A.B., 2005. Crystal structure of nicotinic acetylcholine receptor homolog AChBP in complex with an alpha-conotoxin PnIA variant. Nat. Struct. Mol. Biol. 12 (7), 582-588]. In complex with the receptor, NTII was located at about 30 A from the membrane surface, the tip of its loop II plunges into the ligand-binding pocket between the alpha/gamma or alpha/delta nAChR subunits, while the loops I and III contact nAChR by their tips only in a 'surface-touch' manner. The toxin structure undergoes some changes during the final complex formation (for 1.45 rmsd in 15-25 ps according to AMBER'99 molecular dynamics simulation), which correlates with NMR data. The data on the mobility and accessibility of spin- and fluorescence labels in free and bound NTII were used in MD simulations. The binding process is dependent on spontaneous outward movement of the C-loop earlier found in the AChBP complexes with alpha-cobratoxin and alpha-conotoxin. Among common features in binding of short- and long alpha-neurotoxins is the rearrangement of aromatic residues in the binding pocket not observed for alpha-conotoxin binding. Being in general very similar, the binding modes of short- and long alpha-neurotoxins differ in the ways of loop II entry into nAChR.

Amino Acid Sequence↗

Structure-based 3D QSAR and design of novel acetylcholinesterase inhibitors.

The paper describes the construction, validation and application of a structure-based 3D QSAR model of novel acetylcholinesterase (AChE) inhibitors. Initial use was made of four X-ray structures of AChE complexed with small, non-specific inhibitors to create a model of the binding of recently developed aminopyridazine derivatives. Combined automated and manual docking methods were applied to dock the co-crystallized inhibitors into the binding pocket. Validation of the modelling process was achieved by comparing the predicted enzyme-bound conformation with the known conformation in the X-ray structure. The successful prediction of the binding conformation of the known inhibitors gave confidence that we could use our model to evaluate the binding conformation of the aminopyridazine compounds. The alignment of 42 aminopyridazine compounds derived by the docking procedure was taken as the basis for a 3D QSAR analysis applying the GRID/GOLPE method. A model of high quality was obtained using the GRID water probe, as confirmed by the cross-validation method (q2LOO = 0.937, q2L50%O = 0.910). The validated model, together with the information obtained from the calculated AChE-inhibitor complexes, were considered for the design of novel compounds. Seven designed inhibitors which were synthesized and tested were shown to be highly active. After performing our modelling study the X-ray structure of AChE complexed with donepezil, an inhibitor structurally related to the developed aminopyirdazines, has been made available. The good agreement found between the predicted binding conformation of the aminopyridazines and the one observed for donepezil in the crystal structure further supports our developed model.

Acetylcholinesterase↗

Arginine 336 and asparagine 333 of the human cholecystokinin-A receptor binding site interact with the penultimate aspartic acid and the C-terminal amide of cholecystokinin.

The cholecystokinin-A receptor (CCK-AR) is a G protein-coupled receptor that mediates important central and peripheral cholecystokinin actions. Residues of the CCK-AR binding site that interact with the C-terminal part of CCK that is endowed with biological activity are still unknown. Here we report on the identification of Arg-336 and Asn-333 of CCK-AR, which interact with the Asp-8 carboxylate and the C-terminal amide of CCK-9, respectively. Identification of the two amino acids was achieved by dynamics-based docking of CCK in a refined three-dimensional model of CCK-AR using, as constraints, previous results that demonstrated that Trp-39/Gln-40 and Met-195/Arg-197 interact with the N terminus and the sulfated tyrosine of CCK, respectively. Arg-336-Asp-8 and Asn-333-amide interactions were pharmacologically assessed by mutational exchange of Arg-336 and Asn-333 in the receptor or reciprocal elimination of the partner chemical functions in CCK. This study also allowed us to demonstrate that (i) the identified interactions are crucial for stabilizing the high affinity phospholipase C-coupled state of the CCK-AR.CCK complex, (ii) Arg-336 and Asn-333 are directly involved in interactions with nonpeptide antagonists SR-27,897 and L-364,718, and (iii) Arg-336 but not Asn-333 is directly involved in the binding of the peptide antagonist JMV 179 and the peptide partial agonist JMV 180. These data will be used to obtain an integrated dynamic view of the molecular processes that link agonist binding to receptor activation.

Amides↗

Docking and small angle X-ray scattering studies of purine nucleoside phosphorylase.

Docking simulations have been used to assess protein complexes with some success. Small angle X-ray scattering (SAXS) is a well-established technique to investigate protein spatial configuration. This work describes the integration of geometric docking with SAXS to investigate the quaternary structure of recombinant human purine nucleoside phosphorylase (PNP). This enzyme catalyzes the reversible phosphorolysis of N-ribosidic bonds of purine nucleosides and deoxynucleosides. A genetic deficiency due to mutations in the gene encoding for PNP causes gradual decrease in T-cell immunity. Inappropriate activation of T-cells has been implicated in several clinically relevant human conditions such as transplant rejection, rheumatoid arthritis, lupus, and T-cell lymphomas. PNP is therefore a target for inhibitor development aiming at T-cell immune response modulation and has been submitted to extensive structure-based drug design. The present analysis confirms the trimeric structure observed in the crystal. The potential application of the present procedure to other systems is discussed.

Humans↗

Origins and predictions of stereoselective antibody catalysis: theoretical analysis of Diels-Alder catalysis by 39A11 and its germ-line antibody.

The binding of enantiomeric haptens and transition states by the Schultz Diels-Alderase antibody 39A11 and its germ-line antibody were studied theoretically. The mechanisms by which one hapten and one transition state stereoisomer is recognized selectively are explored with docking simulations and quantum mechanical models. Transition states of the relevant Diels-Alder reaction were located with density functional theory. A prediction is made that the stereoselectivity of 39A11 will be achieved by two strategically placed hydrogen bonds and pi-stacking interactions of the maleimide with a binding-site tryptophan, arranged so as to coordinate one enantiomeric transition state. Binding of other ligands by antibody 39A11 and the germ-line antibody has also been investigated. The polyspecific nature of 39A11 and its germ-line precursor was found to originate from the general ability of the binding pockets to achieve hydrophobic binding of small organic substrates. Comparison of the highly homologous progesterone and Diels-Alderase antibodies (DB3, 1E9, and 39A11) highlights the fact that differences of several key residues in the binding pockets are sufficient to confer selectivity for different antigens.

Amino Acid Sequence↗

Enantioselective substrate binding in a monooxygenase protein model by molecular dynamics and docking.

The two-component flavoenzyme styrene monooxygenase (SMO) is an efficient alternative to several chemical epoxidation catalysts on a preparative scale. A first homology model of the catalytic domain (StyA) of SMO was constructed (Protein Data Bank ID 2HD8) based on the structure of para-hydroxybenzoate hydroxylase. The StyA protein structure was optimized by restrained molecular dynamics to reproduce specific pre-S binding orientations of styrene. Effects of all 10 point mutations examined were explained by the distance of the site to the styrene and FAD binding sites. Thirteen of 20 ligands could be accommodated in a catalytically active binding orientation, and predicted affinities correlated well with experimental turnover and inhibition. The binding cavity is almost completely hydrophobic except for a hydrogen-bonded network formed by three water molecules, the backbone of residues 300-302, and the flavin ribityl, similar to P293, and three crystal waters in para-hydroxybenzoate hydroxylase suggest that P302, T47, and the waters in StyA are a vital component of the catalytic mechanism. The current optimized and validated StyA model provides a good starting point for elucidation of the structural basis of StyA ligand binding and catalysis. Novel insights in the binding of ligands to SMO/StyA, provided by the current protein model, will aid the rational design of mutants with specific, altered enantioselective properties.

Binding Sites↗

From synthetic coiled coils to functional proteins: automated design of a receptor for the calmodulin-binding domain of calcineurin.

A series of synthetic receptors capable of binding to the calmodulin-binding domain of calcineurin (CN393-414) was designed, synthesized and characterized. The design was accomplished by docking CN393-414 against a two-helix receptor, using an idealized three-stranded coiled coil as a starting geometry. The sequence of the receptor was chosen using a side-chain re-packing program, which employed a genetic algorithm to select potential binders from a total of 7.5x10(6) possible sequences. A total of 25 receptors were prepared, representing 13 sequences predicted by the algorithm as well as 12 related sequences that were not predicted. The receptors were characterized by CD spectroscopy, analytical ultracentrifugation, and binding assays. The receptors predicted by the algorithm bound CN393-414 with apparent dissociation constants ranging from 0.2 microM to >50 microM. Many of the receptors that were not predicted by the algorithm also bound to CN393-414. Methods to circumvent this problem and to improve the automated design of functional proteins are discussed.

Algorithms↗

The binding mode of calyculin A to protein phosphatase-1. A novel spiroketal vector model.

The catalytic subunits of serine/threonine protein phosphatases 1 and 2A are subject to inhibition by various toxins such as the microcystins, the nodularins, okadaic acid, tautomycin, and the calyculins. A recent paper (Bagu, J. R., Sykes, B. D, Craig, M. M., and Holmes, C. F. B. (1997) J. Biol. Chem. 272, 5087-5097) reported the successful docking of the crystal structure of calyculin A to the crystal structure of protein phosphatase-1. Unfortunately, the model presented there is based on the structure of the unnatural enantiomer of calyculin A and must therefore be incorrect. We have developed a spiroketal vector model which appears to account for the spatial orientation of the hydrophobic and basic chains extending from the spiroketal-phosphate core of calyculin A. The model also clearly demonstrates why the unnatural enantiomer of calyculin A does not fit properly into the pocket of the active site. Based on our model, we present a possible open binding mode for calyculin A in the enzyme. This open structure is conceptually similar to the predicted binding mode of the peptide inhibitor DARPP-32 to the enzyme; the hydrophobic, metal-binding, and electrostatic interactions are all retained in this model.

Binding Sites↗

Improved prediction of protein-protein binding sites using a support vector machines approach.

MOTIVATION: Structural genomics projects are beginning to produce protein structures with unknown function, therefore, accurate, automated predictors of protein function are required if all these structures are to be properly annotated in reasonable time. Identifying the interface between two interacting proteins provides important clues to the function of a protein and can reduce the search space required by docking algorithms to predict the structures of complexes. RESULTS: We have combined a support vector machine (SVM) approach with surface patch analysis to predict protein-protein binding sites. Using a leave-one-out cross-validation procedure, we were able to successfully predict the location of the binding site on 76% of our dataset made up of proteins with both transient and obligate interfaces. With heterogeneous cross-validation, where we trained the SVM on transient complexes to predict on obligate complexes (and vice versa), we still achieved comparable success rates to the leave-one-out cross-validation suggesting that sufficient properties are shared between transient and obligate interfaces. AVAILABILITY: A web application based on the method can be found at http://www.bioinformatics.leeds.ac.uk/ppi_pred. The dataset of 180 proteins used in this study is also available via the same web site. CONTACT: westhead@bmb.leeds.ac.uk SUPPLEMENTARY INFORMATION: http://www.bioinformatics.leeds.ac.uk/ppi-pred/supp-material.

Algorithms↗

Solution structure of a brodimoprim analogue in its complex with Lactobacillus casei dihydrofolate reductase.

Two-dimensional (2D) double-quantum-filtered correlation spectroscopy (DQF-COSY), total correlation spectroscopy (TOCSY), nuclear Overhauser effect spectroscopy (NOESY), and rotating-frame NOESY (ROESY) spectra were used to assign essentially all the protons in a 1:1 complex of Lactobacillus casei dihydrofolate reductase formed with an analogue of the antibacterial drug brodimoprim [2,4-diamino-5-(3',5'-dimethoxy-4'-bromobenzyl)pyrimidine]. The analogue has a 4,6-dicarboxylic acid side chain substituted on the 3'-O position designed to interact with the Arg 57 and His 28 residues in L. casei dihydrofolate reductase; it binds a factor of 10(3) more tightly to the enzyme than does the parent compound. Thirty-eight intermolecular and 11 intramolecular NOEs were measured involving the bound brodimoprim-4,6-dicarboxylic acid analogue. These provided the distance constraints used in conjunction with an energy minimization and simulated annealing protocol (using Discover from Biosym Ltd.) to dock the brodimoprim analogue into dihydrofolate reductase. In calculations where side chains and backbone fragments for binding-site residues were allowed flexibility, 90% of the 40 calculated structures had reasonable covalent geometry and none of them had NOE distance violations of greater than 0.36 A. The conformations of the aromatic rings in the bound ligand were well-defined in all the structures, with torsion angles tau 1 = -153 degrees +/- 4 degrees (C4-C5-C7-C1') and tau 2 = 53 degrees +/- 4 degrees (C5-C7-C1'-C2'): the aromatic rings of the ligand occupied essentially the same space in all the calculated structures (root mean square deviation value 1.83 A). Inclusion of the electrostatic interactions into the energy minimizations indicated that structures in which the 4,6-dicarboxylate group of the ligand interacts with the side chains of Arg 57 and His 28 are of low energy. Significant differences in side-chain and backbone conformations were detected between binding-site residues in the enzyme complexes with the brodimorpim analogue and methotrexate.

Amino Acid Sequence↗

Cobatoxin 1 from Centruroides noxius scorpion venom: chemical synthesis, three-dimensional structure in solution, pharmacology and docking on K+ channels.

CoTX1 (cobatoxin 1) is a 32-residue toxin with three disulphide bridges that has been isolated from the venom of the Mexican scorpion Centruroides noxius Hoffmann. Here we report the chemical synthesis, disulphide bridge organization, 3-D (three-dimensional) solution structure determination, pharmacology on K+ channel subtypes (voltage-gated and Ca2+-activated) and docking-simulation experiments. An enzyme-based cleavage of the synthetic folded/oxidized CoTX1 indicated half-cystine pairs between Cys3-Cys22, Cys8-Cys27 and Cys12-Cys29. The 3-D structure of CoTX1 (solved by 1H-NMR) showed that it folds according to the common alpha/beta scaffold of scorpion toxins. In vivo, CoTX1 was lethal after intracerebroventricular injection to mice (LD50 value of 0.5 microg/mouse). In vitro, CoTX1 tested on cells expressing various voltage-gated or Ca2+-activated (IKCa1) K+ channels showed potent inhibition of currents from rat K(v)1.2 ( K(d) value of 27 nM). CoTX1 also weakly competed with 125I-labelled apamin for binding to SKCa channels (small-conductance Ca2+-activated K+ channels) on rat brain synaptosomes (IC50 value of 7.2 microM). The 3-D structure of CoTX1 was used in docking experiments which suggests a key role of Arg6 or Lys10, Arg14, Arg18, Lys21 (dyad), Ile23, Asn24, Lys28 and Tyr30 (dyad) residues of CoTX1 in its interaction with the rat K(v)1.2 channel. In addition, a [Pro7,Gln9]-CoTX1 analogue (ACoTX1) was synthesized. The two residue replacements were selected aiming to restore the RPCQ motif in order to increase peptide affinity towards SKCa channels, and to alter the CoTX1 dipole moment such that it is expected to decrease peptide activity on K(v) channels. Unexpectedly, ACoTX1 exhibited an activity similar to that of CoTX1 towards SKCa channels, while it was markedly more potent on IKCa1 and several voltage-gated K+ channels.

Amino Acid Sequence↗

PROXIMO--a new docking algorithm to model protein complexes using data from radical probe mass spectrometry (RP-MS).

The design and implementation of a new algorithm, known as PROXIMO for protein oxidation interface modeller, is described to predict the structure of protein complexes using data generated in radical probe mass spectrometry (RP-MS) experiments. Photochemical radiolysis and discharge sources can be used to effect RP-MS in which hydroxyl radicals are formed directly from the bulk solvent on millisecond timescales and react with surface accessible residues in footprinting-like experiments. The algorithm utilizes a geometric surface fitting routine to predict likely structures for protein complexes. These structures are scored based on a correlation between the measured solvent accessibility of oxidizable residue side chains and oxidation shielding data obtained by RP-MS. The algorithm has been implemented to predict structures for the ribonuclease S-protein-peptide and calmodulin-melittin complexes using RP-MS data generated in this laboratory. The former is in close agreement with the high-resolution experimental structure available.

Amino Acid Sequence↗

Discovery of new inhibitors of aldose reductase from molecular docking and database screening.

Aldose reductase (ALR2) is a target enzyme for the treatment of diabetic complications. Owing to the limited number of currently available drugs for the treatment of diabetic complications, the discovery of new inhibitors of ALR2 that can potentially be optimized as drugs appears highly desirable. In this study, a molecular docking analysis of the structures of more than 127,000 organic compounds contained in the National Cancer Institute database was performed to find and score molecules that are complementary to ALR2. Besides retrieving several carboxylic acid derivatives, which are known to generally inhibit aldose reductase, docking proposed other families of putative inhibitors such as sulfonic acids, nitro-derivatives, sulfonamides and carbonyl derivatives. Twenty-five compounds, chosen as the highest-scoring representatives of each of these families, were tested as aldose reductase inhibitors. Five of them were found to inhibit aldose reductase in the micromolar range. For these active compounds, selectivity with respect to the closely-related aldehyde reductase was determined by measuring the corresponding inhibitory activities. The structures of the complexes between the new lead inhibitors and aldose reductase, here refined with molecular mechanics and molecular dynamics calculations, suggest that new pharmacophoric groups can bind aldose reductase very efficiently. In the case of the family of the nitro-derivative inhibitors, a class of particularly interesting compounds, a round of optimizations was performed with the synthesis and biological evaluation of a series of derivatives aimed at testing the proposed binding mode and at improving interaction with active site residues. Starting from a hit compound having an IC(50) of 42 microM, the most potent compound synthesized showed a 10-fold increase in inhibitory activity and 10-fold selectivity with respect to ALR1, and structure--activity relationships of the designed compounds were in agreement with the proposed mode of binding at the active site.

Aldehyde Reductase↗

Identification of multiple allosteric sites on the M1 muscarinic acetylcholine receptor.

Staurosporine and four staurosporine derivatives were docked on the rhodopsin-based homology model of the M1 muscarinic acetylcholine receptor in order to localize the possible allosteric sites of this receptor. It was found that there were three major allosteric sites, two of which are located at the extracellular face of the receptor, and one in the intracellular domain of the receptor. In the present study, the localization of these binding sites is described for the first time. The present study confirms the existence of multiple allosteric sites on the M1 muscarinic receptor, and lays the ground for further experimental and computational analysis to better understand how muscarinic receptors are modulated via their allosteric sites. These findings will also help to design and develop novel drugs acting as allosteric modulators of the M1 receptor, which can be used in the treatment of the Alzheimer's disease.

Algorithms↗

Aldolase A Ins(1,4,5)P3-binding domains as determined by site-directed mutagenesis.

We substituted neutral amino acids for some positively charged residues (R42, K107, K146, R148 and K229) that line the active site of aldolase A in an effort to determine binding sites for inositol 1, 4,5-trisphosphate. In addition, D33 (involved in carbon-carbon bond cleavage) was mutated. K229A and D33S aldolases showed almost no catalytic activity, but Ins(1,4,5)P(3) binding was similar to that determined with the use of wild-type aldolase A. R42A, K107A, K146R and R148A had markedly decreased affinities for Ins(1,4,5)P(3) binding, increased EC(50) values for Fru(1,6)P(2)-evoked release of bound Ins(1,4,5)P(3) and increased K(i) values for Ins(1,4, 5)P(3)-evoked inhibition of aldolase activity. K146Q (positive charge removal) had essentially no catalytic activity and could not bind Ins(1,4,5)P(3). Computer-simulated docking of Ins(1,4,5)P(3) in the aldolase A structure was consistent with electrostatic binding of Ins(1,4,5)P(3) to K107, K146, R148, R42, R303 and backbone nitrogens, as has been reported for Fru(1,6)P(2) binding. Results indicate that Ins(1,4,5)P(3) binding occurs at the active site and is not dependent on having a catalytically active enzyme; they also suggest that there is competition between Ins(1,4,5)P(3) and Fru(1, 6)P(2) for binding. Although Ins(1,4,5)P(3) binding to aldolase involved electrostatic interactions, the aldolase A Ins(1,4, 5)P(3)-binding domain did not show other similarities to pleckstrin homology domains or phosphotyrosine-binding domains known to bind Ins(1,4,5)P(3) in other proteins.

Amino Acid Sequence↗

Structure prediction of protein complexes by an NMR-based protein docking algorithm.

Protein docking algorithms can be used to study the driving forces and reaction mechanisms of docking processes. They are also able to speed up the lengthy process of experimental structure elucidation of protein complexes by proposing potential structures. In this paper, we are discussing a variant of the protein-protein docking problem, where the input consists of the tertiary structures of proteins A and B plus an unassigned one-dimensional 1H-NMR spectrum of the complex AB. We present a new scoring function for evaluating and ranking potential complex structures produced by a docking algorithm. The scoring function computes a 'theoretical' 1H-NMR spectrum for each tentative complex structure and subtracts the calculated spectrum from the experimental one. The absolute areas of the difference spectra are then used to rank the potential complex structures. In contrast to formerly published approaches (e.g. [Morelli et al. (2000) Biochemistry, 39, 2530-2537]) we do not use distance constraints (intermolecular NOE constraints). We have tested the approach with four protein complexes whose three-dimensional structures are stored in the PDB data bank (Bernstein et al., 1977) and whose 1H-NMR shift assignments are available from the BMRB database. The best result was obtained for an example, where all standard scoring functions failed completely. Here, our new scoring function achieved an almost perfect separation between good approximations of the true complex structure and false positives.

Algorithms↗