Showing papers on "Docking (molecular) published in 2000"

A conserved docking motif in MAP kinases common to substrates, activators and regulators

Abstract: Mitogen-activated protein kinases (MAPKs) are specifically phosphorylated and activated by the MAPK kinases, phosphorylate various targets such as MAPK-activated protein kinases and transcription factors, and are inactivated by specific phosphatases. Recently, docking interactions via the non-catalytic regions of MAPKs have been suggested to be important in regulating these reactions. Here we identify docking sites in MAPKs and in MAPK-interacting enzymes. A docking domain in extracellular-signal-regulated kinase (ERK), a MAPK, serves as a common site for binding to the MAPK kinase MEK1, the MAPK-activated protein kinase MNK1 and the MAPK phosphatase MKP3. Two aspartic acids in this domain are essential for docking, one of which is mutated in the sevenmaker mutant of Drosophila ERK/Rolled. A corresponding domain in the MAPKs p38 and JNK/SAPK also serves as a common docking site for their MEKs, MAPK-activated protein kinases and MKPs. These docking interactions increase the efficiency of the enzymatic reactions. These findings reveal a hitherto unidentified docking motif in MAPKs that is used in common for recognition of their activators, substrates and regulators.

Protein-based virtual screening of chemical databases. 1. Evaluation of different docking/scoring combinations.

Abstract: Three different database docking programs (Dock, FlexX, Gold) have been used in combination with seven scoring functions (Chemscore, Dock, FlexX, Fresno, Gold, Pmf, Score) to assess the accuracy of virtual screening methods against two protein targets (thymidine kinase, estrogen receptor) of known three-dimensional structure. For both targets, it was generally possible to discriminate about 7 out of 10 true hits from a random database of 990 ligands. The use of consensus lists common to two or three scoring functions clearly enhances hit rates among the top 5% scorers from 10% (single scoring) to 25−40% (double scoring) and up to 65−70% (triple scoring). However, in all tested cases, no clear relationships could be found between docking and ranking accuracies. Moreover, predicting the absolute binding free energy of true hits was not possible whatever docking accuracy was achieved and scoring function used. As the best docking/consensus scoring combination varies with the selected target and the physicoch...

BiGGER: a new (soft) docking algorithm for predicting protein interactions.

Abstract: A new computationally efficient and automated “soft docking” algorithm is described to assist the prediction of the mode of binding between two proteins, using the three-dimensional structures of the unbound molecules. The method is implemented in a software package called BiGGER (Bimolecular Complex Generation with Global Evaluation and Ranking) and works in two sequential steps: first, the complete 6-dimensional binding spaces of both molecules is systematically searched. A population of candidate protein-protein docked geometries is thus generated and selected on the basis of the geometric complementarity and amino acid pairwise affinities between the two molecular surfaces. Most of the conformational changes observed during protein association are treated in an implicit way and test results are equally satisfactory, regardless of starting from the bound or the unbound forms of known structures of the interacting proteins. In contrast to other methods, the entire molecular surfaces are searched during the simulation, using absolutely no additional information regarding the binding sites. In a second step, an interaction scoring function is used to rank the putative docked structures. The function incorporates interaction terms that are thought to be relevant to the stabilization of protein complexes. These include: geometric complementarity of the surfaces, explicit electrostatic interactions, desolvation energy, and pairwise propensities of the amino acid side chains to contact across the molecular interface. The relative functional contribution of each of these interaction terms to the global scoring function has been empirically adjusted through a neural network optimizer using a learning set of 25 protein-protein complexes of known crystallographic structures. In 22 out of 25 protein-protein complexes tested, near-native docked geometries were found with Cα RMS deviations ≤ 4.0 A from the experimental structures, of which 14 were found within the 20 top ranking solutions. The program works on widely available personal computers and takes 2 to 8 hours of CPU time to run any of the docking tests herein presented. Finally, the value and limitations of the method for the study of macromolecular interactions, not yet revealed by experimental techniques, are discussed. Proteins 2000;39:372–384. © 2000 Wiley-Liss, Inc.

Structure of a ligand-binding intermediate in wild-type carbonmonoxy myoglobin.

Abstract: Small molecules such as NO, O2, CO or H2 are important biological ligands that bind to metalloproteins to function crucially in processes such as signal transduction, respiration and catalysis. A key issue for understanding the regulation of reaction mechanisms in these systems is whether ligands gain access to the binding sites through specific channels and docking sites, or by random diffusion through the protein matrix. A model system for studying this issue is myoglobin, a simple haem protein. Myoglobin has been studied extensively by spectroscopy, crystallography, computation and theory1,2,3,4,5,6,7,8,9,10,11. It serves as an aid to oxygen diffusion but also binds carbon monoxide, a byproduct of endogenous haem catabolism. Molecular dynamics simulations3,4,5, random mutagenesis6 and flash photolysis studies7,8,9,10 indicate that ligand migration occurs through a limited number of pathways involving docking sites. Here we report the 1.4 A resolution crystal structure of a ligand-binding intermediate in carbonmonoxy myoglobin that may have far-reaching implications for understanding the dynamics of ligand binding and catalysis.

Scoring docked conformations generated by rigid-body protein-protein docking.

Abstract: Rigid-body methods, particularly Fourier correlation techniques, are very efficient for docking bound (co-crystallized) protein confor- mations using measures of surface complementar- ity as the target function. However, when docking unbound (separately crystallized) conformations, the method generally yields hundreds of false posi- tive structures with good scores but high root mean square deviations (RMSDs). This paper describes a two-step scoring algorithm that can discriminate near-native conformations (with less than 5 A RMSD) from other structures. The first step includes two rigid-body filters that use the desolvation free en- ergy and the electrostatic energy to select a manage- able number of conformations for further process- ing, but are unable to eliminate all false positives. Complete discrimination is achieved in the second step that minimizes the molecular mechanics en- ergy of the retained structures, and re-ranks them with a combined free-energy function which in- cludes electrostatic, solvation, and van der Waals energy terms. After minimization, the improved fit in near-native complex conformations provides the free-energy gap required for discrimination. The algorithm has been developed and tested using docking decoys, i.e., docked conformations gener- ated by Fourier correlation techniques. The decoy sets are available on the web for testing other discrimination procedures. Proteins 2000;40:525-537.

Structure-based design, synthesis, and biological evaluation of novel pyrrolyl aryl sulfones: HIV-1 non-nucleoside reverse transcriptase inhibitors active at nanomolar concentrations.

Abstract: Pyrrolyl aryl sulfones (PASs) have been recently reported as a new class of human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT) inhibitors acting at the non-nucleoside binding site of this enzyme (Artico, M.; et al. J. Med. Chem. 1996, 39, 522−530). Compound 3, the most potent inhibitor within the series (EC50 = 0.14 μM, IC50 = 0.4 μM, and SI > 1429), was then selected as a lead compound for a synthetic project based on molecular modeling studies. Using the three-dimensional structure of RT cocrystallized with the α-APA derivative R95845, we derived a model of the RT/3 complex by taking into account previously developed structure−activity relationships. Inspection of this model and docking calculations on virtual compounds prompted the design of novel PAS derivatives and related analogues. Our computational approach proved to be effective in making qualitative predictions, that is in discriminating active versus inactive compounds. Among the compounds synthesized and tested, 20 was the ...

Analysis of Binding Affinities for Celecoxib Analogues with COX-1 and COX-2 from Combined Docking and Monte Carlo Simulations and Insight into the COX-2/COX-1 Selectivity

Abstract: The origins of binding affinity and COX-2/COX-1 selectivity for analogues of celecoxib have been explored using an approach that combines docking with Monte Carlo (MC) simulations. These inhibitors are COX-2-selective nonsteroidal antiinflammatory drugs (NSAIDs) that are of current interest because the gastrointestinal irritation they cause is reduced compared to that caused by traditional NSAIDs. We report a novel docking method, based on a combined Tabu and Monte Carlo protocol, that determines starting conformations for MC simulations. Using the docking-predicted starting conformations, relative changes in binding free energies were computed for methyl, ethyl, hydroxymethyl, hydroxyl, thiomethyl, methoxy, trifluoromethyl, chloro, fluoro, and unsubstituted derivatives with the MC free energy perturbation (FEP) method. The computed free energies are in good accord with IC50 values, and the structural information from the simulations can be used to explain the experimentally observed binding trends. In ad...

Identification of ligands for RNA targets via structure-based virtual screening: HIV-1 TAR

Abstract: Binding of the Tat protein to TAR RNA is necessary for viral replication of HIV-1. We screened the Available Chemicals Directory (ACD) to identify ligands to bind to a TAR RNA structure using a four-step docking procedure: rigid docking first, followed by three steps of flexible docking using a pseudobrownian Monte Carlo minimization in torsion angle space with progressively more detailed conformational sampling on a progressively smaller list of top-ranking compounds. To validate the procedure, we successfully docked ligands for five RNA complexes of known structure. For ranking ligands according to binding avidity, an empirical binding free energy function was developed which accounts, in particular, for solvation, isomerization free energy, and changes in conformational entropy. System-specific parameters for the function were derived on a training set of RNA/ligand complexes with known structure and affinity. To validate the free energy function, we screened the entire ACD for ligands for an RNA aptamer which binds l-arginine tightly. The native ligand ranked 17 out of ca. 153,000 compounds screened, i.e., the procedure is able to filter out >99.98% of the database and still retain the native ligand. Screening of the ACD for TAR ligands yielded a high rank for all known TAR ligands contained in the ACD and suggested several other potential TAR ligands. Eight of the highest ranking compounds not previously known to be ligands were assayed for inhibition of the Tat-TAR interaction, and two exhibited a CD50 of ca. 1 μM.

The structural basis for molecular recognition by the vitamin B 12 RNA aptamer

Abstract: Previous solution structures of ligand-binding RNA aptamers have shown that molecular recognition is achieved by the folding of an initially unstructured RNA around its cognate ligand, coupling the processes of RNA folding and binding. The 3 A crystal structure of the cyanocobalamin (vitamin B12) aptamer reported here suggests a different approach to molecular recognition in which elements of RNA secondary structure combine to create a solvent-accessible docking surface for a large, complex ligand. Central to this structure is a locally folding RNA triplex, stabilized by a novel three-stranded zipper. Perpendicular stacking of a duplex on this triplex creates a cleft that functions as the vitamin B12 binding site. Complementary packing of hydrophobic surfaces, direct hydrogen bonding and dipolar interactions between the ligand and the RNA appear to contribute to binding. The nature of the interactions that stabilize complex formation and the possible uncoupling of folding and binding for this RNA suggest a strong mechanistic similarity to typical protein-ligand complexes.

Active site binding modes of HIV-1 integrase inhibitors.

Abstract: 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.

Molecular Docking Reveals a Novel Binding Site Model for Fentanyl at the μ-Opioid Receptor

Abstract: The ligand binding modes of a series of fentanyl derivatives are examined using a combination of conformational analysis and molecular docking to the μ-opioid receptor. Condensed-phase molecular dynamics simulations are applied to evaluate potential relationships between ligand conformation and fentanyl substitution and to generate probable “bioactive” structures for the ligand series. Automated docking of the largely populated solution conformers identified a common binding site orientation that places the N-phenethyl group of fentanyl deep in a crevice between transmembrane (TM) helices II and III while the N-phenylpropanamide group projected toward a pocket formed by TM-III, -VI, and -VII domains. An analysis of the binding modes indicates the most potent fentanyl derivatives adopt an extended conformation both in solution and in the bound state, suggesting binding affinity may depend on the conformational preferences of the ligands. The results are consistent with ligand binding data derived from chim...

Similarity-driven flexible ligand docking.

Abstract: A similarity-driven approach to flexible ligand docking is presented. Given a reference ligand or a pharmacophore positioned in the protein active site, the method allows inclusion of a similarity term during docking. Two different algorithms have been implemented, namely, a similarity-penalized docking (SP-DOCK) and a similarity-guided docking (SG-DOCK). The basic idea is to maximally exploit the structural information about the ligand binding mode present in cases where ligand-bound protein structures are available, information that is usually ignored in standard docking procedures. SP-DOCK and SG-DOCK have been derived as modified versions of the program DOCK 4.0, where the similarity program MIMIC acts as a module for the calculation of similarity indices that correct docking energy scores at certain steps of the calculation. SP-DOCK applies similarity corrections to the set of ligand orientations at the end of the ligand incremental construction process, penalizing the docking energy and, thus, having only an effect on the relative ordering of the final solutions. SG-DOCK applies similarity corrections throughout the entire ligand incremental construction process, thus affecting not only the relative ordering of solutions but also actively guiding the ligand docking. The performance of SP-DOCK and SG-DOCK for binding mode assessment and molecular database screening is discussed. When applied to a set of 32 thrombin ligands for which crystal structures are available, SG-DOCK improves the average RMSD by ca. 1 A when compared with DOCK. When those 32 thrombin ligands are included into a set of 1,000 diverse molecules from the ACD, DIV, and WDI databases, SP-DOCK significantly improves the retrieval of thrombin ligands within the first 10% of each of the three databases with respect to DOCK, with minimal additional computational cost. In all cases, comparison of SP-DOCK and SG-DOCK results with those obtained by DOCK and MIMIC is performed. Proteins 2000;40:623–636. © 2000 Wiley-Liss, Inc.

Molecular imprinting of biotin derivatives and its application to competitive binding assay using nonisotopic labeled ligands.

Abstract: Synthetic biotin-binding polymers were prepared by molecular imprinting. Methacrylic acid (MAA) was copolymerized with ethylene glycol dimethacrylate in the presence of biotin methyl ester (B-Me) in chloroform. Hydrogen-bonding-based complexation of B-Me with MAA generates the binding sites complementary to B-Me after extracting B-Me from the resulting copolymers. Data from NMR titration suggest a one-to-one prepolymerization complex formation of B-Me with MAA in chloroform. A possible complex structure was estimated by docking of the most stable conformers by intermolecular Monte Carlo conformational search under the assumption of a one-to-one association. The selectivity of the imprinted polymers was investigated and an imprinted polymer-based competitive binding assay for B-Me was demonstrated using biotin p-nitrophenyl ester as a nonisotopic-labeled ligand.

A method for including protein flexibility in protein-ligand docking: improving tools for database mining and virtual screening.

Abstract: Second-generation methods for docking ligands into their biological receptors, such as FLOG, provide for flexibility of the ligand but not of the receptor Molecular dynamics based methods, such as free energy perturbation, account for flexibility, solvent effects, etc, but are very time consuming We combined the use of statistical analysis of conformational samples from short-run protein molecular dynamics with grid-based docking protocols and demonstrated improved performance in two test cases Our statistical analysis explores the importance of the average strength of a potential interaction with the biological target and optionally applies a weighting depending on the variability in the strength of the interaction seen during dynamics simulation Using these methods, we improved the number of known dihydrofolate reductase ligands found in the top-ranked 10% of a database of drug-like molecules, in searches based on the three-dimensional structure of the protein These methods are able to match the ability of manual docking to assess likely inactivity on steric grounds and indeed to rank order ligands from a homologous series of cyclooxygenase-2 inhibitors with good correlation to their true activity Furthermore, these methods reduce the need for human intervention in setting up molecular docking experiments

Rational design of femtomolar inhibitors of isoleucyl tRNA synthetase from a binding model for pseudomonic acid-A.

Abstract: This paper describes the design and characterization of novel inhibitors of IleRS, whose binding affinity approaches the tightest reported for noncovalent inhibition. Compounds were designed from a binding model for the natural product pseudomonic acid-A (PS-A) together with a detailed understanding of the reaction cycle of IleRS and characterization of the mode of binding of the reaction intermediate IleAMP. The interactions of the compounds with IleRS were characterized by inhibition of aminoacylation of tRNA or PP(i)/ATP exchange at supersaturating substrate concentration and by transient kinetics and calorimetry methods. A detailed understanding of the interaction of a comprehensive series of compounds with IleRS allowed the identification of key features and hence the design of exquisitely potent inhibitors. Predictions based on these results have been recently supported by a docking model based on the crystal structure of IleRS with PS-A [Silvian, L. F., Wang J. M., and Steitz T. A. (1999) Science 285 1074-1077].

New approach to molecular docking and its application to virtual screening of chemical databases

Abstract: This paper describes the validation of a molecular docking method and its application to virtual database screening. The code flexibly docks ligand molecules into rigid receptor structures using a tabu search methodology driven by an empirically derived function for estimating the binding affinity of a protein−ligand complex. The docking method has been tested on 70 ligand−receptor complexes for which the experimental binding affinity and binding geometry are known. The lowest energy geometry produced by the docking protocol is within 2.0 A root mean square of the experimental binding mode for 79% of the complexes. The method has been applied to the problem of virtual database screening to identify known ligands for thrombin, factor Xa, and the estrogen receptor. A database of 10 000 randomly chosen “druglike” molecules has been docked into the three receptor structures. In each case known receptor ligands were included in the study. The results showed good separation between the predicted binding affinit...

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

Abstract: 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 β-turn site. The ligand 1α,25-dihydroxyvitamin D3 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 β-turn (site 1). Three residues forming hydrogen bonds with the functionally important 1α- and 25-hydroxyl groups of 1α,25-dihydroxyvitamin D3 were identified and confirmed by mutational analysis: the 1α-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.

Adenosine Deaminase Prefers a Distinct Sugar Ring Conformation for Binding and Catalysis: Kinetic and Structural Studies

Abstract: Several recent X-ray crystal structures of adenosine deaminase (ADA) in complex with various adenosine surrogates have illustrated the preferred mode of substrate binding for this enzyme. To define more specific structural details of substrate preferences for binding and catalysis, we have studied the ADA binding efficiencies and deamination kinetics of several synthetic adenosine analogues in which the furanosyl ring is biased toward a particular conformation. NMR solution studies and pseudorotational analyses were used to ascertain the preferred furanose ring puckers (P, nu(MAX)) and rotamer distributions (chi and gamma) of the nucleoside analogues. It was shown that derivatives which are biased toward a "Northern" (3'-endo, N) sugar ring pucker were deaminated up to 65-fold faster and bound more tightly to the enzyme than those that preferred a "Southern" (2'-endo, S) conformation. This behavior, however, could be modulated by other structural factors. Similarly, purine riboside inhibitors of ADA that prefer the N hemisphere were more potent inhibitors than S analogues. These binding propensities were corroborated by detailed molecular modeling studies. Docking of both N- and S-type analogues into the ADA crystal structure coordinates showed that N-type substrates formed a stable complex with ADA, whereas for S-type substrates, it was necessary for the sugar pucker to adjust to a 3'-endo (N-type) conformation to remain in the ADA substrate binding site. These data outline the intricate structural details for optimum binding in the catalytic cleft of ADA.

DoMCoSAR: a novel approach for establishing the docking mode that is consistent with the structure-activity relationship. Application to HIV-1 protease inhibitors and VEGF receptor tyrosine kinase inhibitors.

Abstract: DoMCoSAR is a novel approach for statistically determining the docking mode that is consistent with a structure−activity relationship. The approach establishes the binding mode for the compounds in a chemical series with the assumption that all molecules exhibit the same binding mode. It involves three stages. In the first stage all molecules that belong to a given chemical series are docked to the active site of the protein target. The only bias used in the docking at this stage involves the location of the protein binding site. Coordinates of the common substructure (CS) that results from the unbiased docking are then clustered to establish the major substructure docking modes. In the second stage all molecules are docked to the major docking modes (MDMs) with constraints based on the common substructure. The third stage generates, for the major docking modes, interaction-based descriptors that include electrostatic, VDW, strain, and solvation contributions. The problem of docking mode evaluation is now...

A new model for 20-hydroxyecdysone and dibenzoylhydrazine binding: a homology modeling and docking approach.

Abstract: The ecdysone receptor (ECR), a nuclear transcription factor controlling insect development, is a novel target for insecticides such as dibenzoylhydrazines with low environmental and toxicological impacts. To understand the high selectivity of such synthetic molecules toward ECR, two homology models of the Chironomus tentans ECR ligand-binding domain (LDB) have been constructed by taking as templates the known LBD crystal structures of the retinoic acid and vitamin D receptors. Docking of 20-hydroxyecdysone (20E) and dibenzoylhydrazines to the receptor suggests a novel superposition of the natural and synthetic molecules; the N-tert-butyl substituent of the dibenzoylhydrazines extends significantly beyond the 20E volume. Our ECR–LBD protein models rationalize how 20E and dibenzoylhy-drazines interact with the ligand-binding pocket. The homology model complexes provide new insights that can be exploited in the rational design of new environmentally safe insecticides.

Docking-Based Development of Purine-like Inhibitors of Cyclin-Dependent Kinase-2

Abstract: The cell division cycle is controlled by cyclin-dependent kinases (cdk), which consist of a catalytic subunit (cdk1-cdk8) and a regulatory subunit (cyclin A-H). Purine-like inhibitors of cyclin-dependent kinases have recently been found to be of potential use as anticancer drugs. Rigid and flexible docking techniques were used for analysis of binding mode and design of new inhibitors. X-ray structures of three (ATP, olomoucine, roscovitine) cdk2 complexes were available at the beginning of the study and were used to optimize the docking parameters. The new potential inhibitors were then docked into the cdk2 enzyme, and the enzyme/inhibitor interaction energies were calculated and tested against the assayed activities of cdk1 (37 compounds) and cdk2 (9 compounds). A significant rank correlation between the activity and the rigid docking interaction energy has been found. This implies that (i) the rigid docking can be used as a tool for qualitative prediction of activity and (ii) values obtained by the rigid docking technique into the cdk2 active site can also be used for the prediction of cdk1 activity. While the resulting geometries obtained by the rigid docking are in good agreement with the X-ray data, the flexible docking did not always produce the same inhibitor conformation as that found in the crystal.

Modeling of the inhibition of retroviral integrases by styrylquinoline derivatives.

Abstract: Styrylquinoline derivatives, known to be potent inhibitors of HIV-1 integrase, have been experimentally tested for their inhibitory effect on the disintegration reaction catalyzed by catalytic cores of HIV-1 and Rous sarcoma virus (RSV) integrases. A modified docking protocol, consisting of coupling a grid search method with full energy minimization, has been specially designed to study the interaction between the inhibitors and the integrases. The inhibitors consist of two moieties that have hydroxyl and/or carboxyl substituents: the first moiety is either benzene, phenol, catechol, resorcinol, or salicycilic acid; the hydroxyl substituents on the second (quinoline) moiety may be in the keto or in the enol forms. Several tautomeric forms of the drugs have been docked to the crystallographic structure of the RSV catalytic core. The computed binding energy of the keto forms correlates best with the measured inhibitory effect. The docking procedure shows that the inhibitors bind closely to the crystallogra...

Evaluation of docking/scoring approaches: a comparative study based on MMP3 inhibitors.

Abstract: An increasing number of docking/scoring programs are available that use different sampling and scoring algorithms. A reliable scoring function is the crucial element of such approaches. Comparative studies are needed to evaluate their current capabilities. DOCK4 with force field and PMF scoring as well as FlexX were used to evaluate the predictive power of these docking/scoring approaches to identify the correct binding mode of 61 MMP-3 inhibitors in a crystal structure of stromelysin and also to rank them according to their different binding affinities. It was found that DOCK4/PMF scoring performs significantly better than FlexX and DOCK4/FF in both ranking ligands and predicting their binding modes. Most notably, DOCK4/PMF was the only scoring/docking approach that found a significant correlation between binding affinity and predicted score of the docked inhibitors. However, comparing only those cases where the correct binding mode was identified (scoring highest among sampled poses), FlexX showed the best `fine tuning' (lowest rmsd) in predicted binding modes. The results suggest that not so much the sampling procedure but rather the scoring function is the crucial element of a docking program.

Heteronuclear NMR and soft docking: an experimental approach for a structural model of the cytochrome c553-ferredoxin complex.

Abstract: The combination of docking algorithms with NMR data has been developed extensively for the studies of protein−ligand interactions. However, to extend this development for the studies of protein−protein interactions, the intermolecular NOE constraints, which are needed, are more difficult to access. In the present work, we describe a new approach that combines an ab initio docking calculation and the mapping of an interaction site using chemical shift variation analysis. The cytochrome c553−ferredoxin complex is used as a model of numerous electron-transfer complexes. The 15N-labeling of both molecules has been obtained, and the mapping of the interacting site on each partner, respectively, has been done using HSQC experiments. 1H and 15N chemical shift analysis defines the area of both molecules involved in the recognition interface. Models of the complex were generated by an ab initio docking software, the BiGGER program (bimolecular complex generation with global evaluation and ranking). This program ge...

Application of a Molecular Dynamics Simulation Method with a Generalized Effective Potential to the Flexible Molecular Docking Problems

Abstract: We present a new molecular dynamics method using Tsallis effective potential for the flexible docking problems of streptavidin/biotin and protein kinase C/phorbol-13-acetate. With a full flexibility of the ligands and a partial flexibility of the receptor active sites included, the new MD scheme accelerates the docking process significantly, by way of infrequent q-jumping and q-relaxation procedures between a normal potential energy surface and its transformed one by the Tsallis scheme. In the transformation, only the nonbonding interaction terms of an empirical potential energy function were employed. It has been found that the current method can predict the correctly docked structures in quite an effective way. Current results strongly indicate that this new MD method can be a very promising strategy for flexible ligand/flexible receptor docking.