PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “molecular docking”

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 829 records · Page 46Linked to original sources

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↗

Molecular interactions of 3',5'-cyclic purine analogues with the binding site of retinal rod ion channels.

Photoreceptor outer segments transduce information about incoming light levels through a class of ion channels that respond directly to changes in cytosolic 3',5'-cyclic guanosine monophosphate levels. A series of 3',5'-cyclic purine analogues with alterations at N1, C2, C6, or C8 positions was used to examine molecular interactions between the nucleotide and the channel. The maximal current activated by C2-altered analogues in excised membrane patches was less than the current activated by cGMP, and the K0.5, the concentration which activates 50% of the current in a patch, was increased. Nonpolar C8-substituted cAMP analogues activated more current than the parent cAMP with lower K0.5 values. This was in contrast to 8-amino-cAMP, which exhibited greatly reduced activity. The rank order of activity, based on K0.5 values, for C8-cAMP substituents was as follows: 8-azido- > 8-methylamino- > 8-benzylamino- > cAMP > 8-bromo- > 8-hydroxy- >> 8-amino-cAMP. 1,N6-Etheno-cAMP and N6-monobutyryl-cAMP activated a small fraction of the total possible current with high K0.5 values. Other analogues with alterations at N1 or C6 positions including N1-oxide-cAMP, 2-aminopurine riboside 3',5'-monophosphate, and N6-monosuccinyl-cAMP do not bind to the channel, suggesting that interactions with the channel in this region are essential for binding. In order to help interpret the changes in maximal current and K0.5 values compared to cGMP, molecular models of the active analogues were constructed and then docked into a molecular model of the cyclic nucleotide binding site of the retinal channel. This model, proposed by Kumar and Weber [(1992) Biochemistry 31, 4643-4649], was based on the crystal structure of cAMP bound to catabolite activator protein. Our modeling showed that the analogues were sterically accommodated within the binding site. No hydrogen bonds were predicted between the purine rings of cAMP and the pocket; however, Phe 533 on the beta 5 strand was predicted to form weak electrostatic interactions with C6 substituents on both cAMP and cGMP. The importance of contacts in this region of the binding pocket is further emphasized by the inactive analogues, all of which are altered at N1 or C6.

Ambystoma↗

Metal-dependent inhibition of HIV-1 integrase by beta-diketo acids and resistance of the soluble double-mutant (F185K/C280S).

The beta-diketo acids (DKAs) represent a major advance for anti-HIV-1 integrase drug development. We compared the inhibition of HIV-1 integrase by six DKA derivatives using the wild-type enzyme or the double-mutant F185K/C280S, which has been previously used for crystal structure determinations. With the wild-type enzyme, we found that DKAs could be classified into two groups: those similarly potent in the presence of magnesium and manganese and those potent in manganese and relatively ineffective in the presence of magnesium. Both the aromatic and the carboxylic or tetrazole functions of DKAs determined their metal selectivity. The F185K/C280S enzyme was markedly more active in the presence of manganese than magnesium. The F185K/C280S integrase was also relatively resistant to the same group of DKAs that were potent in the presence of magnesium with the wild-type enzyme. Resistance was caused by a synergistic effect from both the F185K and C280S mutations. Molecular modeling and docking suggested metal-dependent differences for binding of DKAs. Molecular modeling also indicated that the tetrazole or the azido groups of some derivatives could directly chelate magnesium or manganese in the integrase catalytic site. Together, these experiments suggest that DKAs recognize conformational differences between wild-type and the double-mutant HIV-1 integrase, because they chelate the magnesium or manganese in the enzyme active site and compete for DNA binding.

Acetoacetates↗

Endothelial caveolae have the molecular transport machinery for vesicle budding, docking, and fusion including VAMP, NSF, SNAP, annexins, and GTPases.

Transport by discrete vesicular carriers is well established at least in part because of recent discoveries identifying key protein mediators of vesicle formation, docking, and fusion. A general mechanism sensitive to N-ethylmaleimide (NEM) is required for the transport of a divergent group of vesicular carriers in all eukaryotes. Many endothelia have an abundant population of non-coated plasmalemmal vesicles or caveolae, which have been reported with considerable controversy to function in transport. We recently have shown that like other vesicular transport systems, caveolae-mediated endocytosis and transcytosis are inhibited by NEM (Schnitzer, J. E., Allard, J., and Oh, P. (1995) Am. J. Physiol. 268, H48-H55). Here, we continue this work by utilizing our recently developed method for purifying endothelial caveolae from rat lung tissue (Schnitzer, J. E., Oh, P., Jacobson, B. S., and Dvorak, A. M. (1995) Proc. Natl. Acad. Sci. U. S. A. 92, 1759-1763) to show that these caveolae contain key proteins known to mediate different aspects of vesicle formation, docking, and/or fusion including the vSNARE VAMP-2, monomeric and trimeric GTPases, annexins II and VI, and the NEM-sensitive fusion factor NSF along with its attachment protein SNAP. Like neuronal VAMPs, this endothelial VAMP is sensitive to cleavage by botulinum B and tetanus neurotoxins. Caveolae in endothelium are indeed like other carrier vesicles and contain similar NEM-sensitive molecular machinery for transport.

Animals↗

Predicting molecular interactions and inducible complementarity: fragment docking of Fab-peptide complexes.

Antibody-antigen interactions are representative of a broad class of receptor-ligand interactions involving both specificity and potential inducible complementarity. To test possible mechanisms of antigen-antibody recognition and specificity computationally, we have used a Metropolis Monte Carlo algorithm to dock fragments of the epitope Glu-Val-Val-Pro-His-Lys-Lys to the X-ray structures of both the free and the complexed Fab of the antibody B13I2 (raised against the C-helix of myohemerythrin). The fragments Pro-His and Val-Pro-His, which contain residues experimentally identified as important for binding, docked correctly to both structures, but all tetrapeptide and larger fragments docked correctly only to the complexed Fab, even when torsional flexibility was added to the ligand. However, only tetrapeptide and larger fragments showed significantly more favorable energies when docked to the complexed Fab coordinates than when docked to either the free Fab or a non-specific site remote from the combining site. Comparison of the free and complexed B13I2 structures revealed that atoms within 5 A of Val-Pro-His showed little movement upon peptide binding, but atoms within 5 A of the other four epitope residues showed greater movements. These results computationally distinguished recognition and binding processes with practical implications for drug design strategies. Overall, this new fragment docking approach establishes distinct roles for the "lock-and-key" (recognition) and the "handshake" (binding) paradigms in antibody-antigen interaction, suggests an incremental approach to incorporating flexibility in computational docking, and identifies critical regions within receptor binding sites for ligand recognition.

Algorithms↗

Putative hAPN receptor binding sites in SARS_CoV spike protein.

AIM: To obtain the information of ligand-receptor binding between the S protein of SARS-CoV and CD13, identify the possible interacting domains or motifs related to binding sites, and provide clues for studying the functions of SARS proteins and designing anti-SARS drugs and vaccines. METHODS: On the basis of comparative genomics, the homology search, phylogenetic analyses, and multi-sequence alignment were used to predict CD13 related interacting domains and binding sites in the S protein of SARS-CoV. Molecular modeling and docking simulation methods were employed to address the interaction feature between CD13 and S protein of SARS-CoV in validating the bioinformatics predictions. RESULTS: Possible binding sites in the SARS-CoV S protein to CD13 have been mapped out by using bioinformatics analysis tools. The binding for one protein-protein interaction pair (D757-R761 motif of the SARS-CoV S protein to P585-A653 domain of CD13) has been simulated by molecular modeling and docking simulation methods. CONCLUSION: CD13 may be a possible receptor of the SARS-CoV S protein, which may be associated with the SARS infection. This study also provides a possible strategy for mapping the possible binding receptors of the proteins in a genome.

Amino Acid Sequence↗

Docking by least-squares fitting of molecular surface patterns.

Molecular surfaces are fitted to each other by a new solution to the problem of docking a ligand into the active site of a protein molecule. The procedure constructs patterns of points on the surfaces and superimposes them upon each other using a least-squares best-fit algorithm. This brings the surfaces into contact and provides a direct measure of their local complementarity. The search over the ligand surface produces a large number of dockings, of which a small fraction having the best complementarity and the least steric hindrance are evaluated for electrostatic interaction energy. When applied to molecules taken from crystallographically observed complexes, this procedure consistently assigns the lowest electrostatic energies to correct dockings. On independently determined structures, the ability of the method to discern correct dockings depends on how much conformational difference there is between the free and complexed forms of the molecules. The procedure is found to be fast enough on contemporary workstation computers to permit many conformations to be considered, and tolerant enough to make rather coarse bond dihedral sampling a practicable way to overcome the problem of structural flexibility.

Binding Sites↗

A model of scorpion toxin binding to voltage-gated K+ channels.

Mutational studies have identified part of the S5-S6 loop of voltage-dependent K+ channels (P region) responsible for tetraethylammonium (TEA) block and permeation properties. Several scorpion peptide toxins-charybdotoxin (ChTX), kaliotoxin (KITX), and agitoxin (AgTX)-also block the channel with high affinity and specificity. Here, we examine the interaction predicted when the toxins are docked onto the molecular model of the K+ channel pore that we recently proposed. Docking with the model of the Kv1.3 channel started by location of Lys-27 side chain into the central axis of the pore, followed by energy minimization. In the optimal arrangement, Arg-24 of KITX or AgTX forms a hydrogen bond with the Asp-386 carboxyl of one subunit, and Asn-30 is in immediate contact with Asp-386 of the opposing subunit in the tetramer. Toxin residues in proximity to the side chain of Lys-27 (Phe-25, Thr-36, Met-29, and Ser-11 in KITX) interact with the four C-end His-404s. For ChTX the interaction with Asp-386 is reduced, but this is compensated by additional nonbonded interactions formed by Tyr-36 and Arg-34. Comparison of calculated energy of interaction of these specific toxin-channel residues with experimental studies reveals good agreement. The similar total calculated energy of interaction is consistent with the similar IC50 for Kv1.3 block by KITX and AgTX. Steric contacts of residues in position 380 of the S5-P linker with residues on the upper part of toxins permit reconstruction of the K+ channel outer vestibule walls, which are about 30 A apart and about 9 A high. Molecular modeling shows complementarity of the pore model to toxin spacial structures, and supports the proposal that the N-terminal borders of the P regions surround residues of their C-terminal halves.

Amino Acid Sequence↗

Molecular dynamics, database screening, density functional and docking studies of novel RAR ligands in cancer chemotherapy.

Due to the major challenge which cancer treatment and cure still imposes after many decades to the international scientific community, there is actually considerable interest in new ligands with increased bioactivity. We have focused on the retinoid acid receptor, which is considered an interesting target for drug design. In this work, we have carried out density functional geometry optimizations and different docking procedures. We have performed screening in a large database (hundreds of thousands of molecules which we optimized at the AM1 level) yielding a set of potential bioactive ligands. Two new ligands were selected and optimized at B3LYP/6-31G* level. A flexible docking program was used to investigate the interactions between the receptor and the new ligands. Molecular dynamics were performed in order to investigate the stability of the two ligands as well as the crystallographic RAR ligand inside the RAR active site. We also investigated the stability of all the main protein-ligand contacts. The parameters of the Rule of Five were investigated. The result of this work is compared with a crystallographic ligand of RAR. One of our theoretically bioactive new ligands indicates stronger and more polar and hydrophobic interactions with the receptor.

Algorithms↗

Molecular modelling of protein-carbohydrate interactions. Docking of monosaccharides in the binding site of concanavalin A.

A general procedure is described for addressing the computer simulation of protein-carbohydrate interactions. First, a molecular mechanical force field capable of performing conformational analysis of oligosaccharides has been derived by the addition of new parameters to the Tripos force field; it is also compatible with the simulation of protein. Second, a docking procedure which allows for a systematic exploration of the orientations and positions of a ligand into a protein cavity has been designed. This so-called 'crankshaft' method uses rotations and variations about/of virtual bonds connecting, via dummy atoms, the ligand to the protein binding site. Third, calculation of the relative stability of protein ligand complexes is performed. This strategy has been applied to search for all favourable interactions occurring between a lectin [concanavalin A (ConA)] and methyl alpha-D-mannopyranoside or methyl alpha-D-glucopyranoside. For each monosaccharide, different stable orientations and positions within the binding site can be distinguished. Among them, one corresponds to very favourable interactions, not only in terms of hydrogen bonding, but also in terms of van der Waals interactions. It corresponds precisely to the binding mode of methyl alpha-D-mannopyranoside into ConA as revealed by the 2.9 A resolution of the crystalline complex (Derewenda et al., 1989). Some implications of the present modelling study with respect to the molecular basis of the specificity of the interaction of lectins with various monosaccharides are presented.

Binding Sites↗

Automated docking of 82 N-benzylpiperidine derivatives to mouse acetylcholinesterase and comparative molecular field analysis with 'natural' alignment.

Automated docking and three-dimensional Quantitative Structure-Activity Relationship studies (3D QSAR) were performed for a series of 82 reversible, competitive and selective acetylcholinesterase (AChE) inhibitors. The suggested automated docking technique, making use of constraints taken from experimental crystallographic data, allowed to dock all the 82 substituted N-benzylpiperidines to the crystal structure of mouse AChE, because of short computational times. A 3D QSAR model was then established using the CoMFA method. In contrast to conventional CoMFA studies, the compounds were not fitted to a reference molecule but taken in their 'natural' alignment obtained by the docking study. The established and validated CoMFA model was then applied to another series of 29 N-benzylpiperidine derivatives whose AChE inhibitory activity data were measured under different experimental conditions. A good correlation between predicted and experimental activity data shows that the model can be extended to AChE inhibitory activity data measured on another acetylcholinesterase and/or at different incubation times and pH level.

Acetylcholinesterase↗

Protein-ligand docking: current status and future challenges.

Understanding the ruling principles whereby protein receptors recognize, interact, and associate with molecular substrates and inhibitors is of paramount importance in drug discovery efforts. Protein-ligand docking aims to predict and rank the structure(s) arising from the association between a given ligand and a target protein of known 3D structure. Despite the breathtaking advances in the field over the last decades and the widespread application of docking methods, several downsides still exist. In particular, protein flexibility-a critical aspect for a thorough understanding of the principles that guide ligand binding in proteins-is a major hurdle in current protein-ligand docking efforts that needs to be more efficiently accounted for. In this review the key concepts of protein-ligand docking methods are outlined, with major emphasis being given to the general strengths and weaknesses that presently characterize this methodology. Despite the size of the field, the principal types of search algorithms and scoring functions are reviewed and the most popular docking tools are briefly depicted. Recent advances that aim to address some of the traditional limitations associated with molecular docking are also described. A selection of hand-picked examples is used to illustrate these features.

Algorithms↗

Molecular insight into pseudolysin inhibition using the MM-PBSA and LIE methods.

Pseudolysin, the extracellullar elastase of Pseudomonas aeruginosa (EC: 3.4.24.26) plays an important role in the pathogenesis of P. aeruginosa infections. In the present study, molecular dynamics simulations and theoretical affinity predictions were used to gain molecular insight into pseudolysin inhibition. Four low molecular weight inhibitors were docked at their putative binding sites and molecular dynamics (MD) simulations were performed for 5.0 ns, and the free energy of binding was calculated by the linear interaction energy method. The number and the contact surface area of stabilizing hydrophobic, aromatic, and hydrogen bonding interactions appears to reflect the affinity differences between the inhibitors. The proteinaceous inhibitor, Streptomyces metalloproteinase inhibitor (SMPI) was docked in three different binding positions and MD simulations were performed for 3.0 ns. The MD trajectories were used for molecular mechanics-Poisson-Boltzmann surface area analysis of the three binding positions. Computational alanine scanning of the average pseudolysin-SMPI complexes after MD revealed residues at the pseudolysin-SMPI interface giving the main contribution to the free energy of binding. The calculations indicated that SMPI interacts with pseudolysin via the rigid active site loop, but that also contact sites outside this loop contribute significantly to the free energy of association.

Amino Acid Sequence↗

Oligomerization states of Bowman-Birk inhibitor by atomic force microscopy and computational approaches.

Several methods have been applied to study protein-protein interaction from structural and thermodynamic point of view. The present study reveals that atomic force microscopy (AFM), molecular modeling, and docking approaches represent alternative methods offering new strategy to investigate structural aspects in oligomerization process of proteinase inhibitors. The topography of the black-eyed pea trypsin/chymotrypsin inhibitor (BTCI) was recorded by AFM and compared with computational rigid-bodies docking approaches. Multimeric states of BTCI identified from AFM analysis showed globular-ellipsoidal shapes. Monomers, dimers, trimers, and hexamers were the most prominent molecular arrays observed in AFM images as evaluated by molecular volume calculations and corroborated by in silico docking and theoretical approaches. We therefore propose that BTCI adopts stable and well-packed self-assembled states in monomer-dimer-trimer-hexamer equilibrium. Although there are no correlation between specificity and packing efficiency among proteinases and proteinase inhibitors, the AFM and docked BTCI analyses suggest that these assemblies may exist in situ to play their potential function in oligomerization process.

Computational Biology↗

Targeting aurora2 kinase in oncogenesis: a structural bioinformatics approach to target validation and rational drug design.

The aurora kinases are a novel oncogenic family of mitotic serine/threonine kinases (S/T kinases) that are overexpressed in a number of solid tumors, including pancreas and colorectal cancer. A PSI-BLAST search [National Center for Biotechnology Information (NCBI)] with the sequence of the S/T kinase domain of human aurora1 kinase [also known as AUR1, ARK2, AIk2, AIM-1, and STK12] and human aurora2 kinase (also known as AUR2, ARK1, AIK, BTAK, and STK15) showed a high sequence similarity to the three-dimensional structures of bovine cAMP-dependent kinase [Brookhaven Protein Data Bank code 1CDK], murine cAMP-dependent kinase (1APM), and Caenorhabditis elegans twitchin kinase (1KOA). When the aurora1 or aurora2 sequence was input into the tertiary structure prediction programs THREADER and 3D-PSSM (three-dimensional position-sensitive scoring matrix), the top structural matches were 1CDK, 1APM, and 1KOA, confirming that these domains are structurally conserved. The structural models of aurora1 and aurora2 were built using 1CDK as the template structure. Molecular dynamics and docking simulations, targeting the ATP binding site of aurora2 with adenylyl imidodiphosphate (AMP-PNP), staurosporine, and six small molecular S/T kinase inhibitors, identified active-site residues that interact with these inhibitors differentially. The docked structures of the aurora2-AMP-PNP and aurora2-staurosporine complexes indicated that the adenine ring of AMP-PNP and the indolocarbazole moiety of staurosporine have similar positions and orientations and provided the basis for the docking of the other S/T kinase inhibitors. Inhibitors with isoquinoline and quinazoline moieties were recognized by aurora2 in which H-89 and 6,7-dimethoxyquinazoline compounds exhibited high binding energies compared with that of staurosporine. The calculated binding energies for the docked small-molecule inhibitors were qualitatively consistent with the IC(50) values generated using an in vitro kinase assay. The aurora2 structural model provides a rational basis for site-directed mutagenesis of the active site; design of novel H-89, staurosporine, and quinazoline analogues; and the screening of the available chemical database for the identification of other novel, small-molecular entities.

Amino Acid Sequence↗

Surrogate docking: structure-based virtual screening at high throughput speed.

Structure-based screening using fully flexible docking is still too slow for large molecular libraries. High quality docking of a million molecule library can take days even on a cluster with hundreds of CPUs. This performance issue prohibits the use of fully flexible docking in the design of large combinatorial libraries. We have developed a fast structure-based screening method, which utilizes docking of a limited number of compounds to build a 2D QSAR model used to rapidly score the rest of the database. We compare here a model based on radial basis functions and a Bayesian categorization model. The number of compounds that need to be actually docked depends on the number of docking hits found. In our case studies reasonable quality models are built after docking of the number of molecules containing approximately 50 docking hits. The rest of the library is screened by the QSAR model. Optionally a fraction of the QSAR-prioritized library can be docked in order to find the true docking hits. The quality of the model only depends on the training set size - not on the size of the library to be screened. Therefore, for larger libraries the method yields higher gain in speed no change in performance. Prioritizing a large library with these models provides a significant enrichment with docking hits: it attains the values of approximately 13 and approximately 35 at the beginning of the score-sorted libraries in our two case studies: screening of the NCI collection and a combinatorial libraries on CDK2 kinase structure. With such enrichments, only a fraction of the database must actually be docked to find many of the true hits. The throughput of the method allows its use in screening of large compound collections and in the design of large combinatorial libraries. The strategy proposed has an important effect on efficiency but does not affect retrieval of actives, the latter being determined by the quality of the docking method itself.

Bayes Theorem↗

Docking of flexible ligands to flexible receptors in solution by molecular dynamics simulation.

In this paper, a method of simulating the docking of small flexible ligands to flexible receptors in water is reported. The method is based on molecular dynamics simulations and is an extension of an algorithm previously reported by Di Nola et al. (Di Nola et al., Proteins 1994;19:174-182). The method allows a fast exploration of the receptor surface, using a high temperature of the center of mass translational motion, while the ligand internal motions, the solvent, and the receptor are simulated at room temperature. In addition, the method allows a fast center of mass motion of the ligand, even in solution. The dampening effect of the solvent can be overcome by applying different weights to the interactions between system subsets (solvent, receptor, and ligand). Specific ligand-receptor distances have been used to compare the results of the simulations with the crystal structure. The method is applied, as a test system, to the docking of the phosphocholine to the immunoglobulin McPC603. The results show the similarity of structure between the complex in solution and in the crystal.

Algorithms↗

Exploring the energy landscapes of molecular recognition by a genetic algorithm: analysis of the requirements for robust docking of HIV-1 protease and FKBP-12 complexes.

Energy landscapes of molecular recognition are explored by performing "semi-rigid" docking of FK-506 and rapamycin with the Fukisawa binding protein (FKBP-12), and flexible docking simulations of the Ro-31-8959 and AG-1284 inhibitors with HIV-1 protease by a genetic algorithm. The requirements of a molecular recognition model to meet thermodynamic and kinetic criteria of ligand-protein docking simultaneously are investigated using a family of simple molecular recognition energy functions. The critical factor that determines the success rate in predicting the structure of ligand-protein complexes is found to be the roughness of the binding energy landscape, in accordance with a minimal frustration principle. The results suggest that further progress in structure prediction of ligand-protein complexes can be achieved by designing molecular recognition energy functions that generate binding landscapes with reduced frustration.

Algorithms↗