Toward prediction of functional protein pockets using blind docking and pocket search algorithms

Location of functional binding pockets of bioactive ligands on protein molecules is essential in structural genomics and drug design projects. If the experimental determination of ligand-protein complex structures is complicated, blind docking (BD) and pocket search (PS) calculations can help in the prediction of atomic resolution binding mode and the location of the pocket of a ligand on the entire protein surface. Whereas the number of successful predictions by these methods is increasing even for the complicated cases of exosites or allosteric binding sites, their reliability has not been fully established. For a critical assessment of reliability, we use a set of ligand-protein complexes, which were found to be problematic in previous studies. The robustness of BD and PS methods is addressed in terms of success of the selection of truly functional pockets from among the many putative ones identified on the surfaces of ligand-bound and ligand-free (holo and apo) protein forms. Issues related to BD such as effect of hydration, existence of multiple pockets, and competition of subsidiary ligands are considered. Practical cases of PS are discussed, categorized and strategies are recommended for handling the different situations. PS can be used in conjunction with BD, as we find that a consensus approach combining the techniques improves predictive power.     

Automatic identification and representation of protein binding sites for molecular docking.

Molecular docking is a popular way to screen for novel drug compounds. The method involves aligning small molecules to a protein structure and estimating their binding affinity. To do this rapidly for tens of thousands of molecules requires an effective representation of the binding region of the target protein. This paper presents an algorithm for representing a protein's binding site in a way that is specifically suited to molecular docking applications. Initially the protein's surface is coated with a collection of molecular fragments that could potentially interact with the protein. Each fragment, or probe, serves as a potential alignment point for atoms in a ligand, and is scored to represent that probe's affinity for the protein. Probes are then clustered by accumulating their affinities, where high affinity clusters are identified as being the "stickiest" portions of the protein surface. The stickiest cluster is used as a computational binding "pocket" for docking. This method of site identification was tested on a number of ligand-protein complexes; in each case the pocket constructed by the algorithm coincided with the known ligand binding site. Successful docking experiments demonstrated the effectiveness of the probe representation.     

Towards understanding the mechanisms of molecular recognition by computer simulations of ligand-protein interactions

The thermodynamic and kinetic aspects of molecular recognition for the methotrexate (MTX)-dihydrofolate reductase (DHFR) ligand-protein system are investigated by the binding energy landscape approach. The impact of 'hot' and 'cold' errors in ligand mutations on the thermodynamic stability of the native MTX-DHFR complex is analyzed, and relationships between the molecular recognition mechanism and the degree of ligand optimization are discussed. The nature and relative stability of intermediates and thermodynamic phases on the ligand-protein association pathway are studied, providing new insights into connections between protein folding and molecular recognition mechanisms, and cooperativity of ligand-protein binding. The results of kinetic docking simulations are rationalized based on the thermodynamic properties determined from equilibrium simulations and the shape of the underlying binding energy landscape. We show how evolutionary ligand selection for a receptor active site can produce well-optimized ligand-protein systems such as MTX-DHFR complex with the thermodynamically stable native structure and a direct transition mechanism of binding from unbound conformations to the unique native structure.     

Molecular dynamics simulations suggest Thiosemicarbazones can bind p53 cancer mutant R175H

Cancer pathologies are associated with the unfolding and aggregation of most recurring mutations in the DNA Binding Domain (DBD) of p53 that coordinate the destabilization of protein. Substitution at the 175th codon with arginine to histidine (R175H, a mutation of large to small side-chain amino acid) destabilizes the DBD by 3 kcal/mol and triggers breasts, lung cancer, etc. Stabilizing the p53 mutant by small molecules offers an attractive drug-targeted anti-cancer therapy. The thiosemicarbazone (TSC) molecules NPC and DPT are known to act as zinc-metallochaperones to reactivate p53R175H. Here, a combination of LESMD simulations for 10 TSC conformations with a p53R175H receptor, single ligand-protein conformation MD, and ensemble docking with multiple p53R175H conformations observed during simulations is suggested to identify the potential binding site of the target protein in light of their importance for the direct TSC – p53R175H binding. NPC binds mutant R175H in the loop region L2-L3, forming pivotal hydrogen bonds with HIS175, pi?sulfur bonds with TYR163, and pi-alkyl linkages with ARG174 and PRO190, all of which are contiguous to the zinc-binding native site on p53DBD. DPT, on the other hand, was primarily targeting alternative binding sites such as the loop-helix L1/H2 region and the S8 strand. The similar structural characteristics of TSC-bound p53R175H complexes with wild-type p53DBD are thought to be attributable to involved interactions that favour binding free energy contributions of TSC ligands. Our findings may be useful in the identification of novel pockets with druggable properties.     

Screening and docking chemical ligands onto pocket cavities of a protease for designing a biocatalyst

A detailed study of the trypsin surface has been carried out to gain insight into its biological functions and interactions which helped to determine the binding specificity. Twenty-four cavity pockets were automatically identified on trypsin from PDB file entry 1AUJ using CASTp (Computed Atlas of Surface Topography of proteins). Molecular docking was exploited as an efficient in silico screening tool for studying protein–ligand interactions. A systematic docking study using Autodock 3.05 has been performed on the five largest binding pockets in trypsin. A set of ten putative chemical ligands was used to dock into selected binding pockets. Docking of ligands into the five largest pockets in trypsin showed that 1,10-phenanthroline and ethanolamine preferentially bound at pocket 24 and benzamidine at pocket 22. Thermodynamically, we also found that ethanol, propanol, propandiol and phosphoethanolamine preferentially bound at pocket 21 whereas p-aminobenzamidine, phenylacetic acid and phenylalanine interacted mainly at pocket 20 based on their lowest interaction free energy.     

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. © 1996 Wiley-Liss, Inc.     

In silico assessment of potential druggable pockets on the surface of α1-antitrypsin conformers

The search for druggable pockets on the surface of a protein is often performed on a single conformer, treated as a rigid body. Transient druggable pockets may be missed in this approach. Here, we describe a methodology for systematic in silico analysis of surface clefts across multiple conformers of the metastable protein α(1)-antitrypsin (A1AT). Pathological mutations disturb the conformational landscape of A1AT, triggering polymerisation that leads to emphysema and hepatic cirrhosis. Computational screens for small molecule inhibitors of polymerisation have generally focused on one major druggable site visible in all crystal structures of native A1AT. In an alternative approach, we scan all surface clefts observed in crystal structures of A1AT and in 100 computationally produced conformers, mimicking the native solution ensemble. We assess the persistence, variability and druggability of these pockets. Finally, we employ molecular docking using publicly available libraries of small molecules to explore scaffold preferences for each site. Our approach identifies a number of novel target sites for drug design. In particular one transient site shows favourable characteristics for druggability due to high enclosure and hydrophobicity. Hits against this and other druggable sites achieve docking scores corresponding to a K(d) in the μM-nM range, comparing favourably with a recently identified promising lead. Preliminary ThermoFluor studies support the docking predictions. In conclusion, our strategy shows considerable promise compared with the conventional single pocket/single conformer approach to in silico screening. Our best-scoring ligands warrant further experimental investigation.     

Efficient molecular docking of NMR structures: application to HIV-1 protease

Docking ligands into an ensemble of NMR conformers is essential to structure-based drug discovery if only NMR structures are available for the target. However, sequentially docking ligands into each NMR conformer through standard single-receptor-structure docking, referred to as sequential docking, is computationally expensive for large-scale database screening because of the large number of NMR conformers involved. Recently, we developed an efficient ensemble docking algorithm to consider protein structural variations in ligand binding. The algorithm simultaneously docks ligands into an ensemble of protein structures and achieves comparable performance to sequential docking without significant increase in computational time over single-structure docking. Here, we applied this algorithm to docking with NMR structures. The HIV-1 protease was used for validation in terms of docking accuracy and virtual screening. Ensemble docking of the NMR structures identified 91% of the known inhibitors under the criterion of RMSD < 2.0 A for the best-scored conformation, higher than the average success rate of single docking of individual crystal structures (66%). In the virtual screening test, on average, ensemble docking of the NMR structures obtained higher enrichments than single-structure docking of the crystal structures. In contrast, docking of either the NMR minimized average structure or a single NMR conformer performed less satisfactorily on both binding mode prediction and virtual screening, indicating that a single NMR structure may not be suitable for docking calculations. The success of ensemble docking of the NMR structures suggests an efficient alternative method for standard single docking of crystal structures and for considering protein flexibility.     

Comparative surface geometry of the protein kinase family

Identifying conserved pockets on the surfaces of a family of proteins can provide insight into conserved geometric features and sites of protein–protein interaction. Here we describe mapping and comparison of the surfaces of aligned crystallographic structures, using the protein kinase family as a model. Pockets are rapidly computed using two computer programs, FADE and Crevasse. FADE uses gradients of atomic density to locate grooves and pockets on the molecular surface. Crevasse, a new piece of software, splits the FADE output into distinct pockets. The computation was run on 10 kinase catalytic cores aligned on the αF‐helix, and the resulting pockets spatially clustered. The active site cleft appears as a large, contiguous site that can be subdivided into nucleotide and substrate docking sites. Substrate specificity determinants in the active site cleft between serine/threonine and tyrosine kinases are visible and distinct. The active site clefts cluster tightly, showing a conserved spatial relationship between the active site and αF‐helix in the C‐lobe. When the αC‐helix is examined, there are multiple mechanisms for anchoring the helix using spatially conserved docking sites. A novel site at the top of the N‐lobe is present in all the kinases, and there is a large conserved pocket over the hinge and the αC‐β4 loop. Other pockets on the kinase core are strongly conserved but have not yet been mapped to a protein–protein interaction. Sites identified by this algorithm have revealed structural and spatially conserved features of the kinase family and potential conserved intermolecular and intramolecular binding sites.     

Detection of multiscale pockets on protein surfaces using mathematical morphology

Detection of pockets on protein surfaces is an important step toward finding the binding sites of small molecules. In a previous study, we defined a pocket as a space into which a small spherical probe can enter, but a large probe cannot. The radius of the large probes corresponds to the shallowness of pockets. We showed that each type of binding molecule has a characteristic shallowness distribution. In this study, we introduced fundamental changes to our previous algorithm by using a 3D grid representation of proteins and probes, and the theory of mathematical morphology. We invented an efficient algorithm for calculating deep and shallow pockets (multiscale pockets) simultaneously, using several different sizes of spherical probes (multiscale probes). We implemented our algorithm as a new program, ghecom (grid‐based HECOMi finder). The statistics of calculated pockets for the structural dataset showed that our program had a higher performance of detecting binding pockets, than four other popular pocket‐finding programs proposed previously. The ghecom also calculates the shallowness of binding ligands, R_inaccess (minimum radius of inaccessible spherical probes) that can be obtained from the multiscale molecular volume. We showed that each part of the binding molecule had a bias toward a specific range of shallowness. These findings will be useful for predicting the types of molecules that will be most likely to bind putative binding pockets, as well as the configurations of binding molecules. The program ghecom is available through the Web server ( http://biunit.naist.jp/ghecom ). Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Screening and docking chemical ligands onto pocket cavities of a protease for designing a biocatalyst

A detailed study of the trypsin surface has been carried out to gain insight into its biological functions and interactions which helped to determine the binding specificity. Twenty-four cavity pockets were automatically identified on trypsin from PDB file entry 1AUJ using CASTp (Computed Atlas of Surface Topography of proteins). Molecular docking was exploited as an efficient in silico screening tool for studying protein-ligand interactions. A systematic docking study using Autodock 3.05 has been performed on the five largest binding pockets in trypsin. A set of ten putative chemical ligands was used to dock into selected binding pockets. Docking of ligands into the five largest pockets in trypsin showed that 1,10-phenanthroline and ethanolamine preferentially bound at pocket 24 and benzamidine at pocket 22. Thermodynamically, we also found that ethanol, propanol, propandiol and phosphoethanolamine preferentially bound at pocket 21 whereas p-aminobenzamidine, phenylacetic acid and phenylalanine interacted mainly at pocket 20 based on their lowest interaction free energy.     

Assessment and Challenges of Ligand Docking into Comparative Models of G-Protein Coupled Receptors

The rapidly increasing number of high-resolution X-ray structures of G-protein coupled receptors (GPCRs) creates a unique opportunity to employ comparative modeling and docking to provide valuable insight into the function and ligand binding determinants of novel receptors, to assist in virtual screening and to design and optimize drug candidates. However, low sequence identity between receptors, conformational flexibility, and chemical diversity of ligands present an enormous challenge to molecular modeling approaches. It is our hypothesis that rapid Monte-Carlo sampling of protein backbone and side-chain conformational space with Rosetta can be leveraged to meet this challenge. This study performs unbiased comparative modeling and docking methodologies using 14 distinct high-resolution GPCRs and proposes knowledge-based filtering methods for improvement of sampling performance and identification of correct ligand-receptor interactions. On average, top ranked receptor models built on template structures over 50% sequence identity are within 2.9 Å of the experimental structure, with an average root mean square deviation (RMSD) of 2.2 Å for the transmembrane region and 5 Å for the second extracellular loop. Furthermore, these models are consistently correlated with low Rosetta energy score. To predict their binding modes, ligand conformers of the 14 ligands co-crystalized with the GPCRs were docked against the top ranked comparative models. In contrast to the comparative models themselves, however, it remains difficult to unambiguously identify correct binding modes by score alone. On average, sampling performance was improved by 10³ fold over random using knowledge-based and energy-based filters. In assessing the applicability of experimental constraints, we found that sampling performance is increased by one order of magnitude for every 10 residues known to contact the ligand. Additionally, in the case of DOR, knowledge of a single specific ligand-protein contact improved sampling efficiency 7 fold. These findings offer specific guidelines which may lead to increased success in determining receptor-ligand complexes.     

Transient pockets on XIAP-BIR2: toward the characterization of putative binding sites of small-molecule XIAP inhibitors

Protein-protein interactions are abundant in signal transduction pathways and thus of crucial importance in the regulation of apoptosis. However, designing small-molecule inhibitors for these potential drug targets is very challenging as such proteins often lack well-defined binding pockets. An example for such an interaction is the binding of the anti-apoptotic BIR2 domain of XIAP to the pro-apoptotic caspase-3 that results in the survival of damaged cells. Although small-molecule inhibitors of this interaction have been identified, their exact binding sites on XIAP are not known as its crystal structures reveal no suitable pockets. Here, we apply our previously developed protocol for identifying transient binding pockets to XIAP-BIR2. Transient pockets were identified in snapshots taken during four different molecular dynamics simulations that started from the caspase-3:BIR2 complex or from the unbound BIR2 structure and used water or methanol as solvent. Clustering of these pockets revealed that surprisingly many pockets opened in the flexible linker region that is involved in caspase-3 binding. We docked three known inhibitors into these transient pockets and so determined five putative binding sites. In addition, by docking two inactive compounds of the same series, we show that this protocol is also able to distinguish between binders and nonbinders which was not possible when docking to the crystal structures. These findings represent a first step toward the understanding of the binding of small-molecule XIAP-BIR2 inhibitors on a molecular level and further highlight the importance of considering protein flexibility when designing small-molecule protein-protein interaction inhibitors.     

Ensemble docking of multiple protein structures: considering protein structural variations in molecular docking

One approach to incorporate protein flexibility in molecular docking is the use of an ensemble consisting of multiple protein structures. Sequentially docking each ligand into a large number of protein structures is computationally too expensive to allow large-scale database screening. It is challenging to achieve a good balance between docking accuracy and computational efficiency. In this work, we have developed a fast, novel docking algorithm utilizing multiple protein structures, referred to as ensemble docking, to account for protein structural variations. The algorithm can simultaneously dock a ligand into an ensemble of protein structures and automatically select an optimal protein structure that best fits the ligand by optimizing both ligand coordinates and the conformational variable m, where m represents the m-th structure in the protein ensemble. The docking algorithm was validated on 10 protein ensembles containing 105 crystal structures and 87 ligands in terms of binding mode and energy score predictions. A success rate of 93% was obtained with the criterion of root-mean-square deviation <2.5 A if the top five orientations for each ligand were considered, comparable to that of sequential docking in which scores for individual docking are merged into one list by re-ranking, and significantly better than that of single rigid-receptor docking (75% on average). Similar trends were also observed in binding score predictions and enrichment tests of virtual database screening. The ensemble docking algorithm is computationally efficient, with a computational time comparable to that for docking a ligand into a single protein structure. In contrast, the computational time for the sequential docking method increases linearly with the number of protein structures in the ensemble. The algorithm was further evaluated using a more realistic ensemble in which the corresponding bound protein structures of inhibitors were excluded. The results show that ensemble docking successfully predicts the binding modes of the inhibitors, and discriminates the inhibitors from a set of noninhibitors with similar chemical properties. Although multiple experimental structures were used in the present work, our algorithm can be easily applied to multiple protein conformations generated by computational methods, and helps improve the efficiency of other existing multiple protein structure(MPS)-based methods to accommodate protein flexibility.     

Predicting structural effects in HIV-1 protease mutant complexes with flexible ligand docking and protein side-chain optimization

We present a computational approach for predicting structures of ligand-protein complexes and analyzing binding energy landscapes that combines Monte Carlo simulated annealing technique to determine the ligand bound conformation with the dead-end elimination algorithm for side-chain optimization of the protein active site residues. Flexible ligand docking and optimization of mobile protein side-chains have been performed to predict structural effects in the V32I/I47V/V82I HIV-1 protease mutant bound with the SB203386 ligand and in the V82A HIV-1 protease mutant bound with the A77003 ligand. The computational structure predictions are consistent with the crystal structures of these ligand-protein complexes. The emerging relationships between ligand docking and side-chain optimization of the active site residues are rationalized based on the analysis of the ligand-protein binding energy landscape. Proteins 33:295–310, 1998. © 1998 Wiley-Liss, Inc.     

Protein-protein interactions: hot spots and structurally conserved residues often locate in complemented pockets that pre-organized in the unbound states: implications for docking

Energetic hot spots account for a significant portion of the total binding free energy and correlate with structurally conserved interface residues. Here, we map experimentally determined hot spots and structurally conserved residues to investigate their geometrical organization. Unfilled pockets are pockets that remain unfilled after protein-protein complexation, while complemented pockets are pockets that disappear upon binding, representing tightly fit regions. We find that structurally conserved residues and energetic hot spots are strongly favored to be located in complemented pockets, and are disfavored in unfilled pockets. For the three available protein-protein complexes with complemented pockets where both members of the complex were alanine-scanned, 62% of all hot spots (DeltaDeltaG>2kcal/mol) are within these pockets, and 60% of the residues in the complemented pockets are hot spots. 93% of all red-hot residues (DeltaDeltaG>/=4kcal/mol) either protrude into or are located in complemented pockets. The occurrence of hot spots and conserved residues in complemented pockets highlights the role of local tight packing in protein associations, and rationalizes their energetic contribution and conservation. Complemented pockets and their corresponding protruding residues emerge among the most important geometric features in protein-protein interactions. By screening the solvent, this organization shields backbone hydrogen bonds and charge-charge interactions. Complemented pockets often pre-exist binding. For 18 protein-protein complexes with complemented pockets whose unbound structures are available, in 16 the pockets are identified to pre-exist in the unbound structures. The root-mean-squared deviations of the atoms lining the pockets between the bound and unbound states is as small as 0.9A, suggesting that such pockets constitute features of the populated native state that may be used in docking.