Computer-aided drug design at Boehringer Ingelheim

From Proline to Chlorantraniliprole Mimics: Computer-Aided Design, Simple Preparation, and Excellent Insecticidal Profiles

Chlorantraniliprole (CHL), a favored agricultural insecticide, is renowned for its high efficiency and broad-spectrum effectiveness against lepidoptera insects. However, the urgency for new insecticide development is underscored by the intricate multistep preparation process and modest overall yields of CHL, along with the escalating challenge of insect resistance. In response, we have crafted CHL mimics from proline employing computer-aided drug design. Molecular docking analysis of CHL's interactions with the ryanodine receptor (RyR) revealed that the nitrogen atom within the pyrazole moiety does not engage in pivotal interactions. Its removal may not abolish bioactivity entirely but could substantially simplify the synthetic process, thereby enhancing atom economy. This revelation prompted the exclusion of nitrogen and the subsequent formation of a pyrrole ring, enabling the meticulous design of synthetic pathways characterized by cost-effective precursors, streamlined synthesis, the avoidance of toxic reagents, minimal instrumentation, and high yields in the pursuit of innovative RyR modulators. Among these modulators, A1 and B1, obtained with yields exceeding 60%, showcased exceptional insecticidal potency, with LC50 values spanning from 0.12 to 1.47 mg L-1 against P. xylostella and M. separate. The inhibitory effects of these two compounds on insect detoxification enzymes imply a reduced likelihood of eliciting resistance in comparison to CHL, a finding further corroborated by their insecticidal potency against resistant pests. Moreover, molecular docking, MD simulations, and DFT calculations provided valuable structural insights, potentially unraveling the superior insecticidal activity of these two molecules, and thus paving the way for developing more potent insecticides.

Development and Prospects of Furin Inhibitors for Therapeutic Applications

Furin, a serine protease enzyme located in the Golgi apparatus of animal cells, plays a crucial role in cleaving precursor proteins into their mature, active forms. It is ubiquitously expressed across various tissues, including the brain, lungs, gastrointestinal tract, liver, pancreas, and reproductive organs. Since its discovery in 1990, furin has been recognized as a significant therapeutic target, leading to the active development of furin inhibitors for potential use in antiviral, antibacterial, anticancer, and other therapeutic applications. This review provides a comprehensive overview of the progress in the development and characterization of furin inhibitors, encompassing peptides, linear and macrocyclic peptidomimetics, and non-peptide compounds, highlighting their potential in the treatment of both infectious and non-infectious diseases.

Current perspectives and trend of computer-aided drug design: a review and bibliometric analysis

AIM

Computer-aided drug design (CADD) is a drug design technique for computing ligand-receptor interactions and is involved in various stages of drug development. To better grasp the frontiers and hotspots of CADD, we conducted a review analysis through bibliometrics.

METHODS

A systematic review of studies published between 2000 and 20 July 2023 was conducted following the PRISMA guidelines. Literature on CADD was selected from the Web of Science Core Collection. General information, publications, output trends, countries/regions, institutions, journals, keywords, and influential authors were visually analyzed using software such as Excel, VOSviewer, RStudio, and CiteSpace.

RESULTS

A total of 2031 publications were included. These publications primarily originated from 99 countries or regions led by the U.S. and China. Among the contributors, MacKerell AD had the highest number of articles and the greatest influence. The Journal of Medicinal Chemistry was the most cited journal, whereas the Journal of Chemical Information and Modeling had the highest number of publications.

CONCLUSIONS

Influential authors in the field were identified. Current research shows active collaboration between countries, institutions, and companies. CADD technologies such as homology modeling, pharmacophore modeling, quantitative conformational relationships, molecular docking, molecular dynamics simulation, binding free energy prediction, and high-throughput virtual screening can effectively improve the efficiency of new drug discovery. Artificial intelligence-assisted drug design and screening based on CADD represent key topics that will influence future development. Furthermore, this paper will be helpful in better understanding the frontiers and hotspots of CADD.

DockOpt: A Tool for Automatic Optimization of Docking Models

Molecular docking is a widely used technique for leveraging protein structure for ligand discovery, but it remains difficult to utilize due to limitations that have not been adequately addressed. Despite some progress toward automation, docking still requires expert guidance, hindering its adoption by a broader range of investigators. To make docking more accessible, we developed a new utility called DockOpt, which automates the creation, evaluation, and optimization of docking models prior to their deployment in large-scale prospective screens. DockOpt outperforms our previous automated pipeline across all 43 targets in the DUDE-Z benchmark data set, and the generated models for 84% of targets demonstrate sufficient enrichment to warrant their use in prospective screens, with normalized LogAUC values of at least 15%. DockOpt is available as part of the Python package Pydock3 included in the UCSF DOCK 3.8 distribution, which is available for free to academic researchers at https://dock.compbio.ucsf.edu and free for everyone upon registration at https://tldr.docking.org.

Predictive Models Based on Molecular Images and Molecular Descriptors for Drug Screening

Various toxicity and pharmacokinetic evaluations as screening experiments are needed at the drug discovery stage. Currently, to reduce the use of animal experiments and developmental expenses, the development of high-performance predictive models based on quantitative structure-activity relationship analysis is desired. From these evaluation targets, we selected 50% lethal dose (LD50), blood-brain barrier penetration (BBBP), and the clearance (CL) pathway for this investigation and constructed predictive models for each target using 636-11,886 compounds. First, we constructed predictive models using the DeepSnap-deep learning (DL) method and images of compounds as features. The calculated area under the curve (AUC) and balanced accuracy (BAC) were, respectively, 0.887 and 0.818 for LD50, 0.893 and 0.824 for BBBP, and 0.883 and 0.763 for the CL pathway. Next, molecular descriptors (MDs) of compounds were calculated using Molecular Operating Environment, alvaDesc, and ADMET Predictor to construct predictive models using the MD-based method. Using these MDs, we constructed predictive models using DataRobot. The calculated AUC and BAC were, respectively, 0.931 and 0.805 for LD50, 0.919 and 0.849 for BBBP, and 0.900 and 0.807 for the CL pathway. In this investigation, we constructed predictive models combining the DeepSnap-DL and MD-based methods. In ensemble models using the mean predictive probability of the DeepSnap-DL and MD-based methods, the calculated AUC and BAC were, respectively, 0.942 and 0.842 for LD50, 0.936 and 0.853 for BBBP, and 0.908 and 0.832 for the CL pathway, with improved predictive performance observed for all variables compared with either single method alone. Moreover, in consensus models that adopted only compounds for which the results of the two methods agreed, the calculated BAC for LD50, BBBP, and the CL pathway were 0.916, 0.918, and 0.847, respectively, indicating higher predictive performance than the ensemble models for all three variables. The predictive models combining the DeepSnap-DL and MD-based methods displayed high predictive performance for LD50, BBBP, and the CL pathway. Therefore, the application of this approach to prediction targets in various drug discovery screenings is expected to accelerate drug discovery.

Exploring the Molecular Targets for the Antidepressant and Antisuicidal Effects of Ketamine Enantiomers by Using Network Pharmacology and Molecular Docking

Ketamine, a racemic mixture of esketamine (S-ketamine) and arketamine (R-ketamine), has received particular attention for its rapid antidepressant and antisuicidal effects. NMDA receptor inhibition has been indicated as one of the main mechanisms of action of the racemic mixture, but other pharmacological targets have also been proposed. This study aimed to explore the possible multiple targets of ketamine enantiomers related to their antidepressant and antisuicidal effects. To this end, targets were predicted using Swiss Target Prediction software for each ketamine enantiomer. Targets related to depression and suicide were collected by the Gene Cards database. The intersections of targets were analyzed using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). Network pharmacology analysis was performed using Gene Mania and Cytoscape software. Molecular docking was used to predict the main targets of the network. The results indicated that esketamine and arketamine share some biological targets, particularly NMDA receptor and phosphodiesterases 3A, 7A, and 5A but have specific molecular targets. While esketamine is predicted to interact with the GABAergic system, arketamine may interact with macrophage migration inhibitory factor (MIF). Both ketamine enantiomers activate neuroplasticity-related signaling pathways and show addiction potential. Our results identified novel, poorly explored molecular targets that may be related to the beneficial effects of esketamine and arketamine against depression and suicide.

CADD, AI and ML in drug discovery: A comprehensive review

Computer-aided drug design (CADD) is an emerging field that has drawn a lot of interest because of its potential to expedite and lower the cost of the drug development process. Drug discovery research is expensive and time-consuming, and it frequently took 10-15 years for a drug to be commercially available. CADD has significantly impacted this area of research. Further, the combination of CADD with Artificial Intelligence (AI), Machine Learning (ML), and Deep Learning (DL) technologies to handle enormous amounts of biological data has reduced the time and cost associated with the drug development process. This review will discuss how CADD, AI, ML, and DL approaches help identify drug candidates and various other steps of the drug discovery process. It will also provide a detailed overview of the different in silico tools used and how these approaches interact.

Prioritizing Small Sets of Molecules for Synthesis through in-silico Tools: A Comparison of Common Ranking Methods

Prioritizing molecules for synthesis is a key role of computational methods within medicinal chemistry. Multiple tools exist for ranking molecules, from the cheap and popular molecular docking methods to more computationally expensive molecular-dynamics (MD)-based methods. It is often questioned whether the accuracy of the more rigorous methods justifies the higher computational cost and associated calculation time. Here, we compared the performance on ranking the binding of small molecules for seven scoring functions from five docking programs, one end-point method (MM/GBSA), and two MD-based free energy methods (PMX, FEP+). We investigated 16 pharmaceutically relevant targets with a total of 423 known binders. The performance of docking methods for ligand ranking was strongly system dependent. We observed that MD-based methods predominantly outperformed docking algorithms and MM/GBSA calculations. Based on our results, we recommend the application of MD-based free energy methods for prioritization of molecules for synthesis in lead optimization, whenever feasible.

Plant Metabolites as SARS-CoV-2 Inhibitors Candidates: In Silico and In Vitro Studies

Since it acquired pandemic status, SARS-CoV-2 has been causing all kinds of damage all over the world. More than 6.3 million people have died, and many cases of sequelae are in survivors. Currently, the only products available to most of the world’s population to fight the pandemic are vaccines, which still need improvement since the number of new cases, admissions into intensive care units, and deaths are again reaching worrying rates, which makes it essential to compounds that can be used during infection, reducing the impacts of the disease. Plant metabolites are recognized sources of diverse biological activities and are the safest way to research anti-SARS-CoV-2 compounds. The present study computationally evaluated 55 plant compounds in five SARS-CoV-2 targets such Main Protease (Mpro or 3CL or MainPro), RNA-dependent RNA polymerase (RdRp), Papain-Like Protease (PLpro), NSP15 Endoribonuclease, Spike Protein (Protein S or Spro) and human Angiotensin-converting enzyme 2 (ACE-2) followed by in vitro evaluation of their potential for the inhibition of the interaction of the SARS-CoV-2 Spro with human ACE-2. The in silico results indicated that, in general, amentoflavone, 7-O-galloylquercetin, kaempferitrin, and gallagic acid were the compounds with the strongest electronic interaction parameters with the selected targets. Through the data obtained, we can demonstrate that although the indication of individual interaction of plant metabolites with both Spro and ACE-2, the metabolites evaluated were not able to inhibit the interaction between these two structures in the in vitro test. Despite this, these molecules still must be considered in the research of therapeutic agents for treatment of patients affected by COVID-19 since the activity on other targets and influence on the dynamics of viral infection during the interaction Spro x ACE-2 should be investigated.

Revisiting the therapeutic potential of tocotrienol

The therapeutic potential of the tocotrienol group stems from its nutraceutical properties as a dietary supplement. It is largely considered to be safe when consumed at low doses for attenuating pathophysiology as shown by animal models, in vitro assays, and ongoing human trials. Medical researchers and the allied sciences have experimented with tocotrienols for many decades, but its therapeutic potential was limited to adjuvant or concurrent treatment regimens. Recent studies have focused on targeted drug delivery by enhancing the bioavailability through carriers, self-sustained emulsions, nanoparticles, and ethosomes. Epigenetic modulation and computer remodeling are other means that will help increase chemosensitivity. This review will focus on the systemic intracellular anti-cancer, antioxidant, and anti-inflammatory mechanisms that are stimulated and/or regulated by tocotrienols while highlighting its potent therapeutic properties in a diverse group of clinical diseases.

Novel QSAR Approach for a Regression Model of Clearance That Combines DeepSnap-Deep Learning and Conventional Machine Learning

The toxicity, absorption, distribution, metabolism, and excretion properties of some targets are difficult to predict by quantitative structure-activity relationship analysis. Therefore, there is a need for a new prediction method that performs well for these targets. The aim of this study was to develop a new regression model of rat clearance (CL). We constructed a regression model using 1545 in-house compounds for which we had rat CL data. Molecular descriptors were calculated using molecular operating environment, alvaDesc, and ADMET Predictor software. The classification model of DeepSnap and Deep Learning (DeepSnap-DL) with images of the three-dimensional chemical structures of compounds as features was constructed, and the prediction probabilities for each compound were calculated. For molecular descriptor-based methods that use molecular descriptors and conventional machine learning algorithms selected by DataRobot, the correlation coefficient (R ²) and root mean square error (RMSE) were 0.625-0.669 and 0.295-0.318, respectively. We combined molecular descriptors and prediction probability of DeepSnap-DL as features and developed a novel regression method we called the combination model. In the combination model with these two types of features and conventional algorithms selected by DataRobot, R ² and RMSE were 0.710-0.769 and 0.247-0.278, respectively. This finding shows that the combination model performed better than molecular descriptor-based methods. Our combination model will contribute to the design of more rational compounds for drug discovery. This method may be applicable not only to rat CL but also to other pharmacokinetic and pharmacological activity and toxicity parameters; therefore, applying it to other parameters may help to accelerate drug discovery.

Molecular Modeling Techniques Applied to the Design of Multitarget Drugs: Methods and Applications

Multifactorial diseases, such as cancer and diabetes present a challenge for the traditional "one-target, one disease" paradigm due to their complex pathogenic mechanisms. Although a combination of drugs can be used, a multitarget drug may be a better choice face of its efficacy, lower adverse effects and lower chance of resistance development. The computer-based design of these multitarget drugs can explore the same techniques used for single-target drug design, but the difficulties associated to the obtention of drugs that are capable of modulating two or more targets with similar efficacy impose new challenges, whose solutions involve the adaptation of known techniques and also to the development of new ones, including machine-learning approaches. In this review, some SBDD and LBDD techniques for the multitarget drug design are discussed, together with some cases where the application of such techniques led to effective multitarget ligands.

SeamDock: An Interactive and Collaborative Online Docking Resource to Assist Small Compound Molecular Docking

In silico assessment of protein receptor interactions with small ligands is now part of the standard pipeline for drug discovery, and numerous tools and protocols have been developed for this purpose. With the SeamDock web server, we propose a new approach to facilitate access to small molecule docking for nonspecialists, including students. The SeamDock online service integrates different docking tools in a common framework that allows ligand global and/or local docking and a hierarchical approach combining the two for easy interaction site identification. This service does not require advanced computer knowledge, and it works without the installation of any programs with the exception of a common web browser. The use of the Seamless framework linking the RPBS calculation server to the user’s browser allows the user to navigate smoothly and interactively on the SeamDock web page. A major effort has been put into the 3D visualization of ligand, receptor, and docking poses and their interactions with the receptor. The advanced visualization features combined with the seamless library allow a user to share with an unlimited number of collaborators, a docking session, and its full visualization states. As a result, SeamDock can be seen as a free, simple, didactic, evolving online docking resource best suited for education and training.

Shedding Light on the Molecular Recognition of Sub-Kilodalton Macrocyclic Peptides on Thrombin by Supervised Molecular Dynamics

Macrocycles are attractive structures for drug development due to their favorable structural features, potential in binding to targets with flat featureless surfaces, and their ability to disrupt protein–protein interactions. Moreover, large novel highly diverse libraries of low-molecular-weight macrocycles with therapeutically favorable characteristics have been recently established. Considering the mentioned facts, having a validated, fast, and accurate computational protocol for studying the molecular recognition and binding mode of this interesting new class of macrocyclic peptides deemed to be helpful as well as insightful in the quest of accelerating drug discovery. To that end, the ability of the in-house supervised molecular dynamics protocol called SuMD in the reproduction of the X-ray crystallography final binding state of a macrocyclic non-canonical tetrapeptide—from a novel library of 8,988 sub-kilodalton macrocyclic peptides—in the thrombin active site was successfully validated. A comparable binding mode with the minimum root-mean-square deviation (RMSD) of 1.4 Å at simulation time point 71.6 ns was achieved. This method validation study extended the application domain of the SuMD sampling method for computationally cheap, fast but accurate, and insightful macrocycle–protein molecular recognition studies.

Prediction Model of Clearance by a Novel Quantitative Structure-Activity Relationship Approach, Combination DeepSnap-Deep Learning and Conventional Machine Learning

Some targets predicted by machine learning (ML) in drug discovery remain a challenge because of poor prediction. In this study, a new prediction model was developed and rat clearance (CL) was selected as a target because it is difficult to predict. A classification model was constructed using 1545 in-house compounds with rat CL data. The molecular descriptors calculated by Molecular Operating Environment (MOE), alvaDesc, and ADMET Predictor software were used to construct the prediction model. In conventional ML using 100 descriptors and random forest selected by DataRobot, the area under the curve (AUC) and accuracy (ACC) were 0.883 and 0.825, respectively. Conversely, the prediction model using DeepSnap and Deep Learning (DeepSnap-DL) with compound features as images had AUC and ACC of 0.905 and 0.832, respectively. We combined the two models (conventional ML and DeepSnap-DL) to develop a novel prediction model. Using the ensemble model with the mean of the predicted probabilities from each model improved the evaluation metrics (AUC = 0.943 and ACC = 0.874). In addition, a consensus model using the results of the agreement between classifications had an increased ACC (0.959). These combination models with a high level of predictive performance can be applied to rat CL as well as other pharmacokinetic parameters, pharmacological activity, and toxicity prediction. Therefore, these models will aid in the design of more rational compounds for the development of drugs.

Therapeutic Potential of Curcumin as an Antimycobacterial Agent

Curcumin is the principal curcuminoid obtained from the plant Curcuma longa and has been extensively studied for its biological and chemical properties. Curcumin displays a vast range of pharmacological properties, including antimicrobial, anti-inflammatory, antioxidant, and antitumor activity. Specifically, curcumin has been linked to the improvement of the outcome of tuberculosis. There are many reviews on the pharmacological effects of curcumin; however, reviews of the antitubercular activity are comparatively scarcer. In this review, we attempt to discuss the different aspects of the research on the antitubercular activity of curcumin. These include antimycobacterial activity, modulation of the host immune response, and enhancement of BCG vaccine efficacy. Recent advances in the antimycobacterial activity of curcumin synthetic derivatives, the role of computer aided drug design in identifying curcumin targets, the hepatoprotective role of curcumin, and the dosage and toxicology of curcumin will be discussed. While growing evidence supports the use of curcumin and its derivatives for tuberculosis therapy, further preclinical and clinical investigations are of pivotal importance before recommending the use of curcumin formulations in public health.

Improving de novo Molecule Generation by Embedding LSTM and Attention Mechanism in CycleGAN

The application of deep learning in the field of drug discovery brings the development and expansion of molecular generative models along with new challenges in this field. One of challenges in de novo molecular generation is how to produce new reasonable molecules with desired pharmacological, physical, and chemical properties. To improve the similarity between the generated molecule and the starting molecule, we propose a new molecule generation model by embedding Long Short-Term Memory (LSTM) and Attention mechanism in CycleGAN architecture, LA-CycleGAN. The network layer of the generator in CycleGAN is fused head and tail to improve the similarity of the generated structure. The embedded LSTM and Attention mechanism can overcome long-term dependency problems in treating the normally used SMILES input. From our quantitative evaluation, we present that LA-CycleGAN expands the chemical space of the molecules and improves the ability of structure conversion. The generated molecules are highly similar to the starting compound structures while obtaining expected molecular properties during cycle generative adversarial network learning, which comprehensively improves the performance of the generative model.

Mapping the intellectual structure of the coronavirus field (2000-2020): a co-word analysis

Over the two last decades, coronaviruses have affected human life in different ways, especially in terms of health and economy. Due to the profound effects of novel coronaviruses, growing tides of research are emerging in various research fields. This paper employs a co-word analysis approach to map the intellectual structure of the coronavirus literature for a better understanding of how coronavirus research and the disease itself have developed during the target timeframe. A strategic diagram has been drawn to depict the coronavirus domain's structure and development. A detailed picture of coronavirus literature has been extracted from a huge number of papers to provide a quick overview of the coronavirus literature. The main themes of past coronavirus-related publications are (a) "Antibody-Virus Interactions," (b) "Emerging Infectious Diseases," (c) "Protein Structure-based Drug Design and Antiviral Drug Discovery," (d) "Coronavirus Detection Methods," (e) "Viral Pathogenesis and Immunity," and (f) "Animal Coronaviruses." The emerging infectious diseases are mostly related to fatal diseases (such as Middle East respiratory syndrome, severe acute respiratory syndrome, and COVID-19) and animal coronaviruses (including porcine, turkey, feline, canine, equine, and bovine coronaviruses and infectious bronchitis virus), which are capable of placing animal-dependent industries such as the swine and poultry industries under strong economic pressure. Although considerable research into coronavirus has been done, this unique field has not yet matured sufficiently. Therefore, "Antibody-virus Interactions," "Emerging Infectious Diseases," and "Coronavirus Detection Methods" hold interesting, promising research gaps to be both explored and filled in the future.

Experimental and Computational Approaches to Improve Binding Affinity in Chemical Biology and Drug Discovery

Drug discovery is one of the most complicated processes and establishment of a single drug may require multidisciplinary attempts to design efficient and commercially viable drugs. The main purpose of drug design is to identify a chemical compound or inhibitor that can bind to an active site of a specific cavity on a target protein. The traditional drug design methods involved various experimental based approaches including random screening of chemicals found in nature or can be synthesized directly in chemical laboratories. Except for the long cycle design and time, high cost is also the major issue of concern. Modernized computer-based algorithm including structure-based drug design has accelerated the drug design and discovery process adequately. Surprisingly from the past decade remarkable progress has been made concerned with all area of drug design and discovery. CADD (Computer Aided Drug Designing) based tools shorten the conventional cycle size and also generate chemically more stable and worthy compounds and hence reduce the drug discovery cost. This special edition of editorial comprises the combination of seven research and review articles set emphasis especially on the computational approaches along with the experimental approaches using a chemical synthesizing for the binding affinity in chemical biology and discovery as a salient used in de-novo drug designing. This set of articles exfoliates the role that systems biology and the evaluation of ligand affinity in drug design and discovery for the future.

Deep Learning-Based Ligand Design Using Shared Latent Implicit Fingerprints from Collaborative Filtering

In their previous work, Srinivas et al. [ J. Cheminf. 2018, 10, 56] have shown that implicit fingerprints capture ligands and proteins in a shared latent space, typically for the purposes of virtual screening with collaborative filtering models applied on known bioactivity data. In this work, we extend these implicit fingerprints/descriptors using deep learning techniques to translate latent descriptors into discrete representations of molecules (SMILES), without explicitly optimizing for chemical properties. This allows the design of new compounds based upon the latent representation of nearby proteins, thereby encoding druglike properties including binding affinities to known proteins. The implicit descriptor method does not require any fingerprint similarity search, which makes the method free of any bias arising from the empirical nature of the fingerprint models [Srinivas, R.; J. Cheminf. 2018, 10, 56]. We evaluate the properties of the potentially novel drugs generated by our approach using physical properties of druglike molecules and chemical complexity. Additionally, we analyze the reliability of the biological activity of the new compounds generated using this method by employing models of protein-ligand interaction, which assists in assessing the potential binding affinity of the designed compounds. We find that the generated compounds exhibit properties of chemically feasible compounds and are predicted to be excellent binders to known proteins. Furthermore, we also analyze the diversity of compounds created using the Tanimoto distance and conclude that there is a wide diversity in the generated compounds.

Artificial applicability labels for improving policies in retrosynthesis prediction

Automated retrosynthetic planning algorithms are a research area of increasing importance. Automated reaction-template extraction from large datasets, in conjunction with neural-network-enhanced tree-search algorithms, can find plausible routes to target compounds in seconds. However, the current method for training neural networks to predict suitable templates for a given target product leads to many predictions that are not applicable in silico. Most templates in the top 50 suggested templates cannot be applied to the target molecule to perform the virtual reaction. Here, we describe how to generate data and train a neural network policy that predicts whether templates are applicable or not. First, we generate a massive training dataset by applying each retrosynthetic template to each product from our reaction database. Second, we train a neural network to perform near-perfect prediction of the applicability labels on a held-out test set. The trained network is then joined with a policy model trained to predict and prioritize templates using the labels from the original dataset. The combined model was found to outperform the policy model in a route-finding task using 1700 compounds from our internal drug-discovery projects.

Improved 3D-QSAR Prediction by Multiple Conformational Alignments and Molecular Docking Studies to Design and Discover HIV-I Protease Inhibitors

BACKGROUND

Inhibition of HIV-I protease enzyme is a strategic step for providing better treatment in retrovirus infections, which avoids resistance and possesses less toxicity.

OBJECTIVES

In the course of our research to discover new and potent protease inhibitors, 3D-QSAR (CoMFA and CoMSIA) models were generated using 3 different alignment techniques, including multifit alignment, docking based and Distill based alignment for 63 compounds. Novel molecules were designed from the output of this study.

METHODS

A total of 3 alignment methods were used to generate CoMFA and CoMSIA models. A Distill based alignment method was considered a better method according to different validation parameters. A 3D-QSAR model was generated and contour maps were discussed. The biological activity of designed molecules was predicted using the generated QSAR model to validate QSAR. The newly designed molecules were docked to predict binding affinity.

RESULTS

In CoMFA, leave one out cross-validated coefficient (q²), conventional coefficient (r²) and predicted correlation coefficient (r²_Predicted) values were found to be 0.721, 0.991 and 0.780, respectively. The best obtained CoMSIA model also showed significant cross-validated coefficient (q²), conventional coefficient (r²) and predicted correlation coefficient (r²_Predicted) values of 0.714, 0.987 and 0.721, respectively. Steric and electrostatic contour maps generated from CoMFA and hydrophobic and hydrogen bond donor and hydrogen bond acceptor contour maps from CoMSIA models were used to design new and bioactive protease inhibitors by incorporating bioisosterism and knowledge-based structure-activity relationship.

CONCLUSION

The results from both these approaches, ligand-based drug design and structure-based drug design, are adequate and promising to discover protease inhibitors.

Discovery of Hydroxyamidine Derivatives as Highly Potent, Selective Indoleamine-2,3-dioxygenase 1 Inhibitors

In this study, a series of novel hydroxyamidine derivatives were identified as potent and selective IDO1 inhibitors by structure-based drug design. Among them, compounds 13-15 and 18 exhibited favorable enzymatic and cellular activities. Compound 18 showed improved bioavailability in mouse, rat, and dog (F% = 44%, 58.8%, 102.1%, respectively). With reasonable in vivo pharmacokinetic properties, compound 18 was further evaluated in a transgenic MC38 xenograft mouse model. The combination of compound 18 with PD-1 monoclonal antibody showed a synergistic antitumor effect. These data indicated that compound 18 as a potential cancer immunotherapy agent should warrant further investigation.

Progress of tubulin polymerization activity detection methods

Tubulin, an important target in tumor therapy, is one of the hotspots in the field of antineoplastic drugs in recent years, and it is of great significance to design and screen new inhibitors for this target. Natural products and chemical synthetic drugs are the main sources of tubulin inhibitors. However, due to the variety of compound structure types, it has always been difficult for researchers to screen out polymerization inhibitors with simple operation, high efficiency and low cost. A large number of articles have reported the screening methods of tubulin inhibitors and their biological activity. In this article, the biological activity detection methods of tubulin polymerization inhibitors are reviewed. Thus, it provides a theoretical basis for the further study of tubulin polymerization inhibitors and the selection of methods for tubulin inhibitors.

Accounting for the Central Role of Interfacial Water in Protein-Ligand Binding Free Energy Calculations

Rigorous binding free energy methods in drug discovery are growing in popularity because of a combination of methodological advances, improvements in computer hardware, and workflow automation. These calculations typically use molecular dynamics (MD) to sample from the Boltzmann distribution of conformational states. However, when part or all of the binding sites is inaccessible to the bulk solvent, the time needed for water molecules to equilibrate between bulk solvent and the binding site can be well beyond what is practical with standard MD. This sampling limitation is problematic in relative binding free energy calculations, which compute the reversible work of converting ligand 1 to ligand 2 within the binding site. Thus, if ligand 1 is smaller and/or more polar than ligand 2, the perturbation may allow additional water molecules to occupy a region of the binding site. However, this change in hydration may not be captured by standard MD simulations and may therefore lead to errors in the computed free energy. We recently developed a hybrid Monte Carlo/MD (MC/MD) method, which speeds up the equilibration of water between bulk solvent and buried cavities, while sampling from the intended distribution of states. Here, we report on the use of this approach in the context of alchemical binding free energy calculations. We find that using MC/MD markedly improves the accuracy of the calculations and also reduces hysteresis between the forward and reverse perturbations, relative to matched calculations using only MD with or without the crystallographic water molecules. The present method is available for use in AMBER simulation software.

Docking with AutoDock4

AutoDock is one of the most popular receptor-ligand docking simulation programs. It was first released in the early 1990s and is in continuous development and adapted to specific protein targets. AutoDock has been applied to a wide range of biological systems. It has been used not only for protein-ligand docking simulation but also for the prediction of binding affinity with good correlation with experimental binding affinity for several protein systems. The latest version makes use of a semi-empirical force field to evaluate protein-ligand binding affinity and for selecting the lowest energy pose in docking simulation. AutoDock4.2.6 has an arsenal of four search algorithms to carry out docking simulation including simulated annealing, genetic algorithm, and Lamarckian algorithm. In this chapter, we describe a tutorial about how to perform docking with AutoDock4. We focus our simulations on the protein target cyclin-dependent kinase 2.

User-Friendly Quantum Mechanics: Applications for Drug Discovery

Quantum mechanics (QM) methods provide a fine description of receptor-ligand interactions and of chemical reactions. Their use in drug design and drug discovery is increasing, especially for complex systems including metal ions in the binding sites, for the design of highly selective inhibitors, for the optimization of bi-specific compounds, to understand enzymatic reactions, and for the study of covalent ligands and prodrugs. They are also used for generating molecular descriptors for predictive QSAR/QSPR models and for the parameterization of force fields. Thanks to the continuous increase of computational power offered by GPUs and to the development of sophisticated algorithms, QM methods are becoming part of the standard tools used in computer-aided drug design (CADD). We present the most used QM methods and software packages, and we discuss recent representative applications in drug design and drug discovery.

The In Silico Drug Discovery Toolbox: Applications in Lead Discovery and Optimization

BACKGROUND

Discovery and development of a new drug is a long lasting and expensive journey that takes around 20 years from starting idea to approval and marketing of new medication. Despite R&D expenditures have been constantly increasing in the last few years, the number of new drugs introduced into market has been steadily declining. This is mainly due to preclinical and clinical safety issues, which still represent about 40% of drug discontinuation. To cope with this issue, a number of in silico techniques are currently being used for an early stage evaluation/prediction of potential safety issues, allowing to increase the drug-discovery success rate and reduce costs associated with the development of a new drug.

METHODS

In the present review, we will analyse the early steps of the drug-discovery pipeline, describing the sequence of steps from disease selection to lead optimization and focusing on the most common in silico tools used to assess attrition risks and build a mitigation plan.

RESULTS

A comprehensive list of widely used in silico tools, databases, and public initiatives that can be effectively implemented and used in the drug discovery pipeline has been provided. A few examples of how these tools can be problem-solving and how they may increase the success rate of a drug discovery and development program have been also provided. Finally, selected examples where the application of in silico tools had effectively contributed to the development of marketed drugs or clinical candidates will be given.

CONCLUSION

The in silico toolbox finds great application in every step of early drug discovery: (i) target identification and validation; (ii) hit identification; (iii) hit-to-lead; and (iv) lead optimization. Each of these steps has been described in details, providing a useful overview on the role played by in silico tools in the decision-making process to speed-up the discovery of new drugs.

Investigation of Binding Characteristics of Phosphoinositide-dependent Kinase-1 (PDK1) Co-crystallized Ligands Through Virtual Pharmacophore Modeling Leading to Novel Anti-PDK1 Hits

BACKGROUND

3-Phosphoinositide Dependent Protein Kinase-1 (PDK1) is being lately considered as an attractive and forthcoming anticancer target. A Protein Data Bank (PDB) cocrystallized crystal provides not only rigid theoretical data but also a realistic molecular recognition data that can be explored and used to discover new hits.

OBJECTIVE

This incited us to investigate the co-crystallized ligands' contacts inside the PDK1 binding pocket via a structure-based receptor-ligand pharmacophore generation technique in Discovery Studio 4.5 (DS 4.5).

METHODS

Accordingly, 35 crystals for PDK1 were collected and studied. Every single receptorligand interaction was validated and the significant ones were converted into their corresponding pharmacophoric features. The generated pharmacophores were scored by the Receiver Operating Characteristic (ROC) curve analysis.

RESULTS

Consequently, 169 pharmacophores were generated and sorted, 11 pharmacophores acquired good ROC-AUC results of 0.8 and a selectivity value above 8. Pharmacophore 1UU3_2_01 was used in particular as a searching filter to screen NCI database because of its acceptable validity criteria and its distinctive positive ionizable feature. Several low micromolar PDK1 inhibitors were revealed. The most potent hit illustrated anti-PDK1 IC50 values of 200 nM with 70% inhibition against SW480 cell lines.

CONCLUSION

Eventually, the active hits were docked inside the PDK1 binding pocket and the recognition points between the active hits and the receptor were analyzed that led to the discovery of new scaffolds as potential PDK1 inhibitors.

SAnDReS: A Computational Tool for Docking

Since the early 1980s, we have witnessed considerable progress in the development and application of docking programs to assess protein-ligand interactions. Most of these applications had as a goal the identification of potential new binders to protein targets. Another remarkable progress is taking place in the determination of the structures of protein-ligand complexes, mostly using X-ray diffraction crystallography. Considering these developments, we have a favorable scenario for the creation of a computational tool that integrates into one workflow all steps involved in molecular docking simulations. We had these goals in mind when we developed the program SAnDReS. This program allows the integration of all computational features related to modern docking studies into one workflow. SAnDReS not only carries out docking simulations but also evaluates several docking protocols allowing the selection of the best approach for a given protein system. SAnDReS is a free and open-source (GNU General Public License) computational environment for running docking simulations. Here, we describe the combination of SAnDReS and AutoDock4 for protein-ligand docking simulations. AutoDock4 is a free program that has been applied to over a thousand receptor-ligand docking simulations. The dataset described in this chapter is available for downloading at https://github.com/azevedolab/sandres.

Implicit-descriptor ligand-based virtual screening by means of collaborative filtering

Current ligand-based machine learning methods in virtual screening rely heavily on molecular fingerprinting for preprocessing, i.e., explicit description of ligands’ structural and physicochemical properties in a vectorized form. Of particular importance to current methods are the extent to which molecular fingerprints describe a particular ligand and what metric sufficiently captures similarity among ligands. In this work, we propose and evaluate methods that do not require explicit feature vectorization through fingerprinting, but, instead, provide implicit descriptors based only on other known assays. Our methods are based upon well known collaborative filtering algorithms used in recommendation systems. Our implicit descriptor method does not require any fingerprint similarity search, which makes the method free of the bias arising from the empirical nature of the fingerprint models. We show that implicit methods significantly outperform traditional machine learning methods, and the main strengths of implicit methods are their resilience to target-ligand sparsity and high potential for spotting promiscuous ligands.

Protein structure-based drug design: from docking to molecular dynamics

Recent years have witnessed rapid developments of computer-aided drug design methods, which have reached accuracy that allows their routine practical applications in drug discovery campaigns. Protein structure-based methods are useful for the prediction of binding modes of small molecules and their relative affinity. The high-throughput docking of up to 10⁶ small molecules followed by scoring based on implicit-solvent force field can robustly identify micromolar binders using a rigid protein target. Molecular dynamics with explicit solvent is a low-throughput technique for the characterization of flexible binding sites and accurate evaluation of binding pathways, kinetics, and thermodynamics. In this review we highlight recent advancements in applications of ligand docking tools and molecular dynamics simulations to ligand identification and optimization.