Computational Methods Applied to Rational Drug Design



David Ramírez*
Centro de Bioinformática y Simulación Molecular, Universidad de Talca, 2 Norte 685, Casilla, Talca, Chile


Article Metrics

CrossRef Citations:
0
Total Statistics:

Full-Text HTML Views: 927
Abstract HTML Views: 457
PDF Downloads: 175
ePub Downloads: 96
Total Views/Downloads: 1655
Unique Statistics:

Full-Text HTML Views: 485
Abstract HTML Views: 281
PDF Downloads: 142
ePub Downloads: 80
Total Views/Downloads: 988



© David Ramírez; Licensee Bentham Open.

open-access license: This is an open access article licensed under the terms of the Creative Commons Attribution-Non-Commercial 4.0 International Public License (CC BY-NC 4.0) (https://creativecommons.org/licenses/by-nc/4.0/legalcode), which permits unrestricted, non-commercial use, distribution and reproduction in any medium, provided the work is properly cited.

Correspondence: *Address correspondence to this author at the Centro de Bioinformática y Simulación Molecular, Universidad de Talca, 2 Norte 685, Casilla 721, Talca, Chile; Tel: +56972983600; E-mail: davramirez@utalca.cl


Abstract

Due to the synergic relationship between medical chemistry, bioinformatics and molecular simulation, the development of new accurate computational tools for small molecules drug design has been rising over the last years. The main result is the increased number of publications where computational techniques such as molecular docking, de novo design as well as virtual screening have been used to estimate the binding mode, site and energy of novel small molecules. In this work I review some tools, which enable the study of biological systems at the atomistic level, providing relevant information and thereby, enhancing the process of rational drug design.

Keywords: Bioinformatics, de novo design, medical chemistry, molecular docking, molecular modeling, rational drug design, virtual screening.



INTRODUCTION

Synergic relationship between experimental biochemistry sciences and computational theoretical methods has increased in recent years, providing significant benefits in different biological sciences areas, specially biochemistry and molecular biology. It has also brought advantages in medicinal chemistry for the rational design of drugs [1]. Regarding theoretical study of biochemical systems, one of its greatest profits is the molecular knowledge of the currently studied structures; these studies allow an exhaustive system analysis, and find meaningful answers to research questions that scientists face in their experimental laboratories.

Before exploring some of the computational methodologies currently used in drug design, it is important to consider that researchers have to draw from relevant information such as the three-dimensional structure of one or several substrates and molecular targets in order to start the theoretical study of certain biological systems. This information is essential, and in many cases, it determines the success or failure of the ongoing theoretical study. For instance, if a researcher wants to study the system formed by an endogenous or an exogenous inhibitor with its respective molecular target, the three-dimensional structure of the molecules involved in the interaction to be simulated has to be known.

The number of molecular targets with a well known structure is exponentially increasing due to the dramatic progress of spectroscopic techniques such as X-ray crystallography, Nuclear Magnetic Resonance (NMR) [2], and the development of the super-resolved fluorescence microscopy that shows a 3D image of a single molecule [3]. In addition, structures which have been produced from structural genomics have also become valuable tools for the study of systems whose molecular targets have not been solved yet [4, 5]. This growth has allowed a massive and constant use of computational tools in research centers worldwide. Furthermore, the developed methodologies for the production and optimization of small molecules (as most of the substrates studied in biochemistry sciences), are already registered in special databases and have provided vital information to research and find inherent characteristics which have never been described before [6, 7].

Nowadays, there are multiple computational methodologies used as bioinformatics tools for the study of biological

systems and drug discovery. The use of one or another tool is linked to the researcher interests, the calculation level intended, possible technical limitations, how and also which kind of information can be extracted from data processing and analysis. Some of the main methodologies for these studies are: Molecular Docking, de novo design, Virtual Screening (VS), Quantitative structure-activity relationships (QSAR), Molecular Dynamics simulation (MDs) and Molecular Modeling (MM) [8]. These methodologies are wieldy used for rational drug design and discovery processes, where the ligand receptor binding mode is imperative to understand the molecular interaction mechanism and the structural factors related with the bioactivity of each inhibitor in detail [9].

In this review, basic and significant elements of molecular docking, de novo design and virtual screening are shown. These are some of the greatest computational methodologies used for the analyses of biochemical sciences, which also improve the common techniques used to study proteins structures, the design of new molecules biologically relevant for pharmacological application, the structure activity relationship, among others.

Molecular Docking

Molecular Docking dates from the middle of the 19th century. Chemists Archibald Scott Couper, Friedrich August Kekulé and Aleksandr Mikhailovich introduced the valence in organic chemistry, and submitted the first structures with graphical representations of carbon atoms [10]. In 1861 Johann Josef Loschmidt produced the largest molecules collection at that time (368), including the first accurate benzene structure [11]. Due to the technological development in computational sciences, the concept of ‘force field’ was introduced from the vibrational spectroscopy, which consider the forces acting all over the atoms in a molecule. The scientific community did not adopt the concept of ‘force field’ until 1946, when the first premise about the combination of steric interactions and the Newtonian mechanics model of vibrational modes, bonds, and angles emerged.

Many other scientific advances were developed at that time, but it was in 1953 when Metropolis et al. published the study: "Equation of State Calculations by Fast Computing Machines" [12]. For the last century, the molecular modeling foundations were set, and due to the technological advances of our time, it now allows us to perform better theoretical studies and computational calculations and validates the need of using computational methods for science progress.

Molecular docking main goal is to predict as accurately as possible the best conformations (poses) of a ligand - by using a score function - in a conformational area, which is delimited by the molecular target binding site [13]. This is applied in different stages of the drug design process in order to predict the binding mode of the ligands already known [14]; it is also used to identify novel and potent ligands [15] and as a binding affinity predictive tool [16]. The first molecular docking algorithms considered translation and orientation only (degrees of freedom), and used the ligand and protein as rigid bodies. Improvements in computer science allowed the creation of new and accurate algorithms in order to use the ligand as flexible with a rigid (or nearly rigid) receptor [17].

This methodology is described as a multi-step process, where each step has one or several complexity levels. The process starts with the implementation of the molecular docking algorithms responsible for placing the ligands in the binding sites. This process is challenging, since it depends directly on the freedom levels granted to the ligands. The freedom levels sampling must be done accurately in order to identify the best conformation for the receptor-binding site. In addition, it has to be done fast enough to assess several compounds in a relatively short period of time [14].

These algorithms are complemented with scoring functions to sort different poses with the receptor regarding the union force. The scoring functions have some disadvantages because they do not take into account all the interfering forces between the ligand and the receptor; in addition, they also fail in accounting the ligand solvation, the entropic changes when bonding the receptor, the receptor flexibility and the receptor’s conformational changes induced by the same ligand interaction. Based on this, molecular docking is a technique quite fast and effective in terms of time and computational cost [18].

The main software and algorithms developers performing molecular docking are focused on improving aspects such as: receptor’s flexibility, solvation, coupling fragment, post-processing, molecular docking in homology models and comparisons between couplings. Main progress regarding this technique has taken place in the scoring functions area. The use of this methodology has been enhanced over the years; for instance, by doing a single search in PubMed using the word ‘docking’ as a search criterion, it shows the evolution of the publications where molecular docking has been used (Fig. 1).

Many authors have demonstrated that molecular docking software can produce "accurate" binding modes, providing acceptable solutions to the sampling problem; however, the main challenge between the programs and algorithms developers is the poses scoring problem [19]. In this sense, the molecular docking tools are powerful in the prediction of correct binding poses but have critical issues to estimate accurately the corresponding binding affinities.

Fig. (1).

Reported publications number where molecular docking was used. Results from PubMed using the word ‘docking’ as search criteria.


An effort to improve affinity prediction, a rescoring process with other simple functions or solvated-based scoring functions is typically performed. The poses generated by docking program are taken, and methods such as MM/PBSA (Molecular Mechanics/ Poisson-Boltzmann Surface Area) and MM/GBSA (Molecular Mechanics/ Generalized-Born Surface Area) [19-21] include implicit solvent and can be used in order to correct scoring function values and improve docking accuracy. Other strategy is the use of molecular dynamics simulation (MDs) to get conformational sampling of the complex obtained using docking, performing then a subsequent calculation of the binding energy by averaging the score values for different poses extracted from the trajectory [22-24]. Under this approach, the receptor flexibility and the presence of explicit water molecules contribute to a more realistic description of the system, which could have an influence in binding energy calculations. However, the biggest challenge of molecular docking (and also the principal problem) is the accurate prediction of the binding energies, which has major implications for the prediction of novel effective drugs. This process is performed by using scoring functions that score the predicted poses [25].

Fig. (2).

Use of molecular docking programms between 2010 and 2011. Based on a PubMed search. (Taken and adapted from [18]).


Several classes of scoring functions have been developed [25]. These consider the protein as a rigid body [26, 27], or as a soft body [28-30]. They also may consider flexible side chains [31] or certain flexible domains in the target [32-34]. Other empirical scoring functions have shown be useful in the virtual screening and filtering of databases of drug-like compounds at the early stage of drug development processes [34]. Frequently, scoring functions try to reproduce experimental binding affinities, but the ones with popular docking programs do not always yield the best predictions. New programs have arisen in recent years advertising the importance of Molecular Docking in the study of biochemistry sciences and drug design, the most widely used are AutoDock [33, 34] with 29.5% followed by GOLD [37] with 17.5% and Glide [27, 36] with 13.2% (Fig. 2).

Generally speaking, theoretical molecular docking aspects are within the thermodynamics principles, where the main goal is to accurately predict the complex Receptor-Ligand structure [R-L] in equilibrium conditions (Equation 1). As an example in Fig. (3) it was set the Receptor [R] (acetylcholinesterase), the ligand [L] (donepezil), and the complex [R-L] for this particular system.

(1)

Fig. (3) illustrates the ligand Donepezil [39] union (drug used in pharmacological therapy against Alzheimer's disease [40]) to the receptor Acetylcholinesterase described through the Gibbs free energy of binding (∆GBind); this energy is related to the binding affinity by the ratio of Equation 2:

(2)

KA is the affinity constant and Ki is the inhibition constant. The prediction for the accurate structure of the [R-L] complex does not require the experimental information of the constants; however, the prediction of the ligand biological activity does. The terms influencing the scoring function can be classified as follows: Consideration of the [R-L] complex, where aspects such as the ligand and receptor deformation, hydrogen bonds, and electrostatic and steric factors are important. When the equilibrium is considered (Equation 1), the most important factors are: desolvation, rotational and translational entropy.

Molecular docking has been widely used in rational drug design in the past decades with extensive applications in several issues such as:

  • The theoretical study of antitubulin agents with anti-cancer activity [41].
  • The study of the estrogen receptor binding domains [42].
  • Shityakov et al. have studied the blood-barrier permeation by In silico predictive models [43, 44], the dimer complexes of dopamine and levodopa derivatives to assess drug delivery to the central nervous system [45], and the interaction of nanoparticles with multidrug resistance protein [46].
  • The prediction of anti-cancer drugs for the treatment of hypopharyngeal cancer [47].
  • Sabogal et al. established a protocol to identify an interaction pattern between the sea anemone neurotoxins and the potassium channel Kv1.3 [48].
  • Yi et al. studied the metabolic behavior of anisole by using Molecular Docking [49].
Fig. (3).

Ligand donepezil binding to acetylcholinesterase (based on the crystalline structure PDB: 1EVE [39]).


Academically, the amount of research groups who use this technique to study biological systems have increased. In biology, computational techniques have taken great strength, especially in the structural biochemistry field. To prove it, a large number of articles such as the Structural Biochemistry Group at the University of Edinburgh for instance (website: http://www.bch.ed.ac.uk/index.php) have been published. In addition, at multidisciplinary centers of study such as the Department of Molecular and structural Biochemistry of the North Carolina State University - USA (website: http://biochem.ncsu.edu), study programs like biochemistry, molecular biology and biophysics in undergraduate, masters and doctorate levels have been created. This allows students to familiarize with the technique, including it in their academic activities. For these reasons, molecular docking is not only used in the pharmaceutical industry for drug discovery, but also for anyone interested in the study of biological systems and receptor-ligand interactions.

de novo Design

Another computer-based drug design method known as ‘de novo design’ have been implemented to find new potent and selective ligands. de novo design serves as a tool for the discovery of new ligands for biological targets as well as optimization of new ligands [50]. These tools produce (from scratch) novel molecular structures with desired pharmacological properties [51]; subsequently the new structures are docked to find the binding affinities and interacting modes. This new tools and algorithms identify potent ligand-protein interaction characteristics in the binding site, and construct novel molecular structures by assembling atoms and molecular fragments either combined or sequentially. de novo design is a powerful tool, and some authors have reported most potent ligands designed by this tool that their precursors [52], revealing the importance of this tool in the drug design process. Some of these new tools are shown in Table. 1.

Table 1.

de novo design tools (taken and adapted from [53]).


Tool Concept Reference
Builder Combinatory search by recombination of docked molecules [54]
Caveat Database search for fragment fitting [55]
Concerts Fragment-based, stochastic search [56]
Dynamic ligand design Atom-based, structure sampling by simulated annealing [57]
GenStar Atom-Based, molecules growth based on an enzyme contact model [58]
GroupBuild Fragment-based, combinatorial search [59]
Grow Peptide design, sequential growth [60]
Growmol Fragment-based, sequential growth, stochastic search [61]
Hook Linker search for fragments placed by MCSS [62]
Legend Atom-based, stochastic search [63]
LUDI Fragment-based, combinatorial search [64]
MCDNLG Atom-based, stochastic search [65]
MCSS Fragment-based, stochastic sampling [66]
MolMaker Graph-theoretical 3D design [67]
NewLead Fragment-based, builds on 3D pharmacophore-models [68]
Pro-Ligand Fragment-based search [69]
Pro-Select Fragment-based, scaffold-linker approach [70]
Skelgen Small-fragment based, Monte-Carlo search [71]
SME Peptide design, whole-molecule optimization [72]
SMoG Fragment-based, sequential growth, stochastic search [73]
Splice Recombination of ligands retrieved by a 3D database search [74]
Sprout Fragment-based, sequential growth, combinatorial search [75]
Topas Fragment-based, evolutionary search [76]

Automated de novo design is a very useful tool for hit and lead-ligand identification; design molecules, provide ideas for medical chemist and pharmaceutical industries, and to develop novel leads with desired chemical characteristics.

Virtual Screening

A computational technique widely used in drug design is Virtual Screening (VS). Its main objective is to search for specific information in compounds or molecules libraries with similar structural properties that can acceptably interact with a therapeutic target and is an important tool to access novel drug-like compounds [77]. Prior to the VS, massive molecular docking was used but it had a big computational cost and spent a lot of computational time; due to this, a technique - which enables the efficiency improvement to search for compounds in large molecular libraries - was needed. The use of VS (as well as molecular docking) has increased in recent years, rising from an average of 10 publications reported in the 70s - 80s, to an average of more than 1000 publications reported in previous years (Fig. 4).

Fig. (4).

Reported publications number where virtual screening was used. Results from PubMed using the words ‘virtual screening’ as search criteria.


Virtual screening can be performed through two approaches, Ligand-Based Virtual Screening (LBVS) and Structure-Based Virtual Screening (SBVS). LBVS is used when there is not knowledge about the three-dimensional structure of the therapeutic target (receptor), and yet you have a group of molecules with some biological activity against that target; in that sense, structure relationships can be determined. The molecules group with reported activity can also establish the molecular characteristics. These characteristics allow the interaction between the molecules and the molecular target. Based on those (such as acceptor and donors hydrogen bonds, metals, hydrophobic moieties, aromatic rings among others), researchers can screen databases to find other molecules to fit in the profile of the initial group – with the same molecular descriptors – which can possibly present biological activity against the studied therapeutic target [78]. In the LBSV the lead-ligands generated are ranked based on their similarity score obtained by different methods or algorithms [77].

On the other hand, Structure-Based Virtual Screening (SBVS) establishes its action mechanism in the studied molecular target structure and begins with the identification of the potential ligand-binding site on the target. Then, one of its greatest challenges is to determine the target 3D structure by NMR, X-ray crystallography or molecular modeling. Later, an automated and rapid docking of a large number of chemical compounds against the 3D structure of the molecular target is done. This technique has helped with the identification of potential therapeutic molecules for specific pathologies [77].

To use this technique in an effective way, a protocol with several stages should be established:

  1. Ligands Database Preparation
    It is necessary to prepare all the molecules to be studied with the SBVS. Tautomers, isomers, protonation states, ionization states, enantiomers, etc. must be considered. This database must be prepared in accordance with the software to use.
  2. Receptor Preparation
    Next, the molecular target (Receptor) must be provided. It may uses the same ligands database preparation software.
  3. Molecular Docking
    It is necessary to determine the common pharmacophore in the ligand database; after this - and through the appropriate program - massive molecular docking is carried out.
  4. Post-Processing
    The post-processing is performed in accordance with the user requirements through the analysis of the score function, geometric analysis, shape complementarity, solvation corrections, entropic changes, pose clustering, etc.
  5. Compounds Selection
    The final step is the compounds selection for biological tests, or to be further analyzed by other computational techniques [79].

Nowadays, there are several comprehensive programs used to perform VS such as: AutoDock Vina [35], Glide [28, 38, 80], DOCK [81], FlexX [82], GOLD [37], ICM [83], among others. Additionally - and due to the high computational cost for the SBVS implementation - several research groups around the world have developed free software and algorithms to be accessed through virtual platforms. These developers also have the advantage of having extensive updated databases, which allows the virtual screening of millions of molecules. For instance, among the major online servers to perform VS there is the ZINCPharmer, sponsored by the University of Pittsburgh (website: http://zincpharmer.csb.pitt.edu) [84]. This software allows researchers to perform online SBVS for free, mainly based on the established ligand pharmacophore used as input. The screening is basically performed in the ZINC database [7], a non-commercial database with more than 22 million compounds.

In recent years, researchers who use this technique in their computational research have increased, and the establishment of data sets for potential therapeutic agents from VS can be observed, for example: SENP2 inhibitors [85], kinases proteins agonists [86], inhibitors with vasodilator activity [87], inhibitors of HIV-1 reverse transcriptase [88], among others.

The development of these dataset has allowed the access and use of information more efficiently, by evaluating compounds against new molecular targets and finding many potential drug candidates (better than those before the implementation of this technique).

Table 2.

Some active compounds identified by computational methods.


Compound structure Use Computational method References
Antimycobacterial agents Virtual Screening [89]
SENP2 inhibitors Virtual screening, Molecular Docking [85]
TASK-3 channel antagonist Pharmacophore based virtual screening [90]
GAPGH inhibitor Combinatorial docking [91]
Aldose reductase inhibitor LBVS [92]
Ca2+ antagonist Pharmacophore searching [93]
Kv1.5 channel blocker Fragment based, de novo design [76]
Thrombin inhibitor Combinatorial docking, de novo design [94]
Antiretroviral agent Virtual Screening [95]

On the other hand, VS is used as a mechanism for the evaluation of new potential drugs against certain molecular targets in rational drug design. In biochemical terms, this enables the computational evaluation of large molecules quantities in a short time with a relatively low computational cost. It means that after a screening, score, post-processing and visual inspection, it can continue to the laboratory to perform in vitro testing with the best candidates obtained via VS allowing the properties evaluation such as the biological activity, the enzyme kinetics, among others. This contributes with the biological system study in the drug discovery process.

Several bioactive molecules identified, optimized or designed by using molecular docking, de novo design or VS methods (Table. 2) let us understand the importance of the application and use of them as computational tools for drug discovery to explore the chemical space of a binding site and present novel, potent and selective pharmacological alternatives in medical chemistry and drug discovery.

CONCLUSION

Although in appearance, a direct relationship between the computational methodologies - molecular docking and virtual screening - with research in biological and medical sciences does not exist, in recent years most of the authors whom develop biological systems research in areas such as biochemistry, medicinal chemistry, pharmacology, organic chemistry, theoretical chemistry, etc., use computational methods to support, improve and reinforce their research. The structural studies of these systems at theoretical level assist with the construction of a wide and accurate perspective of the issue, and enable the knowledge consolidation.

In addition, there is still a long way to develop methodologies with “thinner” approaches (accurate and precise). Although methodologies such as the molecular docking and virtual screening are effective (in terms of time and computational cost), they have also a poor predictive capability. There are other methods involving a higher calculation level for researchers who study dynamic biological systems (Molecular Dynamics), or bonds and interactions in an electronic level (Quantum Mechanics). These like the other methods are more time-consuming in terms of computational cost, but may provide better results; it all depends on the researcher interests and funding. However, to improve computational techniques such as Molecular Docking and Virtual Screening, a change in the structural paradigm is needed, and this will be achieved when the understanding of the interactions between a drug and a receptor improves.

CONFLICT OF INTEREST

The author confirms that this article content has no conflict of interest.

ACKNOWLEDGEMENTS

I thank to Universidad de Talca for doctoral scholarship.

REFERENCES

[1] Zheng M, Liu X, Xu Y, Li H, Luo C, Jiang H. Computational methods for drug design and discovery: focus on China. Trends Pharmacol Sci 2013; 34(10): 549-59.
[2] Blundell TL, Jhoti H, Abell C. High-throughput crystallography for lead discovery in drug design. Nat Rev Drug Discov 2002; 1(1): 45-54.
[3] Hsiao-lu DL, Sahl SJ, Lew MD, Moerner WE. The double-helix microscope microscope super-resolves extended biological structures by localizing single blinking molecules in three dimensions with nanoscale precision. Appl Phys Lett 2012; 100(15): 153701.
[4] Berman HM, Bhat TN, Bourne PE, et al. The protein data bank and the challenge of structural genomics. Nat Struct Mol Biol 2000; 7: 957-9.
[5] Westbrook J, Feng Z, Chen L, Yang H, Berman HM. The protein data bank and structural genomics. Nucleic Acids Res 2003; 31(1): 489-91.
[6] Seiler KP, George GA, Happ MP, et al. others. Chembank: a small-molecule screening and cheminformatics resource database. Nucleic Acids Res 2008; 36(Suppl. 1): D351-9.
[7] Irwin JJ, Shoichet BK. ZINC-a free database of commercially available compounds for virtual screening. J Chem Inf Model 2005; 45(1): 177-82.
[8] Mandal S, Moudgil M, Mandal SK. Rational drug design. Eur J Pharmacol 2009; 625(1): 90-100.
[9] De-la-Torre P, Treuer A V, Gutierrez M, et al. Synthesis and in silico analysis of the quantitative structure-activity relationship of heteroaryl-acrylonitriles as ache inhibitors J Taiwan Inst Chem Eng 2015.
[10] Brown HC. Foundations of the structural theory. J Chem Educ 1959; 36(3): 104.
[11] Wiswesser WJ. Johann Josef Loschmidt (1821-1895): A forgotten genius. Aldrichim Acta 1989; 22: 17-9.
[12] Metropolis N, Rosenbluth AW, Rosenbluth MN, Teller AH, Teller E. Equation of state calculations by fast computing machines. J Chem Phys 1953; 21(6): 1087-92.
[13] Takahashi O, Masuda Y, Muroya A, Furuya T. Theory of docking scores and its application to a customizable scoring function. SAR QSAR Environ Res 2010; 21(5-6): 547-58.
[14] Kitchen DB, Decornez H, Furr JR, Bajorath J. Docking and scoring in virtual screening for drug discovery: Methods and applications. Nat Rev Drug Discov 2004; 3(11): 935-49.
[15] Mohan V, Gibbs AC, Cummings MD, Jaeger EP, DesJarlais RL. Docking: successes and challenges. Curr Pharm Des 2005; 11(3): 323-33.
[16] Krovat EM, Steindl T, Langer T. Recent advances in docking and scoring. Curr Comput Aided Drug Des 2005; 1(1): 93-102.
[17] Leach AR, Shoichet BK, Peishoff CE. Prediction of protein-ligand interactions. Docking and scoring: successes and gaps.docking and scoring in virtual screening for drug discovery: Methods and applications. J Med Chem 2006; 49(20): 5851-5.
[18] Yuriev E, Ramsland Pa. Latest developments in molecular docking: 2010-2011 in review. J Mol Recognit 2013; 26(5): 215-39.
[19] Yuriev E, Agostino M, Ramsland PA. Challenges and advances in computational docking: 2009 in review. J Mol Recognit 2011; 24(2): 149-64.
[20] Lee MR, Sun Y. Improving docking accuracy through molecular mechanics generalized born optimization and scoring. J Chem Theory Comput 2007; 3(3): 1106-19.
[21] Adasme-Carreño F, Muñoz-Gutierrez C, Caballero J, Alzate-Morales JH. Performance of the mm/gbsa scoring using a binding site hydrogen bond network-based frame selection: the protein kinase case. Phys Chem Chem Phys 2014; 16(27): 14047-58.
[22] Mena-Ulecia K, Vergara-Jaque A, Poblete H, Tiznado W, Caballero J. Study of the affinity between the protein kinase pka and peptide substrates derived from kemptide using molecular dynamics simulations and MM/GBSA. PLoS One 2014; 9(10): e109639.
[23] Pak Y, Wang S. Application of a molecular dynamics simulation method with a generalized effective potential to the flexible molecular docking problems. J Phys Chem B 2000; 104(2): 354-9.
[24] Caballero J, Alzate-Morales JH. Molecular dynamics of protein kinase-inhibitor complexes: a valid structural information. Curr Pharm Des 2012; 18(20): 2946-63.
[25] Elokely KM, Doerksen RJ. Docking challenge: protein sampling and molecular docking performance. J Chem Inf Model 2013; 53(8): 1934-45.
[26] McGann MR, Almond HR, Nicholls A, Grant JA, Brown FK. Gaussian docking functions. Biopolymers 2003; 68(1): 76-90.
[27] McGaughey GB, Sheridan RP, Bayly CI, et al. Comparison of topological, shape, and docking methods in virtual screening. J Chem Inf Model 2007; 47(4): 1504-19.
[28] Friesner RA, Banks JL, Murphy RB, et al. others. Glide: a new approach for rapid, accurate docking and scoring. 1. method and assessment of docking accuracy. J Med Chem 2004; 47(7): 1739-49.
[29] Jones G, Willett P, Glen RC. Molecular recognition of receptor sites using a genetic algorithm with a description of desolvation. J Mol Biol 1995; 245(1): 43-53.
[30] Jones G, Willett P, Glen RC, Leach AR, Taylor R. Development and validation of a genetic algorithm for flexible docking. J Mol Biol 1997; 267(3): 727-48.
[31] Barreca ML, Iraci N, De Luca L, Chimirri A. Induced-fit docking spproach provides insight into the binding mode and mechanism of action of HIV-1 integrase inhibitors. Chem Med Chem 2009; 4(9): 1446-56.
[32] Cavasotto CN, Abagyan RA. Protein flexibility in ligand docking and virtual screening to protein kinases. J Mol Biol 2004; 337(1): 209-25.
[33] Davis IW, Baker D. Rosettaligand docking with full ligand and receptor flexibility. J Mol Biol 2009; 385(2): 381-92.
[34] Shityakov S, Förster C. In silico structure-based screening of versatile p-glycoprotein inhibitors using polynomial empirical scoring functions. Adv Appl Bioinform Chem 2014; 7: 1.
[35] Trott O, Olson AJ. Autodock vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J Comput Chem 2010; 31(2): 455-61.
[36] Goodsell DS, Morris GM, Olson AJ. Automated docking of flexible ligands: applications of auto dock. J Mol Recognit 1996; 9(1): 1-5.
[37] Verdonk ML, Cole JC, Hartshorn MJ, Murray CW, Taylor RD. Improved protein-ligand docking using gold. Proteins Struct Funct Bioinforma 2003; 52(4): 609-23.
[38] Halgren TA, Murphy RB, Friesner RA, et al. Glide: a new approach for rapid, accurate docking and scoring. 2. Enrichment factors in database screening. J Med Chem 2004; 47(7): 1750-9.
[39] Kryger G, Silman I, Sussman JL. Structure of acetylcholinesterase complexed with E2020 (Aricept®): implications for the design of new anti-alzheimer drugs. Structure 1999; 7(3): 297-307.
[40] Krishnan KR, Charles HC, Doraiswamy PM, et al. Randomized. Placebo-Controlled Trial of the Effects of Donepezil on Neuronal Markers and Hippocampal Volumes in Alzheimer’s Disease 2014.
[41] Penthala NR, Zong H, Ketkar A, et al. Synthesis, anticancer activity and molecular docking studies on a series of heterocyclic trans-cyanocombretastatin analogues as antitubulin agents. Eur J Med Chem 2015; 92: 212-20.
[42] Subramanian SR, Zuckerman DM. Docking to the highly flexible estrogen receptor ligand binding domain via mixed-resolution Monte Carlo. Biophys J 2015; 108(2): 319a.
[43] Shityakov S, Neuhaus W, Dandekar T, Förster C. Analysing molecular polar surface descriptors to predict blood-brain barrier permeation. Int J Comput Biol Drug Des 2013; 6(1-2): 146-56.
[44] Shityakov S, Förster C. In silico predictive model to determine vector-mediated transport properties for the blood-brain barrier choline transporter. Adv Appl Bioinform Chem 2014; 7: 23.
[45] Shityakov S, Broscheit J, Foerster C. α-cyclodextrin dimer complexes of dopamine and levodopa derivatives to assess drug delivery to the central nervous system: adme and molecular docking studies.in silico predictive model to determine vector-mediated transport properties for the blood--brain barrier choline transporter Int J Nanomedicine 2012; 7: 3211.
[46] Shityakov S, Förster C. Multidrug resistance protein P-Gp interaction with nanoparticles (fullerenes and carbon nanotube) to assess their drug delivery potential: A theoretical molecular docking study. Int J Comput Biol Drug Des 2013; 6(4): 343-57.
[47] Gao J, Chao LV. Molecular docking to predict hypopharyngeal cancer treatment drugs: tiliroside-noscapine interaction. Lat Am J Pharm 2014; 33(8): 1397-9.
[48] Sabogal-arango A, Barreto GE, Ramírez-sánchez D, et al. Computational insights of the interaction among sea anemones neurotoxins and Kv1.3 channel. Bioinform Biol Insights 2014; 8: 73-81.
[49] Yi Z, Lu M. Deep understanding of metabolic behaviour of alisol a using molecular docking. Lat Am J Pharm 2014; 33(3): 522-4.
[50] Kutchukian PS, Shakhnovich EI. de novo design: balancing novelty and confined chemical space. Expert Opin Drug Discov 2010; 5(8): 789-812.
[51] Schneider G, Fechner U. Computer-based de novo design of drug-like molecules. Nat Rev Drug Discov 2005; 4(8): 649-63.
[52] Chirapu SR, Pachaiyappan B, Nural HF, et al. others. Molecular modeling, synthesis, and activity studies of novel biaryl and fused-ring bace1 inhibitors. Bioorg Med Chem Lett 2009; 19(1): 264-74.
[53] Schneider G, Böhm HJ. Virtual screening and fast automated docking methods. Drug Discov Today 2002; 7(1): 64-70.
[54] Roe DC, Kuntz ID. BUILDER v. 2: Improving the chemistry of a de novo design strategy. J Comput Aided Mol Des 1995; 9(3): 269-82.
[55] Lauri G, Bartlett PA. CAVEAT: A program to facilitate the design of organic moleculesimproving the chemistry of a de novo design strategy.. J Comput Aided Mol Des 1994; 8(1): 51-66.
[56] Pearlman DA, Murcko MA. CONCERTS: Dynamic connection of fragments as an approach to de novo ligand design. J Med Chem 1996; 39(8): 1651-63.
[57] Miranker A, Karplus M. An automated method for dynamic ligand design. Proteins Struct Funct Bioinforma 1995; 23(4): 472-90.
[58] Rotstein SH, Murcko MA. GenStar: a method for de novo drug design. J Comput Aided Mol Des 1993; 7(1): 23-43.
[59] Rotstein SH, Murcko MA. GroupBuild: a fragment-based method for de novo drug design. J Med Chem 1993; 36(12): 1700-10.
[60] Moon JB, Howe WJ. Computer design of bioactive molecules: a method for receptor-based de novo ligand designan automated method for dynamic ligand design.. Proteins Struct Funct Bioinforma 1991; 11(4): 314-28.
[61] Bohacek RS, McMartin C. Multiple highly diverse structures complementary to enzyme binding sites: results of extensive application of a de novo design method incorporating combinatorial growth. J Am Chem Soc 1994; 116(13): 5560-71.
[62] Eisen MB, Wiley DC, Karplus M, Hubbard RE. HOOK: a program for finding novel molecular architectures that satisfy the chemical and steric requirements of a macromolecule binding site. Proteins Struct Funct Bioinforma 1994; 19(3): 199-221.
[63] Nishibata Y, Itai A. Automatic creation of drug candidate structures based on receptor structure. Starting point for artificial lead generation. Tetrahedron 1991; 47(43): 8985-90.
[64] Böhm H-J. The computer program LUDI: s new method for the de novo design of enzyme inhibitors. J Comput Aided Mol Des 1992; 6(1): 61-78.
[65] Gehlhaar DK, Moerder KE, Zichi D, Sherman CJ, Ogden RC, Freer ST. de novo design of enzyme inhibitors by monte carlo ligand generation. J Med Chem 1995; 38(3): 466-72.
[66] Miranker A, Karplus M. Functionality maps of binding sites: A multiple copy simultaneous search method. Proteins Struct Funct Bioinforma 1991; 11(1): 29-34.
[67] Clark DE, Firth MA, Murray CW. MOLMAKER: de novo generation of 3D databases for use in drug design. J Chem Inf Comput Sci 1996; 36(1): 137-45.
[68] Tschinke V, Cohen NC. The NEWLEAD program: A new method for the design of candidate structures from pharmacophoric hypotheses. J Med Chem 1993; 36(24): 3863-70.
[69] Clark DE, Frenkel D, Levy SA, et al. PRO_LIGAND: An approach to de novo molecular design. 1. Application to the design of organic molecules. J Comput Aided Mol Des 1995; 9(1): 13-32.
[70] Murray CW, Clark DE, Auton TR, et al. PRO_SELECT: Combining structure-based drug design and combinatorial chemistry for rapid lead discovery. 1. Technology. J Comput Aided Mol Des 1997; 11(2): 193-207.
[71] Todorov NP, Dean PM. Evaluation of a method for controlling molecular scaffold diversity in de novo ligand design. J Comput Aided Mol Des 1997; 11(2): 175-92.
[72] Schneider G, Schrödl W, Wallukat G, et al. Peptide design by artificial neural networks and computer-based evolutionary search. Proc Natl Acad Sci USA 1998; 95(21): 12179-84.
[73] DeWitte RS, Ishchenko AV, Shakhnovich EI. SMoG: de novo design method based on simple, fast, and accurate free energy estimates. 2. Case studies in molecular design. J Am Chem Soc 1997; 119(20): 4608-17.
[74] Ho CM, Marshall GR. FOUNDATION: A program to retrieve all possible structures containing a user-defined minimum number of matching query elements from three-dimensional databases. J Comput Aided Mol Des 1993; 7(1): 3-22.
[75] Mata P, Gillet VJ, Johnson AP, et al. SPROUT: 3D structure generation using templates. J Chem Inf Comput Sci 1995; 35(3): 479-93.
[76] Schneider G, Lee M-L, Stahl M, Schneider P. de novo design of molecular architectures by evolutionary assembly of drug-derived building blocks. J Comput Aided Mol Des 2000; 14(5): 487-94.
[77] Reddy, a S, Pati S P, Kumar P P, Pradeep H N, Sastry G N. Virtual screening in drug discovery - a computational perspective. Curr Protein Pept Sci 2007; 8(4): 329-51.
[78] Koeppen H, Kriegl J, Lessel U, Tautermann C S, Wellenzohn B. Structure-based virtual screening: an overview. Drug Discov 2011; 7(20): 1047-55.
[79] Lyne PD. Structure-based virtual screening: an overview. Drug Discov Today 2002; 7(20): 1047-55.
[80] Friesner RA, Murphy RB, Repasky MP, et al. Extra precision glide: docking and scoring incorporating a model of hydrophobic enclosure for protein-ligand complexes. J Med Chem 2006; 49(21): 6177-96.
[81] Ewing TJ, Makino S, Skillman AG, Kuntz ID. DOCK 4.0: search strategies for automated molecular docking of flexible molecule databases. J Comput Aided Mol Des 2001; 15(5): 411-28.
[82] Schellhammer I, Rarey M. FlexX-Scan: fast, structure-based virtual screening. PROTEINS Struct Funct Bioinforma 2004; 57(3): 504-17.
[83] Abagyan R, Totrov M, Kuznetsov D. ICM-a new method for protein modeling and design: applications to docking and structure prediction from the distorted native conformation. J Comput Chem 1994; 15(5): 488-506.
[84] Koes D R, Camacho C J. ZINC Pharmer: Pharmacophore search of the ZINC database. Nucleic Acids Res 2012; 40(Web Server issue): W409-14.
[85] Kumar A, Ito A, Takemoto M, Yoshida M, Zhang KY. Identification of 1, 2, 5-oxadiazoles as a new class of SENP2 inhibitors using structure based virtual screening. J Chem Inf Model 2014; 54(3): 870-80.
[86] Jain A, Trivedi V. Docking and virtual screening to identify pkc agonists: potentials in anticancer therapeutics. Curr Comput Aided Drug Des 2014; 10(1): 50-8.
[87] Mittal A, Paliwal S, Sharma M, Singh A, Sharma S, Yadav D. Pharmacophore based virtual screening, molecular docking and biological evaluation to identify novel PDE5 inhibitors with vasodilatory activity. Bioorg Med Chem Lett 2014; 24(14): 3137-41.
[88] Hecht FM, Grant RM, Petropoulos CJ, et al. others. Sexual transmission of an HIV-1 variant resistant to multiple reverse-transcriptase and protease inhibitors. N Engl J Med 1998; 339(5): 307-11.
[89] Kincaid VA, London N, Wangkanont K, et al. others. Virtual screening for UDP-Galactopyranose mutase ligands identifies a new class of antimycobacterial agents. ACS Chem Biol 2015.
[90] Coburn C, Luo Y, Cui M, et al. Discovery of a pharmacologically active antagonist of the two-pore-domain potassium channel K2P9.1 (TASK-3). Chem Med Chem 2012; 7(1): 123-33.
[91] Bressi JC, Verlinde CL, Aronov AM, et al. others. Adenosine analogues as selective inhibitors of glyceraldehyde-3-phosphate dehydrogenase of trypanosomatidae via structure-based drug design. J Med Chem 2001; 44(13): 2080-93.
[92] Iwata Y, Arisawa M, Hamada R, et al. Discovery of novel aldose reductase inhibitors using a protein structure-based approach: 3D-database search followed by design and synthesis. J Med Chem 2001; 44(11): 1718-28.
[93] Schneider G, Clément-Chomienne O, Hilfiger L, et al. Virtual screening for bioactive molecules by evolutionary de novo design. Angew Chem Int Ed 2000; 39(22): 4130-3.
[94] Böhm H-J, Banner DW, Weber L. Combinatorial docking and combinatorial chemistry: design of potent non-peptide thrombin inhibitors. J Comput Aided Mol Des 1999; 13(1): 51-6.
[95] Shityakov S, Dandekar T. Lead Expansion and virtual screening of indinavir derivate HIV-1 protease inhibitors using pharmacophoric-shape similarity scoring function. Bioinformation 2010; 4(7): 295.