LEA3D: a computer-aided ligand design for structure-based drug design

Computational Exploration of Potential Pharmacological Inhibitors Targeting the Envelope Protein of the Kyasanur Forest Disease Virus

The limitations of the current vaccination strategy for the Kyasanur Forest Disease virus (KFDV) underscore the critical need for effective antiviral treatments, highlighting the crucial importance of exploring novel therapeutic approaches through in silico drug design. Kyasanur Forest Disease, caused by KFDV, is a tick-borne disease with a mortality of 3-5% and an annual incidence of 400 to 500 cases. In the early stage of infection, the envelope protein plays a crucial role by facilitating host-virus interactions. The objective of this research is to develop effective antivirals targeting the envelope protein to disrupt the virus-host interaction. In line with this, the 3D structure of the envelope protein was modeled and refined through molecular modeling techniques, and subsequently, ligands were designed via de novo design and pharmacophore screening, yielding 12 potential hits followed by ADMET analysis. The top five candidates underwent geometry optimization and molecular docking. Notably, compounds L4 (SA28) and L3 (CNP0247967) are predicted to have significant binding affinities of -8.91 and -7.58 kcal/mol, respectively, toward the envelope protein, based on computational models. Both compounds demonstrated stability during 200 ns molecular dynamics simulations, and the MM-GBSA binding free-energy values were -85.26 ± 4.63 kcal/mol and -66.60 ± 2.92 kcal/mol for the envelope protein L3 and L4 complexes, respectively. Based on the computational prediction, it is suggested that both compounds have potential as drug candidates for controlling host-virus interactions by targeting the envelope protein. Further validation through in-vitro assays would complement the findings of the present in silico investigations.

A hybrid energy-based and AI-based screening approach for the discovery of novel inhibitors of JAK3

The JAKs protein family is composed of four isoforms, and JAK3 has been regarded as a druggable target for the development of drugs to treat various diseases, including hematologic tumors, cancer, and neuronal death. Therefore, the discovery of JAK3 inhibitors with novel scaffolds possesses the potential to provide additional options for drug development. This article presents a structure-based hybrid high-throughput virtual screening (HTVS) protocol as well as the DeepDock algorithm, which is based on geometric deep learning. These techniques were used to identify inhibitors of JAK3 with a novel sketch from a specific "In-house" database. Using molecular docking with varying precision, MM/GBSA, geometric deep learning scoring, and manual selection, 10 compounds were obtained for subsequent biological evaluation. One of these 10 compounds, compound 8, was found to have inhibitory potency against JAK3 and the MOLM-16 cell line, providing a valuable lead compound for further development of JAK3 inhibitors. To gain a better understanding of the interaction between compound 8 and JAK3, molecular dynamics (MD) simulations were conducted to provide more details on the binding conformation of compound 8 with JAK3 to guide the subsequent structure optimization. In this article, we achieved compound 8 with a novel sketch possessing inhibitory bioactivity against JAK3, and it would provide an acceptable "hit" for further structure optimization and modification to develop JAK3 inhibitors.

Drug repurposing and sequence analysis in S-glycoprotein variants reveals critical signature patterns and destabilization of receptor-binding domain in omicron variant

The evolution of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus since its emergence in 2019 has yielded several new viral variants with varied infectivity, disease severity, and antigenicity. Although most mutations are expected to be relatively neutral, mutations at the Spike region of the genome have shown to have a major impact on the viral transmission and infection in humans. Therefore, it is crucial to survey the structures of spike protein across the global virus population to contextualize the rate of therapeutic success against these variants. In this study, high-frequency mutational variants from different geographic regions were pooled in order to study the structural evolution of the spike protein through drug docking and MD simulations. We investigated the mutational burden in the spike subregions and have observed that the different variants harbour unique signature patterns in the spike subregions, with certain domains being highly prone to mutations. Further, the MD simulations and docking study revealed that different variants show differential stability when docked for the same set of drug targets. This work sheds light on the mutational burden and the stability landscape of the spike protein across the variants from different geographical regions.Communicated by Ramaswamy H. Sarma.

BMaps: A Web Application for Fragment-Based Drug Design and Compound Binding Evaluation

Fragment-based drug design uses data about where, and how strongly, small chemical fragments bind to proteins, to assemble new drug molecules. Over the past decade, we have been successfully using fragment data, derived from thermodynamically rigorous Monte Carlo fragment-protein binding simulations, in dozens of preclinical drug programs. However, this approach has not been available to the broader research community because of the cost and complexity of doing simulations and using design tools. We have developed a web application, called BMaps, to make fragment-based drug design widely available with greatly simplified user interfaces. BMaps provides access to a large repository (>550) of proteins with 100s of precomputed fragment maps, druggable hot spots, and high-quality water maps. Users can also employ their own structures or those from the Protein Data Bank and AlphaFold DB. Multigigabyte data sets are searched to find fragments in bondable orientations, ranked by a binding-free energy metric. The designers use this to select modifications that improve affinity and other properties. BMaps is unique in combining conventional tools such as docking and energy minimization with fragment-based design, in a very easy to use and automated web application. The service is available at https://www.boltzmannmaps.com.

Artificial intelligence based virtual screening study for competitive and allosteric inhibitors of the SARS-CoV-2 main protease

SARS-CoV-2 is a highly contagious and dangerous coronavirus that first appeared in late 2019 causing COVID-19, a pandemic of acute respiratory illnesses that is still a threat to health and the general public safety. We performed deep docking studies of 800 M unique compounds in both the active and allosteric sites of the SARS-COV-2 Main Protease (Mpro) and the 15 M and 13 M virtual hits obtained were further taken for conventional docking and molecular dynamic (MD) studies. The best XP Glide docking scores obtained were -14.242 and -12.059 kcal/mol by CHEMBL591669 and the highest binding affinities were -10.5 kcal/mol (from 444215) and -11.2 kcal/mol (from NPC95421) for active and allosteric sites, respectively. Some hits can bind both sites making them a great area of concern. Re-docking of 8 random allosteric complexes in the active site shows a significant reduction in docking scores with a t-test P value of 2.532 × 10-11 at 95% confidence. Some specific interactions have higher elevations in docking scores. MD studies on 15 complexes show that single-ligand systems are stable as compared to double-ligand systems, and the allosteric binders identified are shown to modulate the active site binding as evidenced by the changes in the interaction patterns and stability of ligands and active site residues. When an allosteric complex was docked to the second monomer to check for homodimer formation, the validated homodimer could not be re-established, further supporting the potential of the identified allosteric binders. These findings could be important in developing novel therapeutics against SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

A molecular evolution algorithm for ligand design in DOCK

As a complement to virtual screening, de novo design of small molecules is an alternative approach for identifying potential drug candidates. Here, we present a new 3D genetic algorithm to evolve molecules through breeding, mutation, fitness pressure, and selection. The method, termed DOCK_GA, builds upon and leverages powerful sampling, scoring, and searching routines previously implemented into DOCK6. Three primary experiments were used during development: Single-molecule evolution evaluated three selection methods (elitism, tournament, and roulette), in four clinically relevant systems, in terms of mutation type and crossover success, chemical properties, ensemble diversity, and fitness convergence, among others. Large scale benchmarking assessed performance across 651 different protein-ligand systems. Ensemble-based evolution demonstrated using multiple inhibitors simultaneously to seed growth in a SARS-CoV-2 target. Key takeaways include: (1) The algorithm is robust as demonstrated by the successful evolution of molecules across a large diverse dataset. (2) Users have flexibility with regards to parent input, selection method, fitness function, and molecular descriptors. (3) The program is straightforward to run and only requires a single executable and input file at run-time. (4) The elitism selection method yields more tightly clustered molecules in terms of 2D/3D similarity, with more favorable fitness, followed by tournament and roulette.

Design and Study of In Silico Binding Dynamics of Certain Isoxazole Bearing Leads Against Aβ-42 and BACE-1 Loop in Protein Fibrillation

Aims:

Design isoxazole bearing leads as dual inhibitors against Amyloid β and BACE-1 loop in protein fibrillation.

Background:

Protein fibrillation is one of the key reasons for several diseases, namely Alzheimer’s, Parkinson’s, and many others. One of the key strategies of preventing protein fibrillation is destabilizing the protein fibrils themselves or inhibiting the amyloid fibril-forming pathway in the initial stage.

Introduction:

Attempts have been taken to design newer leads to inhibit protein fibrillation by targeting the β-amyloidogenesis pathway in the brain. To exploit interfenestration between Amyloid β -42 protein and BACE-1 (β-site amyloid precursor protein cleaving enzyme) for amyloidogenesis, studies are undertaken to design dual inhibitors against the same.

Method:

In vitro binding interactions were found using docking, de novo ligand design, and MD simulation study.

Results:

Three compounds bearing an isoxazole heterocyclic nucleus were designed which could successfully bind to the hydrophobic raft and salt bridge residues Asp 23-Lys-26 of Amyloid β, destabilizing the growing fibril. Additionally, one of our candidate compounds exhibited force of interaction with Thr232 at the S3 pocket of BACE-1, interacted with key residue Asp228, Tyr71, and Thr72 of the β-hairpin flap and hydrogen bonding with Gly11 at loop 10s.

Conclusion:

Protein flexibility dynamics of the Aβ-42 protein revealed that there is a considerable conformational change of the same with or without ligand binding. The lower RMSF of the bound region and reprogramming residual contacts within the Aβ-42 protein suggested successful binding of the ligand with the protein, lowering the access for further β-β dimerization.

Targeting the BspC-vimentin interaction to develop anti-virulence therapies during Group B streptococcal meningitis

Bacterial infections are a major cause of morbidity and mortality worldwide and the rise of antibiotic resistance necessitates development of alternative treatments. Pathogen adhesins that bind to host cells initiate disease pathogenesis and represent potential therapeutic targets. We have shown previously that the BspC adhesin in Group B Streptococcus (GBS), the leading cause of bacterial neonatal meningitis, interacts with host vimentin to promote attachment to brain endothelium and disease development. Here we determined that the BspC variable (V-) domain contains the vimentin binding site and promotes GBS adherence to brain endothelium. Site directed mutagenesis identified a binding pocket necessary for GBS host cell interaction and development of meningitis. Using a virtual structure-based drug screen we identified compounds that targeted the V-domain binding pocket, which blocked GBS adherence and entry into the brain in vivo. These data indicate the utility of targeting the pathogen-host interface to develop anti-virulence therapeutics.

LEADD: Lamarckian evolutionary algorithm for de novo drug design

Given an objective function that predicts key properties of a molecule, goal-directed de novo molecular design is a useful tool to identify molecules that maximize or minimize said objective function. Nonetheless, a common drawback of these methods is that they tend to design synthetically unfeasible molecules. In this paper we describe a Lamarckian evolutionary algorithm for de novo drug design (LEADD). LEADD attempts to strike a balance between optimization power, synthetic accessibility of designed molecules and computational efficiency. To increase the likelihood of designing synthetically accessible molecules, LEADD represents molecules as graphs of molecular fragments, and limits the bonds that can be formed between them through knowledge-based pairwise atom type compatibility rules. A reference library of drug-like molecules is used to extract fragments, fragment preferences and compatibility rules. A novel set of genetic operators that enforce these rules in a computationally efficient manner is presented. To sample chemical space more efficiently we also explore a Lamarckian evolutionary mechanism that adapts the reproductive behavior of molecules. LEADD has been compared to both standard virtual screening and a comparable evolutionary algorithm using a standardized benchmark suite and was shown to be able to identify fitter molecules more efficiently. Moreover, the designed molecules are predicted to be easier to synthesize than those designed by other evolutionary algorithms.

Fragment-to-lead tailored in silico design

Fragment-based drug discovery (FBDD) emerged as a disruptive technology and became established during the last two decades. Its rationality and low entry costs make it appealing, and the numerous examples of approved drugs discovered through FBDD validate the approach. However, FBDD still faces numerous challenges. Perhaps the most important one is the transformation of the initial fragment hits into viable leads. Fragment-to-lead (F2L) optimization is resource-intensive and is therefore limited in the possibilities that can be actively pursued. In silico strategies play an important role in F2L, as they can perform a deeper exploration of chemical space, prioritize molecules with high probabilities of being active and generate non-obvious ideas. Here we provide a critical overview of current in silico strategies in F2L optimization and highlight their remarkable impact. While very effective, most solutions are target- or fragment- specific. We propose that fully integrated in silico strategies, capable of automatically and systematically exploring the fast-growing available chemical space can have a significant impact on accelerating the release of fragment originated drugs.

CHARMM-GUI Ligand Designer for Template-Based Virtual Ligand Design in a Binding Site

Rational drug design involves a task of finding ligands that would bind to a specific target protein. This work presents CHARMM-GUI Ligand Designer that is an intuitive and interactive web-based tool to design virtual ligands that match the shape and chemical features of a given protein binding site. Ligand Designer provides ligand modification capabilities with 3D visualization that allow researchers to modify and redesign virtual ligands while viewing how the protein-ligand interactions are affected. Virtual ligands can also be parameterized for further molecular dynamics (MD) simulations and free energy calculations. Using 8 targets from 8 different protein classes in the directory of useful decoys, enhanced (DUD-E) data set, we show that Ligand Designer can produce similar ligands to the known active ligands in the crystal structures. Ligand Designer also produces stable protein-ligand complex structures when tested using short MD simulations. We expect that Ligand Designer can be a useful and user-friendly tool to design small molecules in any given potential ligand binding site on a protein of interest.

Advances in the Exploration of the Epigenetic Relevant Chemical Space

Epigenetic drug discovery is a promising avenue to find therapeutic agents for treating several diseases and developing novel chemical probes for research. In order to identify hit and lead compounds, the chemical space has been explored and screened, generating valuable bioactivity information that can be used for multiple purposes such as prediction of the activity of existing chemicals, e.g., small molecules, guiding the design or optimization of compounds, and expanding the epigenetic relevant chemical space. Herein, we review the chemical spaces explored for epigenetic drug discovery and discuss the advances in using structure-activity relationships stored in public chemogenomic databases. We also review current efforts to chart and identify novel regions of the epigenetic relevant chemical space. In particular, we discuss the development and accessibility of two significant types of compound libraries focused on epigenetic targets: commercially available libraries for screening and targeted chemical libraries using de novo design. In this mini-review, we emphasize inhibitors of DNA methyltransferases.

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.

A comprehensive SARS-CoV-2 genomic analysis identifies potential targets for drug repurposing

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) which is a novel human coronavirus strain (HCoV) was initially reported in December 2019 in Wuhan City, China. This acute infection caused pneumonia-like symptoms and other respiratory tract illness. Its higher transmission and infection rate has successfully enabled it to have a global spread over a matter of small time. One of the major concerns involving the SARS-COV-2 is the mutation rate, which enhances the virus evolution and genome variability, thereby making the design of therapeutics difficult. In this study, we identified the most common haplotypes from the haplotype network. The conserved genes and population level variants were analysed. Non-Structural Protein 10 (NSP10), Nucleoprotein, Papain-like protease (Plpro or NSP3) and 3-Chymotrypsin like protease (3CLpro or NSP5), which were conserved at the highest threshold, were used as drug targets for molecular dynamics simulations. Darifenacin, Nebivolol, Bictegravir, Alvimopan and Irbesartan are among the potential drugs, which are suggested for further pre-clinical and clinical trials. This particular study provides a comprehensive targeting of the conserved genes. We also identified the mutation frequencies across the viral genome.

Advances in de Novo Drug Design: From Conventional to Machine Learning Methods

. De novo drug design is a computational approach that generates novel molecular structures from atomic building blocks with no a priori relationships. Conventional methods include structure-based and ligand-based design, which depend on the properties of the active site of a biological target or its known active binders, respectively. Artificial intelligence, including machine learning, is an emerging field that has positively impacted the drug discovery process. Deep reinforcement learning is a subdivision of machine learning that combines artificial neural networks with reinforcement-learning architectures. This method has successfully been employed to develop novel de novo drug design approaches using a variety of artificial networks including recurrent neural networks, convolutional neural networks, generative adversarial networks, and autoencoders. This review article summarizes advances in de novo drug design, from conventional growth algorithms to advanced machine-learning methodologies and highlights hot topics for further development.

Computational Approaches for De Novo Drug Design: Past, Present, and Future

Drug discovery is time- and resource-consuming. To this end, computational approaches that are applied in de novo drug design play an important role to improve the efficiency and decrease costs to develop novel drugs. Over several decades, a variety of methods have been proposed and applied in practice. Traditionally, drug design problems are always taken as combinational optimization in discrete chemical space. Hence optimization methods were exploited to search for new drug molecules to meet multiple objectives. With the accumulation of data and the development of machine learning methods, computational drug design methods have gradually shifted to a new paradigm. There has been particular interest in the potential application of deep learning methods to drug design. In this chapter, we will give a brief description of these two different de novo methods, compare their application scopes and discuss their possible development in the future.

Carfilzomib as a potential inhibitor of NADH-dependent enoyl-acyl carrier protein reductases of Klebsiella pneumoniae and Mycobacterium tuberculosis as a drug target enzyme: insights from molecular docking and molecular dynamics

Multiple antibiotic-resistant strains of Klebsiella pneumoniae can cause life-threatening infections. Bacterial enoyl-acyl carrier protein (ACP) reductases (ENRs) are considered critical targets for developing antibiotics. Our current study aims to identify inhibitors of K. pneumoniae ENRs (FabI and FabV). Due to the unavailability of experimental structures, protein models of FabI and FabV were predicted and validated in this study. Virtual screening of the 1930 FDA-approved drug database was conducted against the active site of the FabI protein with the help of the LEA3D server, and carfilzomib was chosen among the screened drugs for further docking studies. Carfilzomib, a proteasome inhibitor used in the treatment of multiple myeloma, was among the best-suited compounds obtained from the virtual screening and was found to be bactericidal in the in vitro experiment. Carfilzomib was docked against the active sites of the FabI and FabV proteins, and the ENR of Mycobacterium tuberculosis, InhA. Carfilzomib showed a high binding affinity with all three proteins. Molecular dynamics (MD) simulations were conducted following the docking studies. MD simulations revealed that carfilzomib binds strongly to the active sites of the above mentioned ENRs. Our study found that carfilzomib is a potential inhibitor of the ENRs of K. pneumoniae and M. tuberculosis. This is a possible mechanism of its bactericidal property against M. tuberculosis observed in vitro in addition to its predicted actions on zinc-dependent metalloprotease-1 and peptide deformylase, two other drug target enzymes of M. tuberculosis. Our study suggests that this drug could be used as a lead compound to develop antibiotics that can selectively act against ENRs of bacteria, without interfering with the activities of human proteasome. Communicated by Ramaswamy H. Sarma.

In silico analysis and identification of promising hits against 2019 novel coronavirus 3C-like main protease enzyme

The recent outbreak of the 2019 novel coronavirus disease (COVID-19) has been proved as a global threat. No particular drug or vaccine has not yet been discovered which may act specifically against severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) and causes COVID-19. For this highly infectious virus, 3CL-like main protease (3CL^pro) plays a key role in the virus life cycle and can be considered as a pivotal drug target. Structure-based virtual screening of DrugBank database resulted in 20 hits against 3CL^pro. Atomistic 100 ns molecular dynamics of five top hits and binding energy calculation analyses were performed for main protease-hit complexes. Among the top five hits, Nafarelin and Icatibant affirmed the binding energy (g_MMPBSA) of -712.94 kJ/mol and -851.74 kJ/mol, respectively. Based on binding energy and stability of protein-ligand complex; the present work reports these two drug-like hits against SARS-CoV-2 main protease.Communicated by Ramaswamy H. Sarma.

Lead Molecules as Novel Aromatase Inhibitors: In Silico De Novo Designing and Binding Affinity Studies

Background:

Aromatase inhibitors emerged as a pivotal moiety to selectively block estrogen production, prevention and treatment of tumour growth in breast cancer. De novo drug design is an alternative approach to blind virtual screening for successful designing of the novel molecule against various therapeutic targets.

Objective:

In the present study, we have explored the de novo approach to design novel aromatase inhibitors.

Method:

The e-LEA3D, a computational-aided drug design web server was used to design novel drug-like candidates against the target aromatase. For drug-likeness ADME parameters (molecular weight, H-bond acceptors, H-bond donors, LogP and number of rotatable bonds) of designed molecules were calculated in TSAR software package, geometry optimization and energy minimization was accomplished using Chem Office. Further, molecular docking study was performed in Molegro Virtual Docker (MVD).

Results:

Among 17 generated molecules using the de novo pathway, 13 molecules passed the Lipinski filter pertaining to their bioavailability characteristics. De novo designed molecules with drug-likeness were further docked into the mapped active site of aromatase to scale up their affinity and binding fitness with the target. Among de novo fabricated drug like candidates (1-13), two molecules (5, 6) exhibited higher affinity with aromatase in terms of MolDock score (-150.650, -172.680 Kcal/mol, respectively) while molecule 8 showed lowest target affinity (-85.588 Kcal/mol).

Conclusion:

The binding patterns of lead molecules (5, 6) could be used as a pharmacophore for medicinal chemists to explore these molecules for their aromatase inhibitory potential.

Virtual screening web servers: designing chemical probes and drug candidates in the cyberspace

The interplay between life sciences and advancing technology drives a continuous cycle of chemical data growth; these data are most often stored in open or partially open databases. In parallel, many different types of algorithms are being developed to manipulate these chemical objects and associated bioactivity data. Virtual screening methods are among the most popular computational approaches in pharmaceutical research. Today, user-friendly web-based tools are available to help scientists perform virtual screening experiments. This article provides an overview of internet resources enabling and supporting chemical biology and early drug discovery with a main emphasis on web servers dedicated to virtual ligand screening and small-molecule docking. This survey first introduces some key concepts and then presents recent and easily accessible virtual screening and related target-fishing tools as well as briefly discusses case studies enabled by some of these web services. Notwithstanding further improvements, already available web-based tools not only contribute to the design of bioactive molecules and assist drug repositioning but also help to generate new ideas and explore different hypotheses in a timely fashion while contributing to teaching in the field of drug development.

PETRA: Drug Engineering via Rigidity Analysis

Rational drug design aims to develop pharmaceutical agents that impart maximal therapeutic benefits via their interaction with their intended biological targets. In the past several decades, advances in computational tools that inform wet-lab techniques have aided the development of a wide variety of new medicines with high efficacies. Nonetheless, drug development remains a time and cost intensive process. In this work, we have developed a computational pipeline for assessing how individual atoms contribute to a ligand’s effect on the structural stability of a biological target. Our approach takes as input a protein-ligand resolved PDB structure file and systematically generates all possible ligand variants. We assess how the atomic-level edits to the ligand alter the drug’s effect via a graph theoretic rigidity analysis approach. We demonstrate, via four case studies of common drugs, the utility of our pipeline and corroborate our analyses with known biophysical properties of the medicines, as reported in the literature.

FragPELE: Dynamic Ligand Growing within a Binding Site. A Novel Tool for Hit-To-Lead Drug Design

The early stages of drug discovery rely on hit-to-lead programs, where initial hits undergo partial optimization to improve binding affinities for their biological target. This is an expensive and time-consuming process, requiring multiple iterations of trial and error designs, an ideal scenario for applying computer simulation. However, most state-of-the-art modeling techniques fail to provide a fast and reliable answer to the Induced-Fit protein-ligand problem. To aid in this matter, we present FragPELE, a new tool for in silico hit-to-lead drug design, capable of growing a fragment from a bound core while exploring the protein-ligand conformational space. We tested the ability of FragPELE to predict crystallographic data, even in cases where cryptic sub-pockets open because of the presence of particular R-groups. Additionally, we evaluated the potential of the software on growing and scoring five congeneric series from the 2015 FEP+ dataset, comparing them to FEP+, SP and Induced-Fit Glide, and MMGBSA simulations. Results show that FragPELE could be useful not only for finding new cavities and novel binding modes in cases where standard docking tools cannot but also to rank ligand activities in a reasonable amount of time and with acceptable precision.

New Data Structure for Computational Molecular Design with Atomic or Fragment Resolution

A new molecular data structure and molecular structure operation algorithms are proposed for general purpose molecular design. The data structure allows for a variety of molecular operations for creating new molecules. Two types of molecular operations were developed, unimolecular and bimolecular operations. In unimolecular operations, a child molecule can be created from a parent via addition of a functional group, deletion of a fragment, mutation of an atom, etc. In bimolecular operations, children molecules are generated from two parent molecules through combination or crossover (hybridization). These molecular operations are essential for the creation and modification of molecules for the purpose of molecular design. The data structure is capable of representing linear, branched, multifunctional, and multivalent compounds. Algorithms are developed for deriving the molecular data structure of a molecule from its atomic coordinates and vice versa. We show that this new molecular data structure and the developed algorithms, referred to as Molecular Assembling and Representation Suite, allow one to generate a comprehensive library of new molecules via performing every possible molecular structure modification.

Identification of promising anti-DNA gyrase antibacterial compounds using de novo design, molecular docking and molecular dynamics studies

The rapidly increasing rate of antibiotic resistance is of great concern. Approximately two million deaths result annually from bacterial infections worldwide. Therefore, there is a paramount requirement to develop innovative and novel antibacterial agents with new mechanisms of action and activity against resistant bacterial strains. For this purpose, a set of benzothiazole and N-phenylpyrrolamides derivatives reported as DNA Gyrase B (GyrB) inhibitors were collected from the literature and docked inside the receptor cavity of DNA Gyrase B (PDB ID: 5L3J). The best 10 docked complexes were used to identify novel antibacterial chemical agents through a de novo design approach. Out of initial 300 chemical analogues, the best six analogues were identified using screening with a set of criteria followed by pharmacokinetic analysis. The binding interactions of the best six analogues revealed that all molecules formed a number of critical interactions with catalytic amino residues of DNA Gyrase B with high binding energy. The predicted inhibitory constant biological activity based on binding energy supported the potential of the molecules as DNA Gyrase B ligands. The RMSD, RMSF, and radius of gyration parameters obtained from the 100 ns molecular dynamics simulation study clearly demonstrated that all six analogues were efficient enough to form stable complexes with DNA Gyrase B. High negative binding energy of all ligands obtained from MM-GBSA approach undoubtedly explained the strong affinity toward the DNA Gyrase B. Therefore, the proposed de novo designed molecules can be considered as promising antibacterial chemical agents subject to experimental validation, in vitro.Communicated by Ramaswamy H. Sarma.

MoleGear: A Java-Based Platform for Evolutionary De Novo Molecular Design

A Java-based platform, MoleGear, is developed for de novo molecular design based on the chemistry development kit (CDK) and other Java packages. MoleGear uses evolutionary algorithm (EA) to explore chemical space, and a suite of fragment-based operators of growing, crossover, and mutation for assembling novel molecules that can be scored by prediction of binding free energy or a weighted-sum multi-objective fitness function. The EA can be conducted in parallel over multiple nodes to support large-scale molecular optimizations. Some complementary utilities such as fragment library design, chemical space analysis, and graphical user interface are also integrated into MoleGear. The candidate molecules as inhibitors for the human immunodeficiency virus 1 (HIV-1) protease were designed by MoleGear, which validates the potential capability for de novo molecular design.

Identifying de-NEDDylation inhibitors: Virtual high-throughput screens targeting SENP8

Protein modification can have far-reaching effects. NEDDylation, a protein modification process with the protein NEDD8, stabilizes and modifies how the targeted protein interacts with other proteins. Its role in system regulation makes it a prime therapeutic target, and virtual high-throughput screening has already identified new NEDD8 inhibitors. SENP8 matures the NEDD8 proenzyme into the active form and regulates NEDDylation by removing NEDD8 from over-NEDDylated proteins. In this work, SENP8 inhibitor candidates were identified in two rounds of virtual high-throughput screening. Of the ten candidates identified in the first round of screening, four were active in validation experiments to yield an experimental hit rate of 40%. Of the five candidates identified in the second round of screening, one was active in validation experiments to yield an experimental hit rate of 20%. Results indicate virtual high-throughput screening improved hit rates over traditional high-throughput screening. The SENP8 inhibitor candidates can be used to interrogate the NEDDylation regulation mechanism.

Design, synthesis and algicides activities of thiourea derivatives as the novel scaffold aldolase inhibitors

By using a new Fragment-Based Virtual Screen strategy, two series of novel FBA-II inhibitors (thiourea derivatives) were de novo discovered based on the active site of fructose-1, 6-bisphosphate aldolase from Cyanobacterial (CyFBA). In comparison, most of the N-(2-benzoylhydrazine-1-carbonothioyl) benzamide derivatives (L14∼L22) exhibit higher CyFBA-II inhibitory activities compared to N-(phenylcarbamothioyl) benzamide derivatives (L1∼L13). Especially, compound L14 not only shows higher CyFBA-II activity (K_i = 0.65 μM), but also exhibits most potent in vivo activity against Synechocystis sp. PCC 6803 (EC₅₀ = 0.09 ppm), higher (7-fold) than that of our previous inhibitor (EC₅₀ = 0.6 ppm). The binding modes of compound L14 and CyFBA-II were further elucidated by jointly using DOX computational protocol, MM-PBSA and site-directed mutagenesis assays. The positive results suggest that strategy adopted in this study was promising to rapidly discovery the potent inhibitors with novel scaffolds. The satisfactory algicide activities suggest that the thiourea derivatives is very likely to be a promising lead for the development of novel specific algicides to solve Cyanobacterial harmful algal blooms (CHABs).

Virtual high-throughput screens identifying hPK-M2 inhibitors: Exploration of model extrapolation

Glycolysis with PK-M2 occurs typically in anaerobic conditions and atypically in aerobic conditions, which is known as the Warburg effect. The Warburg effect is found in many oncogenic situations and is believed to provide energy and biomass for oncogenesis to persist. The work presented targets human PK-M2 (hPK-M2) in a virtual high-throughput screen to identify new inhibitors and leads for further study. In the initial screen, one of the 12 candidates selected for experimental validation showed biological activity (hit-rate = 8.13%). In the second screen with retrained models, six of 11 candidates selected for experimental validation showed biological activity (hit-rate: 54.5%). Additionally, four different scaffolds were identified for further analysis when examining the tested candidates and compounds in the training data. Finally, extrapolation was necessary to identify a sufficient number of candidates to test in the second screen. Examination of the results suggested stepwise extrapolation to maximize efficiency.

Identifying new clotting factor XIa inhibitors in virtual high-throughput screens using PCA-GA-SVM models and signature

Blood Clotting Factor XI is an important actor in the clotting mechanism: it activates downstream zymogen involved in the clotting process. It can be targeted for activation or inhibition depending on treatment goals to enhance or inhibit clotting. In terms of antithrombosis treatment, Factor XI has emerged as a promising target to focus on. In this work, an iterative virtual high-throughput screening pipeline was proposed that can supplement current efforts to find inhibitors. The first iteration identified 11 compounds to test with 3 active for a hit-rate of 27.3%. The second iteration of the pipeline identified another 11 compounds to test with 7 active for a hit-rate of 63.6%. © 2018 American Institute of Chemical Engineers Biotechnol. Prog., 34:1553-1565, 2018.

Pharmaceutical Machine Learning: Virtual High-Throughput Screens Identifying Promising and Economical Small Molecule Inhibitors of Complement Factor C1s

When excessively activated, C1 is insufficiently regulated, which results in tissue damage. Such tissue damage causes the complement system to become further activated to remove the resulting tissue damage, and a vicious cycle of activation/tissue damage occurs. Current Food and Drug Administration approved treatments include supplemental recombinant C1 inhibitor, but these are extremely costly and a more economical solution is desired. In our work, we have utilized an existing data set of 136 compounds that have been previously tested for activity against C1. Using these compounds and the activity data, we have created models using principal component analysis, genetic algorithm, and support vector machine approaches to characterize activity. The models were then utilized to virtually screen the 72 million compound PubChem repository. This first round of virtual high-throughput screening identified many economical and promising inhibitor candidates, a subset of which was tested to validate their biological activity. These results were used to retrain the models and rescreen PubChem in a second round vHTS. Hit rates for the first round vHTS were 57%, while hit rates for the second round vHTS were 50%. Additional structure⁻property analysis was performed on the active and inactive compounds to identify interesting scaffolds for further investigation.

Ribonucleotide reductase as a drug target against drug resistance Mycobacterium leprae: A molecular docking study

Leprosy is a chronic infection of skin and nerve caused by Mycobacterium leprae. The treatment is based on standard multi drug therapy consisting of dapsone, rifampicin and clofazamine. The use of rifampicin alone or with dapsone led to the emergence of rifampicin-resistant Mycobacterium leprae strains. The emergence of drug-resistant leprosy put a hurdle in the leprosy eradication programme. The present study aimed to predict the molecular model of ribonucleotide reductase (RNR), the enzyme responsible for biosynthesis of nucleotides, to screen new drugs for treatment of drug-resistant leprosy. The study was conducted by retrieving RNR of M. leprae from GenBank. A molecular 3D model of M. leprae was predicted using homology modelling and validated. A total of 325 characters were included in the analysis. The predicted 3D model of RNR showed that the ϕ and φ angles of 251 (96.9%) residues were positioned in the most favoured regions. It was also conferred that 18 α-helices, 6 β turns, 2 γ turns and 48 helix-helix interactions contributed to the predicted 3D structure. Virtual screening of Food and Drug Administration approved drug molecules recovered 1829 drugs of which three molecules, viz., lincomycin, novobiocin and telithromycin, were taken for the docking study. It was observed that the selected drug molecules had a strong affinity towards the modelled protein RNR. This was evident from the binding energy of the drug molecules towards the modelled protein RNR (-6.10, -6.25 and -7.10). Three FDA-approved drugs, viz., lincomycin, novobiocin and telithromycin, could be taken for further clinical studies to find their efficacy against drug resistant leprosy.

De Novo Design of Ligands Using Computational Methods

De novo design technique is complementary to high-throughput virtual screening and is believed to contribute in pharmaceutical development of novel drugs with desired properties at a very low cost and time-efficient manner. In this chapter, we outline the basic de novo design concepts based on computational methods with an example.

Prediction of Medicinal Properties Using Mathematical Models and Computation, and Selection of Plant Materials

In any phytochemical drug discovery programme, one of the major issues is the appropriate selection of target plant species that may provide lead for new drug discovery and development. Conducting research without any working hypotheses may produce serendipitous discoveries, but the chances of success are much slimmer than any information-based targeted approach. Therefore, the plant selection process is extremely important for ensuring success. In recent years, there have been significant amounts of work involving applications of various mathematical modelling and computational techniques to predict medicinal properties of plants, and thus to provide information-based selection of plant materials for further studies aiming at potential drug discovery and development. This chapter presents an overview of methods and processes involved in plant selection by utilizing various mathematical modelling and computational techniques.

Customizable de novo design strategies for DOCK: Application to HIVgp41 and other therapeutic targets

De novo design can be used to explore vast areas of chemical space in computational lead discovery. As a complement to virtual screening, from-scratch construction of molecules is not limited to compounds in pre-existing vendor catalogs. Here, we present an iterative fragment growth method, integrated into the program DOCK, in which new molecules are built using rules for allowable connections based on known molecules. The method leverages DOCK's advanced scoring and pruning approaches and users can define very specific criteria in terms of properties or features to customize growth toward a particular region of chemical space. The code was validated using three increasingly difficult classes of calculations: (1) Rebuilding known X-ray ligands taken from 663 complexes using only their component parts (focused libraries), (2) construction of new ligands in 57 drug target sites using a library derived from ∼13M drug-like compounds (generic libraries), and (3) application to a challenging protein-protein interface on the viral drug target HIVgp41. The computational testing confirms that the de novo DOCK routines are robust and working as envisioned, and the compelling results highlight the potential utility for designing new molecules against a wide variety of important protein targets. © 2017 Wiley Periodicals, Inc.

Fragment-based drug discovery as alternative strategy to the drug development for neglected diseases

Neglected diseases (NDs) affect large populations and almost whole continents, representing 12% of the global health burden. In contrast, the treatment available today is limited and sometimes ineffective. Under this scenery, the Fragment-Based Drug Discovery emerged as one of the most promising alternatives to the traditional methods of drug development. This method allows achieving new lead compounds with smaller size of fragment libraries. Even with the wide Fragment-Based Drug Discovery success resulting in new effective therapeutic agents against different diseases, until this moment few studies have been applied this approach for NDs area. In this article, we discuss the basic Fragment-Based Drug Discovery process, brief successful ideas of general applications and show a landscape of its use in NDs, encouraging the implementation of this strategy as an interesting way to optimize the development of new drugs to NDs.

Break Down in Order To Build Up: Decomposing Small Molecules for Fragment-Based Drug Design with eMolFrag

Constructing high-quality libraries of molecular building blocks is essential for successful fragment-based drug discovery. In this communication, we describe eMolFrag, a new open-source software to decompose organic compounds into nonredundant fragments retaining molecular connectivity information. Given a collection of molecules, eMolFrag generates a set of unique fragments comprising larger moieties, bricks, and smaller linkers connecting bricks. These building blocks can subsequently be used to construct virtual screening libraries for targeted drug discovery. The robustness and computational performance of eMolFrag is assessed against the Directory of Useful Decoys, Enhanced database conducted in serial and parallel modes with up to 16 computing cores. Further, the application of eMolFrag in de novo drug design is illustrated using the adenosine receptor. eMolFrag is implemented in Python, and it is available as stand-alone software and a web server at www.brylinski.org/emolfrag and https://github.com/liutairan/eMolFrag.