Virtual screening and structure-based 3D pharmacophore approach to identify small-molecule inhibitors of SARS-CoV-2 Mpro

Exploring binding mechanisms of omicron spike protein with dolutegravir and etravirine by molecular dynamics simulation, principal component analysis, and free binding energy calculations

The COVID-19 pandemic was caused by the SARS-CoV-2 virus, frequent mutations occurred to the wild-type virus resulting in evolved new variants. The WHO classified the new variants as 'Variants of Concern'. SARS-CoV-2 omicron evolved as the dominating variant at the end of 2021. Dolutegravir and etravirine were identified as inhibitors of SARS-CoV-2 entry to host cells in Omicron variants. In this study, combined in silico methods such as molecular docking, molecular dynamics, Principal component analysis, binding-free energy calculations, and Per Residues calculations were applied to investigate the mechanism of the bindings of the two inhibitors. The molecular dynamics results revealed the stability of dolutegravir-spike and etravirine-spike complexes in a similar manner to apo-protein. Dolutegravir and etravirine formed H-bonds and salt bridges with Omicron spike protein. The 2,4-difluoro phenyl moiety of dolutegravir plays an important role in binding the ligand. The binding mode and interactions of the two compounds indicated that Arg403, Tyr449, Tyr453, Arg493, Ser496, Arg498, Thr500, Tyr501, Gln502 and His505 are the key residues. The Principal Component Analyses suggested that no significant conformational changes happened for the two complexes during the simulations. Binding-free energy calculations showed that van der Waals interactions were the most important interactions for ligands' binding. Per-residue free energy decomposition revealed residues Arg493, Arg498, and Tyr501 contributed to the binding of the ligands through H-bonds and salt bridges formation while His505 contributed to H-bonds and Pi-Pi stacking and Phe497 contributed to hydrophobic interactions between ligands and Omicron spike protein.Communicated by Ramaswamy H. Sarma.

Identification of mitoxantrone as a potent inhibitor of CDK7/Cyclin H via structure-based virtual screening and In-Vitro validation by ADP-Glo kinase assay

Cyclin-dependent kinases, CDK7 and CDK9 play critical roles in cancer by regulating transcriptional processes essential for cell proliferation and survival. Their dysregulation leads to aberrant gene expression, promoting oncogenic pathways and contributing to tumor growth and progression. This study aimed to identify a new chemotype for CDK7/9 inhibitors using a structure-based virtual screening approach. Our research led to the discovery of mitoxantrone as an inhibitor of CDK7/H and CDK9/T1 from a library of FDA-approved small molecule drugs. Mitoxantrone, a chemotherapy agent used to treat acute nonlymphocytic leukemia, works by disrupting DNA synthesis and repair, thus inhibiting cancer cell growth. The study found that mitoxantrone effectively inhibits both CDK7/H and CDK9/T1 with IC50 values of 0.675 µM and 5.15 µM, respectively, while showing no inhibition of CDK2/E1 (IC50 > 100 µM) in in-vitro ADP-Glo kinase assay. It binds to the ATP pocket of CDK7 and CDK9, forming crucial H-bonds with MET 94 and CYS 106, respectively. It achieves dock scores of - 12.93 and - 12.59 kcal/mol, and MMGBSA binding energies of - 82.87 and - 81.59 kcal/mol, respectively. Molecular dynamics simulations over 100 ns confirmed stable interactions with MET 94 and CYS 106 in the hinge region of CDK7 and CDK9. The active site sequence alignment helped to understand the differential activity of mitoxantrone for CDK7, 9 and 2 inhibitions. The findings of the paper reveal a novel mechanism of mitoxantrone action that may contribute to its anticancer efficacy.

Identifying and characterising promising small molecule inhibitors of kinesin spindle protein using ligand-based virtual screening, molecular docking, molecular dynamics and MM‑GBSA calculations

The kinesin spindle protein (Eg5) is a mitotic protein that plays an essential role in the formation of the bipolar spindles during the mitotic phase. Eg5 protein controls the segregation of the chromosomes in mitosis which renders it a vital target for cancer treatment. In this study our approach to identifying novel scaffold for Eg5 inhibitors is based on targeting the novel allosteric pocket (α4/α6/L11). Extensive computational techniques were applied using ligand-based virtual screening and molecular docking by two approaches, MOE and AutoDock, to screen a library of commercial compounds. We identified compound 8-(3-(1H-imidazol-1-ylpropylamino)-3-methyl-7-((naphthalen-3-yl)methyl)-1H-purine-2, 6 (3H,7H)-dione (compound 5) as a novel scaffold for Eg5 inhibitors. This compound inhibited cancer cell Eg5 ATPase at 2.37 ± 0.15 µM. The molecular dynamics simulations revealed that the identified compound formed stable interactions in the allosteric pocket (α4/α6/L11) of the receptor, indicating its potential as a novel Eg5 inhibitor.

Consensus Pharmacophore Strategy For Identifying Novel SARS-Cov-2 Mpro Inhibitors from Large Chemical Libraries

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main Protease (Mpro) is an enzyme that cleaves viral polyproteins translated from the viral genome and is critical for viral replication. Mpro is a target for anti-SARS-CoV-2 drug development, and multiple Mpro crystals complexed with competitive inhibitors have been reported. In this study, we aimed to develop an Mpro consensus pharmacophore as a tool to expand the search for inhibitors. We generated a consensus model by aligning and summarizing pharmacophoric points from 152 bioactive conformers of SARS-CoV-2 Mpro inhibitors. Validation against a library of conformers from a subset of ligands showed that our model retrieved poses that reproduced the crystal-binding mode in 77% of the cases. Using models derived from a consensus pharmacophore, we screened >340 million compounds. Pharmacophore-matching and chemoinformatics analyses identified new potential Mpro inhibitors. The candidate compounds were chemically dissimilar to the reference set, and among them, demonstrating the relevance of our model. We evaluated the effect of 16 candidates on Mpro enzymatic activity finding that seven have inhibitory activity. Three compounds (1, 4, and 5) had IC50 values in the midmicromolar range. The Mpro consensus pharmacophore reported herein can be used to identify compounds with improved activity and novel chemical scaffolds against Mpro. The method developed for its generation is provided as an open-access code (https://github.com/AngelRuizMoreno/ConcensusPharmacophore) and can be applied to other pharmacological targets.

Discovery of novel benzofuran-based derivatives as acetylcholinesterase inhibitors for the treatment of Alzheimer's disease: Design, synthesis, biological evaluation, molecular docking and 3D-QSAR investigation

A series of novel benzofuran-based compounds 7a-s were designed, synthesized, and investigated in vitro as acetylcholinesterase inhibitors (AChEIs). Compounds 7c and 7e displayed promising inhibitory activity with IC50 values of 0.058 and 0.086 μM in comparison to donepezil with an IC50 value of 0.049 μM. The new molecules' antioxidant evaluation revealed that 7c, 7e, 7j, 7n, and 7q produced the strongest DPPH scavenging activity when compared to vitamin C. As it was the most promising AChEI, compound 7c was selected for further biological evaluation. Acute and chronic toxicity studies exhibited that 7c showed no signs of toxicity or adverse events, no significant differences in the blood profile, and an insignificant difference in hepatic enzymes, glucose, urea, creatinine, and albumin levels in the experimental rat group. Furthermore, 7c did not produce histopathological damage to normal liver, kidney, heart, and brain tissues, ameliorated tissue malonaldehyde (MDA) and glutathione (GSH) levels and reduced the expression levels of the APP and Tau genes in AD rats. Molecular docking results of compounds 7c and 7e showed good binding modes in the active site of the acetylcholinesterase enzyme, which are similar to the native ligand donepezil. 3D-QSAR analysis revealed the importance of the alkyl group in positions 2 and 3 of the phenyl moiety for the activity. Overall, these findings suggested that compound 7c could be deemed a promising candidate for the management of Alzheimer's disease.

Development of a New Benzofuran-Pyrazole-Pyridine-Based Molecule for the Management of Osteoarthritis

Osteoarthritis is a substantial burden for patients with the disease. The known medications for the disease target the mitigation of the disease's symptoms. So, drug development for the management of osteoarthritis represents an important challenge in the medical field. This work is based on the development of a new benzofuran-pyrazole-pyridine-based compound 8 with potential anti-inflammatory and anti-osteoarthritis properties. Microanalytical and spectral data confirmed the chemical structure of compound 8. The biological assays indicated that compound 8 produces multifunctional activity as an anti-osteoarthritic candidate via inhibition of pro-inflammatory mediators, including RANTES, CRP, COMP, CK, and LPO in OA rats. Histopathological and pharmacokinetic studies confirmed the safety profile of the latter molecule. Accordingly, compound 8 is considered a promising anti-osteoarthritis agent and deserves deeper investigation in future trials.

Aurones: A Promising Scaffold to Inhibit SARS-CoV-2 Replication

Aurones are a small subgroup of flavonoids in which the basic C6-C3-C6 skeleton is arranged as (Z)-2-benzylidenebenzofuran-3(2H)-one. These compounds are structural isomers of flavones and flavonols, natural products reported as potent inhibitors of SARS-CoV-2 replication. Herein, we report the design, synthesis, and anti-SARS-CoV-2 activity of a series of 25 aurones bearing different oxygenated groups (OH, OCH3, OCH2OCH3, OCH2O, OCF2H, and OCH2C6H4R) at the A- and/or B-rings using cell-based screening assays. We observed that 12 of the 25 compounds exhibit EC50 < 3 μM (8e, 8h, 8j, 8k, 8l, 8m, 8p, 8q, 8r, 8w, 8x, and 8y), of which five presented EC50 < 1 μM (8h, 8m, 8p, 8q, and 8w) without evident cytotoxic effect in Calu-3 cells. The substitution of the A- and/or B-ring with OCH3, OCH2OCH3, and OCF2H groups seems beneficial for the antiviral activity, while the corresponding phenolic derivatives showed a significant decrease in the anti-SARS-CoV-2 activity. The most potent compound of the series, aurone 8q (EC50 = 0.4 μM, SI = 2441.3), is 2 to 3 times more effective than the polyphenolic flavonoids myricetin (2) and baicalein (1), respectively. Investigation of the five more active compounds as inhibitors of SARS-CoV-2 3CLpro based on molecular dynamic calculations suggested that these aurones should detach from the active site of 3CLpro, and, probably, they could bind to another SARS-CoV-2 protein target (either receptor or enzyme).

Anti-Alzheimer activity of new coumarin-based derivatives targeting acetylcholinesterase inhibition

New 2-oxo-chromene-7-oxymethylene acetohydrazide derivatives 4a-d were designed and synthesized with a variety of bioactive chemical fragments. The newly synthesized compounds were evaluated as acetylcholinesterase (AChE) inhibitors and antioxidant agents in comparison to donepezil and ascorbic acid, respectively. Compound 4c exhibited a promising inhibitory impact with an IC₅₀ value of 0.802 μM and DPPH scavenging activity of 57.14 ± 2.77%. Furthermore, biochemical and haematological studies revealed that compound 4c had no effect on the blood profile, hepatic enzyme levels (AST, ALT, and ALP), or total urea in 4c-treated rats compared to the controls. Moreover, the histopathological studies of 4c-treated rats revealed the normal architecture of the hepatic lobules and renal parenchyma, as well as no histopathological damage in the examined hepatic, kidney, heart, and brain tissues. In addition, an in vivo study investigated the amelioration in the cognitive function of AD-rats treated with 4c through the T-maze and beam balance behavioural tests. Also, 4c detectably ameliorated MDA and GSH, reaching 90.64 and 27.17%, respectively, in comparison to the standard drug (90.64% and 35.03% for MDA and GSH, respectively). The molecular docking study exhibited a good fitting of compound 4c in the active site of the AChE enzyme and a promising safety profile. Compound 4c exhibited a promising anti-Alzheimer's disease efficiency compared to the standard drug donepezil.

Pan-Genomics of Escherichia albertii for Antibiotic Resistance Profiling in Different Genome Fractions and Natural Product Mediated Intervention: In Silico Approach

Escherichia albertii is an emerging, enteric pathogen of significance. It was first isolated in 2003 from a pediatric diarrheal sample from Bangladesh. In this study, a comprehensive in silico strategy was followed to first list out antibiotic-resistant genes from core, accessory and unique genome fractions of 95 available genomes of E. albertii. Then, 56 drug targets were identified from the core essential genome. Finally, ZipA, an essential cell division protein that stabilizes the FtsZ protofilaments by cross-linking them and serves as a cytoplasmic membrane anchor for the Z ring, was selected for further downstream processing. It was computationally modeled using a threading approach, followed by virtual screening of two phytochemical libraries, Ayurvedic (n = 2103 compounds) and Traditional Chinese Medicine (n = 36,043 compounds). ADMET profiling, followed by PBPK modeling in the central body compartment, in a population of 250 non-diseased, 250 cirrhotic and 250 renally impaired people was attempted. ZINC85624912 from Chinese medicinal library showed the highest bioavailability and plasma retention. This is the first attempt to simulate the fate of natural products in the body through PBPK. Dynamics simulation of 20 ns for the top three compounds from both libraries was also performed to validate the stability of the compounds. The obtained information from the current study could aid wet-lab scientists to work on the scaffold of screened drug-like compounds from natural resources and could be useful in our quest for therapy against antibiotic-resistant E. albertii.

Small molecules in the treatment of COVID-19

The outbreak of COVID-19 has become a global crisis, and brought severe disruptions to societies and economies. Until now, effective therapeutics against COVID-19 are in high demand. Along with our improved understanding of the structure, function, and pathogenic process of SARS-CoV-2, many small molecules with potential anti-COVID-19 effects have been developed. So far, several antiviral strategies were explored. Besides directly inhibition of viral proteins such as RdRp and Mpro, interference of host enzymes including ACE2 and proteases, and blocking relevant immunoregulatory pathways represented by JAK/STAT, BTK, NF-κB, and NLRP3 pathways, are regarded feasible in drug development. The development of small molecules to treat COVID-19 has been achieved by several strategies, including computer-aided lead compound design and screening, natural product discovery, drug repurposing, and combination therapy. Several small molecules representative by remdesivir and paxlovid have been proved or authorized emergency use in many countries. And many candidates have entered clinical-trial stage. Nevertheless, due to the epidemiological features and variability issues of SARS-CoV-2, it is necessary to continue exploring novel strategies against COVID-19. This review discusses the current findings in the development of small molecules for COVID-19 treatment. Moreover, their detailed mechanism of action, chemical structures, and preclinical and clinical efficacies are discussed.

Repositioning of anti-dengue compounds against SARS-CoV-2 as viral polyprotein processing inhibitor

A therapy for COVID-19 (Coronavirus Disease 19) caused by Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) remains elusive due to the lack of an effective antiviral therapeutic molecule. The SARS-CoV-2 main protease (Mpro), which plays a vital role in the viral life cycle, is one of the most studied and validated drug targets. In Several prior studies, numerous possible chemical entities were proposed as potential Mpro inhibitors; however, most failed at various stages of drug discovery. Repositioning of existing antiviral compounds accelerates the discovery and development of potent therapeutic molecules. Hence, this study examines the applicability of anti-dengue compounds against the substrate binding site of Mpro for disrupting its polyprotein processing mechanism. An in-silico structure-based virtual screening approach is applied to screen 330 experimentally validated anti-dengue compounds to determine their affinity to the substrate binding site of Mpro. This study identified the top five compounds (CHEMBL1940602, CHEMBL2036486, CHEMBL3628485, CHEMBL200972, CHEMBL2036488) that showed a high affinity to Mpro with a docking score > -10.0 kcal/mol. The best-docked pose of these compounds with Mpro was subjected to 100 ns molecular dynamic (MD) simulation followed by MM/GBSA binding energy. This showed the maximum stability and comparable ΔG binding energy against the reference compound (X77 inhibitor). Overall, we repurposed the reported anti-dengue compounds against SARS-CoV-2-Mpro to impede its polyprotein processing for inhibiting SARS-CoV-2 infection.

Pharmacophore-Based Virtual Screening and Molecular Dynamics Simulation for Identification of a Novel DNA Gyrase B Inhibitor with Benzoxazine Acetamide Scaffold

DNA gyrase B is one of the enzyme targets for antimicrobial drug development, and its absence in mammals makes it a suitable target for the creation of safe antibacterial drugs. We identified six novel hits as DNA gyrase B inhibitors in the present study by employing 3D-pharmacophore structure-based virtual screening. The lead compounds complied with drug-likeness rules and lacked toxicity. Compound 4 (ZINC32858011) showed the highest inhibitory activity with an IC₅₀ value of 6.3 ± 0.1 μM against the DNA gyrase enzyme. In contrast, the positive controls ciprofloxacin and novobiocin used in enzyme inhibition assay had IC₅₀ values of 14.4 ± 0.2 and 12.4 ± 0.2 μM, respectively. The molecular docking of the six hits demonstrated that compounds 1, 2, 4, and 6 had suitable fitting modes inside the binding pocket. Molecular dynamics simulations were carried out for the six hits and the rmsd, rmsf, radius of gyration, and solvent accessible surface area parameters obtained from 100 ns molecular dynamics simulations for the six compounds complexed with a DNA gyrase B protein indicated that compound 4 (ZINC32858011) formed the most stable complex with DNA gyrase B. The binding free energy calculation with the MM-PBSA method suggested that the van der Waals interaction, followed by electrostatic force, played a significant role in the binding. Per-residue free binding energy decomposition showed that Ile78 contributed the most for the binding energy followed by Asn46, Asp49, Glu50, Asp73, Ile78, Pro79, Ala86, Ile90, Val120, Thr165, and Val167.

A Structural Approach to Anti-Virulence: A Discovery Pipeline

The anti-virulence strategy is designed to prevent bacterial virulence factors produced by pathogenic bacteria from initiating and sustaining an infection. One family of bacterial virulence factors is the mono-ADP-ribosyltransferase toxins, which are produced by pathogens as tools to compromise the target host cell. These toxins are bacterial enzymes that exploit host cellular NAD+ as the donor substrate to modify an essential macromolecule acceptor target in the host cell. This biochemical reaction modifies the target macromolecule (often protein or DNA) and functions in a binary fashion to turn the target activity on or off by blocking or impairing a critical process or pathway in the host. A structural biology approach to the anti-virulence method to neutralize the cytotoxic effect of these factors requires the search and design of small molecules that bind tightly to the enzyme active site and prevent catalytic function essentially disarming the pathogen. This method requires a high-resolution structure to serve as the model for small molecule inhibitor development, which illuminates the path to drug development. This alternative strategy to antibiotic therapy represents a paradigm shift that may circumvent multi-drug resistance in the offending microbe through anti-virulence therapy. In this report, the rationale for the anti-virulence structural approach will be discussed along with recent efforts to apply this method to treat honey bee diseases using natural products.

In silico analysis of SARS-CoV-2 papain-like protease potential inhibitors

The emergent outbreak caused by severe acute respiratory syndrome coronavirus 2 continues spreading and causing huge social and economic disruption. Papain-like protease (PLpro) has a crucial role in the cleavage of viral polyproteins, and disruption of host responses. PLpro is considered an important goal for the development of SARS-CoV-2 inhibitors. ZINC101291108 (lead 1) and ZINC16449029 (lead 2) were identified as potent SARS-CoV-2 PLpro inhibitors with IC₅₀ values of 0.085 μM and 0.063 μM, respectively. Molecular dynamics simulations (MD) were carried out for lead 1, 2 and several reported SARS-CoV-2 inhibitors. Analysis results of the simulations confirmed the stability of both compounds and showed that they adopted two confirmations along the simulation period. The per-residue decomposition results revealed that the key residues involved in inhibitor binding were E167, P247, P248, Y264, Y268 and Q269. H-bond analyses showed H-bonds with G266 and N267 and salt bridges with G209 and Y273, which are essential for strengthening the substrate-binding pocket. Both inhibitors showed hydrophobic interactions with the S4 site and BL2 loop residues. The RMSD of the BL2 loop with the two inhibitors was investigated, and the results showed that the Y268 and Q269 BL2 loop residues moved outward to accommodate the large size of lead 2. The van der Waals interaction was the main energy contribution that stabilized lead 2, while van der Waals and electrostatic interactions were the main energy contributions stabilizing lead 1. Rational design strategies were suggested to replace the 2-(2-hydroxybenzylidene) hydrazine moiety with naphthalene or nitrobenzene at the P4 position of lead 2 and introduce polar substituents as aniline and benzoate groups at position P1 to enhance hydrophobic interactions and H-bonds, respectively.

The emergent outbreak caused by severe acute respiratory syndrome coronavirus 2 continues spreading and causing huge social and economic disruption.