Structure-Based Virtual Screening: From Classical to Artificial Intelligence

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.

Inverse Mixed-Solvent Molecular Dynamics for Visualization of the Residue Interaction Profile of Molecular Probes

To ensure efficiency in discovery and development, the application of computational technology is essential. Although virtual screening techniques are widely applied in the early stages of drug discovery research, the computational methods used in lead optimization to improve activity and reduce the toxicity of compounds are still evolving. In this study, we propose a method to construct the residue interaction profile of the chemical structure used in the lead optimization by performing “inverse” mixed-solvent molecular dynamics (MSMD) simulation. Contrary to constructing a protein-based, atom interaction profile, we constructed a probe-based, protein residue interaction profile using MSMD trajectories. It provides us the profile of the preferred protein environments of probes without co-crystallized structures. We assessed the method using three probes: benzamidine, catechol, and benzene. As a result, the residue interaction profile of each probe obtained by MSMD was a reasonable physicochemical description of the general non-covalent interaction. Moreover, comparison with the X-ray structure containing each probe as a ligand shows that the map of the interaction profile matches the arrangement of amino acid residues in the X-ray structure.

Development of New Inhibitors of HDAC1-3 Enzymes Aided by In Silico Design Strategies

Histone deacetylases (HDACs) are overexpressed in cancer, and their inhibition shows promising results in cancer therapy. In particular, selective class I HDAC inhibitors such as entinostat are proposed to be more beneficial in breast cancer treatment. Computational drug design is an inevitable part of today's drug discovery projects because of its unequivocal role in saving time and cost. Using three HDAC inhibitors trichostatin, vorinostat, and entinostat as template structures and a diverse fragment library, all synthetically accessible compounds thereof (∼3200) were generated virtually and filtered based on similarity against the templates and PAINS removal. The 298 selected structures were docked into the active site of HDAC I and ranked using a calculated binding affinity. Top-ranking structures were inspected manually, and, considering the ease of synthesis and drug-likeness, two new structures (3a and 3b) were proposed for synthesis and biological evaluation. The synthesized compounds were purified to a degree of more than 95% and structurally verified using various methods. The designed compounds 3a and 3b showed 65-80 and 5% inhibition on HDAC 1, 2, and 3 isoforms at a concentration of 10 μM, respectively. The novel compound 3a may be used as a lead structure for designing new HDAC inhibitors.

Target-Based Virtual Screening of Natural Compounds Identifies a Potent Antimalarial With Selective Falcipain-2 Inhibitory Activity

We employed a comprehensive approach of target-based virtual high-throughput screening to find potential hits from the ZINC database of natural compounds against cysteine proteases falcipain-2 and falcipain-3 (FP2 and FP3). Molecular docking studies showed the initial hits showing high binding affinity and specificity toward FP2 were selected. Furthermore, the enzyme inhibition and surface plasmon resonance assays were performed which resulted in a compound ZINC12900664 (ST72) with potent inhibitory effects on purified FP2. ST72 exhibited strong growth inhibition of chloroquine-sensitive (3D7; EC₅₀ = 2.8 µM) and chloroquine-resistant (RKL-9; EC₅₀ = 6.7 µM) strains of Plasmodium falciparum. Stage-specific inhibition assays revealed a delayed and growth defect during parasite growth and development in parasites treated with ST72. Furthermore, ST72 significantly reduced parasite load and increased host survival in a murine model infected with Plasmodium berghei ANKA. No Evans blue staining in ST72 treatment indicated that ST72 mediated protection of blood–brain barrier integrity in mice infected with P. berghei. ST72 did not show any significant hemolysis or cytotoxicity against human HepG2 cells suggesting a good safety profile. Importantly, ST72 with CQ resulted in improved growth inhibitory activity than individual drugs in both in vitro and in vivo studies.

In Silico Screening of Cinnamon (Cinnamomum burmannii) Bioactive Compounds as Acetylcholinesterase Inhibitors

Alzheimer’s is a progressive and neurodegenerative disease that mainly affects people aged 65 years and older. The pathophysiology of Alzheimer’s is possibly related to the depletion of the neurotransmitter acetylcholine (ACh) due to beta-amyloid plaques and neurofibrillary tangles. Secondary metabolites found in cinnamon bark (Cinnamomum burmannii) have the potential as anticholinesterases to treat Alzheimer’s symptoms. This study aimed to identify the potency of bioactive compounds from cinnamon bark as AChE inhibitors in silico through analysis of binding energy, inhibition constants, and types of interactions. The research was conducted by screening virtually 60 test ligands using the PyRx program and molecular docking using the Autodock Tools program. The results of the ligand-receptor interaction analysis showed that 12 of the 15 tested ligands had potential as AChE inhibitors. Epicatechin and medioresinol are the ligands with the best potential for AChE inhibition with affinity close to the natural ligand or donepezil. Epicatechin has a binding energy of −10.0 kcal/mol and inhibition constant of 0.0459 M, with four hydrogen bonds and seven hydrophobic bonds. Meanwhile, medioresinol has −9.9 kcal/mol binding energy and inhibition constant of 0.0543 M, with one hydrogen bond and thirteen hydrophobic bonds.

PIGNet: a physics-informed deep learning model toward generalized drug-target interaction predictions

Recently, deep neural network (DNN)-based drug–target interaction (DTI) models were highlighted for their high accuracy with affordable computational costs. Yet, the models' insufficient generalization remains a challenging problem in the practice of in silico drug discovery. We propose two key strategies to enhance generalization in the DTI model. The first is to predict the atom–atom pairwise interactions via physics-informed equations parameterized with neural networks and provides the total binding affinity of a protein–ligand complex as their sum. We further improved the model generalization by augmenting a broader range of binding poses and ligands to training data. We validated our model, PIGNet, in the comparative assessment of scoring functions (CASF) 2016, demonstrating the outperforming docking and screening powers than previous methods. Our physics-informing strategy also enables the interpretation of predicted affinities by visualizing the contribution of ligand substructures, providing insights for further ligand optimization.

PIGNet, a deep neural network-based drug–target interaction model guided by physics and extensive data augmentation, shows significantly improved generalization ability and model performance.

A chemical probe inhibitor targeting STAT1 restricts cancer stem cell traits and angiogenesis in colorectal cancer

Background

Colorectal cancer (CRC) is a worldwide cancer with rising annual incidence. New medications for patients with CRC are still needed. Recently, fluorescent chemical probes have been developed for cancer imaging and therapy. Signal transducer and activator of transcription 1 (STAT1) has complex functions in tumorigenesis and its role in CRC still needs further investigation.

Methods

RNA sequencing datasets in the NCBI GEO repository were analyzed to investigate the expression of STAT1 in patients with CRC. Xenograft mouse models, tail vein injection mouse models, and azoxymethane/dextran sodium sulfate (AOM/DSS) mouse models were generated to study the roles of STAT1 in CRC. A ligand-based high-throughput virtual screening approach combined with SWEETLEAD chemical database analysis was used to discover new STAT1 inhibitors. A newly designed and synthesized fluorescently labeled 4’,5,7-trihydroxyisoflavone (THIF) probe (BODIPY-THIF) elucidated the mechanistic actions of STAT1 and THIF in vitro and in vivo. Colonosphere formation assay and chick chorioallantoic membrane assay were used to evaluate stemness and angiogenesis, respectively.

Results

Upregulation of STAT1 was observed in patients with CRC and in mouse models of AOM/DSS-induced CRC and metastatic CRC. Knockout of STAT1 in CRC cells reduced tumor growth in vivo. We then combined a high-throughput virtual screening approach and analysis of the SWEETLEAD chemical database and found that THIF, a flavonoid abundant in soybeans, was a novel STAT1 inhibitor. THIF inhibited STAT1 phosphorylation and might bind to the STAT1 SH2 domain, leading to blockade of STAT1-STAT1 dimerization. The results of in vitro and in vivo binding studies of THIF and STAT1 were validated. The pharmacological treatment with BODIPY-THIF or ablation of STAT1 via a CRISPR/Cas9-based strategy abolished stemness and angiogenesis in CRC. Oral administration of BODIPY-THIF attenuated colitis symptoms and tumor growth in the mouse model of AOM/DSS-induced CRC.

Conclusions

This study demonstrates that STAT1 plays an oncogenic role in CRC. BODIPY-THIF is a new chemical probe inhibitor of STAT1 that reduces stemness and angiogenesis in CRC. BODIPY-THIF can be a potential tool for CRC therapy as well as cancer cell imaging.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12929-022-00803-4.

Structural insights into inhibition of PRRSV Nsp4 revealed by structure-based virtual screening, molecular dynamics, and MM-PBSA studies

BACKGROUND

Porcine reproductive and respiratory syndrome respiratory sickness in weaned and growing pigs, as well as sow reproductive failure, and its infection is regarded as one of the most serious swine illnesses worldwide. Given the current lack of an effective treatment, in this study, we identified natural compounds capable of inhibiting non-structural protein 4 (Nsp4) of the virus, which is involved in their replication and pathogenesis.

RESULTS

We screened natural compounds (n = 97,999) obtained from the ZINC database against Nsp4 and selected the top 10 compounds for analysing protein-ligand interactions and physicochemical properties. The five compounds demonstrating strong binding affinity were then subjected to molecular dynamics simulations (100 ns) and binding free energy calculations. Based on analysis, we identified four possible lead compounds that represent potentially effective drug-like inhibitors.

CONCLUSIONS

These methods identified that these natural compounds are capable of inhibiting Nsp4 and possibly effective as antiviral therapeutics against PRRSV.

In Vitro and In Vivo Pharmacological Characterization of a Novel TRPM8 Inhibitor Chemotype Identified by Small-Scale Preclinical Screening

Transient receptor potential melastatin type 8 (TRPM8) is a target for the treatment of different physio-pathological processes. While TRPM8 antagonists are reported as potential drugs for pain, cancer, and inflammation, to date only a limited number of chemotypes have been investigated and thus a limited number of compounds have reached clinical trials. Hence there is high value in searching for new TRPM8 antagonistic to broaden clues to structure-activity relationships, improve pharmacological properties and explore underlying molecular mechanisms. To address this, the EDASA Scientific in-house molecular library has been screened in silico, leading to identifying twenty-one potentially antagonist compounds of TRPM8. Calcium fluorometric assays were used to validate the in-silico hypothesis and assess compound selectivity. Four compounds were identified as selective TRPM8 antagonists, of which two were dual-acting TRPM8/TRPV1 modulators. The most potent TRPM8 antagonists (BB 0322703 and BB 0322720) underwent molecular modelling studies to highlight key structural features responsible for drug–protein interaction. The two compounds were also investigated by patch-clamp assays, confirming low micromolar potencies. The most potent compound (BB 0322703, IC₅₀ 1.25 ± 0.26 μM) was then profiled in vivo in a cold allodinya model, showing pharmacological efficacy at 30 μM dose. The new chemotypes identified showed remarkable pharmacological properties paving the way to further investigations for drug discovery and pharmacological purposes.

Structural analysis, virtual screening and molecular simulation to identify potential inhibitors targeting 2'-O-ribose methyltransferase of SARS-CoV-2 coronavirus

SARS-CoV-2, an emerging coronavirus, has spread rapidly around the world, resulting in over ten million cases and more than half a million deaths as of July 1, 2020. Effective treatments and vaccines for SARS-CoV-2 infection do not currently exist. Previous studies demonstrated that nonstructural protein 16 (nsp16) of coronavirus is an S-adenosyl methionine (SAM)-dependent 2'-O-methyltransferase (2'-O-MTase) that has an important role in viral replication and prevents recognition by the host innate immune system. In the present study, we employed structural analysis, virtual screening, and molecular simulation approaches to identify clinically investigated and approved drugs which can act as promising inhibitors against nsp16 2'-O-MTase of SARS-CoV-2. Comparative analysis of primary amino acid sequences and crystal structures of seven human CoVs defined the key residues for nsp16 2-O'-MTase functions. Virtual screening and docking analysis ranked the potential inhibitors of nsp16 from more than 4,500 clinically investigated and approved drugs. Furthermore, molecular dynamics simulations were carried out on eight top candidates, including Hesperidin, Rimegepant, Gs-9667, and Sonedenoson, to calculate various structural parameters and understand the dynamic behavior of the drug-protein complexes. Our studies provided the foundation to further test and repurpose these candidate drugs experimentally and/or clinically for COVID-19 treatment.Communicated by Ramaswamy H. Sarma.

Virtual Combinatorial Chemistry and Pharmacological Screening: A Short Guide to Drug Design

Traditionally, drug development involved the individual synthesis and biological evaluation of hundreds to thousands of compounds with the intention of highlighting their biological activity, selectivity, and bioavailability, as well as their low toxicity. On average, this process of new drug development involved, in addition to high economic costs, a period of several years before hopefully finding a drug with suitable characteristics to drive its commercialization. Therefore, the chemical synthesis of new compounds became the limiting step in the process of searching for or optimizing leads for new drug development. This need for large chemical libraries led to the birth of high-throughput synthesis methods and combinatorial chemistry. Virtual combinatorial chemistry is based on the same principle as real chemistry—many different compounds can be generated from a few building blocks at once. The difference lies in its speed, as millions of compounds can be produced in a few seconds. On the other hand, many virtual screening methods, such as QSAR (Quantitative Sturcture-Activity Relationship), pharmacophore models, and molecular docking, have been developed to study these libraries. These models allow for the selection of molecules to be synthesized and tested with a high probability of success. The virtual combinatorial chemistry–virtual screening tandem has become a fundamental tool in the process of searching for and developing a drug, as it allows the process to be accelerated with extraordinary economic savings.

Structure-based discovery of a specific SHP2 inhibitor with enhanced blood-brain barrier penetration from PubChem database

The thiophene [2,3-d]pyrimidine structure-like small molecules were discovered from structure-based virtual screening of 1 billion compounds. Base on enzyme activity assay results, a SHP2-specific molecule inhibitor Comp#2 with IC₅₀ of 1.174 μM, 85-fold more selective for SHP2 than the highly related SHP1 (IC₅₀ > 100 μM). The compound can effectively inhibit SHP2-mediated cell signaling and cancer cell proliferation, including cervix cancer, human pancreatic cancer, large cell lung cancer, and mouse glioma cell. Moreover, the in vivo assay indicated that Comp#2 could inhibit cervix cancer tumors growth in BABL/c mice. This work has shown the specific SHP2 inhibitor can inhibit glioblastoma growth in vivo.

In Silico Drug Repurposing Approach: Investigation of Mycobacterium tuberculosis FadD32 Targeted by FDA-Approved Drugs

Background: Despite the enormous efforts made towards combating tuberculosis (TB), the disease remains a major global threat. Hence, new drugs with novel mechanisms against TB are urgently needed. Fatty acid degradation protein D32 (FadD32) has been identified as a promising drug target against TB, the protein is required for the biosynthesis of mycolic acids, hence, essential for the growth and multiplication of the mycobacterium. However, the FadD32 mechanism upon the binding of FDA-approved drugs is not well established. Herein, we applied virtual screening (VS), molecular docking, and molecular dynamic (MD) simulation to identify potential FDA-approved drugs against FadD32. Methodology/Results: VS technique was found promising to identify four FDA-approved drugs (accolate, sorafenib, mefloquine, and loperamide) with higher molecular docking scores, ranging from -8.0 to -10.0 kcal/mol. Post-MD analysis showed that the accolate hit displayed the highest total binding energy of -45.13 kcal/mol. Results also showed that the accolate hit formed more interactions with FadD32 active site residues and all active site residues displayed an increase in total binding contribution. RMSD, RMSF, Rg, and DCCM analysis further supported that the presence of accolate exhibited more structural stability, lower bimolecular flexibility, and more compactness into the FadD32 protein. Conclusions: Our study revealed accolate as the best potential drug against FadD32, hence a prospective anti-TB drug in TB therapy. In addition, we believe that the approach presented in the current study will serve as a cornerstone to identifying new potential inhibitors against a wide range of biological targets.

Graph Neural Networks as a Potential Tool in Improving Virtual Screening Programs

Despite the increasing number of pharmaceutical companies, university laboratories and funding, less than one percent of initially researched drugs enter the commercial market. In this context, virtual screening (VS) has gained much attention due to several advantages, including timesaving, reduced reagent and consumable costs and the performance of selective analyses regarding the affinity between test molecules and pharmacological targets. Currently, VS is based mainly on algorithms that apply physical and chemistry principles and quantum mechanics to estimate molecule affinities and conformations, among others. Nevertheless, VS has not reached the expected results concerning the improvement of market-approved drugs, comprising less than twenty drugs that have reached this goal to date. In this context, graph neural networks (GNN), a recent deep-learning subtype, may comprise a powerful tool to improve VS results concerning natural products that may be used both simultaneously with standard algorithms or isolated. This review discusses the pros and cons of GNN applied to VS and the future perspectives of this learnable algorithm, which may revolutionize drug discovery if certain obstacles concerning spatial coordinates and adequate datasets, among others, can be overcome.

Virtual Screening Inhibitors of Ubiquitin-specific Protease 7 combining Pharmacophore Modeling and Molecular Docking

Ubiquitin-specific protease 7 (USP7) is one of the most extensively studied deubiquitinases. USP7 exhibits a high expression signature in various malignant tumors, suggesting that it is a marker of tumor prognosis and a potential drug target for anti-tumor therapy. In this study, virtual screening based on pharmacophore model and biological evaluation have been applied for the discovery of novel USP7 inhibitors targeting the catalytic active site. The TS-4 was screened from 215,480 small molecules and was found to have USP7 inhibitory activity. Preliminary in vitro studies disclosed its antiproliferative activity on human colon cancer cell lines (HCT-116 and RKO), compared with normal colon cell line (CCD841CoN). Molecular dynamics (MD) simulation revealed the combine mechanism between USP7 with the TS-4. The TS-4 formed stable interactions with Asp295, Phe409 and Tyr514, which were critical to enhance its biological activity. This compound will serve as a promising hit compound for facilitating the further design of novel USP7 inhibitors.

HyperAttentionDTI: improving drug-protein interaction prediction by sequence-based deep learning with attention mechanism

MOTIVATION

Identifying drug-target interactions (DTIs) is a crucial step in drug repurposing and drug discovery. Accurately identifying DTIs in silico can significantly shorten development time and reduce costs. Recently, many sequence-based methods are proposed for DTI prediction and improve performance by introducing the attention mechanism. However, these methods only model single non-covalent inter-molecular interactions among drugs and proteins and ignore the complex interaction between atoms and amino acids.

RESULTS

In this article, we propose an end-to-end bio-inspired model based on the convolutional neural network (CNN) and attention mechanism, named HyperAttentionDTI, for predicting DTIs. We use deep CNNs to learn the feature matrices of drugs and proteins. To model complex non-covalent inter-molecular interactions among atoms and amino acids, we utilize the attention mechanism on the feature matrices and assign an attention vector to each atom or amino acid. We evaluate HpyerAttentionDTI on three benchmark datasets and the results show that our model achieves significantly improved performance compared with the state-of-the-art baselines. Moreover, a case study on the human Gamma-aminobutyric acid receptors confirm that our model can be used as a powerful tool to predict DTIs.

AVAILABILITY AND IMPLEMENTATION

The codes of our model are available at https://github.com/zhaoqichang/HpyerAttentionDTI and https://zenodo.org/record/5039589.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Fundamental considerations in drug design

The drug discovery paradigm has been very time-consuming, challenging, and expensive; however, the disease conditions originating from bacteria, virus, protozoa, fungus and other microorganisms are steadily shooting up. For instance, COVID-19 is the latest viral infection that affects millions of people and the world’s economy very severely. Therefore, the quest for discovery of novel and potent drug compounds against deadly pathogens is crucial at the moment. Despite a lot of drawbacks in drug discovery and development and its pertaining technology, the advancement must be taken into account so the time duration and cost would be minimized. In this chapter, basic principles in drug design and discovery have been discussed together with advances in drug development.

Present and future challenges in therapeutic designing using computational approaches

Currently, various computational methods are being used for the purpose of therapeutic design. The advent of the Coronavirus disease-2019 (COVID-19) pandemic has created a lot of problems due to which the development of effective treatment options is urgently needed. Computational intelligence is used in the control, prevention, prediction, diagnosis, and treatment of the disease. Several important drug targets have been identified in severe acute respiratory syndrome-Coronavirus-2 using in silico methods. Computer-aided drug design includes a variety of theoretical and computational approaches that are part of modern drug discovery. Advances in machine learning methods and their applications speed up the drug discovery process. Exploration of nucleic acid-based therapeutics is playing an important role in healthcare also. But a lot of challenges have also been seen that complicate the therapeutic design. Therefore, investigation of challenges associated with therapeutic design is important, and the present chapter is aimed to cover various therapeutic design approaches and challenges associated with them. Moreover, the role of computational strategies in the exploration of potential therapeutics against COVID-19 has been investigated.

Parallel tempered genetic algorithm guided by deep neural networks for inverse molecular design

We present JANUS, an evolutionary algorithm for inverse molecular design. It propagates an explorative and an exploitative population exchanging members via parallel tempering and uses active learning via deep neural networks to enhance sampling.

MoleGuLAR: Molecule Generation Using Reinforcement Learning with Alternating Rewards

The design of new inhibitors for novel targets is a very important problem especially in the current scenario with the world being plagued by COVID-19. Conventional approaches such as high-throughput virtual screening require extensive combing through existing data sets in the hope of finding possible matches. In this study, we propose a computational strategy for de novo generation of molecules with high binding affinities to the specified target and other desirable properties for druglike molecules using reinforcement learning. A deep generative model built using a stack-augmented recurrent neural network initially trained to generate druglike molecules is optimized using reinforcement learning to start generating molecules with desirable properties like LogP, quantitative estimate of drug likeliness, topological polar surface area, and hydration free energy along with the binding affinity. For multiobjective optimization, we have devised a novel strategy in which the property being used to calculate the reward is changed periodically. In comparison to the conventional approach of taking a weighted sum of all rewards, this strategy shows an enhanced ability to generate a significantly higher number of molecules with desirable properties.

Overcoming Depression with 5-HT2A Receptor Ligands

Depression is a multifactorial disorder that affects millions of people worldwide, and none of the currently available therapeutics can completely cure it. Thus, there is a need for developing novel, potent, and safer agents. Recent medicinal chemistry findings on the structure and function of the serotonin 2A (5-HT2A) receptor facilitated design and discovery of novel compounds with antidepressant action. Eligible papers highlighting the importance of 5-HT2A receptors in the pathomechanism of the disorder were identified in the content-screening performed on the popular databases (PubMed, Google Scholar). Articles were critically assessed based on their titles and abstracts. The most accurate papers were chosen to be read and presented in the manuscript. The review summarizes current knowledge on the applicability of 5-HT2A receptor signaling modulators in the treatment of depression. It provides an insight into the structural and physiological features of this receptor. Moreover, it presents an overview of recently conducted virtual screening campaigns aiming to identify novel, potent 5-HT2A receptor ligands and additional data on currently synthesized ligands acting through this protein.

Update and New Insights on Future Cancer Drug Candidates From Plant-Based Alkaloids

Cancer is a complex multifactorial disease that results from alterations in many physiological and biochemical functions. Over the last few decades, it has become clear that cancer cells can acquire multidrug resistance to conventional anticancer drugs, resulting in tumor relapse. Thus, there is a continuous need to discover new and effective anticancer drugs. Natural products from plants have served as a primary source of cancer drugs and continue to provide new plant-derived anticancer drugs. The present review describes plant-based alkaloids, which have been reported as active or potentially active in cancer treatment within the past 4 years (2017-2020), both in preclinical research and/or in clinical trials. In addition, recent insights into the possible molecular mechanism of action of alkaloid prodrugs naturally present in plants are also highlighted.

Identification of Phytoconstituents as Potent Inhibitors of Casein Kinase-1 Alpha Using Virtual Screening and Molecular Dynamics Simulations

Casein kinase-1 alpha (CK1α) is a multifunctional protein kinase that belongs to the serine/threonine kinases of the CK1α family. It is involved in various signaling pathways associated with chromosome segregation, cell metabolism, cell cycle progression, apoptosis, autophagy, etc. It has been known to involve in the progression of many diseases, including cancer, neurodegeneration, obesity, and behavioral disorders. The elevated expression of CK1α in diseased conditions facilitates its selective targeting for therapeutic management. Here, we have performed virtual screening of phytoconstituents from the IMPPAT database seeking potential inhibitors of CK1α. First, a cluster of compounds was retrieved based on physicochemical parameters following Lipinski's rules and PAINS filter. Further, high-affinity hits against CK1α were obtained based on their binding affinity score. Furthermore, the ADMET, PAINS, and PASS evaluation was carried out to select more potent hits. Finally, following the interaction analysis, we elucidated three phytoconstituents, Semiglabrinol, Curcusone_A, and Liriodenine, posturing considerable affinity and specificity towards the CK1α binding pocket. The result was further evaluated by molecular dynamics (MD) simulations, dynamical cross-correlation matrix (DCCM), and principal components analysis (PCA), which revealed that binding of the selected compounds, especially Semiglabrinol, stabilizes CK1α and leads to fewer conformational fluctuations. The MM-PBSA analysis suggested an appreciable binding affinity of all three compounds toward CK1α.

Identification of Novel Potential Anti-Diabetic Candidates targeting Human Pancreatic ɑ-Amylase and Human ɑ-Glycosidase: An Exhaustive Structure-based Screening

Diabetes is a major health issue that half a billion people affected worldwide. It is a serious, long-term medical condition majorly impacting the lives and well-being of individuals, families, and societies at large. It is amongst the top 10 diseases responsible for the death amongst adults with an expected rise to 10.2% (578 million) by 2030 and 10.9% (700 million) by 2045. The carbohydrates absorbed into the body are hydrolyzed by pancreatic α-amylase and other enzymes, human α-glucosidase. The α-amylase and α-glucosidase are validated therapeutic targets in the treatment of Type II diabetes (T2DM) as they play a vital role in modulating the blood glucose post meal. Herein, we report novel and diverse molecules as potential candidates, with predicted affinity for α-amylase and α-glucosidase. These molecules have been identified via hierarchical multistep docking of small molecules database with the estimated binding free energies. A Glide XP Score cutoff −8.00 kcal/mol was implemented to filter out non potential molecules. Four molecules viz. amb22034702, amb18105639, amb17153304, and amb9760832 have been identified after an exhaustive computational study involving, evaluation of binding interactions and assessment of the pharmacokinetics and toxicity profiles. The in-depth analysis of protein– ligand interactions was performed using a 100ns molecular dynamics (MD) simulation to establish the dynamic stability. Furthermore MM-GBSA based binding free energies were computed for 1000 trajectory snapshots to ascertain the strong binding affinity of these molecules for α-amylase and αglucosidase. The identified molecules can be considered as promising candidates for further drug development through necessary experimental assessments.

Drug Discovery for Mycobacterium tuberculosis Using Structure-Based Computer-Aided Drug Design Approach

Developing new, more effective antibiotics against resistant Mycobacterium tuberculosis that inhibit its essential proteins is an appealing strategy for combating the global tuberculosis (TB) epidemic. Finding a compound that can target a particular cavity in a protein and interrupt its enzymatic activity is the crucial objective of drug design and discovery. Such a compound is then subjected to different tests, including clinical trials, to study its effectiveness against the pathogen in the host. In recent times, new techniques, which involve computational and analytical methods, enhanced the chances of drug development, as opposed to traditional drug design methods, which are laborious and time-consuming. The computational techniques in drug design have been improved with a new generation of software used to develop and optimize active compounds that can be used in future chemotherapeutic development to combat global tuberculosis resistance. This review provides an overview of the evolution of tuberculosis resistance, existing drug management, and the design of new anti-tuberculosis drugs developed based on the contributions of computational techniques. Also, we show an appraisal of available software and databases on computational drug design with an insight into the application of this software and databases in the development of anti-tubercular drugs. The review features a perspective involving machine learning, artificial intelligence, quantum computing, and CRISPR combination with available computational techniques as a prospective pathway to design new anti-tubercular drugs to combat resistant tuberculosis.

Bacterial Quorum-Sensing Molecules as Promising Natural Inhibitors of Candida albicans Virulence Dimorphism: An In Silico and In Vitro Study

Farnesol, a self-secreted quorum-sensing molecule (QSM) of Candida albicans, has been known to limit yeast-to-hyphal transition by blocking the RAS1-cAMP-PKA pathway. In a similar fashion, certain bacterial QSMs have also been reported to be successful in attenuating C. albicans biofilm and hyphal formation at relatively high cell density. This prompted us to investigate the antihyphal efficacy of certain bacterial QSMs through virtual docking against seminal drug targets, viz., CYCc and RAS1, that have been reported to be the hallmark players in C. albicans dimorphic virulence cascade. Against this backdrop, 64 QSMs belonging to five different bacterial QS signaling systems were subjected to initial virtual screening with farnesol as reference. Data of the virtual screening unveiled QSMs belonging to diketopiperazines (DKPs), i.e., 3-benzyl-6-isobutylidene-2,5-piperazinedione (QSSM 1157) and cyclo(l-Pro-l-Leu) (QSSM 1112), as potential inhibitors of CYCc and RAS1 with binding energies of -8.2 and -7.3 kcal mol-1, respectively. Further, the molecular dynamics simulations (for 50 ns) of CYCc-QSSM 1157 and RAS1-QSSM 1112 complexes revealed the mean ligand root mean square deviation (RMSD) values of 0.35 and 0.27 Å, respectively, which endorsed the rigid nature, less fluctuation in binding stiffness, and conformation of binding complexes. Furthermore, the identified two QSMs were found to be good in solubility, absorption, and permeation and less toxic in nature, as revealed by pharmacokinetics and toxicity analyses. In addition, the in vitro antihyphal assays using liquid and solid media, germ-tube experiment, and microscopic analysis strongly validated DKP-QSSM 1112 as a promising inhibitor of hyphal transition. Taken together, the present study unequivocally proves that DKPs can be used as potent inhibitors of C. albicans virulence dimorphism.

Applications of Artificial Intelligence in Drug Design: Opportunities and Challenges

Artificial intelligence (AI) has undergone rapid development in recent years and has been successfully applied to real-world problems such as drug design. In this chapter, we review recent applications of AI to problems in drug design including virtual screening, computer-aided synthesis planning, and de novo molecule generation, with a focus on the limitations of the application of AI therein and opportunities for improvement. Furthermore, we discuss the broader challenges imposed by AI in translating theoretical practice to real-world drug design; including quantifying prediction uncertainty and explaining model behavior.

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

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

Fine tuning for success in structure-based virtual screening

Structure-based virtual screening plays a significant role in drug-discovery. The method virtually docks millions of compounds from corporate or public libraries into a binding site of a disease-related protein structure, allowing for the selection of a small list of potential ligands for experimental testing. Many algorithms are available for docking and assessing the affinity of compounds for a targeted protein site. The performance of affinity estimation calculations is highly dependent on the size and nature of the site, therefore a rationale for selecting the best protocol is required. To address this issue, we have developed an automated calibration process, implemented in a Knime workflow. It consists of four steps: preparation of a protein test set with structures and models of the target, preparation of a compound test set with target-related ligands and decoys, automatic test of 24 scoring/rescoring protocols for each target structure and model, and graphical display of results. The automation of the process combined with execution on high performance computing resources greatly reduces the duration of the calibration phase, and the test of many combinations of algorithms on various target conformations results in a rational and optimal choice of the best protocol. Here, we present this tool and exemplify its application in setting-up an optimal protocol for SBVS against Retinoid X Receptor alpha.

In Silico End-to-End Protein-Ligand Interaction Characterization Pipeline: The Case of SARS-CoV-2

SARS-CoV-2 triggered a worldwide pandemic disease, COVID-19, for which an effective treatment has not yet been settled. Among the most promising targets to fight this disease is SARS-CoV-2 main protease (Mpro), which has been extensively studied in the last few months. There is an urgency for developing effective computational protocols that can help us tackle these key viral proteins. Hence, we have put together a robust and thorough pipeline of in silico protein-ligand characterization methods to address one of the biggest biological problems currently plaguing our world. These methodologies were used to characterize the interaction of SARS-CoV-2 Mpro with an α-ketoamide inhibitor and include details on how to upload, visualize, and manage the three-dimensional structure of the complex and acquire high-quality figures for scientific publications using PyMOL (Protocol 1); perform homology modeling with MODELLER (Protocol 2); perform protein-ligand docking calculations using HADDOCK (Protocol 3); run a virtual screening protocol of a small compound database of SARS-CoV-2 candidate inhibitors with AutoDock 4 and AutoDock Vina (Protocol 4); and, finally, sample the conformational space at the atomic level between SARS-CoV-2 Mpro and the α-ketoamide inhibitor with Molecular Dynamics simulations using GROMACS (Protocol 5). Guidelines for careful data analysis and interpretation are also provided for each Protocol.

Screening of β1- and β2-Adrenergic Receptor Modulators through Advanced Pharmacoinformatics and Machine Learning Approaches

Cardiovascular diseases (CDs) are a major concern in the human race and one of the leading causes of death worldwide. β-Adrenergic receptors (β1-AR and β2-AR) play a crucial role in the overall regulation of cardiac function. In the present study, structure-based virtual screening, machine learning (ML), and a ligand-based similarity search were conducted for the PubChem database against both β1- and β2-AR. Initially, all docked molecules were screened using the threshold binding energy value. Molecules with a better binding affinity were further used for segregation as active and inactive through ML. The pharmacokinetic assessment was carried out on molecules retained in the above step. Further, similarity searching of the ChEMBL and DrugBank databases was performed. From detailed analysis of the above data, four compounds for each of β1- and β2-AR were found to be promising in nature. A number of critical ligand-binding amino acids formed potential hydrogen bonds and hydrophobic interactions. Finally, a molecular dynamics (MD) simulation study of each molecule bound with the respective target was performed. A number of parameters obtained from the MD simulation trajectories were calculated and substantiated the stability between the protein-ligand complex. Hence, it can be postulated that the final molecules might be crucial for CDs subjected to experimental validation.

Effect of drug metabolism in the treatment of SARS-CoV-2 from an entirely computational perspective

Understanding the effects of metabolism on the rational design of novel and more effective drugs is still a considerable challenge. To the best of our knowledge, there are no entirely computational strategies that make it possible to predict these effects. From this perspective, the development of such methodologies could contribute to significantly reduce the side effects of medicines, leading to the emergence of more effective and safer drugs. Thereby, in this study, our strategy is based on simulating the electron ionization mass spectrometry (EI-MS) fragmentation of the drug molecules and combined with molecular docking and ADMET models in two different situations. In the first model, the drug is docked without considering the possible metabolic effects. In the second model, each of the intermediates from the EI-MS results is docked, and metabolism occurs before the drug accesses the biological target. As a proof of concept, in this work, we investigate the main antiviral drugs used in clinical research to treat COVID-19. As a result, our strategy made it possible to assess the biological activity and toxicity of all potential by-products. We believed that our findings provide new chemical insights that can benefit the rational development of novel drugs in the future.

Molecular dynamics: a powerful tool for studying the medicinal chemistry of ion channel modulators

Molecular dynamics (MD) simulations allow researchers to investigate the behavior of desired biological targets at ever-decreasing costs with ever-increasing precision. Among the biological macromolecules, ion channels are remarkable transmembrane proteins, capable of performing special biological processes and revealing a complex regulatory matrix, including modulation by small molecules, either endogenous or exogenous. Recently, given the developments in ion channel structure determination and accessibility of bio-computational techniques, MD and related tools are becoming increasingly popular in the intense research area regarding ligand-channel interactions. This review synthesizes and presents the most important fields of MD involvement in investigating channel-molecule interactions, including, but not limited to, deciphering the binding modes of ligands to their ion channel targets and the mechanisms through which chemical compounds exert their effect on channel function. Special attention is devoted to the importance of more elaborate methods, such as free energy calculations, while principles regarding drug design and discovery are highlighted. Several technical aspects involving the creation and simulation of channel-molecule MD systems (ligand parameterization, proper membrane setup, system building, etc.) are also presented.

Structure-based virtual screening, in silico docking, ADME properties prediction and molecular dynamics studies for the identification of potential inhibitors against SARS-CoV-2 Mpro

COVID-19 is a viral pandemic caused by SARS-CoV-2. Due to its highly contagious nature, millions of people are getting affected worldwide knocking down the delicate global socio-economic equilibrium. According to the World Health Organization, COVID-19 has affected over 186 million people with a mortality of around 4 million as of July 09, 2021. Currently, there are few therapeutic options available for COVID-19 control. The rapid mutations in SARS-CoV-2 genome and development of new virulent strains with increased infection and mortality among COVID-19 patients, there is a great need to discover more potential drugs for SARS-CoV-2 on a priority basis. One of the key viral enzymes responsible for the replication and maturation of SARS-CoV-2 is M^pro protein. In the current study, structure-based virtual screening was used to identify four potential ligands against SARS-CoV-2 M^pro from a set of 8,722 ASINEX library compounds. These four compounds were evaluated using ADME filter to check their ADME profile and druggability, and all the four compounds were found to be within the current pharmacological acceptable range. They were individually docked to SARS-CoV-2 M^pro protein to assess their molecular interactions. Further, molecular dynamics (MD) simulations was carried out on protein–ligand complex using Desmond at 100 ns to explore their binding conformational stability. Based on RMSD, RMSF and hydrogen bond interactions, it was found that the stability of protein–ligand complex was maintained throughout the entire 100 ns simulations for all the four compounds. Some of the key ligand amino acid residues participated in stabilizing the protein–ligand interactions includes GLN 189, SER 10, GLU 166, ASN 142 with PHE 66 and TRP 132 of SARS-CoV-2 M^pro. Further optimization of these compounds could lead to promising drug candidates for SARS-CoV-2 M^pro target.

Strategies and Approaches for Discovery of Small Molecule Disruptors of Biofilm Physiology

Biofilms, the predominant growth mode of microorganisms, pose a significant risk to human health. The protective biofilm matrix, typically composed of exopolysaccharides, proteins, nucleic acids, and lipids, combined with biofilm-grown bacteria’s heterogenous physiology, leads to enhanced fitness and tolerance to traditional methods for treatment. There is a need to identify biofilm inhibitors using diverse approaches and targeting different stages of biofilm formation. This review discusses discovery strategies that successfully identified a wide range of inhibitors and the processes used to characterize their inhibition mechanism and further improvement. Additionally, we examine the structure–activity relationship (SAR) for some of these inhibitors to optimize inhibitor activity.

Structure-based virtual screening suggests inhibitors of 3-Chymotrypsin-Like Protease of SARS-CoV-2 from Vernonia amygdalina and Occinum gratissimum

Antiviral culinary plants are potential bioresources for preventive nutraceuticals and/or antiviral drugs in COVID-19. Structure-based virtual screening was undertaken to screen 173 compounds previously reported from Vernonia amygdalina and Occinum gratissimum for direct interaction with the active site of the 3-Chymotrypsin-Like Protease (3CL^pro) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Based on docking scores and comparison with reference inhibitors, a hit-list of 10 top phytocompounds was defined, which also had strong interactions with the catalytic centre of 3CL^pro from three related strains of coronavirus (SARS-CoV, MERS-CoV, HKU4). Among these, six compounds (neoandrographolide, vernolide, isorhamnetin, chicoric acid, luteolin, and myricetin) exhibited the highest binding tendencies to the equilibrated conformers of SARS-CoV-2 3CL^pro in an in-depth docking analysis to 5 different representative conformations from the cluster analysis of the molecular dynamics simulation (MDS) trajectories of the protein. In silico drug-likeness analyses revealed two drug-like terpenoids viz: neoandrographolide and vernolide as promising inhibitors of SARS-CoV-2 3CL^pro. These structures were accommodated within the substrate-binding pocket; and interacted with the catalytic dyad (Cys¹⁴⁵ and His⁴¹), the oxyanion loop (residues 138-145), and the S1/S2 sub-sites of the enzyme active site through the formation of an array of hydrogen bonds and hydrophobic interactions. Molecular dynamics simulation and binding free energy calculation revealed that the terpenoid-enzyme complexes exhibit strong interactions and structural stability. Therefore, these compounds may stabilize the conformation of the flexible oxyanion loop; and thereby interfere with the tetrahedral oxyanion intermediate formation during the proteolytic activity of the enzyme.

Current trends in computer aided drug design and a highlight of drugs discovered via computational techniques: A review

Computer-aided drug design (CADD) is one of the pivotal approaches to contemporary pre-clinical drug discovery, and various computational techniques and software programs are typically used in combination, in a bid to achieve the desired outcome. Several approved drugs have been developed with the aid of CADD. On SciFinder®, we evaluated more than 600 publications through systematic searching and refining, using the terms, virtual screening; software methods; computational studies and publication year, in order to obtain data concerning particular aspects of CADD. The primary focus of this review was on the databases screened, virtual screening and/or molecular docking software program used. Furthermore, we evaluated the studies that subsequently performed molecular dynamics (MD) simulations and we reviewed the software programs applied, the application of density functional theory (DFT) calculations and experimental assays. To represent the latest trends, the most recent data obtained was between 2015 and 2020, consequently the most frequently employed techniques and software programs were recorded. Among these, the ZINC database was the most widely preferred with an average use of 31.2%. Structure-based virtual screening (SBVS) was the most prominently used type of virtual screening and it accounted for an average of 57.6%, with AutoDock being the preferred virtual screening/molecular docking program with 41.8% usage. Following the screening process, 38.5% of the studies performed MD simulations to complement the virtual screening and GROMACS with 39.3% usage, was the popular MD software program. Among the computational techniques, DFT was the least applied whereby it only accounts for 0.02% average use. An average of 36.5% of the studies included reports on experimental evaluations following virtual screening. Ultimately, since the inception and application of CADD in pre-clinical drug discovery, more than 70 approved drugs have been discovered, and this number is steadily increasing over time.

Strategies for discovering resistance-breaking, safe and sustainable commercial herbicides with novel modes of action and chemotypes

Farmers need to manage weeds to grow and harvest crops that are essential to our food and energy supply, and herbicides are the most important tool in the farmers' armory to combat weeds. There is now a crisis in agriculture that has been brought about by herbicides being rendered ineffective by resistant weeds or withdrawn from the market due to safety concerns. Efficacious herbicides with novel modes of action (MoAs) and chemotypes are urgently needed to control resistant weeds and satisfy public and regulators' stringent requirements for safe and sustainable products. This article explores the main strategies being deployed by academic and industrial institutions to discover the next generation of commercial herbicides: phenotypic and in vitro target based approaches. There are early signs that much needed innovation and herbicidal products with novel MoAs are on the horizon from start-ups and established agrochemical companies. © 2021 Society of Chemical Industry. © 2021 Society of Chemical Industry.

Modern computational intelligence based drug repurposing for diabetes epidemic

BACKGROUND AND AIM

Objectives are to explore recent advances in discovery of new antidiabetic agents using repurposing strategies and to discuss modern technologies used for drug repurposing highlighting diabetic specific web portal.

METHODS

Recent literature were studied and analyzed from various sources such as Scopus, PubMed, and IEEE Xplore databases.

RESULTS

Drugs like Niclosamideethanolamine, Methazolamide, Diacerein, Berberine, Clobetasol, etc. with possibility of repurposing to curb diabetes can be potential late-stage clinical candidates, providing access to information on pharmacology, formulation, and probable toxicity if any.

CONCLUSIONS

With collaboration of artificial intelligence (AI) with pharmacology, the efficiency of drug repurposing can improve significantly.

Plant-Based Phytochemical Screening by Targeting Main Protease of SARS-CoV-2 to Design Effective Potent Inhibitors

Currently, a worldwide pandemic has been declared in response to the spread of coronavirus disease 2019 (COVID-19), a fatal and fast-spreading viral infection caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The low availability of efficient vaccines and treatment options has resulted in a high mortality rate, bringing the world economy to its knees. Thus, mechanistic investigations of drugs capable of counteracting this disease are in high demand. The main protease (Mpro) expressed by SARS-CoV-2 has been targeted for the development of potential drug candidates due to the crucial role played by Mpro in viral replication and transcription. We generated a phytochemical library containing 1672 phytochemicals derived from 56 plants, which have been reported as having antiviral, antibacterial, and antifungal activity. A molecular docking program was used to screen the top three candidate compounds: epicatechin-3-O-gallate, psi-taraxasterol, and catechin gallate, which had respective binding affinities of -8.4, -8.5, and -8.8 kcal/mol. Several active sites in the targeted protein, including Cys145, His41, Met49, Glu66, and Met165, were found to interact with the top three candidate compounds. The multiple simulation profile, root-mean-square deviation, root-mean-square fluctuation, radius of gyration, and solvent-accessible surface area values supported the inflexible nature of the docked protein-compound complexes. The toxicity and carcinogenicity profiles were assessed, which showed that epicatechin-3-O-gallate, psi-taraxasterol, and catechin gallate had favorable pharmacological properties with no adverse effects. These findings suggest that these compounds could be developed as part of an effective drug development pathway to treat COVID-19.

The Application and Development of Deep Learning in Radiotherapy: A Systematic Review

With the massive use of computers, the growth and explosion of data has greatly promoted the development of artificial intelligence (AI). The rise of deep learning (DL) algorithms, such as convolutional neural networks (CNN), has provided radiation oncologists with many promising tools that can simplify the complex radiotherapy process in the clinical work of radiation oncology, improve the accuracy and objectivity of diagnosis, and reduce the workload, thus enabling clinicians to spend more time on advanced decision-making tasks. As the development of DL gets closer to clinical practice, radiation oncologists will need to be more familiar with its principles to properly evaluate and use this powerful tool. In this paper, we explain the development and basic concepts of AI and discuss its application in radiation oncology based on different task categories of DL algorithms. This work clarifies the possibility of further development of DL in radiation oncology.

Combination of consensus and ensemble docking strategies for the discovery of human dihydroorotate dehydrogenase inhibitors

The inconsistencies in the performance of the virtual screening (VS) process, depending on the used software and structural conformation of the protein, is a challenging issue in the drug design and discovery field. Varying performance, especially in terms of early recognition of the potential hit compounds, negatively affects the whole process and leads to unnecessary waste of the time and resources. Appropriate application of the ensemble docking and consensus-scoring approaches can significantly increase reliability of the VS results. Dihydroorotate dehydrogenase (DHODH) is a key enzyme in the pyrimidine biosynthesis pathway. It is considered as a valuable therapeutic target in cancer, autoimmune and viral diseases. Based on the conducted benchmark study and analysis of the effect of different combinations of the applied methods and approaches, here we suggested a structure-based virtual screening (SBVS) workflow that can be used to increase the reliability of VS.

Schistosomiasis Drug Discovery in the Era of Automation and Artificial Intelligence

Schistosomiasis is a parasitic disease caused by trematode worms of the genus Schistosoma and affects over 200 million people worldwide. The control and treatment of this neglected tropical disease is based on a single drug, praziquantel, which raises concerns about the development of drug resistance. This, and the lack of efficacy of praziquantel against juvenile worms, highlights the urgency for new antischistosomal therapies. In this review we focus on innovative approaches to the identification of antischistosomal drug candidates, including the use of automated assays, fragment-based screening, computer-aided and artificial intelligence-based computational methods. We highlight the current developments that may contribute to optimizing research outputs and lead to more effective drugs for this highly prevalent disease, in a more cost-effective drug discovery endeavor.

Prospective Role of Peptide-Based Antiviral Therapy Against the Main Protease of SARS-CoV-2

The recently emerged coronavirus (SARS-CoV-2) has created a crisis in world health, and economic sectors as an effective treatment or vaccine candidates are still developing. Besides, negative results in clinical trials and effective cheap solution against this deadly virus have brought new challenges. The viral protein, the main protease from SARS-CoV-2, can be effectively targeted due to its viral replication and pathogenesis role. In this study, we have enlisted 88 peptides from the AVPdb database. The peptide molecules were modeled to carry out the docking interactions. The four peptides molecules, P14, P39, P41, and P74, had more binding energy than the rest of the peptides in multiple docking programs. Interestingly, the active points of the main protease from SARS-CoV-2, Cys145, Leu141, Ser139, Phe140, Leu167, and Gln189, showed nonbonded interaction with the peptide molecules. The molecular dynamics simulation study was carried out for 200 ns to find out the docked complex’s stability where their stability index was proved to be positive compared to the apo and control complex. Our computational works based on peptide molecules may aid the future development of therapeutic options against SARS-CoV-2.