GlamDock: Development and Validation of a New Docking Tool on Several Thousand Protein−Ligand Complexes

In silico docking and molecular dynamics for the discovery of inhibitors of enteric methane production in ruminants - A review

The increase in methane emissions, a major greenhouse gas, threatens human well-being and global ecosystems due to its contribution to global warming. Livestock, particularly ruminants, have been a major research topic in recent decades due to their methane production. Therefore, the objective of the current review was to comprehensively discuss the in silico techniques used to mitigate methane production from ruminants. The review covers the principles of in silico docking and molecular dynamics, which can be used to develop methanogenesis inhibitors. It also discusses specific methanogen enzymes as potential targets for inhibitor development. Furthermore, in silico-based methanogenesis inhibitor development studies have been reviewed with the authors' opinions. The further use of in silico-based research techniques, including artificial intelligence-based systems, is encouraged to help reduce methane production from livestock more efficiently and costeffectively.

Watermelon: setup and validation of an in silico fragment-based approach

We present a new computational approach, named Watermelon, designed for the development of pharmacophore models based on receptor structures. The methodology involves the sampling of potential hotspots for ligand interactions within a protein target's binding site, utilising molecular fragments as probes. By employing docking and molecular dynamics (MD) simulations, the most significant interactions formed by these probes within distinct regions of the binding site are identified. These interactions are subsequently transformed into pharmacophore features that delineates key anchoring sites for potential ligands. The reliability of the approach was experimentally validated using the monoacylglycerol lipase (MAGL) enzyme. The generated pharmacophore model captured features representing ligand-MAGL interactions observed in various X-ray co-crystal structures and was employed to screen a database of commercially available compounds, in combination with consensus docking and MD simulations. The screening successfully identified two new MAGL inhibitors with micromolar potency, thus confirming the reliability of the Watermelon approach.

Enhancing SILCS-MC via GPU Acceleration and Ligand Conformational Optimization with Genetic and Parallel Tempering Algorithms

In the domain of computer-aided drug design, achieving precise and accurate estimates of ligand-protein binding is paramount in the context of screening extensive drug libraries and performing ligand optimization. A fundamental aspect of the SILCS (site identification by ligand competitive saturation) methodology lies in the generation of comprehensive 3D free-energy functional group affinity maps (FragMaps), encompassing the entirety of the target molecule structure. These FragMaps offer an intricate landscape of functional group affinities across the protein, bilayer, or RNA, acting as the basis for subsequent SILCS-Monte Carlo (MC) simulations wherein ligands are docked to the target molecule. To augment the efficiency and breadth of ligand sampling capabilities, we implemented an improved SILCS-MC methodology. By harnessing the parallel computing capability of GPUs, our approach facilitates concurrent calculations over multiple ligands and binding sites, markedly enhancing the computational efficiency. Moreover, the integration of a genetic algorithm (GA) with MC allows us to employ an evolutionary approach to perform ligand sampling, assuring enhanced convergence characteristics. In addition, the potential utility of parallel tempering (PT) to improve sampling was investigated. Implementation of SILCS-MC on GPU architecture is shown to accelerate the speed of SILCS-MC calculations by over 2-orders of magnitude. Use of GA and PT yield improvements over Markov-chain MC, increasing the precision of the resultant docked orientations and binding free energies, though the extent of improvements is relatively small. Accordingly, significant improvements in speed are obtained through the GPU implementation with minor improvements in the precision of the docking obtained via the tested GA and PT algorithms.

Applying molecular docking to pesticides

Pesticide creation is related to the development of sustainable agricultural and ecological safety, and molecular docking technology can effectively help in pesticide innovation. This paper introduces the basic theory behind molecular docking, pesticide databases, and docking software. It also summarizes the application of molecular docking in the pesticide field, including the virtual screening of lead compounds, detection of pesticides and their metabolites in the environment, reverse screening of pesticide targets, and the study of resistance mechanisms. Finally, problems with the use of molecular docking technology in pesticide creation are discussed, and prospects for the future use of molecular docking technology in new pesticide development are discussed. © 2023 Society of Chemical Industry.

Computational Approaches to Designing Antiviral Drugs against COVID-19: A Comprehensive Review

The global impact of the COVID-19 pandemic caused by SARS-CoV-2 necessitates innovative strategies for the rapid development of effective treatments. Computational methodologies, such as molecular modelling, molecular dynamics simulations, and artificial intelligence, have emerged as indispensable tools in the drug discovery process. This review aimed to provide a comprehensive overview of these computational approaches and their application in the design of antiviral agents for COVID-19. Starting with an examination of ligand-based and structure-based drug discovery, the review has delved into the intricate ways through which molecular modelling can accelerate the identification of potential therapies. Additionally, the investigation extends to phytochemicals sourced from nature, which have shown promise as potential antiviral agents. Noteworthy compounds, including gallic acid, naringin, hesperidin, Tinospora cordifolia, curcumin, nimbin, azadironic acid, nimbionone, nimbionol, and nimocinol, have exhibited high affinity for COVID-19 Mpro and favourable binding energy profiles compared to current drugs. Although these compounds hold potential, their further validation through in vitro and in vivo experimentation is imperative. Throughout this exploration, the review has emphasized the pivotal role of computational biologists, bioinformaticians, and biotechnologists in driving rapid advancements in clinical research and therapeutic development. By combining state-of-the-art computational techniques with insights from structural and molecular biology, the search for potent antiviral agents has been accelerated. The collaboration between these disciplines holds immense promise in addressing the transmissibility and virulence of SARS-CoV-2.

Using a genetic algorithm to find molecules with good docking scores

A graph-based genetic algorithm (GA) is used to identify molecules (ligands) with high absolute docking scores as estimated by the Glide software package, starting from randomly chosen molecules from the ZINC database, for four different targets: Bacillus subtilis chorismate mutase (CM), human β2-adrenergic G protein-coupled receptor (β2AR), the DDR1 kinase domain (DDR1), and β-cyclodextrin (BCD). By the combined use of functional group filters and a score modifier based on a heuristic synthetic accessibility (SA) score our approach identifies between ca 500 and 6,000 structurally diverse molecules with scores better than known binders by screening a total of 400,000 molecules starting from 8,000 randomly selected molecules from the ZINC database. Screening 250,000 molecules from the ZINC database identifies significantly more molecules with better docking scores than known binders, with the exception of CM, where the conventional screening approach only identifies 60 compounds compared to 511 with GA+Filter+SA. In the case of β2AR and DDR1, the GA+Filter+SA approach finds significantly more molecules with docking scores lower than −9.0 and −10.0. The GA+Filters+SA docking methodology is thus effective in generating a large and diverse set of synthetically accessible molecules with very good docking scores for a particular target. An early incarnation of the GA+Filter+SA approach was used to identify potential binders to the COVID-19 main protease and submitted to the early stages of the COVID Moonshot project, a crowd-sourced initiative to accelerate the development of a COVID antiviral.

Protein Binding Pocket Optimization for Virtual High-Throughput Screening (vHTS) Drug Discovery

The virtual high-throughput screening (vHTS) approach has been widely used for large database screening to identify potential lead compounds for drug discovery. Due to its high computational demands, docking that allows receptor flexibility has been a challenging problem for virtual screening. Therefore, the selection of protein target conformations is crucial to produce useful vHTS results. Since only a single protein structure is used to screen large databases in most vHTS studies, the main challenge is to reduce false negative rates in selecting compounds for in vitro tests. False negatives are most likely to occur when using apo structures or homology models of protein targets due to the small volume of the binding pocket formed by incorrect side-chain conformations. Even holo protein structures can exhibit high false negative rates due to ligand-induced fit effects, since the shape of the binding pocket highly depends on its bound ligand. To reduce false negative rates and improve success rates for vHTS in drug discovery, we have developed a new Monte Carlo-based approach that optimizes the binding pocket of protein targets. This newly developed Monte Carlo pocket optimization (MCPO) approach was assessed on several datasets showing promising results. The binding pocket optimization approach could be a useful tool for vHTS-based drug discovery, especially in cases when only apo structures or homology models are available.

Structure-Based Virtual Screening: From Classical to Artificial Intelligence

The drug development process is a major challenge in the pharmaceutical industry since it takes a substantial amount of time and money to move through all the phases of developing of a new drug. One extensively used method to minimize the cost and time for the drug development process is computer-aided drug design (CADD). CADD allows better focusing on experiments, which can reduce the time and cost involved in researching new drugs. In this context, structure-based virtual screening (SBVS) is robust and useful and is one of the most promising in silico techniques for drug design. SBVS attempts to predict the best interaction mode between two molecules to form a stable complex, and it uses scoring functions to estimate the force of non-covalent interactions between a ligand and molecular target. Thus, scoring functions are the main reason for the success or failure of SBVS software. Many software programs are used to perform SBVS, and since they use different algorithms, it is possible to obtain different results from different software using the same input. In the last decade, a new technique of SBVS called consensus virtual screening (CVS) has been used in some studies to increase the accuracy of SBVS and to reduce the false positives obtained in these experiments. An indispensable condition to be able to utilize SBVS is the availability of a 3D structure of the target protein. Some virtual databases, such as the Protein Data Bank, have been created to store the 3D structures of molecules. However, sometimes it is not possible to experimentally obtain the 3D structure. In this situation, the homology modeling methodology allows the prediction of the 3D structure of a protein from its amino acid sequence. This review presents an overview of the challenges involved in the use of CADD to perform SBVS, the areas where CADD tools support SBVS, a comparison between the most commonly used tools, and the techniques currently used in an attempt to reduce the time and cost in the drug development process. Finally, the final considerations demonstrate the importance of using SBVS in the drug development process.

Screening of different drug design tools to predict the mode of action of steroidal derivatives as anti-cancer agents

There is a pressing need to discover and develop novel drugs against cancer. With the new era of bioinformatics, which integrates different aspects, drug development has been tremendously improved. Recently, extensive research was directed towards the rational modification of steroid molecules against different disease especially cancer. Moreover, heterocyclic steroid derivatives have shown a lot of different biological activities such as antimicrobial, anti-inflammatory, and anti-cancer activities. Molecular docking methods can be used to explore how the steroid derivatives conformations can adopt within the binding sites of specific macromolecular targets involved in cancer progression. We conducted this study to investigate the accuracy of different molecular docking calculations using different steroidal molecular targets, and to define the most accurate one to study the mode of action of steroid derivatives as potential anti-cancer drugs. Our results revealed that the Dock6, PLANTS, AutoDock, GLIDE (SP and XP), and GOLD (ASP, Chemscore, and PLP) software were able to maintain the binding mode of the co-crystallized ligands inside their proteins by achieving RMSD values lower than two. Moreover, molecular docking study revealed that compound 4, and 5 are promising steroidal derivatives as anti-cancer drugs. Further on, the cytotoxic activity of the selected steroidal derivatives were tested against leukemia cell line using MTT assay. The results revealed that compound 4, and 5 were potential cytotoxic agents against THP-1 cells (IC₅₀s were 44.67 µM, and 46.77 µM, respectively), these results are in agreement with the molecular docking study.

A Molecular Docking Approach to Evaluate the Pharmacological Properties of Natural and Synthetic Treatment Candidates for Use against Hypertension

Cardiovascular diseases (CVDs) have become the leading cause of disability and death worldwide, particularly in low- and middle-income countries. Hypertension, a major cause of CVD progression, is widely attributable to genetic, behavioral, and environmental risk factors. Among the genetic reasons, angiotensin II enzyme, produced as a result of abnormal functioning of the renin⁻angiotensin system, is reported as the foremost cause of hypertension. A cascade of genes, including those encoding for WNK kinases (WNK1 and WNK4), Bp1, Bp2, angiotensinogen, and other enzymes, is involved in the conversion of angiotensin I to angiotensin II. However, the angiotensin-converting enzyme (ACE) plays a crucial role in this pathway. Therefore, ACE could be a potential therapeutic target in regulating the conversion of angiotensin I to angiotensin II and eventually controlling hypertension. In this study, a molecular docking-based approach was utilized for identifying and evaluating potential inhibitors of ACE present in herbs, other natural sources, and synthetic sources, on the basis of these compounds' binding affinities and other physicochemical features. In addition, the suitability of these inhibitors as drugs for biological systems, considering their adsorption, distribution, metabolism, and excretion (ADME), was predicted using Lipinski's rule. In conclusion, our study provides a novel and clearer insight into the interaction properties of known putative inhibitors of ACE.

Extensive Reliability Evaluation of Docking-Based Target-Fishing Strategies

The development of target-fishing approaches, aimed at identifying the possible protein targets of a small molecule, represents a hot topic in medicinal chemistry. A successful target-fishing approach would allow for the elucidation of the mechanism of action of all therapeutically interesting compounds for which the actual target is still unknown. Moreover, target-fishing would be essential for preventing adverse effects of drug candidates, by predicting their potential off-targets, and it would speed up drug repurposing campaigns. However, due to the huge number of possible protein targets that a small-molecule might interact with, experimental target-fishing approaches are out of reach. In silico target-fishing represents a valuable alternative, and examples of receptor-based approaches, exploiting the large number of crystallographic protein structures determined to date, have been reported in the literature. To the best of our knowledge, no proper evaluation of such approaches is, however, reported yet. In the present work, we extensively assessed the reliability of docking-based target-fishing strategies. For this purpose, a set of X-ray structures belonging to different targets was selected, and a dataset of compounds, including 10 experimentally active ligands for each target, was created. A target-fishing benchmark database was then obtained, and used to assess the performance of 13 different docking procedures, in identifying the correct target of the dataset ligands. Moreover, a consensus docking-based target-fishing strategy was developed and evaluated. The analysis highlighted that specific features of the target proteins could affect the reliability of the protocol, which however, proved to represent a valuable tool in the proper applicability domain. Our study represents the first extensive performance assessment of docking-based target-fishing approaches, paving the way for the development of novel efficient receptor-based target fishing strategies.

Docking-based inverse virtual screening: methods, applications, and challenges

Identifying potential protein targets for a small-compound ligand query is crucial to the process of drug development. However, there are tens of thousands of proteins in human alone, and it is almost impossible to scan all the existing proteins for a query ligand using current experimental methods. Recently, a computational technology called docking-based inverse virtual screening (IVS) has attracted much attention. In docking-based IVS, a panel of proteins is screened by a molecular docking program to identify potential targets for a query ligand. Ever since the first paper describing a docking-based IVS program was published about a decade ago, the approach has been gradually improved and utilized for a variety of purposes in the field of drug discovery. In this article, the methods employed in docking-based IVS are reviewed in detail, including target databases, docking engines, and scoring function methodologies. Several web servers developed for non-expert users are also reviewed. Then, a number of applications are presented according to different research purposes, such as target identification, side effects/toxicity, drug repositioning, drug-target network development, and receptor design. The review concludes by discussing the challenges that docking-based IVS needs to overcome to become a robust tool for pharmaceutical engineering.

Computational Approaches for the Discovery of Novel Hepatitis C Virus NS3/4A and NS5B Inhibitors

Nonstructural 5B (NS5B) polymerase and Nonstructural 3/4A (NS3/4A) protease have proven to be promising targets for the development of anti-HCV (Hepatitis C Virus) agents. The NS5B polymerase is of paramount importance in HCV viral replication; therefore, employing NS5B inhibitors was considered an effective way for the treatment of HCV. Identifying inhibitors against NS3/4A serine protease represents another attractive approach applied in anti-HCV drug discovery, which is evidenced by its crucial role of in the biogenesis of the viral replication activity. In this chapter, many different computational approaches including Quantitative Structure-Activity Relationship (QSAR) and virtual screening in anti-HCV drug discovery were considered and discussed in detail. Virtual Screening (VS) techniques, including ligand-based and structure-based, and QSAR have been utilized for the discovery of NS5B inhibitors. Moreover, using various in silico protocols and workflows, a number of studies have been conducted with an aim of identifying potential NS3/4A blockage agents.

Molecular Docking at a Glance

The current chapter introduces different aspects of molecular docking technique in order to give an overview to the readers about the topics which will be dealt with throughout this volume. Like many other fields of science, molecular docking studies has experienced a lagging period of slow and steady increase in terms of acquiring attention of scientific community as well as its frequency of application, followed by a pronounced era of exponential expansion in theory, methodology, areas of application and performance due to developments in related technologies such as computational resources and theoretical as well as experimental biophysical methods. In the following sections the evolution of molecular docking will be reviewed and its different components including methods, search algorithms, scoring functions, validation of the methods, and area of applications plus few case studies will be touched briefly.

Reliability analysis and optimization of the consensus docking approach for the development of virtual screening studies

Ligand-protein docking is one of the most common techniques used in virtual screening campaigns. Despite the large number of docking software available, there is still the need of improving the efficacy of docking-based virtual screenings. To date, only very few studies evaluated the possibility of combining the results of different docking methods to achieve higher success rates in virtual screening studies (consensus docking). In order to better understand the range of applicability of this approach, we carried out an extensive enriched database analysis using the DUD dataset. The consensus docking protocol was then refined by applying modifications concerning the calculation of pose consensus and the combination of docking methods included in the procedure. The results obtained suggest that this approach performs as well as the best available methods found in literature, confirming the idea that this procedure can be profitably used for the identification of new hit compounds.

Using reverse docking for target identification and its applications for drug discovery

INTRODUCTION

In contrast to traditional molecular docking, inverse or reverse docking is used for identifying receptors for a given ligand among a large number of receptors. Reverse docking can be used to discover new targets for existing drugs and natural compounds, explain polypharmacology and the molecular mechanism of a substance, find alternative indications of drugs through drug repositioning, and detecting adverse drug reactions and drug toxicity.

AREAS COVERED

In this review, the authors examine how reverse docking methods have evolved over the past fifteen years and how they have been used for target identification and related applications for drug discovery. They discuss various aspects of target databases, reverse docking tools and servers.

EXPERT OPINION

There are several issues related to reverse docking methods such as target structure dataset construction, computational efficiency, how to include receptor flexibility, and most importantly, how to properly normalize the docking scores. In order for reverse docking to become a truly useful tool for the drug discovery, these issues need to be adequately resolved.

Building a virtual ligand screening pipeline using free software: a survey

Virtual screening, the search for bioactive compounds via computational methods, provides a wide range of opportunities to speed up drug development and reduce the associated risks and costs. While virtual screening is already a standard practice in pharmaceutical companies, its applications in preclinical academic research still remain under-exploited, in spite of an increasing availability of dedicated free databases and software tools. In this survey, an overview of recent developments in this field is presented, focusing on free software and data repositories for screening as alternatives to their commercial counterparts, and outlining how available resources can be interlinked into a comprehensive virtual screening pipeline using typical academic computing facilities. Finally, to facilitate the set-up of corresponding pipelines, a downloadable software system is provided, using platform virtualization to integrate pre-installed screening tools and scripts for reproducible application across different operating systems.

Molecular docking and structure-based drug design strategies

Pharmaceutical research has successfully incorporated a wealth of molecular modeling methods, within a variety of drug discovery programs, to study complex biological and chemical systems. The integration of computational and experimental strategies has been of great value in the identification and development of novel promising compounds. Broadly used in modern drug design, molecular docking methods explore the ligand conformations adopted within the binding sites of macromolecular targets. This approach also estimates the ligand-receptor binding free energy by evaluating critical phenomena involved in the intermolecular recognition process. Today, as a variety of docking algorithms are available, an understanding of the advantages and limitations of each method is of fundamental importance in the development of effective strategies and the generation of relevant results. The purpose of this review is to examine current molecular docking strategies used in drug discovery and medicinal chemistry, exploring the advances in the field and the role played by the integration of structure- and ligand-based methods.

Crystallographic analysis and structure-guided engineering of NADPH-dependent Ralstonia sp. alcohol dehydrogenase toward NADH cosubstrate specificity

The NADP⁺-dependent alcohol dehydrogenase from Ralstonia sp. (RasADH) belongs to the protein superfamily of short-chain dehydrogenases/reductases (SDRs). As an enzyme that accepts different types of substrates--including bulky-bulky as well as small-bulky secondary alcohols or ketones--with high stereoselectivity, it offers potential as a biocatalyst for industrial biotechnology. To understand substrate and cosubstrate specificities of RasADH we determined the crystal structure of the apo-enzyme as well as its NADP⁺-bound state with resolutions down to 2.8 Å. RasADH displays a homotetrameric quaternary structure that can be described as a dimer of homodimers while in each subunit a seven-stranded parallel β-sheet, flanked by three α-helices on each side, forms a Rossmann fold-type dinucleotide binding domain. Docking of the well-known substrate (S)-1-phenylethanol clearly revealed the structural determinants of stereospecificity. To favor practical RasADH application in the context of established cofactor recycling systems, for example, those involving an NADH-dependent amino acid dehydrogenase, we attempted to rationally change its cosubstrate specificity from NADP⁺ to NAD⁺ utilizing the structural information that NADP⁺ specificity is largely governed by the residues Asn15, Gly37, Arg38, and Arg39. Furthermore, an extensive sequence alignment with homologous dehydrogenases that have different cosubstrate specificities revealed a modified general SDR motif ASNG (instead of NNAG) at positions 86-89 of RasADH. Consequently, we constructed mutant enzymes with one (G37D), four (N15G/G37D/R38V/R39S), and six (N15G/G37D/R38V/R39S/A86N/S88A) amino acid exchanges. RasADH (N15G/G37D/R38V/R39S) was better able to accept NAD⁺ while showing much reduced catalytic efficiency with NADP⁺, leading to a change in NADH/NADPH specificity by a factor of ∼3.6 million.

Sanjeevini: a freely accessible web-server for target directed lead molecule discovery

Background

Computational methods utilizing the structural and functional information help to understand specific molecular recognition events between the target biomolecule and candidate hits and make it possible to design improved lead molecules for the target.

Results

Sanjeevini represents a massive on-going scientific endeavor to provide to the user, a freely accessible state of the art software suite for protein and DNA targeted lead molecule discovery. It builds in several features, including automated detection of active sites, scanning against a million compound library for identifying hit molecules, all atom based docking and scoring and various other utilities to design molecules with desired affinity and specificity against biomolecular targets. Each of the modules is thoroughly validated on a large dataset of protein/DNA drug targets.

Conclusions

The article presents Sanjeevini, a freely accessible user friendly web-server, to aid in drug discovery. It is implemented on a tera flop cluster and made accessible via a web-interface at http://www.scfbio-iitd.res.in/sanjeevini/sanjeevini.jsp. A brief description of various modules, their scientific basis, validation, and how to use the server to develop in silico suggestions of lead molecules is provided.

AADS--an automated active site identification, docking, and scoring protocol for protein targets based on physicochemical descriptors

We report here a robust automated active site detection, docking, and scoring (AADS) protocol for proteins with known structures. The active site finder identifies all cavities in a protein and scores them based on the physicochemical properties of functional groups lining the cavities in the protein. The accuracy realized on 620 proteins with sizes ranging from 100 to 600 amino acids with known drug active sites is 100% when the top ten cavity points are considered. These top ten cavity points identified are then submitted for an automated docking of an input ligand/candidate molecule. The docking protocol uses an all atom energy based Monte Carlo method. Eight low energy docked structures corresponding to different locations and orientations of the candidate molecule are stored at each cavity point giving 80 docked structures overall which are then ranked using an effective free energy function and top five structures are selected. The predicted structure and energetics of the complexes agree quite well with experiment when tested on a data set of 170 protein-ligand complexes with known structures and binding affinities. The AADS methodology is implemented on an 80 processor cluster and presented as a freely accessible, easy to use tool at http://www.scfbio-iitd.res.in/dock/ActiveSite_new.jsp .

Molecular library design using multi-objective optimization methods

Advancements in combinatorial chemistry and high-throughput screening technology have enabled the synthesis and screening of large molecular libraries for the purposes of drug discovery. Contrary to initial expectations, the increase in screening library size, typically combined with an emphasis on compound structural diversity, did not result in a comparable increase in the number of promising hits found. In an effort to improve the likelihood of discovering hits with greater optimization potential, more recent approaches attempt to incorporate additional knowledge to the library design process to effectively guide the search. Multi-objective optimization methods capable of taking into account several chemical and biological criteria have been used to design collections of compounds satisfying simultaneously multiple pharmaceutically relevant objectives. In this chapter, we present our efforts to implement a multi-objective optimization method, MEGALib, custom-designed to the library design problem. The method exploits existing knowledge, e.g. from previous biological screening experiments, to identify and profile molecular fragments used subsequently to design compounds compromising the various objectives.

A new Lamarckian genetic algorithm for flexible ligand-receptor docking

We present a Lamarckian genetic algorithm (LGA) variant for flexible ligand-receptor docking which allows to handle a large number of degrees of freedom. Our hybrid method combines a multi-deme LGA with a recently published gradient-based method for local optimization of molecular complexes. We compared the performance of our new hybrid method to two non gradient-based search heuristics on the Astex diverse set for flexible ligand-receptor docking. Our results show that the novel approach is clearly superior to other LGAs employing a stochastic optimization method. The new algorithm features a shorter run time and gives substantially better results, especially with increasing complexity of the ligands. Thus, it may be used to dock ligands with many rotatable bonds with high efficiency.

Structure-Based Approaches to Target Fishing and Ligand Profiling

Chemogenomics is an emerging interdisciplinary field aiming at identifying all possible ligands of all possible targets. If one groups targets in columns and ligands in rows, chemogenomic approaches to drug discovery just fill the interaction matrix. Since experimental data do not suffice, several computational methods are currently actively developed to supplement time-consuming and costly experiments. They are either designed to fill rows and thus profile a ligand towards a heterogeneous set of targets (target profiling) or to fill columns and thus identify novel ligands for an existing target (standard virtual screening). At the interface of both strategies are now true chemogenomic computational methods filling well defined areas in the matrix. The present review will focus on (protein) structure-based approaches and illustrates major advances in this novel exciting field which is supposed to massively impact rational drug design in the next decade.

Further analysis and comparative study of intermolecular interactions using dimers from the S22 database

Accurate MP2 and CCSD(T) complete basis set (CBS) interaction energy curves (14 points for each curve) have been obtained for 20 of the dimers reported in the S22 set and analytical Morse curves have been fitted that can be used in developing updated density functional theory (DFT) and force field models. The magnitude and the effect of the basis set superposition error (BSSE) were carefully investigated. We found that going up to aug-cc-pVDZ and aug-cc-pVTZ basis sets is enough to obtain accurate CBS MP2 energies when BSSE corrected values are used but aug-cc-pVTZ and aug-cc-pVQZ basis sets are needed when the BSSE uncorrected total energies are used in CBS extrapolations. MP2 interaction energies with smaller basis sets such as 6-31G* are found to have very little dispersion energy and that the true source of dispersion attributed attractive interactions is almost entirely due to BSSE. MP2 and CCSD(T) CBS interaction energies are found to be very close to one another if aromatic systems are not involved. Comparative analyses have been performed with semiempirical and ab initio methods utilizing the moderate in size but affordable 6-31G* basis set both of which can be readily applied to macromolecular systems. The new M06-2X and M06-L DFT functionals were found to be more accurate than all methods tested herein. Interaction energy curves using the SG1 grid showed discontinuities for several of the dimer systems but this problem disappeared when finer DFT numerical grids were used.

De novo drug design using multiobjective evolutionary graphs

Drug discovery and development is a complex, lengthy process, and failure of a candidate molecule can occur as a result of a combination of reasons, such as poor pharmacokinetics, lack of efficacy, or toxicity. Successful drug candidates necessarily represent a compromise between the numerous, sometimes competing objectives so that the benefits to patients outweigh potential drawbacks and risks. De novo drug design involves searching an immense space of feasible, druglike molecules to select those with the highest chances of becoming drugs using computational technology. Traditionally, de novo design has focused on designing molecules satisfying a single objective, such as similarity to a known ligand or an interaction score, and ignored the presence of the multiple objectives required for druglike behavior. Recently, methods have appeared in the literature that attempt to design molecules satisfying multiple predefined objectives and thereby produce candidate solutions with a higher chance of serving as viable drug leads. This paper describes the Multiobjective Evolutionary Graph Algorithm (MEGA), a new multiobjective optimization de novo design algorithmic framework that can be used to design structurally diverse molecules satisfying one or more objectives. The algorithm combines evolutionary techniques with graph-theory to directly manipulate graphs and perform an efficient global search for promising solutions. In the Experimental Section we present results from the application of MEGA for designing molecules that selectively bind to a known pharmaceutical target using the ChillScore interaction score family. The primary constraints applied to the design are based on the identified structure of the protein target and a known ligand currently marketed as a drug. A detailed explanation of the key elements of the specific implementation of the algorithm is given, including the methods for obtaining molecular building blocks, evolving the chemical graphs, and scoring the designed molecules. Our findings demonstrate that MEGA can produce structurally diverse candidate molecules representing a wide range of compromises of the supplied constraints and thus can be used as an "idea generator" to support expert chemists assigned with the task of molecular design.

Lead finder: an approach to improve accuracy of protein-ligand docking, binding energy estimation, and virtual screening

An innovative molecular docking algorithm and three specialized high accuracy scoring functions are introduced in the Lead Finder docking software. Lead Finder's algorithm for ligand docking combines the classical genetic algorithm with various local optimization procedures and resourceful exploitation of the knowledge generated during docking process. Lead Finder's scoring functions are based on a molecular mechanics functional which explicitly accounts for different types of energy contributions scaled with empiric coefficients to produce three scoring functions tailored for (a) accurate binding energy predictions; (b) correct energy-ranking of docked ligand poses; and (c) correct rank-ordering of active and inactive compounds in virtual screening experiments. The predicted values of the free energy of protein-ligand binding were benchmarked against a set of experimentally measured binding energies for 330 diverse protein-ligand complexes yielding rmsd of 1.50 kcal/mol. The accuracy of ligand docking was assessed on a set of 407 structures, which included almost all published test sets of the following programs: FlexX, Glide SP, Glide XP, Gold, LigandFit, MolDock, and Surflex. rmsd of 2 A or less was observed for 80-96% of the structures in the test sets (80.0% on the Glide XP and FlexX test sets, 96.0% on the Surflex and MolDock test sets). The ability of Lead Finder to distinguish between active and inactive compounds during virtual screening experiments was benchmarked against 34 therapeutically relevant protein targets. Impressive enrichment factors were obtained for almost all of the targets with the average area under receiver operator curve being equal to 0.92.

Towards the development of universal, fast and highly accurate docking/scoring methods: a long way to go

Accelerating the drug discovery process requires predictive computational protocols capable of reducing or simplifying the synthetic and/or combinatorial challenge. Docking-based virtual screening methods have been developed and successfully applied to a number of pharmaceutical targets. In this review, we first present the current status of docking and scoring methods, with exhaustive lists of these. We next discuss reported comparative studies, outlining criteria for their interpretation. In the final section, we describe some of the remaining developments that would potentially lead to a universally applicable docking/scoring method.

Similarity based docking

We have recently introduced GMA, a highly efficient method for flexible molecular alignment. Here we show how this approach can be used to improve docking accuracy and efficiency, in cases where a complex structure of a ligand with the target protein is known. In cases where a known ligand exists, yet the complex structure is unknown it is possible to make use of the advantages offered by this approach, by combining it with standard ligand docking.