JSM Chemistry

Challenges in Docking: Mini Review

Mini Review | Open Access

  • 1. Department of Pharmaceutical Chemistry, Jamia Hamdard, India
+ Show More - Show Less
Corresponding Authors
Mymoona Akhter, Drug Department of Pharmaceutical Chemistry, Jamia Hamdard, New Delhi

Docking is a computational technique that helps in understanding of ligand macromolecule interactions by placing the ligand in the binding site of macromolecular target. Understanding the ruling principles whereby protein receptors recognize, interact, and associate with molecular substrates and inhibitors is of paramount importance in drug discovery efforts. Despite enormous advances in the field of computational science and biology over the last few decades and the widespread application of docking methods, several pitfalls still exist. In this review the key concepts of protein-ligand docking methods are outlined, with emphasis on challenges faced currently are highlighted. In particular, ligand and protein flexibility—a critical aspect for a thorough understanding of the principles that guide ligand binding in proteins has been discussed. Developments related to receptor representation like flexibility, water molecules, ligand representation like protonation, tautomerism and stereoisomerism have been touched. This review is strongly focused on docking challenges in the context of drug design, specifically structure based drug design.


•    Computational docking
•    Receptor flexibility
•    Ligand flexibility
•    Protein-ligand complexes 


Akhter M (2016) Challenges in Docking: Mini Review. JSM Chem 4(2): 1025.


RNA: Ribonucleic Acid; DNA: Deoxyribonucleic Acid; GA: Genetic Algorithms; PSO: Particle Swarm Optimisation Algorithms; ACO: Ant Colony Optimisation Algorithms; TS: Tabu Search Algorithms; MC: Monte Carlo Algorithms; NMR: Nuclear Magnetic Resonance


Molecular docking is a key tool in structural molecular biology and computer-assisted drug design. Post-genomic era, docking has becomes an indispensable part of the drug discovery and genome informatics. Docking is a computational method that attempt to find the “best” matching between two molecules a small molecule (ligand) and a macromolecule (Protein, DNA, RNA) or between two macromolecules of different size (Figure 1). From the very first docking program ‘DOCK’ developed by Kuntz in 1982 [1] to dock rigid molecules one by one to rigid proteins to state-of-the-art programs that can dock complete libraries of highly flexible small molecules to flexible proteins in the presence of water and metals molecules, large amount of work has been reported and the complexity of the programs has gradually increased. Broadly two docking procedures are commonly considered 1. Rigid docking and 2. Flexible docking. In Rigid Docking (or Lock and Key) considers essentially only geometric complementarities between a ligand and receptor, it does not takes into account the flexibility of both the molecules in question. The rigidity may have limited the specificity and accuracy of the results but the technique is capable of identifying ligand binding site for several different proteins precisely. This approach is more commonly used for protein –protein docking now because complexity involved in modeling the flexibility of protein molecules. In flexible Docking (or Induced fit) the ligand or protein or both molecules at the same time are kept flexible and the energy for different conformations of the ligand fitting into the protein is calculated. Compared to rigid docking flexible docking is more specific but at the same time demands for computational power and CPU time. Docking protocols includes combination of a search algorithm and a scoring function. Execution of molecular docking requires.

1. Structural data
2. Protein target of interest
3. Procedure to estimate protein-ligand interaction poses and strengths

Ligand sampling algorithms are essentials for generating acceptable ligand poses. The search algorithm should allow the degrees of freedom of the protein–ligand complex to hold the true binding poses. A good search algorithm must have good speed and must effectively cover relevant conformational space. Figure (2) gives the outlay of flexible docking search algorithm. These algorithms are of following types- Search algorithm (Pose generation): The search algorithm is a process where all possible conformations and orientations of the ligand- receptor complex in a space i.e. binding pocket of receptor is being searched [2].

Matching Algorithms

In matching algorithm the shape of the ligand is matched with that of the binding pocket by placing the ligand in the binding site of the protein1 . But before actually performing the docking procedure different conformations of the ligand are generated using other programs (e.g., OMEGA [3,4 ], Corina[5]).

Incremental Construction Algorithm or Fragmentation Based

In fragment based algorithm reconstruction of the ligand in the active site of receptor by using an anchor fragment forms the basis of Incremental construction algorithm. Rotatable bonds are identified in the ligand and then the ligand is dissected down at rotatable bonds to obtain rigid pieces. Then docking is performed with a terminal fragment and its best scoring poses is kept and used as anchors on which the next fragment are progressively added and minimized until the molecule has been fully reconstructed. Docking programs which use incremental construction algorithms include DOCK, FLEXX, FLOG. Glide also uses the same algorithm and number of articles can be seen in literature using these software’s because of this algorithm [6].

Stochastic Algorithms

Stochastic search methods modify the conformation of the small molecule in the receptor site and assess it on the fly. Stochastic algorithms represent another class of search algorithms like Genetic Algorithms (GA), Particle Swarm Optimisation Algorithms (PSO), Ant Colony Optimisation Algorithms (ACO), Tabu Search Algorithms (TS) and Monte Carlo Algorithms (MC). Several docking programs based on GA or other evolutionary techniques are Auto Dock, FITTED and GOLD (Figure 3).techniques are Auto Dock, FITTED and GOLD (Figure 3).

Docking Optimization (Scoring Function or Pose Selection)

Scoring functions help to estimate the binding affinity of ligand poses. Good scoring function should represent the thermodynamics of interaction of the protein–ligand system, so as to assist in distinguishing the true binding modes from all the others explored poses, and to rank them accordingly. The accuracy of docking depends on the quality of scoring functions which are mathematical approximating methods for estimating binding affinity i.e. finding the highest-affinity ligand against a target. In general the score is represented by following equation:

Scoring functions help to estimate the binding affinity of ligand poses. Good scoring function should represent the thermodynamics of interaction of the protein–ligand system, so as to assist in distinguishing the true binding modes from all the others explored poses, and to rank them accordingly. The accuracy of docking depends on the quality of scoring functions which are mathematical approximating methods for estimating binding affinity i.e. finding the highest-affinity ligand against a target. In general the score is represented by following equation:

Score = S (target-ligand) + S (ligand)

The term S (target-ligand) is a sum over contributions from all heavy atom contacts between the ligand and the receptor involved in binding.

The binding between a ligand and its receptor is controlled by several factors like interaction between the ligand and receptor, the entropic factors that occur upon binding, the desolvation and solvation energies associated with the interacting molecules. The final free energy of binding (G) will depend on the overall balance of these factors. The interaction forces between the ligand and receptor are Electrostatic, Hydrogen bond, Van der Waals interactions and hydrophobic forces

Docking score is generally calculated by formula

ΔGbinding= ΔGvdW + ΔGelec + ΔGhbond + ΔGdesolv + ΔGtors

where ΔGvdW: 12-6 Lennard-Jones potential; ΔGelec: Coulombic with Solmajer-dielectric; ΔGhbond :12-10 Potential with Goodford Directionality; ΔGdesolv: Stouten Pairwise Atomic Solvation Parameters;

ΔGtors: Number of rotatable bonds

Various Scoring functions include

1. Emperical Scoring

2. Shape and chemical complementary scoring

3. Force field Scoring

4. Knowledge based Scoring

5. Clustering and entropy based Scoring

6. Consensus Scoring.

Challenges in Docking

Docking is currently in a mature stage of development, but it is still far from perfect. Most docking programs available now are normally able to predict known protein bound poses with good accuracies of about 1.5-2 Å with reported success rates in the range of 70–80%. Selecting the docking program that will give the best result for any given target is not straight forward. The most stringent test of docking is the accurate prediction of the binding affinities of a series of related compounds.

Challenges in Docking Docking is currently in a mature stage of development, but it is still far from perfect. Most docking programs available now are normally able to predict known protein bound poses with good accuracies of about 1.5-2 Å with reported success rates in the range of 70–80%. Selecting the docking program that will give the best result for any given target is not straight forward. The most stringent test of docking is the accurate prediction of the binding affinities of a series of related compounds

Ligand representations: Ligand representation and preparation has potential effect on the results of docking as the ligand recognition by a protein depends on shape i.e. 3D structure and electrostatic complementarities. Also ligand conformational sampling is as important as correct ligand preparation. The question of appropriate representation of molecules in databases has been addressed recently [8]. The tautomeric and protomeric states of the small molecules to be docked are user-defined in most of the docking programs. Typically, the structure most likely to be dominant at neutral pH is generated. The structures can be further optimized by removing or adding hydrogens provided approximate pKa values are known a priori.

The accuracy of atom typing is extremely important as the wrong definition of donor and acceptor properties of heteroatoms may lead to serious docking errors. In cases unknown stereochemistry of a synthesized compound, it would be better to generate all possible diastereoisomers of the structure and dock them individually to the receptor. Commercial software programs which use enumeration of all possible diastereoisomers of a given compound include: Stereoplex [9], Stergen [10], and Pipeline Pilot [11]. Accounting various tautomeric and protomeric states of the molecules is challenging at times during docking procedures. Many databases stores molecules such as acids or amines in their neutral forms. They are considered ionized under physiological conditions, so it is necessary to ionize them prior docking. However, while standard ionization is easy to achieve, the problem of generation of tautomer is already much more challenging: which tautomer should one use? Or should one use all possible tautomer’s of less for a given molecule? Not only for tautomers, but also for different ionization states balanced equilibria between the various ionization forms provide real challenges in docking

Baker and co-workers [12] have also discussed the issue of problem involved in sampling tautomeric and protonation states, given the possible difference of free and bound ligand states in these respects. Enumeration of tautomeric and protonation states is a possible solution as suggested by them but have warned about the potentially prohibitive computational cost. Another alternative suggested included incremental and segmentation construction of the docked ligand, whereby the protonation and tautomerism ‘‘decisions’’ are independent and hence decrease the problem size.

Some docking programs (like GOLD and Protein-Ligand ANT System (PLANTS)) have been found having problem in identifying the correct stereoisomer therefore affecting the final outcome. Ten Brink and Exner [13] have recommended that, a preselection of plausible protomers / tautomers should be routinely performed. They have developed Structure Protonation and Recognition System (SPORES) – a tool for preprocessing of protein and ligand-protein complexes and for the setup of 3D ligand databases. Spores perform rule-based assignment of atom types and generate tautomer and protonation as well as stereoisomer states, based on these assignments

Receptor representations: The common source of 3D structures of receptor for docking is X-ray crystallography and NMR (for smaller proteins). However, the growing gap between the sequence and structure availability has to some extent filled by homology modeling, threading, and de novo methods. But the quality of such models for the purposes of docking generally and virtual screening specifically has to be evaluated before use.

The quality of the receptor structure employed plays an important role in determining the success of docking protocol [14-17]. In general, the higher the resolution of the 3D crystal structure employed for studies, the better the observed docking results. The reproducibility of the program increases with increase in resolution of the co-crystal structures (less than 2.0 Å) employed for study [16]. A recent review of the accuracy, limitations and challenges of the structure refinement protocols of protein ligand complexes in general provided a critical assessment of the available structures [7]. Regardless of the possible ambiguities, success has been reported for large number of high throughput docking studies using X-ray receptor structures. Recent examples of this type of study include: kinesin [17], HIV protease [18], phosphoribosyl transferase [19], FKBP12 [20], farnesyl transferase [21], beta-lactamase [22], and PTP1B [23].

Ligand binding to protein often results in conformational changes of protein, so ignoring protein flexibility during molecular docking may give results that are incorrect [24].

One of the major challenges faced in the field of docking is handling of flexible protein receptor. A protein can adopts different conformations depending upon the ligand to which it binds. As a result, docking performed using a rigid receptor will correspond to a single receptor conformation. However, certain ligands require different receptor conformation in order to bind, where we need to keep receptor flexible. Proteins exists in constant motion between different conformational states having similar energies, which is usually neglected in docking studies, although it is known that protein flexibility accounts for increased affinity to be achieved between a given drug and its target. The number of degrees of freedom included in the conformational search is an important aspect that determines the searching efficiency.

A biological system usually consists of a ligand, the macromolecular receptor and solvent molecules. Large numbers of degrees of freedom are associated with the solvent molecules, which is normally excluded from the problem and sometimes they are implicitly considered in the scoring functions to understand the solvent effect. Rest of the degrees of freedom involved with ligand and receptor, can be reduced through the use of different approximations, allowing the search space to be more effectively sampled.

Some approaches which have been proposed to deal with flexible receptor include:

•  Letting the receptor or parts thereof move during docking
• Docking the compounds into several different  conformations of the same receptor and aggregating the results
• Docking into averaged receptor representations. 

Receptor flexibility has been extensively treated in some software’s like by Monte Carlo (MC) simulations and rotamer libraries Rosetta Ligand offers one of the most extensive receptor flexibility treatments developed to date [12,25]. The binding site side chain rotamers are optimized using a simulated annealing procedure and the backbone is minimized subject to restraints. Auto Dock 4 also fully models the flexibility of selected portions of the protein[26]. Side chains of interest are separated from the protein and treated precisely during the simulation, allowing rotation around torsional degrees of freedom. Induced Fit Docking (IFD) Workflow [27, 28] in Schrodinger involves rigid receptor docking with Glide [29, 30], combined with minimization of protein-ligand complex with the homology modeling module Prime [28]. In the MADAMM procedure, the protein is made flexibile by the side chain rotamer libraries of the Insight II [31]. IFD has been successfully used for studies of HIV-1 Integrase [32], kinases [33,34], monoacylglycerol lipase [35], heat shock protein 90[36].

Receptor ensembles by molecular dynamics (MD) have also been widely used to handle the problem of receptor flexibility [37,38]. A 4D-docking protocol for Internal Coordinate Mechanics (ICM), has been developed by Abagyan and Totrov where the fourth dimension is receptor conformation[39-41]. In this protocol, multiple receptor conformations are represented by multiple grids and each is represented as a variable in the global optimization. By this approach increased accuracy with no loss in effectiveness compared to single grid methods has been achieved.

Active site water molecules are another important aspect of target flexibility. Water molecules should be checked carefully to avoid using artifact waters (i.e. water molecules not essential to the protein structure) in the docking process. Using artifact active site water molecules can have a deleterious effect by providing false energetic stability to the protein-ligand complex.

Receptors bind to their ligands in solution and the solvation aspects are commonly treated implicitly, that is, by the use of implicit solvents, knowledge-based scoring functions or by modification or calibration of other scoring functions. Cincilla et al. [42], modified the solvation treatment in the scoring function of Auto Dock 3 [43] to enhance the interactions/predictions of weak complexes containing ligands with polar atoms in the binding site. Kuntz and co-workers [44] have used two implicit solvent scoring functions AMBER/GBSA and AMBER/PBSA, implemented in DOCK 6, for docking small molecules to RNA. To neutralize the backbone charge sodium ions have been used and to shield the charges a double shell of explicit water. They have found that the quality of pose prediction increased from 70% to 80% for moderately flexible ligands (less than 7 rotatable bonds) and from 26% to 50–60% for highly flexible ligands (7–13 rotatable bonds). The effect of structural water molecules on docking as described in the literature suggested that explicit water molecules improve docking outcomes, both in hit identification and pose prediction in virtual screening [45]. Englebienne and Moitessier [46] have shown that the consideration of displaceable water molecules, implemented in FITTED, improves pose prediction, but does not significantly affect scoring accuracy and suggested that the latter may be the outcome of most scoring functions that having been developed for ‘‘dry’’ proteins

An imperfection in scoring function is another challenge in docking. Since search algorithm is capable of generating right conformations, scoring function should also be able to distinct the true binding modes from all the other alternative modes. A very rigorous scoring function would be computationally too expensive, unfeasible for analyzing several binding modes. Scoring functions make number of simplifications and assumptions for ligand affinity evaluation, but at the cost of accuracy. Certain physical phenomenon such as entropy and electrostatic interactions are disregarded in contemporary scoring schemes. Lack of a suitable scoring function, both in terms of accuracy and speed, is the major bottleneck in docking algorithm.

In 2009 several scientists have explored the potential of QM/MM scoring [47-53]. Fong et al. [47] have tested three functions (AM1d, HF/6-31G* , and PM3) for ligand treatment in combination with GoldScore, AMBER, and Chem Score for prediction of successful poses of six HIV protease inhibitors. Gleeson and Gleeson [48] have used the combination of 6-31G**/ B3LYP and Universal Force Field (UFF) for successful crossdocking and re-scoring of nine kinase ligands. A combination of QPLD, a QM/MM docking program, with Site Map [54] a binding site classification module is explored by Chung et al. [49]. They have used 455 protein-ligand complexes and demonstrated a scoring improvement, over Glide, for three possible binding site types (hydrophobic, hydrophilic, and metalloproteins). Cho and co-workers have tested QM/MM scoring function for different types of binding sites, namely for those with hydrophobic groups, polar groups, and metalloproteins [50,53]

Many scoring functions perform very well for the purpose of pose prediction, but a goal of predicting binding affinities using scoring functions is still unfulfilled. It is also clear from the literature and advances reported in scoring functions that the currently available functions could be used more efficiently in combination [55]. Commonly, a consensus scoring involves multiple rescoring of a docked pose with different scoring functions or a combination thereof [56]

Table 1: Non-exhaustive list of docking programs available and their basic characteristics.

Supported Platform Search Algorithm Scoring Algorithm Protein 
Dock1 Windows, Mac, Sun, 
SGI, Unix, Linux, IBM 
Shape fitting GB/SA Solvation scoring, Chemscore, 
Bump filters, Contact score, Grid 
based score, DOCK score
P Improved prediction of 
binding poses by Forcefield scoring
Auto Dock 57 Mac, OSX, Unix, Linux, 
Lamarckian genetic, 
Genetic and Simulated 
Auto Dock Foce-field methods P Improved prediction of 
binding poses by Forcefield scoring
FlexX58 Windows, Linux, Unix, 
SGI, Sun
Screen Score, Drug Score, PLP, FlexX 
NP Scored is on the basis of 
protein-ligand interaction
FRED 59 Windows, Unix, Linux, 
Mac, SGI, OSX, IBM, 
Shape fitting 
Gaussian shape score, PL, Screen 
score, User defined
  Fastest docking tool, 
suitable for ultrahighthroughput docking
GOLD60 Sun, SGI, Linux, IBM, 
Genetic algorithm ChemScore, GoldScore, User defined 
P GOLD has very high 
docking accuracy.
Glide61 IBM, SGI, Linux, Unix, 
MC sampling GlideScore, GlideComp P High docking accuracy but 
very slow, not good for 
large data sets
LigandFit62 SGI, Linux, AIX- IBM MC sampling LigScore, PLP, PMF P Very fast, suitable for 
virtual high or ultra HTS.
ICM63 Windows, Unix, Mac Stochastic algorithms
(MC techniques
adapted in flexible 
docking algorithms
ICM scoring function P  
ProDock64   Monte Carlo sampling AMBER force-field P Bezier spline energy grid 
has been incorporated to 
speed up optimization 
QXP65   Monte Carlo sampling QXP force-field P Uses consensus scoring and 
ranking protocol for highly 
flexible ligands
Surflex66 Linux, Windows and 
construction search
Hammerhead’s empirical scoring 
P Suitable for virtual HTS
Abbreviations: ICM: Internal Coordinate Mechanics; MC: Monte Carlo; SGI: Silicon Graphics Inc; AIX: Advanced Interactive executive: P: Provided; NP: 
Not Provided; HTS: High-Throughput Screening



It is evident from docking literature, that it has attained a good amount of maturity but accounting for flexibility and successful scoring remain significant challenges. Nevertheless important advances are being made in all aspects of docking programs. The selection of best docking technique should be done after thoroughly studying the target, ligands, and docking method performance. The issue of ligand flexibility seems more or less resolved and does not create much problem however protein flexibility needs more attention and improvement. Careful analysis of active site crystal water molecules is required is the key in good docking. Inclusion of water molecules should be considered after studying the hydrogen bonding with nonwater residues and after studying the relative abundance of water molecules by analysis of multiple crystal structures. In general using a single input conformation and performing a single docking run is likely to decrease docking accuracy, having multiple ligand input does not guarantee said accuracy.


The author acknowledges research associate (Dr. Munnazah Tasleem) and research scholar (Apeksha Shrivastava) for their help in compilation of the material


1. Kuntz ID, Blaney JM, Oatley SJ, Langridge R, Ferrin TE. A geometric approach to macromolecule-ligand interactions. J Mol Biol. 1982; 16: 269-288.

2. Janin J, Cherfils J. Protein docking algorithms: simulating molecular recognition. Curr Opin Struct Biol. 1993; 3: 265-269.

3. Hawkins PC, Skillman AG, Warren GL, Ellingson BA, Stahl MT. Conformer Generation with OMEGA: Algorithm and Validation Using High Quality Structures from the Protein Databank and the Cambridge Structural Database. J Chem Inf Model. 2010; 50: 572-584.

4. Sadowski J, Gasteiger J, Klebe G. Comparison of Automatic ThreeDimensional Model Builders Using 639 X-Ray Structures. J Chem Inf Comp Sci. 1994; 34: 1000-1008.

5. Kirchmair J, Wolber G, Laggner C, Langer T. Comparative Performance Assessment of the Conformational Model Generators Omega and Catalyst: A Large-Scale Survey on the Retrieval of Protein-Bound Ligand Conformations. J Chem Inf Model. 2006; 46: 1848-1861.

6. Lodhi SS, Farmer R, Singh AK, Jaiswal YK, Wadhwa G. 3D structure generation, virtual screening and docking of human Ras-associated binding (Rab3A) protein involved in tumourigenesis. Mol Biol Rep. 2014; 41: 3951-3959.

7. Colotta V, Catarzi D, Varano F, Lenzi O, Filacchioni G, Costagli C, et al. Structural investigation of the 7-chloro-3-hydroxy-1H-quinazoline2,4-dione scaffold to obtain AMPA and kainate receptor selective antagonists. Synthesis, pharmacological, and molecular modeling studies. J Med Chem. 2006; 49: 6015-6026.

8. Elokely KM, Doerksen RJ. Docking Challenge: Protein Sampling and Molecular Docking Performance. J Chem Inf Model. 2013; 53: 1934.

9. Tripos: St. Louis. MO. USA.

10. Molecular Networks GmbH. Erlangen. Germany.

11. SciTegic: San Diego. CA. USA.

12. Davis IW, Raha K, Head MS, Baker D. Blind docking of pharmaceutically relevant compounds using Rosetta Ligand. Protein Sci. 2009; 18: 1998-2002.

13. Ten Brink T, Exner TE. Influence of protonation, tautomeric, and stereoisomeric states on protein-ligand docking results. J Chem Inf Model. 2009; 49: 1535-1546.

14. Muegge I. The effect of small changes in protein structure on predicted binding modes of known inhibitors of influenza virus neuraminidase: PMF-scoring in DOCK4. Med Chem Res. 1999; 9: 490-500.

15. McGovern SL, Shoichet BK. Information decay in molecular docking screens against holo, apo, and modeled conformations of enzymes. J Med Chem. 2003; 46: 2895-2907.

16. Schapira M, Abagyan R, Totrov M. Nuclear hormone receptor targeted virtual screening. J Med Chem. 2003; 46: 3045-59.

17. Hopkins SC, Vale RD, Kuntz ID. Inhibitors of kinesin activity from structure-based computer screening. Biochemistry. 2000; 39: 2805- 2814.

18. Olson AJ, Goodsell DS. Automated docking and the search for HIV protease inhibitors. SAR QSAR Environ Res. 1998; 8: 273-285.

19. Aronov AM, Munagala NR, Ortiz De Montellano PR, Kuntz ID, Wang CC. Rational design of selective submicromolar inhibitors of Tritrichomonas foetus hypoxanthine-guanine-xanthine phosphoribosyl transferase. Biochemistry. 2000; 39: 4684-4691.

20. Burkhard P, Hommel U, Sanner M, Walkinshaw MD. The discovery of steroids and other novel FKBP inhibitors using a molecular docking program. J Mol Biol. 1999; 287: 853-858.

21. Perola E, Xu K, Kollmeyer TM, Kaufmann SH, Prendergast FG, Pang YP. Successful virtual screening of a chemical database for farnesyl transferase inhibitor leads. J Med Chem. 2000; 43: 401-408.

22. Powers RA, Morandi F, Shoichet BK. Structure-based discovery of a novel, noncovalent inhibitor of AmpC beta-lactamase. Structure. 2002; 10: 1013-1023.

23. Doman TN, McGovern SL, Witherbee BJ, Kasten TP, Kurumbail R, Stallings WC, et al. Molecular docking and high-throughput screening for novel inhibitors of protein tyrosine phosphatase-1B. J Med Chem. 2002; 45: 2213-2221.

24. Totrov M. Abagyan R. Flexible ligand docking to multiple receptor conformations: A practical alternative. Curr Opin Struct Biol. 2008; 18: 178-184.

25. Davis IW, Baker D. RosettaLigand docking with full ligand and receptor flexibility. J Mol Biol. 2009; 385: 381-392.

26. Morris GM, Huey R, Lindstrom W, Sanner MF, Belew RK, Goodsell DS, Olson AJ. AutoDock4 and AutoDockTools4: automated docking with selective receptor flexibility. J Comput Chem. 2009; 30: 2785-2791.

27. Sherman W, Beard HS, Farid R. Use of an induced fit receptor structure in virtual screening. Chem Biol Drug Des. 2006; 67: 83-84.

28. Sherman W, Day T, Jacobson MP, Friesner RA, Farid R. Novel procedure for modeling ligand/receptor induced fit effects. J Med Chem. 2006; 49: 534-553.

29. Friesner RA, Banks JL, Murphy RB, Halgren TA, Klicic JJ, Mainz DT, et al. Glide: a new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy. J Med Chem. 2004; 47: 1739-1749.

30. Halgren TA, Murphy RB, Friesner RA, Beard HS, Frye LL, Pollard WT, et al. Glide: a new approach for rapid, accurate docking and scoring. 2. Enrichment factors in database screening. J Med Chem. 2004; 47: 1750-1759.

31. Cerqueira NM, Bras NF, Fernandes PA, Ramos MJ. MADAMM: a multistaged docking with an automated molecular modeling protocol. Proteins. 2009; 74: 192-206.

32. Barreca ML, Iraci N, De Luca L, Chimirri A. Induced-fit docking approach provides insight into the binding mode and mechanism of action of HIV-1 integrase inhibitors. Chem Med Chem. 2009; 4: 1446- 1456.

33. Zhong H, Tran LM, Stang JL. Induced-fit docking studies of the active and inactive states of protein tyrosine kinases. J Mol Graph Model. 2009; 28: 336-346.

34. Cho YY, Yao K, Pugliese A, Malakhova ML, Bode AM, Dong Z. A regulatory mechanism for RSK2 NH(2)-terminal kinase activity. Cancer Res. 2009; 69: 4398-4406.

35. King AR, Dotsey EY, Lodola A, Jung KM, Ghomian A, Qiu Y, et al. Discovery of potent and reversible monoacylglycerol lipase inhibitors. Chem Biol. 2009; 16: 1045-1052.

36. Lauria A, Ippolito M, Almerico AM. Inside the Hsp90 inhibitors binding mode through induced fit docking. J Mol Graph Model. 2009; 27: 712- 722.

37. Huang Z, Wong CF. Conformational selection of protein kinase A revealed by flexible-ligand flexible-protein docking. J Comput Chem. 2009; 30: 631-644.

38. Armen RS, Chen JH, Brooks CL. An evaluation of explicit receptor flexibility in molecular docking using molecular dynamics and torsion angle molecular dynamics. J Chem Theory Comput. 2009; 5: 2909- 2923.

39. Abagyan R, Totrov M. Biased probability Monte Carlo conformational searches and electrostatic calculations for peptides and proteins. J Mol Biol. 1994; 235: 983-1002.

40. Bottegoni G, Kufareva I, Totrov M, Abagyan R. Fourdimensional docking: a fast and accurate account of discrete receptor flexibility in ligand docking. J Med Chem. 2009; 52: 397-406.

41. Abagyan R, Kufareva I. The flexible pocketome engine for structural chemogenomics. Methods Mol Biol. 2009; 575: 249-279.

42. Cincilla G, Vidal D, Pons M. An improved scoring function for suboptimal polar ligand complexes. J Comput Aided Mol Des. 2009; 23: 143-152.

43. Morris GM, Goodsell DS, Halliday RS, Huey R, Hart WE, Belew RK, et al. Automated docking using a Lamarckian genetic algorithm and an empirical binding free energy function. J Comput Chem. 1998; 19: 1639-1662.

44. Lang PT, Brozell SR, Mukherjee S, Pettersen EF,Meng EC, Thomas V, Rizzo RC, Case DA, James TL, Kuntz ID. DOCK 6: combining techniques to model RNA-small molecule complexes. RNA. 2009; 15: 1219-1230.

45. Villacanas O, Madurga S, Giralt E, Belda I. Explicit treatment of water molecules in protein-ligand docking. Curr Comp Aided Drug Des. 2009; 5: 145-154.

46. Englebienne P, Moitessier N. Docking ligands into flexible and solvated macromolecules. 4. Are popular scoring functions accurate for this class of proteins? J Chem Inf Model. 2009; 49: 1568-1580.

47. Fong P, McNamara JP, Hillier IH, Bryce RA. Assessment of QM/ MM scoring functions for molecular docking to HIV-1 protease. J Chem Inf Model. 2009; 49: 913-924.

48. Gleeson MP, Gleeson D. QM/MM as a tool in fragment based drug discovery. A cross-docking, rescoring study of kinase inhibitors. J Chem Inf Model. 2009; 49: 1437-1448.

49. Chung JY, Hah JM, Cho AE. Correlation between performance of QM/ MM docking and simple classification of binding sites. J Chem Inf Model. 2009; 49: 2382-2387.

50. Cho AE, Rinaldo D. Extension of QM/MM docking and its applications to metalloproteins. J Comput Chem. 2009; 30: 2609-2616.

51. Anisimov VM, Bugaenko VL. QM/QM docking method based on the variational finite localized molecular orbital approximation. J Comput Chem. 2009; 30: 784-798.

52. Sgrignani J, Bonaccini C, Grazioso G, Chioccioli M, Cavalli A, Gratteri P. Insights into docking and scoring neuronal alpha4beta2 nicotinic receptor agonists using molecular dynamics simulations and QM/MM calculations. J Comput Chem. 2009; 30: 2443-2454.

53. Cho AE, Chung JY, Kim M, Park K. Quantum mechanical scoring for protein docking. J Chem Phys. 2009; 131: 134108.

54. Halgren T. New method for fast and accurate binding-site identification and analysis. Chem Biol Drug Des. 2007; 69: 146-148.

55. Yurieva E, Agostinoa M, Ramsland PA. Challenges and advances in computational docking. J Mol Recognit. 2011; 24: 149-164.

56. Charifson PS, Corkery JJ, Murcko MA, Walters WP. Consensus scoring: a method for obtaining improved hit rates from docking databases of three-dimensional structures into proteins. J Med Chem. 1999; 42: 5100-5109.

57. Goodsell DS, Olson AJ. Automated docking of substrates to proteins by simulated annealing. Proteins 1990; 8: 195-202.

58. Rarey M, Kramer B, Lengauer T. Multiple automatic base selection: protein-ligand docking based on incremental construction without manual intervention. J Comput Aided Mol Des. 1997; 11: 369-384.

59. Schulz-Gasch, Stahl M. Binding Site Characterstics in Structure-based Virtual Screening: Evaluation of Current Docking Tools. J Mol Model. 2003; 9: 47-57.

60. Jones G, Wilett P, Glein RC, Leach AR, Taylor R. Development and Validation of Genetic Algorithm and an Empirical Binding Free Energy Function. J Mol Biol. 1997: 267: 727-748.

61. Friesner RA, Banks JL, Murphy RB, Halgren TA, Klicic JJ, Mainz DT, et al. Glide: A New Approach for Rapid, Accurate Docking and Scoring. 1. Method and Assessment of Docking Acuracy. J Med Chem. 2004; 47: 1739-1749.

62. Venkatachalam CM, Jiang X, Oldfield T, Waldan M. LigandFit: A Novel Method for the Shape-directed Rapid Docking of Ligands to Protein Active Site. J Mol Graphics Modell. 2003; 21: 289-307.

63. Abagyan RA, Totrov MM, Kuznetsov DA. ICM: A New Method For Protein Modeling and Design: Applications To Docking and Structure Prediction From The Distorted Native Conformation. J Comp Chem. 1994; 15: 488-506.

64. Trosset JY, Scheraga HA. PRODOCK: Software Package for Protein Modeling and Docking. J Comput Chem. 1999; 20: 412-427.

65. McMartin C, Bohacek RS. QXP: Powerful, Rapid Computer Algorithms for Structure-based Drug Design. J Comput Aid Mol Des. 1997; 11: 333-344.

66. Jain AN. Surflex: Fully Automatic Flexible Molecular Docking Using a Molecular Similarity-Based Search Engine. J Med Chem. 2003; 46: 499-451

Received : 13 May 2016
Accepted : 12 Jul 2016
Published : 14 Jul 2016
Annals of Otolaryngology and Rhinology
ISSN : 2379-948X
Launched : 2014
JSM Schizophrenia
Launched : 2016
Journal of Nausea
Launched : 2020
JSM Internal Medicine
Launched : 2016
JSM Hepatitis
Launched : 2016
JSM Oro Facial Surgeries
ISSN : 2578-3211
Launched : 2016
Journal of Human Nutrition and Food Science
ISSN : 2333-6706
Launched : 2013
JSM Regenerative Medicine and Bioengineering
ISSN : 2379-0490
Launched : 2013
JSM Spine
ISSN : 2578-3181
Launched : 2016
Archives of Palliative Care
ISSN : 2573-1165
Launched : 2016
JSM Nutritional Disorders
ISSN : 2578-3203
Launched : 2017
Annals of Neurodegenerative Disorders
ISSN : 2476-2032
Launched : 2016
Journal of Fever
ISSN : 2641-7782
Launched : 2017
JSM Bone Marrow Research
ISSN : 2578-3351
Launched : 2016
JSM Mathematics and Statistics
ISSN : 2578-3173
Launched : 2014
Journal of Autoimmunity and Research
ISSN : 2573-1173
Launched : 2014
JSM Arthritis
ISSN : 2475-9155
Launched : 2016
JSM Head and Neck Cancer-Cases and Reviews
ISSN : 2573-1610
Launched : 2016
JSM General Surgery Cases and Images
ISSN : 2573-1564
Launched : 2016
JSM Anatomy and Physiology
ISSN : 2573-1262
Launched : 2016
JSM Dental Surgery
ISSN : 2573-1548
Launched : 2016
Annals of Emergency Surgery
ISSN : 2573-1017
Launched : 2016
Annals of Mens Health and Wellness
ISSN : 2641-7707
Launched : 2017
Journal of Preventive Medicine and Health Care
ISSN : 2576-0084
Launched : 2018
Journal of Chronic Diseases and Management
ISSN : 2573-1300
Launched : 2016
Annals of Vaccines and Immunization
ISSN : 2378-9379
Launched : 2014
JSM Heart Surgery Cases and Images
ISSN : 2578-3157
Launched : 2016
Annals of Reproductive Medicine and Treatment
ISSN : 2573-1092
Launched : 2016
JSM Brain Science
ISSN : 2573-1289
Launched : 2016
JSM Biomarkers
ISSN : 2578-3815
Launched : 2014
JSM Biology
ISSN : 2475-9392
Launched : 2016
Archives of Stem Cell and Research
ISSN : 2578-3580
Launched : 2014
Annals of Clinical and Medical Microbiology
ISSN : 2578-3629
Launched : 2014
JSM Pediatric Surgery
ISSN : 2578-3149
Launched : 2017
Journal of Memory Disorder and Rehabilitation
ISSN : 2578-319X
Launched : 2016
JSM Tropical Medicine and Research
ISSN : 2578-3165
Launched : 2016
JSM Head and Face Medicine
ISSN : 2578-3793
Launched : 2016
JSM Cardiothoracic Surgery
ISSN : 2573-1297
Launched : 2016
JSM Bone and Joint Diseases
ISSN : 2578-3351
Launched : 2017
JSM Bioavailability and Bioequivalence
ISSN : 2641-7812
Launched : 2017
JSM Atherosclerosis
ISSN : 2573-1270
Launched : 2016
Journal of Genitourinary Disorders
ISSN : 2641-7790
Launched : 2017
Journal of Fractures and Sprains
ISSN : 2578-3831
Launched : 2016
Journal of Autism and Epilepsy
ISSN : 2641-7774
Launched : 2016
Annals of Marine Biology and Research
ISSN : 2573-105X
Launched : 2014
JSM Health Education & Primary Health Care
ISSN : 2578-3777
Launched : 2016
JSM Communication Disorders
ISSN : 2578-3807
Launched : 2016
Annals of Musculoskeletal Disorders
ISSN : 2578-3599
Launched : 2016
Annals of Virology and Research
ISSN : 2573-1122
Launched : 2014
JSM Renal Medicine
ISSN : 2573-1637
Launched : 2016
Journal of Muscle Health
ISSN : 2578-3823
Launched : 2016
JSM Genetics and Genomics
ISSN : 2334-1823
Launched : 2013
JSM Anxiety and Depression
ISSN : 2475-9139
Launched : 2016
Clinical Journal of Heart Diseases
ISSN : 2641-7766
Launched : 2016
Annals of Medicinal Chemistry and Research
ISSN : 2378-9336
Launched : 2014
JSM Pain and Management
ISSN : 2578-3378
Launched : 2016
JSM Women's Health
ISSN : 2578-3696
Launched : 2016
Clinical Research in HIV or AIDS
ISSN : 2374-0094
Launched : 2013
Journal of Endocrinology, Diabetes and Obesity
ISSN : 2333-6692
Launched : 2013
Journal of Substance Abuse and Alcoholism
ISSN : 2373-9363
Launched : 2013
JSM Neurosurgery and Spine
ISSN : 2373-9479
Launched : 2013
Journal of Liver and Clinical Research
ISSN : 2379-0830
Launched : 2014
Journal of Drug Design and Research
ISSN : 2379-089X
Launched : 2014
JSM Clinical Oncology and Research
ISSN : 2373-938X
Launched : 2013
JSM Bioinformatics, Genomics and Proteomics
ISSN : 2576-1102
Launched : 2014
Journal of Trauma and Care
ISSN : 2573-1246
Launched : 2014
JSM Surgical Oncology and Research
ISSN : 2578-3688
Launched : 2016
Annals of Food Processing and Preservation
ISSN : 2573-1033
Launched : 2016
Journal of Radiology and Radiation Therapy
ISSN : 2333-7095
Launched : 2013
JSM Physical Medicine and Rehabilitation
ISSN : 2578-3572
Launched : 2016
Annals of Clinical Pathology
ISSN : 2373-9282
Launched : 2013
Annals of Cardiovascular Diseases
ISSN : 2641-7731
Launched : 2016
Journal of Behavior
ISSN : 2576-0076
Launched : 2016
Annals of Clinical and Experimental Metabolism
ISSN : 2572-2492
Launched : 2016
Clinical Research in Infectious Diseases
ISSN : 2379-0636
Launched : 2013
JSM Microbiology
ISSN : 2333-6455
Launched : 2013
Journal of Urology and Research
ISSN : 2379-951X
Launched : 2014
Journal of Family Medicine and Community Health
ISSN : 2379-0547
Launched : 2013
Annals of Pregnancy and Care
ISSN : 2578-336X
Launched : 2017
JSM Cell and Developmental Biology
ISSN : 2379-061X
Launched : 2013
Annals of Aquaculture and Research
ISSN : 2379-0881
Launched : 2014
Clinical Research in Pulmonology
ISSN : 2333-6625
Launched : 2013
Journal of Immunology and Clinical Research
ISSN : 2333-6714
Launched : 2013
Annals of Forensic Research and Analysis
ISSN : 2378-9476
Launched : 2014
JSM Biochemistry and Molecular Biology
ISSN : 2333-7109
Launched : 2013
Annals of Breast Cancer Research
ISSN : 2641-7685
Launched : 2016
Annals of Gerontology and Geriatric Research
ISSN : 2378-9409
Launched : 2014
Journal of Sleep Medicine and Disorders
ISSN : 2379-0822
Launched : 2014
JSM Burns and Trauma
ISSN : 2475-9406
Launched : 2016
Chemical Engineering and Process Techniques
ISSN : 2333-6633
Launched : 2013
Annals of Clinical Cytology and Pathology
ISSN : 2475-9430
Launched : 2014
JSM Allergy and Asthma
ISSN : 2573-1254
Launched : 2016
Journal of Neurological Disorders and Stroke
ISSN : 2334-2307
Launched : 2013
Annals of Sports Medicine and Research
ISSN : 2379-0571
Launched : 2014
JSM Sexual Medicine
ISSN : 2578-3718
Launched : 2016
Annals of Vascular Medicine and Research
ISSN : 2378-9344
Launched : 2014
JSM Biotechnology and Biomedical Engineering
ISSN : 2333-7117
Launched : 2013
Journal of Hematology and Transfusion
ISSN : 2333-6684
Launched : 2013
JSM Environmental Science and Ecology
ISSN : 2333-7141
Launched : 2013
Journal of Cardiology and Clinical Research
ISSN : 2333-6676
Launched : 2013
JSM Nanotechnology and Nanomedicine
ISSN : 2334-1815
Launched : 2013
Journal of Ear, Nose and Throat Disorders
ISSN : 2475-9473
Launched : 2016
JSM Ophthalmology
ISSN : 2333-6447
Launched : 2013
Journal of Pharmacology and Clinical Toxicology
ISSN : 2333-7079
Launched : 2013
Annals of Psychiatry and Mental Health
ISSN : 2374-0124
Launched : 2013
Medical Journal of Obstetrics and Gynecology
ISSN : 2333-6439
Launched : 2013
Annals of Pediatrics and Child Health
ISSN : 2373-9312
Launched : 2013
JSM Clinical Pharmaceutics
ISSN : 2379-9498
Launched : 2014
JSM Foot and Ankle
ISSN : 2475-9112
Launched : 2016
JSM Alzheimer's Disease and Related Dementia
ISSN : 2378-9565
Launched : 2014
Journal of Addiction Medicine and Therapy
ISSN : 2333-665X
Launched : 2013
Journal of Veterinary Medicine and Research
ISSN : 2378-931X
Launched : 2013
Annals of Public Health and Research
ISSN : 2378-9328
Launched : 2014
Annals of Orthopedics and Rheumatology
ISSN : 2373-9290
Launched : 2013
Journal of Clinical Nephrology and Research
ISSN : 2379-0652
Launched : 2014
Annals of Community Medicine and Practice
ISSN : 2475-9465
Launched : 2014
Annals of Biometrics and Biostatistics
ISSN : 2374-0116
Launched : 2013
JSM Clinical Case Reports
ISSN : 2373-9819
Launched : 2013
Journal of Cancer Biology and Research
ISSN : 2373-9436
Launched : 2013
Journal of Surgery and Transplantation Science
ISSN : 2379-0911
Launched : 2013
Journal of Dermatology and Clinical Research
ISSN : 2373-9371
Launched : 2013
JSM Gastroenterology and Hepatology
ISSN : 2373-9487
Launched : 2013
Annals of Nursing and Practice
ISSN : 2379-9501
Launched : 2014
JSM Dentistry
ISSN : 2333-7133
Launched : 2013
Author Information X