Challenges and advances in computational docking: 2009 in review

Docking is a computational technique that places a small molecule (ligand) in the binding site of its macromolecular target (receptor) and estimates its binding affinity. This review addresses methodological developments that have occurred in the docking field in 2009, with a particular focus on the more difficult, and sometimes controversial, aspects of this promising computational discipline. These developments aim to address the main challenges of docking: receptor representation (such aspects as structural waters, side chain protonation, and, most of all, flexibility (from side chain rotation to domain movement)), ligand representation (protonation, tautomerism and stereoisomerism, and the effect of input conformation), as well as accounting for solvation and entropy of binding. This review is strongly focused on docking advances in the context of drug design, specifically in virtual screening and fragment-based drug design.     

Improvements,trends, and new ideas in molecular docking: 2012–2013 in review

Molecular docking is a computational method for predicting the placement of ligands in the binding sites of their receptor(s). In this review, we discuss the methodological developments that occurred in the docking field in 2012 and 2013, with a particular focus on the more difficult aspects of this computational discipline. The main challenges and therefore focal points for developments in docking, covered in this review, are receptor flexibility, solvation, scoring, and virtual screening. We specifically deal with such aspects of molecular docking and its applications as selection criteria for constructing receptor ensembles, target dependence of scoring functions, integration of higher‐level theory into scoring, implicit and explicit handling of solvation in the binding process, and comparison and evaluation of docking and scoring methods. Copyright © 2015 John Wiley & Sons, Ltd.     

Receptor–ligand molecular docking

Docking methodology aims to predict the experimental binding modes and affinities of small molecules within the binding site of particular receptor targets and is currently used as a standard computational tool in drug design for lead compound optimisation and in virtual screening studies to find novel biologically active molecules. The basic tools of a docking methodology include a search algorithm and an energy scoring function for generating and evaluating ligand poses. In this review, we present the search algorithms and scoring functions most commonly used in current molecular docking methods that focus on protein–ligand applications. We summarise the main topics and recent computational and methodological advances in protein–ligand docking. Protein flexibility, multiple ligand binding modes and the free-energy landscape profile for binding affinity prediction are important and interconnected challenges to be overcome by further methodological developments in the docking field.     

Receptor-specific scoring functions derived from quantum chemical models improve affinity estimates for in-silico drug discovery

The adaptation of forcefield-based scoring function to specific receptors remains an important challenge for in-silico drug discovery. Here we compare binding energies of forcefield-based scoring functions with models that are reparameterized on the basis of large-scale quantum calculations of the receptor. We compute binding energies of eleven ligands to the human estrogen receptor subtype alpha (ERalpha) and four ligands to the human retinoic acid receptor of isotype gamma (RARgamma). Using the FlexScreen all-atom receptor-ligand docking approach, we compare docking simulations parameterized by quantum-mechanical calculation of a large protein fragment with purely forcefield-based models. The use of receptor flexibility in the FlexScreen permits the treatment of all ligands in the same receptor model. We find a high correlation between the classical binding energy obtained in the docking simulation and quantum mechanical binding energies and a good correlation with experimental affinities R=0.81 for ERalpha and R=0.95 for RARgamma using the quantum derived scoring functions. A significant part of this improvement is retained, when only the receptor is treated with quantum-based parameters, while the ligands are parameterized with a purely classical model.     

Molecular docking towards drug discovery

Fueled by advances in molecular structure determination, tools for structure-based drug design are proliferating rapidly. Lead discovery through searching of ligand databases with molecular docking techniques represents an attractive alternative to high-throughput random screening. The size of commercial databases imposes severe computational constraints on molecular docking, compromising the level of calculational detail permitted for each putative ligand. We describe alternative philosophies for docking which effectively address this challenge. With respect to the dynamic aspects of molecular recognition, these strategies lie along a spectrum of models bounded by the Lock-and-Key and Induced-Fit theories for ligand binding. We explore the potential of a rigid model in exploiting species specificity and of a tolerant model in predicting absolute ligand binding affinity. Current molecular docking methods are limited primarily by their ability to rank docked complexes; we therefore place particular emphasis on this aspect of the problem throughout our validation of docking strategies.     

How to choose relevant multiple receptor conformations for virtual screening: a test case of Cdk2 and normal mode analysis

Better treatment of protein flexibility is essential in structure-based drug design projects such as virtual screening and protein-ligand docking. Diversity in ligand-binding mechanisms and receptor conformational changes makes it difficult to treat dynamic features of the receptor during the docking simulation. Thus, the use of pregenerated multiple receptor conformations is applied today in virtual screening studies. However, generation of a small relevant set of receptor conformations remains challenging. To address this problem, we propose a new protocol for the generation of multiple receptor conformations via normal mode analysis and for the selection of several receptor conformations suitable for docking/virtual screening. We validated this protocol on cyclin-dependent kinase 2, which possesses a binding site located at the interface between two subdomains and is known to undergo significant conformational changes in the active site region upon ligand binding. We believe that the suggested rules for the choice of suitable receptor conformations can be applied to other targets when dealing with in silico screening on flexible receptors.     

From G Protein-coupled Receptor Structure Resolution to Rational Drug Design

A number of recent technical solutions have led to significant advances in G protein-coupled receptor (GPCR) structural biology. Apart from a detailed mechanistic view of receptor activation, the new structures have revealed novel ligand binding sites. Together, these insights provide avenues for rational drug design to modulate the activities of these important drug targets. The application of structural data to GPCR drug discovery ushers in an exciting era with the potential to improve existing drugs and discover new ones. In this review, we focus on technical solutions that have accelerated GPCR crystallography as well as some of the salient findings from structures that are relevant to drug discovery. Finally, we outline some of the approaches used in GPCR structure based drug design.     

An NMR-based scoring function improves the accuracy of binding pose predictions by docking by two orders of magnitude

Low-affinity ligands can be efficiently optimized into high-affinity drug leads by structure based drug design when atomic-resolution structural information on the protein/ligand complexes is available. In this work we show that the use of a few, easily obtainable, experimental restraints improves the accuracy of the docking experiments by two orders of magnitude. The experimental data are measured in nuclear magnetic resonance spectra and consist of protein-mediated NOEs between two competitively binding ligands. The methodology can be widely applied as the data are readily obtained for low-affinity ligands in the presence of non-labelled receptor at low concentration. The experimental inter-ligand NOEs are efficiently used to filter and rank complex model structures that have been pre-selected by docking protocols. This approach dramatically reduces the degeneracy and inaccuracy of the chosen model in docking experiments, is robust with respect to inaccuracy of the structural model used to represent the free receptor and is suitable for high-throughput docking campaigns.     

Ensemble Docking in Drug Discovery

Ensemble docking corresponds to the generation of an “ensemble” of drug target conformations in computational structure-based drug discovery, often obtained by using molecular dynamics simulation, that is used in docking candidate ligands. This approach is now well established in the field of early-stage drug discovery. This review gives a historical account of the development of ensemble docking and discusses some pertinent methodological advances in conformational sampling.     

Mapping of ligand‐binding cavities in proteins

The complex interactions between proteins and small organic molecules (ligands) are intensively studied because they play key roles in biological processes and drug activities. Here, we present a novel approach to characterize and map the ligand‐binding cavities of proteins without direct geometric comparison of structures, based on Principal Component Analysis of cavity properties (related mainly to size, polarity, and charge). This approach can provide valuable information on the similarities and dissimilarities, of binding cavities due to mutations, between‐species differences and flexibility upon ligand‐binding. The presented results show that information on ligand‐binding cavity variations can complement information on protein similarity obtained from sequence comparisons. The predictive aspect of the method is exemplified by successful predictions of serine proteases that were not included in the model construction. The presented strategy to compare ligand‐binding cavities of related and unrelated proteins has many potential applications within protein and medicinal chemistry, for example in the characterization and mapping of “orphan structures”, selection of protein structures for docking studies in structure‐based design, and identification of proteins for selectivity screens in drug design programs. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

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.     

Ligand-protein inverse docking and its potential use in the computer search of protein targets of a small molecule

Ligand-protein docking has been developed and used in facilitating new drug discoveries. In this approach, docking single or multiple small molecules to a receptor site is attempted to find putative ligands. A number of studies have shown that docking algorithms are capable of finding ligands and binding conformations at a receptor site close to experimentally determined structures. These algorithms are expected to be equally applicable to the identification of multiple proteins to which a small molecule can bind or weakly bind. We introduce a ligand-protein inverse-docking approach for finding potential protein targets of a small molecule by the computer-automated docking search of a protein cavity database. This database is developed from protein structures in the Protein Data Bank (PDB). Docking is conducted with a procedure involving multiple-conformer shape-matching alignment of a molecule to a cavity followed by molecular-mechanics torsion optimization and energy minimization on both the molecule and the protein residues at the binding region. Scoring is conducted by the evaluation of molecular-mechanics energy and, when applicable, by the further analysis of binding competitiveness against other ligands that bind to the same receptor site in at least one PDB entry. Testing results on two therapeutic agents, 4H-tamoxifen and vitamin E, showed that 50% of the computer-identified potential protein targets were implicated or confirmed by experiments. The application of this approach may facilitate the prediction of unknown and secondary therapeutic target proteins and those related to the side effects and toxicity of a drug or drug candidate. Proteins 2001;43:217-226.     

Discovery of binding proteins for a protein target using protein–protein docking‐based virtual screening

Target structure‐based virtual screening, which employs protein‐small molecule docking to identify potential ligands, has been widely used in small‐molecule drug discovery. In the present study, we used a protein–protein docking program to identify proteins that bind to a specific target protein. In the testing phase, an all‐to‐all protein–protein docking run on a large dataset was performed. The three‐dimensional rigid docking program SDOCK was used to examine protein–protein docking on all protein pairs in the dataset. Both the binding affinity and features of the binding energy landscape were considered in the scoring function in order to distinguish positive binding pairs from negative binding pairs. Thus, the lowest docking score, the average Z‐score, and convergency of the low‐score solutions were incorporated in the analysis. The hybrid scoring function was optimized in the all‐to‐all docking test. The docking method and the hybrid scoring function were then used to screen for proteins that bind to tumor necrosis factor‐α (TNFα), which is a well‐known therapeutic target for rheumatoid arthritis and other autoimmune diseases. A protein library containing 677 proteins was used for the screen. Proteins with scores among the top 20% were further examined. Sixteen proteins from the top‐ranking 67 proteins were selected for experimental study. Two of these proteins showed significant binding to TNFα in an in vitro binding study. The results of the present study demonstrate the power and potential application of protein–protein docking for the discovery of novel binding proteins for specific protein targets. Proteins 2014; 82:2472–2482. © 2014 Wiley Periodicals, Inc.     

in Silico investigation of the structural requirements for the AMPA receptor antagonism by quinoxaline derivatives

Glutamate receptors have been implicated in various neurological disorders and their antagonism offers a suitable approach for the treatment of such disorders. The field of drug design and discovery aims to find best medicines to prevent, treat and cure diseases quickly and efficiently. In this regard, computational tools have helped medicinal chemists modify and optimize molecules to potent drug candidates with better pharmacokinetic profiles, and guiding biologists and pharmacologists to explore new disease genes as well as novel drug targets. In the present study, to understand the structural requirements for AMPA receptor antagonism, molecular docking study was performed on 41 structurally diverse antagonists based on quinoxaline nucleus. Lamarckian genetic algorithm methodology was employed for docking simulations using AutoDock 4.2 program. The results obtained signify that the molecular docking approach is reliable and produces a good correlation coefficient (r² = 0.6) between experimental and docking predicted AMPA receptor antagonistic activity. The aromatic moiety of quinoxaline core has been proved to be vital for hydrophobic contacts exhibiting - interactions in docked conformations. However, polar moieties such as carboxylic group and 1,2,4-triazole moieties were noted to be sites for hydrophilic interactions in terms of hydrogen bonding with the receptor. These analyses can be exploited to design and develop novel AMPA receptor antagonists for the treatment of different neurological disorders.     

Molecular insights into the binding selectivity of a synthetic ligand DAAG to Fc fragment of IgG

Affinity chromatography with synthetic ligands has been focused as the potential alternative to protein A‐based chromatography for antibody capture because of its comparable selectivity and efficiency. Better understanding on the molecular interactions between synthetic ligand and antibody is crucial for improving and designing novel ligands. In this work, the molecular interaction mechanism between Fc fragment of IgG and a synthetic ligand (DAAG) was studied with molecular docking and dynamics simulation. The docking results on the consensus binding site (CBS) indicated that DAAG could bind to the CBS with the favorable orientation like a tripod for the top‐ranked binding complexes. The ligand‐Fc fragment complexes were then tested by molecular dynamics simulation at neutral condition (pH 7.0) for 10 ns. The results indicated that the binding of DAAG on the CBS of Fc fragment was achieved by the multimodal interactions, combining the hydrophobic interaction, electrostatic interaction, hydrogen bond, and so on. It was also found that multiple secondary interactions endowed DAAG with an excellent selectivity to Fc fragment. In addition, molecular dynamics simulation conducted at acidic condition (pH 3.0) showed that the departure of DAAG ligand from the surface of Fc fragment was the result of reduced interaction energies. The binding modes between DAAG and CBS not only shed light on the molecular mechanisms of DAAG for antibody purification but also provide useful information for the improvement of ligand design. Copyright © 2014 John Wiley & Sons, Ltd.     

Hit clustering can improve virtual fragment screening: CDK2 and PARP1 case studies

Virtual fragment screening could be a promising alternative to existing experimental screening techniques. However, reliable methods of in silico fragment screening are yet to be established and validated. In order to develop such an approach we first checked how successful the existing molecular docking methods can be in predicting fragment binding affinities and poses. Using our Lead Finder docking software the RMSD of the binding energy prediction was observed to be 1.35 kcal/mol(-1) on a set of 26 experimentally characterized fragment inhibitors, and the RMSD of the predicted binding pose from the experimental one was <1.5 ?. Then, we explored docking of 68 fragments obtained from 39 drug molecules for which co-crystal structures were available from the PDB. It appeared that fragments that participate in oriented non-covalent interactions, such as hydrogen bonds and metal coordination, could be correctly docked in 70-80% of cases suggesting the potential success of rediscovering of corresponding drugs by in silico fragment approach. Based on these findings we've developed a virtual fragment screening technique which involved structural filtration of protein-ligand complexes for specific interactions and subsequent clustering in order to minimize the number of preferable starting fragment candidates. Application of this method led to 2 millimolar-scale fragment PARP1 inhibitors with a new scaffold.     

Flexible ligand docking to multiple receptor conformations: a practical alternative

State of the art docking algorithms predict an incorrect binding pose for about 50-70% of all ligands when only a single fixed receptor conformation is considered. In many more cases, lack of receptor flexibility results in meaningless ligand binding scores, even when the correct pose is obtained. Incorporating conformational rearrangements of the receptor binding pocket into predictions of both ligand binding pose and binding score is crucial for improving structure-based drug design and virtual ligand screening methodologies. However, direct modeling of protein binding site flexibility remains challenging because of the large conformational space that must be sampled, and difficulties remain in constructing a suitably accurate energy function. Here we show that using multiple fixed receptor conformations, either experimentally determined by crystallography or NMR, or computationally generated, is a practical shortcut that may improve docking calculations. In several cases, such an approach has led to experimentally validated predictions.     

Modeling G protein‐coupled receptors for structure‐based drug discovery using low‐frequency normal modes for refinement of homology models: Application to H3 antagonists

G Protein‐Coupled Receptors (GPCRs) are integral membrane proteins that play important role in regulating key physiological functions, and are targets of about 50% of all recently launched drugs. High‐resolution experimental structures are available only for very few GPCRs. As a result, structure‐based drug design efforts for GPCRs continue to rely on in silico modeling, which is considered to be an extremely difficult task especially for these receptors. Here, we describe Gmodel, a novel approach for building 3D atomic models of GPCRs using a normal mode‐based refinement of homology models. Gmodel uses a small set of relevant low‐frequency vibrational modes derived from Random Elastic Network model to efficiently sample the large‐scale receptor conformation changes and generate an ensemble of alternative models. These are used to assemble receptor–ligand complexes by docking a known active into each of the alternative models. Each of these is next filtered using restraints derived from known mutation and binding affinity data and is refined in the presence of the active ligand. In this study, Gmodel was applied to generate models of the antagonist form of histamine 3 (H3) receptor. The validity of this novel modeling approach is demonstrated by performing virtual screening (using the refined models) that consistently produces highly enriched hit lists. The models are further validated by analyzing the available SAR related to classical H3 antagonists, and are found to be in good agreement with the available experimental data, thus providing novel insights into the receptor–ligand interactions. Proteins 2010. © 2009 Wiley‐Liss, Inc.