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,225 records · Page 68Linked to original sources

PubMed1999

Model of maltose-binding protein/chemoreceptor complex supports intrasubunit signaling mechanism.

The Tar protein of Escherichia coli is unique among known bacterial chemoreceptors in that it generates additive responses to two very disparate ligands, aspartate and maltose. Aspartate binds directly to the periplasmic (extracytoplasmic) domain of Tar. Maltose first binds to maltose-binding protein (MBP). MBP then assumes a closed conformation in which it can interact with the periplasmic domain of Tar. MBP residues critical for binding Tar were identified in a screen of mutations that cause specific defects in maltose chemotaxis. Mutations were introduced into a plasmid-borne malE gene that encodes a mutant form of MBP in which two engineered Cys residues spontaneously generate a disulfide bond in the oxidizing environment of the periplasmic space. This disulfide covalently crosslinks the NH3-terminal and COOH-terminal domains of MBP and locks the protein into a closed conformation. Double-Cys MBP confers a dominant-negative phenotype for maltose taxis, and we reasoned that third mutations that relieve this negative dominance probably alter residues that are important for the initial interaction of MBP with Tar. The published three-dimensional structures of MBP and the periplasmic domain of E. coli Tar were docked in a computer simulation that juxtaposed the residues in MBP identified in this way with residues in Tar that have been implicated in maltose taxis. The resulting model of the MBP-Tar complex exhibits good complementarity between the surfaces of the two proteins and supports the idea that aspartate and MBP may each initiate an attractant signal through Tar by inducing similar conformational changes in the chemoreceptor.

ATP-Binding Cassette Transporters↗

PubMed2006

Protein farnesyltransferase: flexible docking studies on inhibitors using computational modeling.

Using MacroModel, peptide, peptidomimetic and non-peptidomimetic inhibitors of the zinc metalloenzyme, farnesyltransferase (FTase), were docked into the enzyme binding site. Inhibitor flexibility, farnesyl pyrophosphate substrate flexibility, and partial protein flexibility were taken into account in these docking studies. In addition to CVFM and CVIM, as well as our own inhibitors FTI-276 and FTI-2148, we have docked other farnesyltransferase inhibitors (FTIs) including Zarnestra, which presently is in advanced clinical trials. The AMBER* force field was employed, augmented with parameters that were derived for zinc. A single binding site model that was derived from the crystal structure of CVFM complexed with farnesyltransferase and farnesylpyrophosphate was used for these studies. The docking results using the lowest energy structure from the simulation, or one of the lowest energy structures, were generally in excellent agreement with the X-ray structures. One of the most important findings of this study is that numerous alternative conformations for the methionine side chain can be accommodated by the enzyme suggesting that the methionine pocket can tolerate groups larger than methionine at the C-terminus of the tetrapeptide and suggesting alternative locations for the placement of side chains that may improve potency.

Alkyl and Aryl Transferases↗

PubMed1996

Structure-function relations of antithrombin III-heparin interactions as assessed by biophysical and biological assays and molecular modeling of peptide-pentasaccharide-docked complexes.

The serine proteinase inhibitor antithrombin III (ATIII) is a key regulatory protein of intrinsic blood coagulation. ATIII attains its full biological activity only upon binding polysulfated glycosaminoglycans, such as heparin. Peptide K121-A134, based on the sequence of ATIII in the D helix region, was previously shown by us (Tyler-Cross et al., Protein Sci. 3, 620-627, 1994) to encompass part (or all) of the purported high-affinity heparin binding region of ATIII. A series of peptide analogs has now been prepared whose sequences are identical to K121-A134 except that single cationic residues of K121-A134 have been successively replaced with Ala. In one case, the Arg residue of the reference peptide corresponding to R129 of ATIII has been replaced by Gln (R129deltaQ peptide), thus mimicking the naturally occurring mutant protein, ATIII Geneva. The heparin affinity of all peptides was quantitated by isothermal titration calorimetry and by peptide/ATIII competition binding assays. Replacement of any single cationic residue with Ala had a deleterious effect on heparin binding. The greatest reduction in heparin affinity (more than 30-fold) was observed with the R129deltaQ peptide (KD = 1.5 +/- 0.06 microM vs 51 +/- 2 nM for the reference peptide, K121-A134). Furthermore, each of the Ala-replacement peptides was a less-effective inhibitor of ATIII-heparin complex formation than the reference peptide. The poorest inhibitor was the R129deltaQ peptide which showed nearly 30% decrease in inhibition potency (60% inhibition at 100 microM peptide vs 90% inhibition with the reference peptide). The relative heparin affinities of the peptides measured by biological assay were the same as determined by titration calorimetry. Consequently, we modeled the complexes formed between the pentasaccharide unit structure and the R129deltaQ peptide or the reference peptide, K121-A134. In the "docked" complex, the assumed conformation of K121-A134 permitted juxtaposition of the cationic residues of the peptide with functional anionic groups of the pentasaccharide known to be involved in binding. A docked complex could also be formed between the R129deltaQ peptide and the pentasaccharide, but misalignment of critical peptide and saccharide functional groups was observed. The structure of the R129deltaQ-pentasaccharide complex was highly irregular because F123 and Y131 were completely surface exposed, likely yielding an unfavorable structure in aqueous solution. The observations from molecular modeling allow us to suggest that ATIII Geneva displays decreased heparin binding affinity due to its inability to form a productive binding complex in which essential electrostatic contacts are made between suitably juxtaposed saccharide anionic functional groups and cationic amino acid side chains.

Amino Acid Sequence↗

PubMed2004

Flexible docking of pyridinone derivatives into the non-nucleoside inhibitor binding site of HIV-1 reverse transcriptase.

Potent non-nucleoside reverse transcriptase inhibitors (NNRTIs) of the pyridinone derivative type were docked into nine NNRTIs binding pockets of HIV-1 reverse transcriptase (RT) structures. The docking results indicate that pyridinone analogues adopt a butterfly conformation and share the same binding mode as the crystal inhibitors in the pocket geometries of nevirapine, 1051U91, 9-Cl-TIBO, Cl-alpha-APA, efavirenz, UC-781, and S-1153. The results are in agreement with the data concerning mutational and structure-activity relationships available for pyridinone analogues and aid in the understanding, at the molecular level, of the biological response of published hybrid pyridinone molecules. Strategies to design further pyridinone derivatives active against RT containing mutations are discussed.

Anti-HIV Agents↗

PubMed2000

Active site binding modes of HIV-1 integrase inhibitors.

Using the crystal structure of the first complex of the HIV-1 integrase catalytic core domain with an inhibitor bound to the active site, structural models for the interaction of various inhibitors with integrase were generated by computational docking. For the compound of the crystallographic study, binding modes unaffected by crystal packing have recently been proposed. Although a large search region was used for the docking simulations, the ligands investigated here are found to bind preferably in similar ways close to the active site. The binding site is formed by residues 64-67, 116, 148, 151-152, 155-156, and 159, as well as by residue 92 in case of the largest ligand of the series. The coherent picture of possible interactions of small-molecule inhibitors at the active site provides an improved basis for structure-based ligand design. The recurring motif of tight interaction with the two lysine residues 156 and 159 is suggested to be of prime importance.

Anthraquinones↗

PubMed2005

Discovery of a novel class of reversible non-peptide caspase inhibitors via a structure-based approach.

In this paper, we report a simple structure-based iterative optimizations (SUBITO) strategy to identify and optimize new protein ligands and inhibitors. The approach is based on a combination of NMR-based screening and computational docking methods and enabled the identification of novel chemical leads among hundreds of thousands of commercially available compounds by screening only a few hundred compounds from a scaffold library followed by iterative screening steps where only few dozen compounds are tested. As an application, we report on the discovery of a novel class of non-peptide reversible caspase inhibitors, with IC(50) values in the low micromolar range.

Caspase 3↗

PubMed1998

Conformational changes in conantokin-G induced upon binding of calcium and magnesium as revealed by NMR structural analysis.

The apo- and metal-bound solution conformations of synthetic conantokin-G (con-G, G1Egamma gammaL5Q gamma NQgamma 10LIRgamma K15SN-CONH2, gamma = gamma-carboxyglutamic acid), an antagonist of N-methyl-D-aspartate receptor-derived neuronal ion channels, have been examined by one- and two-dimensional 1H NMR at neutral pH. A complete structure for the Mg2+-loaded peptide was defined by use of distance geometry calculations and was found to exist as an alpha-helix that spans the entire peptide. The alpha-helical nature of Mg2+/con-G was also supported by the small values (<5.5 Hz) of the 3JHNalpha coupling constants measured for amino acid residues 3-5, 8, 9, and 11-16, and the small values (<4 ppb/K) of the temperature coefficients observed for the alphaNH protons of residues 5-17. This conformation contrasted with that obtained for apo-con-G, which was nearly structureless in solution. Docking of Mg2+ into con-G was accomplished by use of the genetic algorithm/molecular dynamics simulation method, employing the NMR-derived Mg2+-loaded structure for initial coordinates in the midpoint calculations. For the 3 Mg2+/con-G model, it was found that binding of one Mg2+ ion is stabilized by oxygen atoms from three gamma-carboxylates of Gla3, Gla4, and Gla7; another Mg2+ is coordinated by two oxygen atoms, one from each of the gamma-carboxylates of Gla7; and a third metal ion through three donor oxygen atoms of gamma-carboxylates from Gla10 and Gla14. As shown from direct metal binding measurements to mutant con-G peptides, these latter two Gla residues probably stabilized the tightest binding Mg2+ ion. Circular dichroism studies of these same peptide variants demonstrated that all Gla residues contribute to the adoption of the Mg2+-dependent alpha-helical conformation in con-G. The data obtained in this investigation provide a molecular basis for the large conformational alteration observed in apo-con-G as a result of divalent cation binding and allow assessment of the roles of individual Gla residues in defining certain of the structure-function properties of con-G.

1-Carboxyglutamic Acid↗

PubMed2000

An unusual carbohydrate binding site revealed by the structures of two Maackia amurensis lectins complexed with sialic acid-containing oligosaccharides.

Seeds from the legume tree Maackia amurensis contain two lectins that can agglutinate different blood cell types. Their specificity toward sialylated oligosaccharides is unique among legume lectins; the leukoagglutinin preferentially binds to sialyllactosamine (alphaNeuAc(2-3)betaGal(1-4)betaGlcNAc), whereas the hemagglutinin displays higher affinity for a disialylated tetrasaccharide (alphaNeuAc(2-3)betaGal(1-3)[alphaNeuAc(2-6)]alphaG alNAc). The three-dimensional structure of the complex between M. amurensis leukoagglutinin and sialyllactose has been determined at 2.75-A resolution using x-ray crystallography. The carbohydrate binding site consists of a deep cleft that accommodates the three carbohydrate residues of the sialyllactose. The central galactose sits in the primary binding site in an orientation that has not been observed previously in other legume lectins. The carboxyl group of sialic acid establishes a salt bridge with a lysine side chain. The glucose residue is very efficiently docked between two tyrosine aromatic rings. The complex between M. amurensis hemagglutinin and a disialylated tetrasaccharide could be modeled from the leukoagglutinin/sialyllactose crystal structure. The substitution of one tyrosine by an alanine residue is responsible for the difference in fine specificity between the two isolectins. Comparison with other legume lectins indicates that oligosaccharide specificity within this family is achieved by the recycling of structural loops in different combinations.

Amino Acid Sequence↗

PubMed2005

Transient homodimer interactions studied using the electron self-exchange reaction.

Transient homodimer protein interactions have been investigated by analyzing the influence of ionic strength (NaCl) on the electron self-exchange (the bimolecular reaction whereby the two oxidation states of a redox protein interconvert) rate constant (k(ese)) of four plastocyanins. The k(ese) values for the plastocyanins from spinach, Dryopteris crassirhizoma (a fern), and the green alga Ulva pertusa, which possess acidic patches of varying size and locations, increase 190-, 29-, and 21-fold, respectively, at elevated ionic strength (I = 2.03 M). In contrast, the k(ese) for the almost neutral cyanobacterial plastocyanin from Anabaena variabilis exhibits very little dependence on ionic strength. The temperature dependence of the k(ese) for spinach plastocyanin (I = 0.28 M) provides evidence for poor packing at the homodimer interface. Representative structures of the transient homodimers involved in electron self-exchange, which are consistent with fits of the ionic strength dependence of k(ese) to van Leeuwen theory, have been obtained from protein modeling and docking simulations. The Coulombic energy of the docked homodimers follows the order spinach > D. crassirhizoma > U. pertusa > A. variabilis, which matches that of the overall influence of ionic strength on k(ese). Analysis of the homodimer structures indicates that poor packing and high planarity are features of the interface that favor transient interactions. The physiologically relevant Mg2+ ion has a much more pronounced influence on the k(ese) of spinach plastocyanin, which along with the known properties of the thylakoid lumen suggests a biological role for electron self-exchange.

Amino Acid Sequence↗

PubMed1997

Transfer RNA docking pair model in the ribosomal pre- and post-translocational states.

A consensus has been reached that the conformation of the anticodon-codon interactions of two adjacent tRNA molecules on the ribosome is a Sundaralingam-type (S-type). Even if it is kept to the S-type, there are still various possibilities. Various experimental data have been supporting an idea that the conformation of A-site tRNA is different from that of P-site tRNA. Those data as well as the recent result of Brimacombe and co-workers that U20:1 of lupin tRNAmMetbound to the A-site was cross-linked to a region, 875-905, of 23S rRNA in combination with the other recent findings of Nierhaus and co-workers about the spin-contrast method of neutron diffraction of the ribosome and the better accessible nucleotide patterns of phosphorothioated tRNAs on the ribosome have led to a new tRNA docking pair model, in which the highly conserved G18 and G19 of D-loop in A-site tRNA and C56 and C61 of TpsiC-loop in P-site tRNA base pair along with the conventional base pairs of adjacent codon-anticodon interactions. This A-P tRNA pair model can be translocated to the P-E tRNA pair model without changing the conformation except the ACCA termini, keeping the position of the growing nascent polypeptide chain.

Base Composition↗

PubMed2001

Can an optimization/scoring procedure in ligand-protein docking be employed to probe drug-resistant mutations in proteins?

A simple ligand-protein structural optimization and binding evaluation procedure has been routinely used in high-speed ligand-protein docking studies. In this work, we examine whether such an optimization/scoring procedure is useful in indicating possible drug-resistant mutations in proteins. Crystal structures of three wild-type enzymes (HIV-1 protease, HIV-1 reverse transcriptase, and Mycobacterium tuberculosis H37Rv enoyl-ACP reductase) complexed to a variety of inhibitors are studied. Mutations are introduced into these structures by using the molecular modeling software, SYBYL. Structural optimization and scoring of a mutant complex is conducted by a procedure similar to that used in a recent docking study (Wang et al., 1999). The computed results are compared with observed drug resistance data and the profile of nonresistant mutations. Most mutations studied show an energy change in the same direction as those indicated by observed resistance data. 50% of the polar to polar or nonpolar to nonpolar mutations are found to correlate qualitatively with observed drug resistance data. Van der Waals interactions account for most of these changes, which is in agreement with conclusions from structural studies. Substantially larger deviations are found between computed results and observed data for most polar to nonpolar or nonpolar to polar mutations, which result from deficiency in modelling and scoring ligand-protein interactions in our procedure. Our results suggest that an optimization/docking scoring procedure is useful for qualitatively probing polar to polar or nonpolar to nonpolar resistant mutations in addition to its application in screening active compounds. More accurate description of ligand-protein interactions and the use of methods such as free energy perturbation and Poisson-Boltzmann may be needed to further improve the quality of prediction.

Computer Simulation↗

PubMed1991

Protein-protein recognition analyzed by docking simulation.

Antibody-lysozyme and protease-inhibitor complexes are reconstituted by docking lysozyme as a rigid body onto the combining site of the antibodies and the inhibitors onto the active site of the proteases. Simplified protein models with one sphere per residue are subjected to simulated annealing using a crude energy function where the attractive component is proportional to the interface area. The procedure finds clusters of orientations in which a steric fit between the two protein components is achieved over a large contact surface. With five out of six complexes, the native structure of the complexes determined by X-ray crystallography is among those retained. Docked complexes are then subjected to conformational energy refinement with full atomic detail. With Fab HyHEL 5 and lysozyme, a native-like complex has the lowest refined energy. It can also be retrieved when starting with the X-ray structure of free lysozyme. However, some non-native complexes cannot be rejected: they form large interfaces, have a large number of H-bonds, and few unpaired polar groups. While these are necessary features of protein-protein recognition, they are not sufficient in determining specificity.

Algorithms↗

PubMed2005

Synthesis, structure-activity relationships and antagonistic effects in human MCF-7 breast cancer cells of flexible estrogen receptor modulators.

Estrogen receptors are therapeutic intervention targets for diseases such as osteoporosis and breast cancer with both tamoxifen and raloxifene established as clinical estrogen receptor antagonists. We report a series of novel selective estrogen receptor modulators (SERMs) whose structures are based on a flexible core scaffold differing from the triphenylethylene of tamoxifen analogues through the insertion of a benzylic methylene group as a flexible spacer between the aryl ring C and the ethylene group. A facile synthesis of the target compounds utilises the titanium tetrachloride/zinc mediated McMurry coupling reaction. Successive introduction onto the parent scaffold of hydroxyl functional groups afforded a series of increased potency ligands for the ER - essentially exploring the predicted in vivo metabolic activation of such aromatic SERM ligands. This second generation compound series demonstrated high antiproliferative potency against the MCF-7 human breast cancer cell line, with low cytotoxicity. High ER binding affinity (IC50 20 nM) together with up to 12 fold ERalpha/beta selectivity was also observed. In addition, the compounds displayed antiestrogenic effects at 40 nM when evaluated in the Ishikawa cell line with little estrogenic stimulation. Representative ligands were shown to be pro-apoptotic in human MCF-7 cells in a FACS based assay. Comparison of the docked structure obtained for the most active compound with the X-ray crystal structure reported for the complex of ERalpha and 4-hydroxytamoxifen, predict that these ligands bind in an antiestrogenic manner with some differences being observed in the benzylic Ring C orientation, as expected. This work further demonstrates the tolerance of the estrogen receptor towards flexible modulators.

Antineoplastic Agents↗

PubMed1997

High resolution fast quantitative docking using Fourier domain correlation techniques.

A 'docking' method based on finite grid forcefield sampling is proposed for fast evaluation of interaction energies between macromolecules and ligands. Forcefield used to calculate interaction energies utilizes a potential energy function composed of a 1/r-dependent electrostatic term and a (6-12) Lennard-Jones term for van der Waals interactions. Fast evaluation makes use of the convolution theorem allowing a point-by-point N-dimensional correlation in direct space to be replaced by a simple multiplication in spatial frequency space. Predictive accuracy was assessed by using seven protein-ligand complexes available from the Brookhaven Data Bank and determined crystallographically to high resolution. Successful prediction of ligand position and determination of ligand-protein interaction enthalpy was dependent on forecefield sampling grid size. Minimum interaction enthalpy calculated for four protein-ligand complexes coincided with crystallographic structures that used sampling grid sizes of 0.25 A resolution and was independent of ligand starting position and orientation. Successful docking was obtained for the remaining complexes at same grid resolution provided ligand starting positions were not randomized. Sensitivity of the docking algorithm to starting orientation was a consequence of tight fit of respective ligand structures with their protein target sites for these three cases and can be circumvented by using finer rotational sampling grids for the ligand. Boltzmann statistics derived from calculated interaction energies successfully extracted the observed ribonuclease A cytidylic acid complex from a manifold of similar interaction energies. The proposed method was able to reproduce the observed crystallographic complex by using a dynamical description of ligand.

Algorithms↗

PubMed2004

Protein conformational transitions coupled to binding in molecular recognition of unstructured proteins: hierarchy of structural loss from all-atom Monte Carlo simulations of p27Kip1 unfolding-unbinding and structural determinants of the binding mechanism.

Conformational transitions coupled to binding are studied for the p27(Kip1) protein which undergoes a functional disorder-to-order folding transition during tertiary complex formation with the phosphorylated cyclin A-cyclin-dependent kinase 2 (Cdk2) binary complex. Temperature-induced Monte Carlo simulations of p27(Kip1) unfolding-unbinding carried out from the crystal structure of the tertiary complex have revealed a systematic trend in the hierarchy of structural loss for p27(Kip1) and a considerable difference in mobility of p27(Kip1) secondary structure elements. The most persistent interactions of p27(Kip1) at the intermolecular interface during unfolding-unbinding simulations are formed by beta-hairpin and beta-strand that on average maintain their structural integrity considerably longer than other p27(Kip1) elements. We have found that the ensemble of unfolded p27(Kip1) conformations is characterized by transitions between mostly unbound, collapsed conformations and entropically favorable p27(Kip1) conformations, which are weakly bound to the cyclin A side of the binary complex. The results of this study are consistent with the experimental evidence pointing to this region of the intermolecular interface as a potential initiation docking site during binding reaction and may reconcile conflicting experimental hypotheses on the recognition of substrate recruitment motifs.

Cell Cycle Proteins↗

PubMed1993

NMR and computer-aided modeling studies of the interactions between a cyclic hexapeptide and the two enantiomers of some Boc- and Fmoc-amino acids.

The cyclic hexapeptide cyclo [-Pro1-Gly2-Glu3 (OBzl)-Pro4-Phe5-Leu6-] (1) was modeled and synthesized to be used for chiral discrimination studies. Total correlated spectroscopy and nuclear Overhauser effect spectroscopy spectra of the cyclic hexapeptide 1 in CDCl3 showed the presence of three stereoisomers: two dominant stereoisomers 1a and 1b that exchanged chemically with each other, and a minor stereoisomer 1c (4%) that exchanged exclusively with the stereoisomer 1b. Of the two dominant stereoisomers, only 1a interacted specifically with t-butyloxycarbonyl (Boc-) and 9-flourenylmethyloxycarbonyl (Fmoc-) amino acids in CDCl3. The interaction site of 1a when complexing with the derivatized amino acids was the chain segment Phe5-Leu6. The Phe5 NH and Leu6 NH protons are contiguous and solvent exposed. Their nmr signals shifted strongly downfield with the addition of Boc- or Fmoc- amino acids to the peptide solution. Thus, both NH protons hydrogen bond to the amino acids, forming a two-point hydrogen-bonding complex. The peptide stereoisomer 1b did not interact specifically with the Boc- and Fmoc-amino acids because of the lack of two contiguous and solvent-exposed peptidic NH protons that seem to be needed for specific interactions of the cyclic hexapeptide 1 with the Boc- and Fmoc-amino acids. A clear difference in the interaction of 1a with D- and L-enantiomers of Boc- Trp and Fmoc-Trp was observed with nmr spectroscopy. Docking models and molecular mechanics calculations together with nmr observations showed that the NH proton of the indole ring of the Boc-L-Trp and the Fmoc-L-Trp hydrogen bonded to the Pro1 carbonyl group. In this three-point hydrogen-bonding complex, the indole ring becomes locked underneath the Leu residue. The nmr signals of all the Leu6 protons (except for Leu NH) shifted strongly upfield owing to the shielding effect of the indole aromatic ring currents. The indole NH of the D-enantiomer did not hydrogen bond to the Pro1 carbonyl group because the formation of such a three-point hydrogen-bonding complex was thermodynamically unfavorable.

Amino Acid Sequence↗

PubMed2006

Fragment library-based representation system for protein conformation.

The representation system for protein conformation has a crucial effect on the speed of various protein-related simulations, including ab initio protein structure prediction and protein-protein docking simulation. Usually, the finer a representation system, the longer is the computational time required to employ the representation system in simulations. On the other hand, very coarse lattice systems cannot be directly applied to the simulation problems with real proteins. We report a new, fragment library-based protein conformation representation system, prepared by clustering amino acid conformations from 154 proteins. This system was composed of 64 most representative fragments per each amino acid, and based on the unified residue approach in which two spheres per amino acid were used. It could represent the conformation of the 82 proteins in an independent test set with the mean and standard deviation RMSD of 1.01 and 0.09 A, respectively, based on the position of alpha carbons and the centers of mass of sidechains.

Models, Molecular↗

PubMed2005

ZDOCK and RDOCK performance in CAPRI rounds 3, 4, and 5.

We present an evaluation of the results of our ZDOCK and RDOCK algorithms in Rounds 3, 4, and 5 of the protein docking challenge CAPRI. ZDOCK is a Fast Fourier Transform (FFT)-based, initial-stage rigid-body docking algorithm, and RDOCK is an energy minimization algorithm for refining and reranking ZDOCK results. Of the 9 targets for which we submitted predictions, we attained at least acceptable accuracy for 7, at least medium accuracy for 6, and high accuracy for 3. These results are evidence that ZDOCK in combination with RDOCK is capable of making accurate predictions on a diverse set of protein complexes.

Algorithms↗