Refinement and Rescoring of Virtual Screening Results

Assessing the Martini 3 protein model: A review of its path and potential

Coarse-grained (CG) protein models have become indispensable tools for studying many biological protein details, from conformational dynamics to the organization of protein macro-complexes, and even the interaction of proteins with other molecules. The Martini force field is one of the most widely used CG models for bio-molecular simulations, partly because of the enormous success of its protein model. With the recent release of a new and improved version of the Martini force field - Martini 3 - a new iteration of its protein model was also made available. The Martini 3 protein force field is an evolution of its Martini 2 counterpart, aimed at improving many of the shortcomings that had been previously identified. In this mini-review, we first provide a general overview of the model and then focus on the successful advances made in the short time since its release, many of which would not have been possible before. Furthermore, we discuss reported limitations, potential directions for model improvement and comment on what the likely future development and application avenues are.

Harnessing Brazilian biodiversity database: identification of flavonoids as potential inhibitors of SARS-CoV-2 main protease using computational approaches and all-atom molecular dynamics simulation

SARS-CoV-2 (Severe Acute Respiratory Syndrome Coronavirus 2) is the etiological agent responsible for the global outbreak of COVID-19 (Coronavirus Disease 2019). The main protease of SARS-CoV-2, Mpro, is a key enzyme that plays a vital role in mediating viral replication and transcription. In this study, a comprehensive computational approach was employed to investigate the binding affinity, selectivity, and stability of natural product candidates as potential new antivirals acting on the viral polyprotein processing mediated by SARS-CoV-2 Mpro. A library of 288 flavonoids extracted from Brazilian biodiversity was screened to select potential Mpro inhibitors. An initial filter based on Lipinski's rule of five was applied, and 204 compounds that did not violate any of the Lipinski rules were selected. The compounds were then docked into the active site of Mpro using the GOLD program, and the poses were subsequently re-scored using MM-GBSA (Molecular Mechanics Generalized Born Surface Area) binding free energy calculations performed by AmberTools23. The top five flavonoids with the best MM-GBSA binding free energy values were selected for analysis of their interactions with the active site residues of the protein. Next, we conducted a toxicity and drug-likeness analysis, and non-toxic compounds were subjected to molecular dynamics simulation and free energy calculation using the MM-PBSA (Molecular Mechanics Poisson-Boltzmann Surface Area) method. It was observed that the five selected flavonoids had lower MM-GBSA binding free energy with Mpro than the co-crystal ligand. Furthermore, these compounds also formed hydrogen bonds with two important residues, Cys145 and Glu166, in the active site of Mpro. Two compounds that passed the drug-likeness filter showed stable conformations during the molecular dynamics simulations. Among these, NuBBE_867 exhibited the best MM-PBSA binding free energy value compared to the crystallographic inhibitor. Therefore, this study suggests that NuBBE_867 could be a potential inhibitor against the main protease of SARS-CoV-2 and may be further examined to confirm our results.

Structure-Based Discovery of Mouse Trace Amine-Associated Receptor 5 Antagonists

Trace amine-associated receptors (TAARs) were discovered in 2001 as new members of class A G protein-coupled receptors (GPCRs). With the only exception of TAAR1, TAAR members (TAAR2-9, also known as noncanonical olfactory receptors) were originally described exclusively in the olfactory epithelium and believed to mediate the innate perception of volatile amines. However, most noncanonical olfactory receptors are still orphan receptors. Given its recently discovered nonolfactory expression and therapeutic potential, TAAR5 has been the focus of deorphanization campaigns that led to the discovery of a few druglike antagonists. Here, we report four novel TAAR5 antagonists identified through high-throughput screening, which, along with the four ligands published in the literature, constituted our starting point to design a computational strategy for the identification of TAAR5 ligands. We developed a structure-based virtual screening protocol that allowed us to identify three new TAAR5 antagonists with a hit rate of 10%. Despite lacking an experimental structure, we accurately modeled the TAAR5 binding site by integrating comparative sequence- and structure-based analyses of serotonin receptors with homology modeling and side-chain optimization. In summary, we have identified seven new TAAR5 antagonists that could serve as lead candidates for the development of new treatments for depression, anxiety, and neurodegenerative diseases.

Insights into molecular docking and dynamics to reveal therapeutic potential of natural compounds against P53 protein

P53 is eminent tumour suppressor protein that plays a prominent role in cell cycle arrest, DNA repair, senescence, differentiation and initiation of apoptosis. P53 is an attractive drug target and the high toxicity of some cancer chemotherapy drugs increase the demand for new anti-cancer drugs from natural products. In this current scenario, identification of promising anticancer compounds from natural sources by repurposing approach is still relevant for the early prevention and effective management of cancer. In present study, we docked natural compounds like podophyllotoxin, quercetin and rutin along standard drugs (MG-132 and Bay 61-3606) against p53 protein. ADME/T analysis predicted toxicity of phytochemicals and drugs. In silico docking analysis of podophyllotoxin, quercetin and rutin gave HDOCK docking scores of -187.87, -148. 97 and -143.85, whereas control drugs MG-132 and Bay 61-3606 showed docking scores of -159.59 and -140.71 against p53 respectively. AutoDock analysis of rutin and MG-132 showed highest binding affinity scores of -7.3 and -6.8 kcal/mol against p53. Molecular dynamic simulation for p53 protein displayed stable conformation and convergence. In this study, P53-rutin complex showed free binding energy score of 11.84 kcal/mol and P53-MG-132 complex reported free energy score of 16.3 kcal/mol. Protein contacts atlas gives non-covalent contacts framework by exploring interfaces of individual subunits and protein-ligand interactions. STRING tool predicts physical and functional interactions between proteins. The results of this study revealed that rutin and MG-132 could be promising inhibitors against targeted p53 protein and this could prove detrimental for molecular therapeutics and drug-designing strategies.Communicated by Ramaswamy H. Sarma.

In Silico Strategies to Predict Anti-aging Features of Whey Peptides

We have analysed the in silico potential of bioactive peptides from cheese whey, the most relevant by-product from the dairy industry, to bind into the active site of collagenase and elastase. The peptides generated from the hydrolysis of bovine β-lactoglobulin with three proteases (trypsin, chymotrypsin, and subtilisin) were docked onto collagenase and elastase by molecular docking. The interaction models were ranked according to their free binding energy using molecular dynamics simulations, which showed that most complexes presented favourable interactions. Interactions with elastase had significantly lower binding energies than those with collagenase. Regarding the interaction site, it was found that four bioactive peptides were positioned in collagenase's active site, while six were found in elastase's active site. Among these, the most we have found one promising collagen-binding peptide produced by chymotrypsin and two for elastase, produced by subtilisin and chymotrypsin. These in silico results can be used as a tool for designing further experiments aiming at testing the in vitro potential of the peptides found in this work.

D3CARP: a comprehensive platform with multiple-conformation based docking, ligand similarity search and deep learning approaches for target prediction and virtual screening

Resource- and time-consuming biological experiments are unavoidable in traditional drug discovery, which have directly driven the evolution of various computational algorithms and tools for drug-target interaction (DTI) prediction. For improving the prediction reliability, a comprehensive platform is highly expected as some previously reported webservers are small in scale, single-method, or even out of service. In this study, we integrated the multiple-conformation based docking, 2D/3D ligand similarity search and deep learning approaches to construct a comprehensive webserver, namely D3CARP, for target prediction and virtual screening. Specifically, 9352 conformations with positive control of 1970 targets were used for molecular docking, and approximately 2 million target-ligand pairs were used for 2D/3D ligand similarity search and deep learning. Besides, the positive compounds were added as references, and related diseases of therapeutic targets were annotated for further disease-based DTI study. The accuracies of the molecular docking and deep learning approaches were 0.44 and 0.89, respectively. And the average accuracy of five ligand similarity searches was 0.94. The strengths of D3CARP encompass the support for multiple computational methods, ensemble docking, utilization of positive controls as references, cross-validation of predicted outcomes, diverse disease types, and broad applicability in drug discovery. The D3CARP is freely accessible at https://www.d3pharma.com/D3CARP/index.php.

In silico analysis and molecular docking studies of natural compounds of Withania somnifera against bovine NLRP9

CONTEXT

NLRP9 is a member of nucleotide-binding domain leucine-rich repeat-containing receptors and is found to be associated with many inflammatory diseases. In the current scenario, the identification of promising anti-inflammatory compounds from natural sources by repurposing approach is still relevant for the early prevention and effective management of the disease.

METHODS

In the present study, we docked bioactives of Ashwagandha (Withanoside IV, Withanoside V, Withanolide A, Withanolide B, and Sitoindoside IX) and two control drugs against bovine NLRP9 protein. ADME/T analysis was used to determine the physiochemical properties of compounds and standard drugs. Molecular modeling was used to evaluate the correctness and quality of protein structures. In silico docking analysis revealed Withanolide B had the highest binding affinity score of -10.5 kcal/mol, whereas, among control drugs, doxycycline hydrochloride was most effective (-10.3 kcal/mol). The results of this study revealed that bioactives of Withania somnifera could be promising inhibitors against bovine NLRP9. In the present study, molecular simulation was used to measure protein conformational changes over time. The Rg value was found to be 34.77A°. RMSD and B-factor were also estimated to provide insights into the flexibility and mobile regions of protein structure. A functional protein network interaction was constructed from information collected from non-curative sources as protein-protein interactions (PPI) that play an important role in determining the function of the target protein and the ability of the drug molecule. Thus, in the present situation, it is important to identify bioactives with the potential to combat inflammatory diseases and provide strength and immunity to the host. However, there is still a need to study in vitro and in vivo to further support these findings.

Past, Present, and Future Perspectives on Computer-Aided Drug Design Methodologies

The application of computational approaches in drug discovery has been consolidated in the last decades. These families of techniques are usually grouped under the common name of "computer-aided drug design" (CADD), and they now constitute one of the pillars in the pharmaceutical discovery pipelines in many academic and industrial environments. Their implementation has been demonstrated to tremendously improve the speed of the early discovery steps, allowing for the proficient and rational choice of proper compounds for a desired therapeutic need among the extreme vastness of the drug-like chemical space. Moreover, the application of CADD approaches allows the rationalization of biochemical and interactive processes of pharmaceutical interest at the molecular level. Because of this, computational tools are now extensively used also in the field of rational 3D design and optimization of chemical entities starting from the structural information of the targets, which can be experimentally resolved or can also be obtained with other computer-based techniques. In this work, we revised the state-of-the-art computer-aided drug design methods, focusing on their application in different scenarios of pharmaceutical and biological interest, not only highlighting their great potential and their benefits, but also discussing their actual limitations and eventual weaknesses. This work can be considered a brief overview of computational methods for drug discovery.

AI-accelerated protein-ligand docking for SARS-CoV-2 is 100-fold faster with no significant change in detection

Protein-ligand docking is a computational method for identifying drug leads. The method is capable of narrowing a vast library of compounds down to a tractable size for downstream simulation or experimental testing and is widely used in drug discovery. While there has been progress in accelerating scoring of compounds with artificial intelligence, few works have bridged these successes back to the virtual screening community in terms of utility and forward-looking development. We demonstrate the power of high-speed ML models by scoring 1 billion molecules in under a day (50 k predictions per GPU seconds). We showcase a workflow for docking utilizing surrogate AI-based models as a pre-filter to a standard docking workflow. Our workflow is ten times faster at screening a library of compounds than the standard technique, with an error rate less than 0.01% of detecting the underlying best scoring 0.1% of compounds. Our analysis of the speedup explains that another order of magnitude speedup must come from model accuracy rather than computing speed. In order to drive another order of magnitude of acceleration, we share a benchmark dataset consisting of 200 million 3D complex structures and 2D structure scores across a consistent set of 13 million "in-stock" molecules over 15 receptors, or binding sites, across the SARS-CoV-2 proteome. We believe this is strong evidence for the community to begin focusing on improving the accuracy of surrogate models to improve the ability to screen massive compound libraries 100 × or even 1000 × faster than current techniques and reduce missing top hits. The technique outlined aims to be a fast drop-in replacement for docking for screening billion-scale molecular libraries.

Fragmented blind docking: a novel protein-ligand binding prediction protocol

In the present paper we propose a novel blind docking protocol based on Autodock-Vina. The developed docking protocol can provide binding site identification and binding pose prediction at the same time, by a systematical exploration of the protein volume performed with several preliminary docking calculations. In our opinion, this protocol can be successfully applied during the first steps of the virtual screening pipeline, because it provides binding site identification and binding pose prediction at the same time without visual evaluation of the binding site. After the binding pose prediction, MM/GBSA re-scoring rescoring procedures has been applied to improve the accuracy of the protein-ligand bound state. The FRAD protocol has been tested on 116 protein-ligand complexes of the Heat Shock Protein 90 - alpha, on 176 of Human Immunodeficiency virus protease 1, and on more than 100 protein-ligand system taken from the PDBbind dataset. Overall, the FRAD approach combined to MM/GBSA re-scoring can be considered as a powerful tool to increase the accuracy and efficiency with respect to other standard docking approaches when the ligand-binding site is unknown.Communicated by Ramaswamy H. Sarma.

Application of Computational Biology and Artificial Intelligence in Drug Design

Traditional drug design requires a great amount of research time and developmental expense. Booming computational approaches, including computational biology, computer-aided drug design, and artificial intelligence, have the potential to expedite the efficiency of drug discovery by minimizing the time and financial cost. In recent years, computational approaches are being widely used to improve the efficacy and effectiveness of drug discovery and pipeline, leading to the approval of plenty of new drugs for marketing. The present review emphasizes on the applications of these indispensable computational approaches in aiding target identification, lead discovery, and lead optimization. Some challenges of using these approaches for drug design are also discussed. Moreover, we propose a methodology for integrating various computational techniques into new drug discovery and design.

Protocol for iterative optimization of modified peptides bound to protein targets

Peptides are commonly used as therapeutic agents. However, they suffer from easy degradation and instability. Replacing natural by non-natural amino acids can avoid these problems, and potentially improve the affinity towards the target protein. Here, we present a computational pipeline to optimize peptides based on adding non-natural amino acids while improving their binding affinity. The workflow is an iterative computational evolution algorithm, inspired by the PARCE protocol, that performs single-point mutations on the peptide sequence using modules from the Rosetta framework. The modifications can be guided based on the structural properties or previous knowledge of the biological system. At each mutation step, the affinity to the protein is estimated by sampling the complex conformations and applying a consensus metric using various open protein-ligand scoring functions. The mutations are accepted based on the score differences, allowing for an iterative optimization of the initial peptide. The sampling/scoring scheme was benchmarked with a set of protein-peptide complexes where experimental affinity values have been reported. In addition, a basic application using a known protein-peptide complex is also provided. The structure- and dynamic-based approach allows users to optimize bound peptides, with the option to personalize the code for further applications. The protocol, called mPARCE, is available at: https://github.com/rochoa85/mPARCE/.

An Electrostatically-steered Conformational Selection Mechanism Promotes SARS-CoV-2 Spike Protein Variation

After two years since the outbreak, the COVID-19 pandemic remains a global public health emergency. SARS-CoV-2 variants with substitutions on the spike (S) protein emerge increasing the risk of immune evasion and cross-species transmission. Here, we analyzed the evolution of the S protein as recorded in 276,712 samples collected before the start of vaccination efforts. Our analysis shows that most variants destabilize the S protein trimer, increase its conformational heterogeneity and improve the odds of the recognition by the host cell receptor. Most frequent substitutions promote overall hydrophobicity by replacing charged amino acids, reducing stabilizing local interactions in the unbound S protein trimer. Moreover, our results identify "forbidden" regions that rarely show any sequence variation, and which are related to conformational changes occurring upon fusion. These results are significant for understanding the structure and function of SARS-CoV-2 related proteins which is a critical step in vaccine development and epidemiological surveillance.

Computational Investigations on the Natural Small Molecule as an Inhibitor of Programmed Death Ligand 1 for Cancer Immunotherapy

Several therapeutic monoclonal antibodies approved by the FDA are available against the PD-1/PD-L1 (programmed death 1/programmed death ligand 1) immune checkpoint axis, which has been an unprecedented success in cancer treatment. However, existing therapeutics against PD-L1, including small molecule inhibitors, have certain drawbacks such as high cost and drug resistance that challenge the currently available anti-PD-L1 therapy. Therefore, this study presents the screening of 32,552 compounds from the Natural Product Atlas database against PD-L1, including three steps of structure-based virtual screening followed by binding free energy to refine the ideal conformation of potent PD-L1 inhibitors. Subsequently, five natural compounds, i.e., Neoenactin B1, Actinofuranone I, Cosmosporin, Ganocapenoid A, and 3-[3-hydroxy-4-(3-methylbut-2-enyl)phenyl]-5-(4-hydroxybenzyl)-4-methyldihydrofuran-2(3H)-one, were collected based on the ADMET (absorption, distribution, metabolism, excretion, and toxicity) profiling and binding free energy (>−60 kcal/mol) for further computational investigation in comparison to co-crystallized ligand, i.e., JQT inhibitor. Based on interaction mapping, explicit 100 ns molecular dynamics simulation, and end-point binding free energy calculations, the selected natural compounds were marked for substantial stability with PD-L1 via intermolecular interactions (hydrogen and hydrophobic) with essential residues in comparison to the JQT inhibitor. Collectively, the calculated results advocate the selected natural compounds as the putative potent inhibitors of PD-L1 and, therefore, can be considered for further development of PD-L1 immune checkpoint inhibitors in cancer immunotherapy.

Computational strategies to identify new drug candidates against neuroinflammation

The even more increasing application of computational approaches in these last decades has deeply modified the process of discovery and commercialization of new therapeutic entities. This is especially true in the field of neuroinflammation, in which both the peculiar anatomical localization and the presence of the blood-brain barrier makeit mandatory to finely tune the candidates' physicochemical properties from the early stages of the discovery pipeline. The aim of this review is therefore to provide a general overview to the readers about the topic of neuroinflammation, together with the most common computational strategies that can be exploited to discover and design small molecules controlling neuroinflammation, especially those based on the knowledge of the three-dimensional structure of the biological targets of therapeutic interest. The techniques used to describe the molecular recognition mechanisms, such as molecular docking and molecular dynamics, will therefore be eviscerated, highlighting their advantages and their limitations. Finally, we report several case studies in which computational methods have been applied in drug discovery on neuroinflammation, focusing on the last decade's research.

Computer-Aided Drug Design (CADD) to De-Orphanize Marine Molecules: Finding Potential Therapeutic Agents for Neurodegenerative and Cardiovascular Diseases

Computer-aided drug design (CADD) techniques allow the identification of compounds capable of modulating protein functions in pathogenesis-related pathways, which is a promising line on drug discovery. Marine natural products (MNPs) are considered a rich source of bioactive compounds, as the oceans are home to much of the planet's biodiversity. Biodiversity is directly related to chemodiversity, which can inspire new drug discoveries. Therefore, natural products (NPs) in general, and MNPs in particular, have been used for decades as a source of inspiration for the design of new drugs. However, NPs present both opportunities and challenges. These difficulties can be technical, such as the need to dive or trawl to collect the organisms possessing the compounds, or biological, due to their particular marine habitats and the fact that they can be uncultivable in the laboratory. For all these difficulties, the contributions of CADD can play a very relevant role in simplifying their study, since, for example, no biological sample is needed to carry out an in-silico analysis. Therefore, the amount of natural product that needs to be used in the entire preclinical and clinical study is significantly reduced. Here, we exemplify how this combination between CADD and MNPs can help unlock their therapeutic potential. In this study, using a set of marine invertebrate molecules, we elucidate their possible molecular targets and associated therapeutic potential, establishing a pipeline that can be replicated in future studies.

Antifungal Activity of N-(4-Halobenzyl)amides against Candida spp. and Molecular Modeling Studies

Fungal infections remain a high-incidence worldwide health problem that is aggravated by limited therapeutic options and the emergence of drug-resistant strains. Cinnamic and benzoic acid amides have previously shown bioactivity against different species belonging to the Candida genus. Here, 20 cinnamic and benzoic acid amides were synthesized and tested for inhibition of C. krusei ATCC 14243 and C. parapsilosis ATCC 22019. Five compounds inhibited the Candida strains tested, with compound 16 (MIC = 7.8 µg/mL) producing stronger antifungal activity than fluconazole (MIC = 16 µg/mL) against C. krusei ATCC 14243. It was also tested against eight Candida strains, including five clinical strains resistant to fluconazole, and showed an inhibitory effect against all strains tested (MIC = 85.3–341.3 µg/mL). The MIC value against C. krusei ATCC 6258 was 85.3 mcg/mL, while against C. krusei ATCC 14243, it was 10.9 times smaller. This strain had greater sensitivity to the antifungal action of compound 16. The inhibition of C. krusei ATCC 14243 and C. parapsilosis ATCC 22019 was also achieved by compounds 2, 9, 12, 14 and 15. Computational experiments combining target fishing, molecular docking and molecular dynamics simulations were performed to study the potential mechanism of action of compound 16 against C. krusei. From these, a multi-target mechanism of action is proposed for this compound that involves proteins related to critical cellular processes such as the redox balance, kinases-mediated signaling, protein folding and cell wall synthesis. The modeling results might guide future experiments focusing on the wet-lab investigation of the mechanism of action of this series of compounds, as well as on the optimization of their inhibitory potency.

Navigating Chemical Space by Interfacing Generative Artificial Intelligence and Molecular Docking

Here, we report the implementation and application of a simple, structure-aware framework to generate target-specific screening libraries. Our approach combines advances in generative artificial intelligence (AI) with conventional molecular docking to explore chemical space conditioned on the unique physicochemical properties of the active site of a biomolecular target. As a demonstration, we used our framework, which we refer to as sample-and-dock, to construct focused libraries for cyclin-dependent kinase type-2 (CDK2) and the active site of the main protease (M^pro) of the SARS-CoV-2 virus. We envision that the sample-and-dock framework could be used to generate theoretical maps of the chemical space specific to a given target and so provide information about its molecular recognition characteristics.

Key criteria for engineering mycotoxin binding aptamers via computational simulations: Aflatoxin B1 as a case study

Due to the difficulties in monoclonal antibody production specific to mycotoxins, aptameric probes have been considered as suitable alternatives. The low efficiency of the SELEX procedure in screening high affinity aptamers for binding mycotoxins as small molecules can be significantly improved through computational techniques. Previously, we designed five new aptamers to aflatoxin B₁ (AFB1) based on a known aptamer sequence (Patent: PCT/CA2010/001 292, Apt1) through a genetic algorithm-based in silico maturation strategy and experimentally measured their affinity to the target toxin. Here, integrated molecular dynamic simulation (MDs) studies with molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) analysis to clarify the binding modes, critical interacting nucleic bases and energy component contributions in the six AFB1-binding aptamers. The aptamer F20, which was selected in the first work, showed the best free binding energy and complex stability compared to other aptamers. The trajectory analysis revealed that AFB1 recognized F20 through the groove binding mode along with precise shape complementarity. The MD simulation results revealed that dynamic water intermediate interactions also play a key role in promoting complex stability. According to the MM-PBSA calculations, van der Waals contacts were identified as dominant energy components in all complexes. Interestingly, a high consistency is observed between the experimentally obtained binding affinities of the six aptamers with their free energy solvation. The computational findings, confirmed via previous experiments, highlighted the binding modes, the dynamic hydration of complex components and the total free interacting energy as the crucial criteria in discovering high functional aptameric probes.

Cytotoxic and Antifungal Amides Derived from Ferulic Acid: Molecular Docking and Mechanism of Action

Amides derived from ferulic acid have a wide spectrum of pharmacological activities, including antitumor and antifungal activity. In the present study, a series of ten amides were obtained by coupling reactions using the reagents (benzotriazol-1-yloxy) tripyrrolidinophosphonium hexafluorophosphate (PyBOP) and N,N′-dicyclohexylcarbodiimide (DCC). All the compounds were identified on the basis of their IR, ¹H- and ¹³C-NMR, HRMS data, and with yields ranging from 43.17% to 91.37%. The compounds were subjected to cytotoxic tests by the alamar blue technique and antifungal screening by the broth microdilution method to determine the minimum inhibitory concentration (MIC). The amides 10 and 11 displayed the best result in both biological evaluations, and compound 10 was the most potent and selective in HL-60 cancer cells, with no cytotoxicity on healthy cells. This amide had antifungal activity in all strains and had the lowest MIC against Candida albicans and Candida tropicalis. The possible mechanism of antifungal action occurs via the fungal cell wall. Molecular modeling suggested that compounds 10 and 11 interact with the enzymes GWT1 and GSC1, which are essential for the development of C. albicans. The findings of the present study demonstrated that compounds 10 and 11 may be used as a platform in drug development in the future.

Identification of Novel CD39 Inhibitors Based on Virtual Screening and Enzymatic Assays

The accumulation of adenosine in the tumor microenvironment mediates immunosuppression and promotes tumor growth and proliferation. Intervention of the adenosine pathway is an important direction of antitumor immunity research. CD39 is an important ecto-nucleotidases for adenosine generation, therefore targeting the CD39-adenosine pathway is an emerging immune checkpoint for anticancer treatment. However, currently no CD39 inhibitor has been approved by the U.S. Food and Drug Administration. The development of CD39 drugs is urgent for clinical application. In this study, we combined homology modeling, virtual screening, and in vitro enzymatic activity to characterize the structural features of the CD39 protein and identify a triazinoindole-based compound as a CD39 inhibitor. The identified inhibitor and one of its analogues could effectively prevent the enzymatic activity of CD39 with IC₅₀ values of 27.42 ± 5.52 and 79.24 ± 12.21 μM, respectively. At the same time, the inhibitor significantly inhibited the adenosine monophosphate production in colorectal cancer cell lines (HT29 and MC38) and thereafter prevented cell proliferation. Molecular docking studies, mutagenesis, and microscale thermophoresis indicated that residues such as R85 could be the main contributor in binding triazinoindole compounds. The binding mode can potentially be utilized for hit-to-lead optimization, and the identified inhibitor can be further tested for its anticancer activity in vivo or may serve as a chemical agent to study CD39-related functions.

In silico identification of noncompetitive inhibitors targeting an uncharacterized allosteric site of falcipain-2

Falcipain-2 (FP-2) is a Plasmodium falciparum hemoglobinase widely targeted in the search for antimalarials. FP-2 can be allosterically modulated by various noncompetitive inhibitors that have been serendipitously identified. Moreover, the crystal structures of two inhibitors bound to an allosteric site, termed site 6, of the homolog enzyme human cathepsin K (hCatK) suggest that the equivalent region in FP-2 might play a similar role. Here, we conduct the rational identification of FP-2 inhibitors through virtual screenings (VS) of compounds into several pocket-like conformations of site 6, sampled during molecular dynamics (MD) simulations of the free enzyme. Two noncompetitive inhibitors, ZINC03225317 and ZINC72290660, were confirmed using in vitro enzymatic assays and their poses into site 6 led to calculated binding free energies matching the experimental ones. Our results provide strong evidence about the allosteric inhibition of FP-2 through binding of small molecules to site 6, thus opening the way toward the discovery of new inhibitors against this enzyme.

Direct Keap1-kelch inhibitors as potential drug candidates for oxidative stress-orchestrated diseases: A review on In silico perspective

The recent outcry in the search for direct keap1 inhibitors requires a quicker and more effective drug discovery process which is an inherent property of the Computer Aided Drug Discovery (CADD) to bring drug candidates into the clinic for patient's use. This Keap1 (negative regulator of ARE master activator) is emerging as a therapeutic strategy to combat oxidative stress-orchestrated diseases. The advances in computer algorithm and compound databases require that we highlight the functionalities that this technology possesses that can be exploited to target Keap1-Nrf2 PPI. Therefore, in this review, we uncover the in silico approaches that had been exploited towards the identification of keap1 inhibition in the light of appropriate fitting with relevant amino acid residues, we found 3 and 16 other compounds that perfectly fit keap1 kelch pocket/domain. Our goal is to harness the parameters that could orchestrate keap1 surface druggability by utilizing hotspot regions for virtual fragment screening and identification of hotspot residues.

Decision Making in Structure-Based Drug Discovery: Visual Inspection of Docking Results

Molecular docking is a computational method widely used in drug discovery. Due to the inherent inaccuracies of molecular docking, visual inspection of binding modes is a crucial routine in the decision making process of computational medicinal chemists. Despite its apparent importance for medicinal chemistry projects, guidelines for the visual docking pose assessment have been hardly discussed in the literature. Here, we review the medicinal chemistry literature with the aim of identifying consistent principles for visual inspection, highlighting cases of its successful application, and discussing its limitations. In this context, we conducted a survey reaching experts in both academia and the pharmaceutical industry, which also included a challenge to distinguish native from incorrect poses. We were able to collect 93 expert opinions that offer valuable insights into visually supported decision-making processes. This perspective shall motivate discussions among experienced computational medicinal chemists and guide young scientists new to the field to stratify their compounds.

Drug Repurposing and Polypharmacology to Fight SARS-CoV-2 Through Inhibition of the Main Protease

The outbreak of a new coronavirus (SARS-CoV-2), which is responsible for the COVID-19 disease and is spreading rapidly around the world, urgently requires effective therapeutic treatments. In this context, drug repurposing represents a valuable strategy, as it enables accelerating the identification of drug candidates with already known safety profiles, possibly aiding in the late stages of clinical evaluation. Moreover, therapeutic treatments based on drugs with beneficial multi-target activities (polypharmacology) may show an increased antiviral activity or help to counteract severe complications concurrently affecting COVID-19 patients. In this study, we present the results of a computational drug repurposing campaign that aimed at identifying potential inhibitors of the main protease (Mpro) of the SARS-CoV-2. The performed in silico screening allowed the identification of 22 candidates with putative SARS-CoV-2 Mpro inhibitory activity. Interestingly, some of the identified compounds have recently entered clinical trials for COVID-19 treatment, albeit not being assayed for their SARS-CoV-2 antiviral activity. Some candidates present a polypharmacology profile that may be beneficial for COVID-19 treatment and, to the best of our knowledge, have never been considered in clinical trials. For each repurposed compound, its therapeutic relevance and potential beneficial polypharmacological effects that may arise due to its original therapeutic indication are thoroughly discussed.

G protein-coupled receptors: structure- and function-based drug discovery

As one of the most successful therapeutic target families, G protein-coupled receptors (GPCRs) have experienced a transformation from random ligand screening to knowledge-driven drug design. We are eye-witnessing tremendous progresses made recently in the understanding of their structure-function relationships that facilitated drug development at an unprecedented pace. This article intends to provide a comprehensive overview of this important field to a broader readership that shares some common interests in drug discovery.

Potential Inhibitors of Fascin From A Database of Marine Natural Products: A Virtual Screening and Molecular Dynamics Study

Marine nature products are unique compounds that are produced by the marine environment including plants, animals, and microorganisms. The wide diversity of marine natural products have great potential and are versatile in terms of drug discovery. In this paper, we use state-of-the-art computational methods to discover inhibitors from marine natural products to block the function of Fascin, an overexpressed protein in various cancers. First, virtual screening (pharmacophore model and molecular docking) was carried out based on a marine natural products database (12015 molecules) and provided eighteen molecules that could potentially inhibit the function of Fascin. Next, molecular mechanics generalized Born surface area (MM/GBSA) calculations were conducted and indicated that four molecules have higher binding affinities than the inhibitor NP-G2-029, which was validated experimentally. ADMET analyses of pharmacokinetics demonstrated that one of the four molecules does not match the criterion. Finally, ligand Gaussian accelerated molecular dynamics (LiGaMD) simulations were carried out to validate the three inhibitors binding to Fascin stably. In addition, dynamic interactions between protein and ligands were analyzed systematically. Our study will accelerate the development of the cancer drugs targeting Fascin.

Demystifying the Molecular Basis of Pyrazoloquinolinones Recognition at the Extracellular α1+/β3- Interface of the GABAA Receptor by Molecular Modeling

GABA_A receptors are pentameric ligand-gated ion channels that serve as major inhibitory neurotransmitter receptors in the mammalian brain and the target of numerous clinically relevant drugs interacting with different ligand binding sites. Here, we report an in silico approach to investigate the binding of pyrazoloquinolinones (PQs) that mediate allosteric effects through the extracellular α+/β- interface of GABA_A receptors. First, we docked a potent prototype of PQs into the α1+/β3- site of a homology model of the human α1β3γ2 subtype of the GABA_A receptor. Next, for each docking pose, we computationally derived protein-ligand complexes for 18 PQ analogs with known experimental potency. Subsequently, binding energy was calculated for all complexes using the molecular mechanics-generalized Born surface area method. Finally, docking poses were quantitatively assessed in the light of experimental data to derive a binding hypothesis. Collectively, the results indicate that PQs at the α1+/β3- site likely exhibit a common binding mode that can be characterized by a hydrogen bond interaction with β3Q64 and hydrophobic interactions involving residues α1F99, β3Y62, β3M115, α1Y159, and α1Y209. Importantly, our results are in good agreement with the recently resolved cryo-Electron Microscopy structures of the human α1β3γ2 and α1β2γ2 subtypes of GABA_A receptors.

Targeting the Protein Tunnels of the Urease Accessory Complex: A Theoretical Investigation

Urease is a nickel-containing enzyme that is essential for the survival of several and often deadly pathogenic bacterial strains, including Helicobacter pylori. Notwithstanding several attempts, the development of direct urease inhibitors without side effects for the human host remains, to date, elusive. The recently solved X-ray structure of the HpUreDFG accessory complex involved in the activation of urease opens new perspectives for structure-based drug discovery. In particular, the quaternary assembly and the presence of internal tunnels for nickel translocation offer an intriguing possibility to target the HpUreDFG complex in the search of indirect urease inhibitors. In this work, we adopted a theoretical framework to investigate such a hypothesis. Specifically, we searched for putative binding sites located at the protein–protein interfaces on the HpUreDFG complex, and we challenged their druggability through structure-based virtual screening. We show that, by virtue of the presence of tunnels, some protein–protein interfaces on the HpUreDFG complex are intrinsically well suited for hosting small molecules, and, as such, they possess good potential for future drug design endeavors.

How Computational Chemistry and Drug Delivery Techniques Can Support the Development of New Anticancer Drugs

The early and late development of new anticancer drugs, small molecules or peptides can be slowed down by some issues such as poor selectivity for the target or poor ADME properties. Computer-aided drug design (CADD) and target drug delivery (TDD) techniques, although apparently far from each other, are two research fields that can give a significant contribution to overcome these problems. Their combination may provide mechanistic understanding resulting in a synergy that makes possible the rational design of novel anticancer based therapies. Herein, we aim to discuss selected applications, some also from our research experience, in the fields of anticancer small organic drugs and peptides.

Integrating Ligand and Target-Driven Based Virtual Screening Approaches With in vitro Human Cell Line Models and Time-Resolved Fluorescence Resonance Energy Transfer Assay to Identify Novel Hit Compounds Against BCL-2

Antiapoptotic members of B-cell leukemia/lymphoma-2 (BCL-2) family proteins are one of the overexpressed proteins in cancer cells that are oncogenic targets. As such, targeting of BCL-2 family proteins raises hopes for new therapeutic discoveries. Thus, we used multistep screening and filtering approaches that combine structure and ligand-based drug design to identify new, effective BCL-2 inhibitors from a small molecule database (Specs SC), which includes more than 210,000 compounds. This database is first filtered based on binary “cancer-QSAR” model constructed with 886 training and 167 test set compounds and common 26 toxicity quantitative structure-activity relationships (QSAR) models. Predicted non-toxic compounds are considered for target-driven studies. Here, we applied two different approaches to filter and select hit compounds for further in vitro biological assays and human cell line experiments. In the first approach, a molecular docking and filtering approach is used to rank compounds based on their docking scores and only a few top-ranked molecules are selected for further long (100-ns) molecular dynamics (MD) simulations and in vitro tests. While docking algorithms are promising in predicting binding poses, they can be less prone to precisely predict ranking of compounds leading to decrease in the success rate of in silico studies. Hence, in the second approach, top-docking poses of each compound filtered through QSAR studies are subjected to initially short (1 ns) MD simulations and their binding energies are calculated via molecular mechanics generalized Born surface area (MM/GBSA) method. Then, the compounds are ranked based on their average MM/GBSA energy values to select hit molecules for further long MD simulations and in vitro studies. Additionally, we have applied text-mining approaches to identify molecules that contain “indol” phrase as many of the approved drugs contain indole and indol derivatives. Around 2700 compounds are filtered based on “cancer-QSAR” model and are then docked into BCL-2. Short MD simulations are performed for the top-docking poses for each compound in complex with BCL-2. The complexes are again ranked based on their MM/GBSA values to select hit molecules for further long MD simulations and in vitro studies. In total, seven molecules are subjected to biological activity tests in various human cancer cell lines as well as Time-Resolved Fluorescence Resonance Energy Transfer (TR-FRET) assay. Inhibitory concentrations are evaluated, and biological activities and apoptotic potentials are assessed by cell culture studies. Four molecules are found to be limiting the proliferation capacity of cancer cells while increasing the apoptotic cell fractions.

Molecular Docking: Shifting Paradigms in Drug Discovery

Molecular docking is an established in silico structure-based method widely used in drug discovery. Docking enables the identification of novel compounds of therapeutic interest, predicting ligand-target interactions at a molecular level, or delineating structure-activity relationships (SAR), without knowing a priori the chemical structure of other target modulators. Although it was originally developed to help understanding the mechanisms of molecular recognition between small and large molecules, uses and applications of docking in drug discovery have heavily changed over the last years. In this review, we describe how molecular docking was firstly applied to assist in drug discovery tasks. Then, we illustrate newer and emergent uses and applications of docking, including prediction of adverse effects, polypharmacology, drug repurposing, and target fishing and profiling, discussing also future applications and further potential of this technique when combined with emergent techniques, such as artificial intelligence.

In Silico Repositioning of Cannabigerol as a Novel Inhibitor of the Enoyl Acyl Carrier Protein (ACP) Reductase (InhA)

Cannabigerol (CBG) and cannabichromene (CBC) are non-psychoactive cannabinoids that have raised increasing interest in recent years. These compounds exhibit good tolerability and low toxicity, representing promising candidates for drug repositioning. To identify novel potential therapeutic targets for CBG and CBC, an integrated ligand-based and structure-based study was performed. The results of the analysis led to the identification of CBG as a low micromolar inhibitor of the Enoyl acyl carrier protein (ACP) reductase (InhA) enzyme.