A practical guide to large-scale docking

Benchmarking AlphaFold-enabled molecular docking predictions for antibiotic discovery

Efficient identification of drug mechanisms of action remains a challenge. Computational docking approaches have been widely used to predict drug binding targets; yet, such approaches depend on existing protein structures, and accurate structural predictions have only recently become available from AlphaFold2. Here, we combine AlphaFold2 with molecular docking simulations to predict protein‐ligand interactions between 296 proteins spanning Escherichia coli's essential proteome, and 218 active antibacterial compounds and 100 inactive compounds, respectively, pointing to widespread compound and protein promiscuity. We benchmark model performance by measuring enzymatic activity for 12 essential proteins treated with each antibacterial compound. We confirm extensive promiscuity, but find that the average area under the receiver operating characteristic curve (auROC) is 0.48, indicating weak model performance. We demonstrate that rescoring of docking poses using machine learning‐based approaches improves model performance, resulting in average auROCs as large as 0.63, and that ensembles of rescoring functions improve prediction accuracy and the ratio of true‐positive rate to false‐positive rate. This work indicates that advances in modeling protein‐ligand interactions, particularly using machine learning‐based approaches, are needed to better harness AlphaFold2 for drug discovery.

Development of an automated screen for Kv7.2 potassium channels and discovery of a new agonist chemotype

To identify pore domain ligands on Kv7.2 potassium ion channels, we compared wild-type (WT) and W236L mutant Kv7.2 channels in a series of assays with previously validated and novel agonist chemotypes. Positive controls were retigabine, flupirtine, and RL-81; i.e. Kv7.2 channel activators that significantly shift voltage-dependent activation to more negative potentials (ΔV50) at 5 µM. We identified 6 new compounds that exhibited differential enhancing activity between WT and W236L mutant channels. Whole cell patch-clamp electrophysiology studies were conducted to identify Kv7.2. Kv7.2/3, Kv7.4, and Kv7.5 selectivity. Our results validate the SyncroPatch platform and establish new structure activity relationships (SAR). Specifically, in addition to selective Kv7.2, Kv7.2/3, Kv7.4. and Kv7.5 agonists, we identified a novel chemotype, ZK-21, a 4-aminotetrahydroquinoline that is distinct from any of the previously described Kv7 channel modifiers. Using flexible receptor docking, ZK-21 was predicted to be stabilized by W236 and bind perpendicular to retigabine, burying the benzyl carbamate group into a tunnel reaching the core of the pore domain.

Protein-Based Virtual Screening Tools Applied for RNA-Ligand Docking Identify New Binders of the preQ1-Riboswitch

Targeting RNA with small molecules is an emerging field. While several ligands for different RNA targets are reported, structure-based virtual screenings (VSs) against RNAs are still rare. Here, we elucidated the general capabilities of protein-based docking programs to reproduce native binding modes of small-molecule RNA ligands and to discriminate known binders from decoys by the scoring function. The programs were found to perform similar compared to the RNA-based docking tool rDOCK, and the challenges faced during docking, namely, protomer and tautomer selection, target dynamics, and explicit solvent, do not largely differ from challenges in conventional protein-ligand docking. A prospective VS with the Bacillus subtilis preQ₁-riboswitch aptamer domain performed with FRED, HYBRID, and FlexX followed by microscale thermophoresis assays identified six active compounds out of 23 tested VS hits with potencies between 29.5 nM and 11.0 μM. The hits were selected not solely based on their docking score but for resembling key interactions of the native ligand. Therefore, this study demonstrates the general feasibility to perform structure-based VSs against RNA targets, while at the same time it highlights pitfalls and their potential solutions when executing RNA-ligand docking.

Identification of Potential Cytochrome P450 3A5 Inhibitors: An Extensive Virtual Screening through Molecular Docking, Negative Image-Based Screening, Machine Learning and Molecular Dynamics Simulation Studies

Cytochrome P450 3A5 (CYP3A5) is one of the crucial CYP family members and has already proven to be an important drug target for cardiovascular diseases. In the current study, the PubChem database was screened through molecular docking and high-affinity molecules were adopted for further assessment. A negative image-based (NIB) model was used for a similarity search by considering the complementary shape and electrostatics of the target and small molecules. Further, the molecules were segregated into active and inactive groups through six machine learning (ML) matrices. The active molecules found in each ML model were used for in silico pharmacokinetics and toxicity assessments. A total of five molecules followed the acceptable pharmacokinetics and toxicity profiles. Several potential binding interactions between the proposed molecules and CYP3A5 were observed. The dynamic behavior of the selected molecules in the CYP3A5 was explored through a molecular dynamics (MD) simulation study. Several parameters obtained from the MD simulation trajectory explained the stability of the protein–ligand complexes in dynamic states. The high binding affinity of each molecule was revealed by the binding free energy calculation through the MM-GBSA methods. Therefore, it can be concluded that the proposed molecules might be potential CYP3A5 molecules for therapeutic application in cardiovascular diseases subjected to in vitro/in vivo validations.

Ligand Binding Prediction Using Protein Structure Graphs and Residual Graph Attention Networks

Computational prediction of ligand–target interactions is a crucial part of modern drug discovery as it helps to bypass high costs and labor demands of in vitro and in vivo screening. As the wealth of bioactivity data accumulates, it provides opportunities for the development of deep learning (DL) models with increasing predictive powers. Conventionally, such models were either limited to the use of very simplified representations of proteins or ineffective voxelization of their 3D structures. Herein, we present the development of the PSG-BAR (Protein Structure Graph-Binding Affinity Regression) approach that utilizes 3D structural information of the proteins along with 2D graph representations of ligands. The method also introduces attention scores to selectively weight protein regions that are most important for ligand binding. Results: The developed approach demonstrates the state-of-the-art performance on several binding affinity benchmarking datasets. The attention-based pooling of protein graphs enables identification of surface residues as critical residues for protein–ligand binding. Finally, we validate our model predictions against an experimental assay on a viral main protease (Mpro)—the hallmark target of SARS-CoV-2 coronavirus.

DOCKSTRING: Easy Molecular Docking Yields Better Benchmarks for Ligand Design

The field of machine learning for drug discovery is witnessing an explosion of novel methods. These methods are often benchmarked on simple physicochemical properties such as solubility or general druglikeness, which can be readily computed. However, these properties are poor representatives of objective functions in drug design, mainly because they do not depend on the candidate compound's interaction with the target. By contrast, molecular docking is a widely applied method in drug discovery to estimate binding affinities. However, docking studies require a significant amount of domain knowledge to set up correctly, which hampers adoption. Here, we present dockstring, a bundle for meaningful and robust comparison of ML models using docking scores. dockstring consists of three components: (1) an open-source Python package for straightforward computation of docking scores, (2) an extensive dataset of docking scores and poses of more than 260,000 molecules for 58 medically relevant targets, and (3) a set of pharmaceutically relevant benchmark tasks such as virtual screening or de novo design of selective kinase inhibitors. The Python package implements a robust ligand and target preparation protocol that allows nonexperts to obtain meaningful docking scores. Our dataset is the first to include docking poses, as well as the first of its size that is a full matrix, thus facilitating experiments in multiobjective optimization and transfer learning. Overall, our results indicate that docking scores are a more realistic evaluation objective than simple physicochemical properties, yielding benchmark tasks that are more challenging and more closely related to real problems in drug discovery.

New insights into geraniol's antihemolytic, anti-inflammatory, antioxidant, and anticoagulant potentials using a combined biological and in silico screening strategy

The study aims to assess the antihemolytic and antioxidant activities of geraniol versus 2, 2'-azobis, 2-amidinopropane dihydro-chloride- (AAPH-) induced oxidative damage and hemolysis to erythrocytes and its anti-inflammatory potential against lipopolysaccharide- (LPS-) induced inflammation in white blood cells (WBCs) with a focus on its integrated computational strategies against different targeted receptors participating in inflammation and coagulation. The rats' erythrocyte suspension was incubated with different geraniol concentrations. Molecular docking and simulation were used to explore the possible interaction patterns of geraniol against the potential targeted proteins for therapeutic screening. The results displayed that geraniol had a prolonged noteworthy effect on activated partial thromboplastin time and thromboplastin time. Geraniol displayed strong antioxidant effects via reduced malondialdehyde (MDA) formation and increased GSH level and SOD activity. We observed dose-dependent prevention of K⁺ ion leakage along with a remarkable decline of hemolysis in erythrocytes pretreated with geraniol. Geraniol 100 µg/mL and diclofenac 100 µM were nontoxic to WBCs. Geraniol significantly reduces the expression and release of cellular pro-inflammatory factors TNF-α, IL-1β, IL-8, and nitric oxide, accompanied by a significant upregulation of gene expression of anti-inflammatory cytokine IL-10 in LPS-induced WBCs compared to nontreated cells. It demonstrates a much stronger inhibition potential than diclofenac in terms of inflammation inhibition. When comparing molecular docking and simulation data, current work showed that geraniol has a good affinity toward apoptosis signal-regulating kinase 1 (ASK1) and human P2Y12 receptors and could be developed as an antioxidant, anti-inflammatory, and anticoagulant medication in the future. Consequently, geraniol is recommended to have a defensive influence against oxidative stress, and hemolysis also could be developed as a promising anti-inflammatory, antioxidant, and anticoagulant medication.

Emerging frontiers in virtual drug discovery: From quantum mechanical methods to deep learning approaches

Virtual screening-based approaches to discover initial hit and lead compounds have the potential to reduce both the cost and time of early drug discovery stages, as well as to find inhibitors for even challenging target sites such as protein-protein interfaces. Here in this review, we provide an overview of the progress that has been made in virtual screening methodology and technology on multiple fronts in recent years. The advent of ultra-large virtual screens, in which hundreds of millions to billions of compounds are screened, has proven to be a powerful approach to discover highly potent hit compounds. However, these developments are just the tip of the iceberg, with new technologies and methods emerging to propel the field forward. Examples include novel machine-learning approaches, which can reduce the computational costs of virtual screening dramatically, while progress in quantum-mechanical approaches can increase the accuracy of predictions of various small molecule properties.

Dithiocarbazate Ligand-Based Cu(II), Ni(II), and Zn(II) Complexes: Synthesis, Structural Investigations, Cytotoxicity, DNA Binding, and Molecular Docking Studies

S-4-methylbenzyl-β-N-(2-methoxybenzylmethylene)dithiocarbazate ligand, 1, prepared from S-(4-methylbenzyl)dithiocarbazate, was used to produce a novel series of transition metal complexes of the type, [M (L)₂] [M = Cu(II) (2), Ni(II) (3), and Zn(II) (4), L = 1]. The ligand and its complexes were investigated by elemental analysis, FTIR, ¹H and ¹³C-NMR, MS spectrometry, and molar conductivity. In addition, single X-ray crystallography was also performed for ligand, 1, and complex 3. The Hirshfeld surface analyses were also performed to know about various bonding interactions in the ligand, 1, and complex 3. The investigated compounds were also tested to evaluate their cytotoxic behaviour. However, complex 2 showed promising results against MCF-7 and MDA-MB-213 cancer cell lines. Furthermore, the interaction of CT-DNA with ligand, 1, and complex 2 was also studied using the electronic absorption method, revealing that the compounds have potential DNA-binding ability via hydrogen bonding and hydrophobic and van der Waals interactions. A molecular docking study of complex 2 was also carried out, which revealed that free binding free energy value was −7.39 kcal mol⁻¹.

Iterative computational design and crystallographic screening identifies potent inhibitors targeting the Nsp3 Macrodomain of SARS-CoV-2

The nonstructural protein 3 (NSP3) of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) contains a conserved macrodomain enzyme (Mac1) that is critical for pathogenesis and lethality. While small molecule inhibitors of Mac1 have great therapeutic potential, at the outset of the COVID-19 pandemic there were no well-validated inhibitors for this protein nor, indeed, the macrodomain enzyme family, making this target a pharmacological orphan. Here, we report the structure-based discovery and development of several different chemical scaffolds exhibiting low- to sub-micromolar affinity for Mac1 through iterations of computer-aided design, structural characterization by ultra-high resolution protein crystallography, and binding evaluation. Potent scaffolds were designed with in silico fragment linkage and by ultra-large library docking of over 450 million molecules. Both techniques leverage the computational exploration of tangible chemical space and are applicable to other pharmacological orphans. Overall, 160 ligands in 119 different scaffolds were discovered, and 152 Mac1-ligand complex crystal structures were determined, typically to 1 Å resolution or better. Our analyses discovered selective and cell-permeable molecules, unexpected ligand-mediated protein dynamics within the active site, and key inhibitor motifs that will template future drug development against Mac1.

3D-QSAR-Based Pharmacophore Modeling, Virtual Screening, and Molecular Dynamics Simulations for the Identification of Spleen Tyrosine Kinase Inhibitors

Spleen tyrosine kinase (SYK) is an essential mediator of immune cell signaling and has been anticipated as a therapeutic target for autoimmune diseases, notably rheumatoid arthritis, allergic rhinitis, asthma, and cancers. Significant attempts have been undertaken in recent years to develop SYK inhibitors; however, limited success has been achieved due to poor pharmacokinetics and adverse effects of inhibitors. The primary goal of this research was to identify potential inhibitors having high affinity, selectivity based on key molecular interactions, and good drug-like properties than the available inhibitor, fostamatinib. In this study, a 3D-QSAR model was built for SYK based on known inhibitor IC₅₀ values. The best pharmacophore model was then used as a 3D query to screen a drug-like database to retrieve hits with novel chemical scaffolds. The obtained compounds were subjected to binding affinity prediction using the molecular docking approach, and the results were subsequently validated using molecular dynamics (MD) simulations. The simulated compounds were ranked according to binding free energy (ΔG), and the binding affinity was compared with fostamatinib. The binding mode analysis of selected compounds revealed that the hit compounds form hydrogen bond interactions with hinge region residue Ala451, glycine-rich loop residue Lys375, Ser379, and DFG motif Asp512. Identified hits were also observed to form a desirable interaction with Pro455 and Asn457, the rare feature observed in SYK inhibitors. Therefore, we argue that identified hit compounds ZINC98363745, ZINC98365358, ZINC98364133, and ZINC08789982 may help in drug design against SYK.

Application of Mathematical Modeling and Computational Tools in the Modern Drug Design and Development Process

The conventional drug discovery approach is an expensive and time-consuming process, but its limitations have been overcome with the help of mathematical modeling and computational drug design approaches. Previously, finding a small molecular candidate as a drug against a disease was very costly and required a long time to screen a compound against a specific target. The development of novel targets and small molecular candidates against different diseases including emerging and reemerging diseases remains a major concern and necessitates the development of novel therapeutic targets as well as drug candidates as early as possible. In this regard, computational and mathematical modeling approaches for drug development are advantageous due to their fastest predictive ability and cost-effectiveness features. Computer-aided drug design (CADD) techniques utilize different computer programs as well as mathematics formulas to comprehend the interaction of a target and drugs. Traditional methods to determine small-molecule candidates as a drug have several limitations, but CADD utilizes novel methods that require little time and accurately predict a compound against a specific disease with minimal cost. Therefore, this review aims to provide a brief insight into the mathematical modeling and computational approaches for identifying a novel target and small molecular candidates for curing a specific disease. The comprehensive review mainly focuses on biological target prediction, structure-based and ligand-based drug design methods, molecular docking, virtual screening, pharmacophore modeling, quantitative structure-activity relationship (QSAR) models, molecular dynamics simulation, and MM-GBSA/MM-PBSA approaches along with valuable database resources and tools for identifying novel targets and therapeutics against a disease. This review will help researchers in a way that may open the road for the development of effective drugs and preventative measures against a disease in the future as early as possible.

Maximizing the integration of virtual and experimental screening in hit discovery

INTRODUCTION

Experimental and virtual screening contributes to the discovery of more than 50% of clinical candidates. Considering the similar concept and goals, early-phase drug discovery would benefit from the effective integration of these approaches.

AREAS COVERED

After reviewing the recent trends in both experimental and virtual screening, the authors discuss different integration strategies from parallel, focused, sequential, and iterative screening. Strategic considerations are demonstrated in a number of real-life case studies.

EXPERT OPINION

Experimental and virtual screening are complementary approaches that should be integrated in lead discovery settings. Virtual screening can access extremely large synthetically feasible chemical space that can be effectively searched on GPU clusters or cloud architectures. Experimental screening provides reliable datasets by quantitative HTS applications, and DNA-encoded libraries (DEL) have enlarged the chemical space covered by these technologies. These developments, together with the use of artificial intelligence methods, represent new options for their efficient integration. The case studies discussed here demonstrate the benefits of complementary strategies, such as focused and iterative screening.

Insights into the Pharmacological Effects of Flavonoids: The Systematic Review of Computer Modeling

Computer modeling is a method that is widely used in scientific investigations to predict the biological activity, toxicity, pharmacokinetics, and synthesis strategy of compounds based on the structure of the molecule. This work is a systematic review of articles performed in accordance with the recommendations of PRISMA and contains information on computer modeling of the interaction of classical flavonoids with different biological targets. The review of used computational approaches is presented. Furthermore, the affinities of flavonoids to different targets that are associated with the infection, cardiovascular, and oncological diseases are discussed. Additionally, the methodology of bias risks in molecular docking research based on principles of evidentiary medicine was suggested and discussed. Based on this data, the most active groups of flavonoids and lead compounds for different targets were determined. It was concluded that flavonoids are a promising object for drug development and further research of pharmacology by in vitro, ex vivo, and in vivo models is required.

Design, Semisynthesis, and Estrogenic Activity of Lignan Derivatives from Natural Dibenzylbutyrolactones

Based on molecular docking studies on the ERα, a series of lignan derivatives (3–16) were designed and semisynthesized from the natural dibenzylbutyrolactones bursehernin (1) and matairesinol dimethyl ether (2). To examine their estrogenic and antiestrogenic potencies, the effects of these compounds on estrogen receptor element (ERE)-driven reporter gene expression and viability in human ER+ breast cancer cells were evaluated. Lignan compounds induced ERE-driven reporter gene expression with very low potency as compared with the pure agonist E2. However, coincubation of 5 μM of lignan derivatives 1, 3, 4, 7, 8, 9, 11, 13, and 14 with increasing concentrations of E2 (from 0.01 pM to 1 nM) reduced both the potency and efficacy of pure agonists. The binding to the rhERα-LBD was validated by TR-FRET competitive binding assay and lignans bound to the rhERα with IC₅₀ values from 0.16 μM (compound 14) to 6 μM (compound 4). Induced fit docking (IFD) and molecular dynamics (MD) simulations for compound 14 were carried out to further investigate the binding mode interactions. Finally, the in silico ADME predictions indicated that the most potent lignan derivatives exhibited good drug-likeness.

Efficient Screening of Target-Specific Selected Compounds in Mixtures by 19F NMR Binding Assay with Predicted 19F NMR Chemical Shifts

Ligand-based 19 F NMR screening is a highly effective and well-established hit-finding approach. The high sensitivity to protein binding makes it particularly suitable for fragment screening. Different criteria can be considered for generating fluorinated fragment libraries. One common strategy is to assemble a large, diverse, well-designed and characterized fragment library which is screened in mixtures, generated based on experimental 19 F NMR chemical shifts. Here, we introduce a complementary knowledge-based 19 F NMR screening approach, named 19 Focused screening, enabling the efficient screening of putative active molecules selected by computational hit finding methodologies, in mixtures assembled and on-the-fly deconvoluted based on predicted 19 F NMR chemical shifts. In this study, we developed a novel approach, named LEFshift , for 19 F NMR chemical shift prediction using rooted topological fluorine torsion fingerprints in combination with a random forest machine learning method. A demonstration of this approach to a real test case is reported.

Drugsniffer: An Open Source Workflow for Virtually Screening Billions of Molecules for Binding Affinity to Protein Targets

The SARS-CoV2 pandemic has highlighted the importance of efficient and effective methods for identification of therapeutic drugs, and in particular has laid bare the need for methods that allow exploration of the full diversity of synthesizable small molecules. While classical high-throughput screening methods may consider up to millions of molecules, virtual screening methods hold the promise of enabling appraisal of billions of candidate molecules, thus expanding the search space while concurrently reducing costs and speeding discovery. Here, we describe a new screening pipeline, called drugsniffer, that is capable of rapidly exploring drug candidates from a library of billions of molecules, and is designed to support distributed computation on cluster and cloud resources. As an example of performance, our pipeline required ∼40,000 total compute hours to screen for potential drugs targeting three SARS-CoV2 proteins among a library of ∼3.7 billion candidate molecules.

Systemic evolutionary chemical space exploration for drug discovery

Chemical space exploration is a major task of the hit-finding process during the pursuit of novel chemical entities. Compared with other screening technologies, computational de novo design has become a popular approach to overcome the limitation of current chemical libraries. Here, we reported a de novo design platform named systemic evolutionary chemical space explorer (SECSE). The platform was conceptually inspired by fragment-based drug design, that miniaturized a “lego-building” process within the pocket of a certain target. The key to virtual hits generation was then turned into a computational search problem. To enhance search and optimization, human intelligence and deep learning were integrated. Application of SECSE against phosphoglycerate dehydrogenase (PHGDH), proved its potential in finding novel and diverse small molecules that are attractive starting points for further validation. This platform is open-sourced and the code is available at http://github.com/KeenThera/SECSE.

Considerations and challenges for sex-aware drug repurposing

Sex differences are essential factors in disease etiology and manifestation in many diseases such as cardiovascular disease, cancer, and neurodegeneration [33]. The biological influence of sex differences (including genomic, epigenetic, hormonal, immunological, and metabolic differences between males and females) and the lack of biomedical studies considering sex differences in their study design has led to several policies. For example, the National Institute of Health's (NIH) sex as a biological variable (SABV) and Sex and Gender Equity in Research (SAGER) policies to motivate researchers to consider sex differences [204]. However, drug repurposing, a promising alternative to traditional drug discovery by identifying novel uses for FDA-approved drugs, lacks sex-aware methods that can improve the identification of drugs that have sex-specific responses [7, 11, 14, 33]. Sex-aware drug repurposing methods either select drug candidates that are more efficacious in one sex or deprioritize drug candidates based on if they are predicted to cause a sex-bias adverse event (SBAE), unintended therapeutic effects that are more likely to occur in one sex. Computational drug repurposing methods are encouraging approaches to develop for sex-aware drug repurposing because they can prioritize sex-specific drug candidates or SBAEs at lower cost and time than traditional drug discovery. Sex-aware methods currently exist for clinical, genomic, and transcriptomic information [1, 7, 155]. They have not expanded to other data types, such as DNA variation, which has been beneficial in other drug repurposing methods that do not consider sex [114]. Additionally, some sex-aware methods suffer from poorer performance because a disproportionate number of male and female samples are available to train computational methods [7]. However, there is development potential for several different categories (i.e., data mining, ligand binding predictions, molecular associations, and networks). Low-dimensional representations of molecular association and network approaches are also especially promising candidates for future sex-aware drug repurposing methodologies because they reduce the multiple hypothesis testing burden and capture sex-specific variation better than the other methods [151, 159]. Here we review how sex influences drug response, the current state of drug repurposing including with respect to sex-bias drug response, and how model organism study design choices influence drug repurposing validation.

Evaluation of the binding performance of flavonoids to estrogen receptor alpha by Autodock, Autodock Vina and Surflex-Dock

Molecular docking is a widely used method to predict the binding modes of small-molecule ligands to the target binding site. However, it remains a challenge to identify the correct binding conformation and the corresponding binding affinity for a series of structurally similar ligands, especially those with weak binding. An understanding of the various relative attributes of popular docking programs is required to ensure a successful docking outcome. In this study, we systematically compared the performance of three popular docking programs, Autodock, Autodock Vina, and Surflex-Dock for a series of structurally similar weekly binding flavonoids (22) binding to the estrogen receptor alpha (ERα). For these flavonoids-ERα interactions, Surflex-Dock showed higher accuracy than Autodock and Autodock Vina. The hydrogen bond overweighting by Autodock and Autodock Vina led to incorrect binding results, while Surflex-Dock effectively balanced both hydrogen bond and hydrophobic interactions. Moreover, the selection of initial receptor structure is critical as it influences the docking conformations of flavonoids-ERα complexes. The flexible docking method failed to further improve the docking accuracy of the semi-flexible docking method for such chemicals. In addition, binding interaction analysis revealed that 8 residues, including Ala350, Glu353, Leu387, Arg394, Phe404, Gly521, His524, and Leu525, are the key residues in ERα-flavonoids complexes. This work provides reference for assessing molecular interactions between ERα and flavonoid-like chemicals and provides instructive information for other environmental chemicals.

Hypothesis-Driven, Structure-Based Design in Photopharmacology: The Case of eDHFR Inhibitors

Photopharmacology uses light to regulate the biological activity of drugs. This precise control is obtained through the incorporation of molecular photoswitches into bioactive molecules. A major challenge for photopharmacology is the rational design of photoswitchable drugs that show light-induced activation. Computer-aided drug design is an attractive approach toward more effective, targeted design. Herein, we critically evaluated different structure-based approaches for photopharmacology with Escherichia coli dihydrofolate reductase (eDHFR) as a case study. Through the iterative examination of our hypotheses, we progressively tuned the design of azobenzene-based, photoswitchable eDHFR inhibitors in five design–make–switch–test–analyze cycles. Targeting a hydrophobic subpocket of the enzyme and a specific salt bridge only with the thermally metastable cis-isomer emerged as the most promising design strategy. We identified three inhibitors that could be activated upon irradiation and reached potencies in the low-nanomolar range. Above all, this systematic study provided valuable insights for future endeavors toward rational photopharmacology.

Next-Generation Molecular Discovery: From Bottom-Up In Vivo and In Vitro Approaches to In Silico Top-Down Approaches for Therapeutics Neogenesis

Protein and drug engineering comprises a major part of the medical and research industries, and yet approaches to discovering and understanding therapeutic molecular interactions in biological systems rely on trial and error. The general approach to molecular discovery involves screening large libraries of compounds, proteins, or antibodies, or in vivo antibody generation, which could be considered “bottom-up” approaches to therapeutic discovery. In these bottom-up approaches, a minimal amount is known about the therapeutics at the start of the process, but through meticulous and exhaustive laboratory work, the molecule is characterised in detail. In contrast, the advent of “big data” and access to extensive online databases and machine learning technologies offers promising new avenues to understanding molecular interactions. Artificial intelligence (AI) now has the potential to predict protein structure at an unprecedented accuracy using only the genetic sequence. This predictive approach to characterising molecular structure—when accompanied by high-quality experimental data for model training—has the capacity to invert the process of molecular discovery and characterisation. The process has potential to be transformed into a top-down approach, where new molecules can be designed directly based on the structure of a target and the desired function, rather than performing screening of large libraries of molecular variants. This paper will provide a brief evaluation of bottom-up approaches to discovering and characterising biological molecules and will discuss recent advances towards developing top-down approaches and the prospects of this.

In Silico study for acyclovir, ganciclovir and its derivatives to fight the COVID-19: Molecular docking, DFT calculations, ADME and td-Molecular dynamics simulations

In the present work, we have designed three molecules, acyclovir (A), ganciclovir (G) and derivative of hydroxymethyl derivative of ganciclovir (CH₂OH of G, that is D) and investigated their biological potential against the Mpro of nCoV via in silico studies. Further, density functional theory (DFT) calculations of A, G and D were performed using Gaussian 16 on applying B3LYP under default condition to collect the information for the delocalization of electron density in their optimized geometry. Authors have also calculated various energies including free energy of A, G and D in Hartree per particle. It can be seen that D has the least free energy. As mentioned, the molecular docking of the A, G and D against the Mpro of nCoV was performed using iGemdock, an acceptable computational tool and the interaction has been studied in the form of physical data, that is, binding energy for A, G and D were calculated in kcal/mol. It can be seen the D showed effective binding, that is, maximum inhibition that A and G. For a better understanding for the inhibition of the Mpro of nCoV by A, G and D, temperature dependent molecular dynamics simulations were performed. Different trajectories like RMSD, RMSF, Rg and hydrogen bond were extracted and analyzed. The results of molecular docking of A, G and D corroborate with the td-MD simulations and hypothesized that D could be a promising candidate to inhibit the activity of Mpro of nCoV.

Benchmarking ensemble docking methods in D3R Grand Challenge 4

The discovery of new drugs is a time consuming and expensive process. Methods such as virtual screening, which can filter out ineffective compounds from drug libraries prior to expensive experimental study, have become popular research topics. As the computational drug discovery community has grown, in order to benchmark the various advances in methodology, organizations such as the Drug Design Data Resource have begun hosting blinded grand challenges seeking to identify the best methods for ligand pose-prediction, ligand affinity ranking, and free energy calculations. Such open challenges offer a unique opportunity for researchers to partner with junior students (e.g., high school and undergraduate) to validate basic yet fundamental hypotheses considered to be uninteresting to domain experts. Here, we, a group of high school-aged students and their mentors, present the results of our participation in Grand Challenge 4 where we predicted ligand affinity rankings for the Cathepsin S protease, an important protein target for autoimmune diseases. To investigate the effect of incorporating receptor dynamics on ligand affinity rankings, we employed the Relaxed Complex Scheme, a molecular docking method paired with molecular dynamics-generated receptor conformations. We found that Cathepsin S is a difficult target for molecular docking and we explore some advanced methods such as distance-restrained docking to try to improve the correlation with experiments. This project has exemplified the capabilities of high school students when supported with a rigorous curriculum, and demonstrates the value of community-driven competitions for beginners in computational drug discovery.

Importance of Binding Site Hydration and Flexibility Revealed When Optimizing a Macrocyclic Inhibitor of the Keap1-Nrf2 Protein-Protein Interaction

Upregulation of the transcription factor Nrf2 by inhibition of the interaction with its negative regulator Keap1 constitutes an opportunity for the treatment of disease caused by oxidative stress. We report a structurally unique series of nanomolar Keap1 inhibitors obtained from a natural product-derived macrocyclic lead. Initial exploration of the structure-activity relationship of the lead, followed by structure-guided optimization, resulted in a 100-fold improvement in inhibitory potency. The macrocyclic core of the nanomolar inhibitors positions three pharmacophore units for productive interactions with key residues of Keap1, including R415, R483, and Y572. Ligand optimization resulted in the displacement of a coordinated water molecule from the Keap1 binding site and a significantly altered thermodynamic profile. In addition, minor reorganizations of R415 and R483 were accompanied by major differences in affinity between ligands. This study therefore indicates the importance of accounting both for the hydration and flexibility of the Keap1 binding site when designing high-affinity ligands.

Design, Synthesis, and Biological Evaluation of 4,4’-Difluorobenzhydrol Carbamates as Selective M1 Antagonists

Due to their important role in mediating a broad range of physiological functions, muscarinic acetylcholine receptors (mAChRs) have been a promising target for therapeutic and diagnostic applications alike; however, the list of truly subtype-selective ligands is scarce. Within this work, we have identified a series of twelve 4,4’-difluorobenzhydrol carbamates through a rigorous docking campaign leveraging commercially available amine databases. After synthesis, these compounds have been evaluated for their physico–chemical property profiles, including characteristics such as HPLC-logD, tPSA, logBB, and logPS. For all the synthesized carbamates, these characteristics indicate the potential for BBB permeation. In competitive radioligand binding experiments using Chinese hamster ovary cell membranes expressing the individual human mAChR subtype hM1-hM5, the most promising compound 2 displayed a high binding affinitiy towards hM1R (1.2 nM) while exhibiting modest-to-excellent selectivity versus the hM2-5R (4–189-fold). All 12 compounds were shown to act in an antagonistic fashion towards hM1R using a dose-dependent calcium mobilization assay. The structural eligibility for radiolabeling and their pharmacological and physico–chemical property profiles render compounds 2, 5, and 7 promising candidates for future position emission tomography (PET) tracer development.

Ultralarge Virtual Screening Identifies SARS-CoV-2 Main Protease Inhibitors with Broad-Spectrum Activity against Coronaviruses

Drugs targeting SARS-CoV-2 could have saved millions of lives during the COVID-19 pandemic, and it is now crucial to develop inhibitors of coronavirus replication in preparation for future outbreaks. We explored two virtual screening strategies to find inhibitors of the SARS-CoV-2 main protease in ultralarge chemical libraries. First, structure-based docking was used to screen a diverse library of 235 million virtual compounds against the active site. One hundred top-ranked compounds were tested in binding and enzymatic assays. Second, a fragment discovered by crystallographic screening was optimized guided by docking of millions of elaborated molecules and experimental testing of 93 compounds. Three inhibitors were identified in the first library screen, and five of the selected fragment elaborations showed inhibitory effects. Crystal structures of target-inhibitor complexes confirmed docking predictions and guided hit-to-lead optimization, resulting in a noncovalent main protease inhibitor with nanomolar affinity, a promising in vitro pharmacokinetic profile, and broad-spectrum antiviral effect in infected cells.

Artificial intelligence-enabled virtual screening of ultra-large chemical libraries with deep docking

With the recent explosion of chemical libraries beyond a billion molecules, more efficient virtual screening approaches are needed. The Deep Docking (DD) platform enables up to 100-fold acceleration of structure-based virtual screening by docking only a subset of a chemical library, iteratively synchronized with a ligand-based prediction of the remaining docking scores. This method results in hundreds- to thousands-fold virtual hit enrichment (without significant loss of potential drug candidates) and hence enables the screening of billion molecule-sized chemical libraries without using extraordinary computational resources. Herein, we present and discuss the generalized DD protocol that has been proven successful in various computer-aided drug discovery (CADD) campaigns and can be applied in conjunction with any conventional docking program. The protocol encompasses eight consecutive stages: molecular library preparation, receptor preparation, random sampling of a library, ligand preparation, molecular docking, model training, model inference and the residual docking. The standard DD workflow enables iterative application of stages 3-7 with continuous augmentation of the training set, and the number of such iterations can be adjusted by the user. A predefined recall value allows for control of the percentage of top-scoring molecules that are retained by DD and can be adjusted to control the library size reduction. The procedure takes 1-2 weeks (depending on the available resources) and can be completely automated on computing clusters managed by job schedulers. This open-source protocol, at https://github.com/jamesgleave/DD_protocol , can be readily deployed by CADD researchers and can significantly accelerate the effective exploration of ultra-large portions of a chemical space.

Virtual Screening Directly Identifies New Fragment-Sized Inhibitors of Carboxylesterase Notum with Nanomolar Activity

Notum is a negative regulator of Wnt signaling acting through the hydrolysis of a palmitoleoylate ester, which is required for Wnt activity. Inhibitors of Notum could be of use in diseases where dysfunctional Notum activity is an underlying cause. A docking-based virtual screen (VS) of a large commercial library was used to shortlist 952 compounds for experimental validation as inhibitors of Notum. The VS was successful with 31 compounds having an IC50 < 500 nM. A critical selection process was then applied with two clusters and two singletons (1-4d) selected for hit validation. Optimization of 4d guided by structural biology identified potent inhibitors of Notum activity that restored Wnt/β-catenin signaling in cell-based models. The [1,2,4]triazolo[4,3-b]pyradizin-3(2H)-one series 4 represent a new chemical class of Notum inhibitors and the first to be discovered by a VS campaign. These results demonstrate the value of VS with well-designed docking models based on X-ray structures.

Anticancer potential of β-sitosterol and oleanolic acid as through inhibition of human estrogenic 17beta-hydroxysteroid dehydrogenase type-1 based on an in silico approach

We presented pharmacokinetic study, molecular docking, and MD simulation to study β-sitosterol and oleanolic acid compounds and potential HSD17B1 inhibitors.

Virtual screening and in vitro validation of natural compound inhibitors against SARS-CoV-2 spike protein

The COVID-19 pandemic caused by the SARS-CoV-2 virus has led to a major public health burden and has resulted in millions of deaths worldwide. As effective treatments are limited, there is a significant requirement for high-throughput, low resource methods for the discovery of novel antivirals. The SARS-CoV-2 spike protein plays a key role in viral entry and has been identified as a therapeutic target. Using the available spike crystal structure, we performed a virtual screen with a library of 527 209 natural compounds against the receptor binding domain of this protein. Top hits from this screen were subjected to a second, more comprehensive molecular docking experiment and filtered for favourable ADMET properties. The in vitro activity of 10 highly ranked compounds was assessed using a virus neutralisation assay designed to facilitate viral entry in a physiologically relevant manner via the plasma membrane route. Subsequently, four compounds ZINC02111387, ZINC02122196, SN00074072 and ZINC04090608 were identified to possess antiviral activity in the µM range. These findings validate the virtual screening method as a tool for identifying novel antivirals and provide a basis for future drug development against SARS-CoV-2.

New Inhibitors of Laccase and Tyrosinase by Examination of Cross-Inhibition between Copper-Containing Enzymes

Coppers play crucial roles in the maintenance homeostasis in living species. Approximately 20 enzyme families of eukaryotes and prokaryotes are known to utilize copper atoms for catalytic activities. However, small-molecule inhibitors directly targeting catalytic centers are rare, except for those that act against tyrosinase and dopamine-β-hydroxylase (DBH). This study tested whether known tyrosinase inhibitors can inhibit the copper-containing enzymes, ceruloplasmin, DBH, and laccase. While most small molecules minimally reduced the activities of ceruloplasmin and DBH, aside from known inhibitors, 5 of 28 tested molecules significantly inhibited the function of laccase, with the K_i values in the range of 15 to 48 µM. Enzyme inhibitory kinetics classified the molecules as competitive inhibitors, whereas differential scanning fluorimetry and fluorescence quenching supported direct bindings. To the best of our knowledge, this is the first report on organic small-molecule inhibitors for laccase. Comparison of tyrosinase and DBH inhibitors using cheminformatics predicted that the presence of thione moiety would suffice to inhibit tyrosinase. Enzyme assays confirmed this prediction, leading to the discovery of two new dual tyrosinase and DBH inhibitors.

Virtual Screening in Search for a Chemical Probe for Angiotensin-Converting Enzyme 2 (ACE2)

We elaborate new models for ACE and ACE2 receptors with an excellent prediction power compared to previous models. We propose promising workflows for working with huge compound collections, thereby enabling us to discover optimized protocols for virtual screening management. The efficacy of elaborated roadmaps is demonstrated through the cost-effective molecular docking of 1.4 billion compounds. Savings of up to 10-fold in CPU time are demonstrated. These developments allowed us to evaluate ACE2/ACE selectivity in silico, which is a crucial checkpoint for developing chemical probes for ACE2.