Accurate Prediction of Inhibitor Binding to HIV-1 Protease Using CANDOCK

The human immunodeficiency virus 1 (HIV-1) protease is an important target for treating HIV infection. Our goal was to benchmark a novel molecular docking protocol and determine its effectiveness as a therapeutic repurposing tool by predicting inhibitor potency to this target. To accomplish this, we predicted the relative binding scores of various inhibitors of the protease using CANDOCK, a hierarchical fragment-based docking protocol with a knowledge-based scoring function. We first used a set of 30 HIV-1 protease complexes as an initial benchmark to optimize the parameters for CANDOCK. We then compared the results from CANDOCK to two other popular molecular docking protocols Autodock Vina and Smina. Our results showed that CANDOCK is superior to both of these protocols in terms of correlating predicted binding scores to experimental binding affinities with a Pearson coefficient of 0.62 compared to 0.48 and 0.49 for Vina and Smina, respectively. We further leveraged the Database of Useful Decoys: Enhanced (DUD-E) HIV protease set to ascertain the effectiveness of each protocol in discriminating active versus decoy ligands for proteases. CANDOCK again displayed better efficacy over the other commonly used molecular docking protocols with area under the receiver operating characteristic curve (AUROC) of 0.94 compared to 0.71 and 0.74 for Vina and Smina. These findings support the utility of CANDOCK to help discover novel therapeutics that effectively inhibit HIV-1 and possibly other retroviral proteases.

(S)-5-Methylmellein Isolated from an Endogenous Lichen Fungus Rosellinia corticium as a Potent Inhibitor of Human Monoamine Oxidase A

In this study, the inhibitory activities against human monoamine oxidases (hMAOs) were evaluated using a library of 195 endogenous lichen fungi from Ukraine. Among them, the extract ELF68 of the endogenous fungus Rosellinia corticium from the lichen Pseudevernia furfuracea (L.) Zopf. exhibited the strongest inhibitory activity against hMAO-A. Using the activity-guided method, (S)-5-methylmellein (5MM) was isolated from the extract and had an IC50 value of 5.31 µM for hMAO-A with a lower potency for hMAO-B (IC50 = 9.15 µM). Compound 5MM also moderately inhibited acetylcholinesterase (IC50 = 27.07 µM) but very weakly inhibited butyrylcholinesterase and β-secretase. Compound 5MM had a Ki value of 2.45 μM and was a reversible competitive inhibitor of hMAO-A. A molecular docking study predicted that (S)-5MM showed higher binding affinity for hMAO-A (−6.8 kcal/mol) than hMAO-B (−6.4 kcal/mol). Its isomer, (R)-5MM, exhibited lower binding affinities for hMAO-A (−6.6 kcal/mol) and hMAO-B (−5.2 kcal/mol), compared to (S)-5MM. The S-form interacted with hMAO-A through hydrogen bonding with the Phe208 residue (distance: 1.972 Å), while the R-form interacted with the Asn181 residue (2.375 Å). The results of an in silico pharmacokinetic analysis indicated that 5MM did not violate Lipinski’s five rules and showed high gastrointestinal absorption and blood–brain barrier permeability. These results suggest that 5MM can be considered a candidate in the treatment of neuropsychiatric disorders, such as depression and cardiovascular disease.

In Silico Screening of Potential Phytocompounds from Several Herbs against SARS-CoV-2 Indian Delta Variant B.1.617.2 to Inhibit the Spike Glycoprotein Trimer

In October 2020, the SARS-CoV-2 B.1.617 lineage was discovered in India. It has since become a prominent variant in several Indian regions and 156 countries, including the United States of America. The lineage B.1.617.2 is termed the delta variant, harboring diverse spike mutations in the N-terminal domain (NTD) and the receptor-binding domain (RBD), which may heighten its immune evasion potentiality and cause it to be more transmissible than other variants. As a result, it has sparked substantial scientific investigation into the development of effective vaccinations and anti-viral drugs. Several efforts have been made to examine ancient medicinal herbs known for their health benefits and immune-boosting action against SARS-CoV-2, including repurposing existing FDA-approved anti-viral drugs. No efficient anti-viral drugs are available against the SARS-CoV-2 Indian delta variant B.1.617.2. In this study, efforts were made to shed light on the potential of 603 phytocompounds from 22 plant species to inhibit the Indian delta variant B.1.617.2. We also compared these compounds with the standard drug ceftriaxone, which was already suggested as a beneficial drug in COVID-19 treatment; these compounds were compared with other FDA-approved drugs: remdesivir, chloroquine, hydroxy-chloroquine, lopinavir, and ritonavir. From the analysis, the identified phytocompounds acteoside (−7.3 kcal/mol) and verbascoside (−7.1 kcal/mol), from the plants Clerodendrum serratum and Houttuynia cordata, evidenced a strong inhibitory effect against the mutated NTD (MT-NTD). In addition, the phytocompounds kanzonol V (−6.8 kcal/mol), progeldanamycin (−6.4 kcal/mol), and rhodoxanthin (−7.5 kcal/mol), from the plant Houttuynia cordata, manifested significant prohibition against RBD. Nevertheless, the standard drug, ceftriaxone, signals less inhibitory effect against MT-NTD and RBD with binding affinities of −6.3 kcal/mol and −6.5 kcal/mol, respectively. In this study, we also emphasized the pharmacological properties of the plants, which contain the screened phytocompounds. Our research could be used as a lead for future drug design to develop anti-viral drugs, as well as for preening the Siddha formulation to control the Indian delta variant B.1.617.2 and other future SARS-CoV-2 variants.

HIGH-THROUGHPUT SCREENING AND DYNAMIC STUDIES OF SELECTED COMPOUNDS AGAINST SARS-COV-2

Objective: This study was aimed to analyze the inhibitory effect of the drugs used in nanocarrier as well as nanoparticles formulation based drug delivery system selected from PubChem database literature against 3CLpro (3C-like protease) receptor of SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) by implementing several in silico analysis techniques. Methods: This paper detailed a molecular docking-based virtual screening of 5240 compounds previously utilized in nanoparticle and nanocarrier drug delivery systems utilizing AutoDock Vina software on 3CL protease to discover potential inhibitors using a molecular docking technique. Results: According to the results of the screening, the top two compounds, PubChem Id 58823276 and PubChem Id 60838 exhibited a high affinity for the 3CL protease binding region. Their binding affinities were-9.6 and-8.5 kJ/mol, indicating that they were tightly bound to the target receptor, respectively. These results outperformed those obtained using the co-crystallized native ligand, which exhibited a binding affinity of-7.4 kJ/mol. PubChem Id 60838, the main hit compound in terms of both binding affinity and ADMET analysis, displayed substantial deformability after MD simulation. As a result of the VS and molecular docking techniques, novel 3CL protease inhibitors from the PubChem database were discovered using the Lipinski rule of five and functional molecular contacts with the target protein, as evidenced by the findings of this work. Conclusion: The findings suggest that the compounds discovered may represent attractive opportunities for the development of COVID-19 3CLpro inhibitors and that they need further evaluation and investigation.

Interactions of Nereistoxin and Its Analogs with Vertebrate Nicotinic Acetylcholine Receptors and Molluscan ACh Binding Proteins

Nereistoxin (NTX) is a marine toxin isolated from an annelid worm that lives along the coasts of Japan. Its insecticidal properties were discovered decades ago and this stimulated the development of a variety of insecticides such as Cartap that are readily transformed into NTX. One unusual feature of NTX is that it is a small cyclic molecule that contains a disulfide bond. In spite of its size, it acts as an antagonist at insect and mammalian nicotinic acetylcholine receptors (nAChRs). The functional importance of the disulfide bond was assessed by determining the effects of inserting a methylene group between the two sulfur atoms, creating dimethylaminodithiane (DMA-DT). We also assessed the effect of methylating the NTX and DMA-DT dimethylamino groups on binding to three vertebrate nAChRs. Radioligand receptor binding experiments were carried out using washed membranes from rat brain and fish (Torpedo) electric organ; [3H]-cytisine displacement was used to assess binding to the predominantly high affinity alpha4beta2 nAChRs and [125I]-alpha-bungarotoxin displacement was used to measure binding of NTX and analogs to the alpha7 and skeletal muscle type nAChRs. While the two quaternary nitrogen analogs, relative to their respective tertiary amines, displayed lower α4β2 nAChR binding affinities, both displayed much higher affinities for the Torpedo muscle nAChR and rat alpha7 brain receptors than their respective tertiary amine forms. The binding affinities of DMA-DT for the three nAChRs were lower than those of NTX and MeNTX. An AChBP mutant lacking the C loop disulfide bond that would potentially react with the NTX disulfide bond displayed an NTX affinity very similar to the parent AChBP. Inhibition of [3H]-epibatidine binding to the AChBPs was not affected by exposure to NTX or MeNTX for up to 24 hr prior to addition of the radioligand. Thus, the disulfide bond of NTX is not required to react with the vicinal disulfide in the AChBP C loop for inhibition of [3H]-epibatidine binding. However, a reversible disulfide interchange reaction of NTX with nAChRs might still occur, especially under reducing conditions. Labeled MeNTX, because it can be readily prepared with high specific radioactivity and possesses relatively high affinity for the nAChR-rich Torpedo nAChR, would be a useful probe to detect and identify any nereistoxin adducts.

Tripeptides from Allium subhirsitum L. extracts: Pharmacokinetics properties, toxicity prediction and in silico study against SARS-CoV-2 enzymes and pro-inflammatory proteins

Developing new prophylactic and therapeutic agents with broad-spectrum antiviral activities is urgently needed to combat emerging human severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Since no available clinically antiviral drugs have been approved to eradicate COVID-19 as of the writing of this report, this study aimed to investigate bioactive short peptides from Allium subhirsutum L. (Hairy garlic) extracts identified through HR-LC/MS analysis that could potentially hinder the multiplication cycle of SARS-CoV-2 via molecular docking study. The obtained promising results showed that the peptides (Asn-Asn-Asn) possess the highest binding affinities of -8.4 kcal/mol against S protein, (His-Phe-Gln) of -9.8 kcal/mol and (Gln-His-Phe) of -9.7 kcal/mol towards hACE2, (Thr-Leu-Trp) of -10.3 kcal/mol and (Gln-Phe-Tyr) of -9.8 kcal/mol against furin. Additionally, the identified peptides show strong interactions with the targeted and pro-inflammatory ranging from -8.1 to -10.5 kcal/mol for NF−κB-inducing kinase (NIK), from -8.2 to -10 kcal/mol for phospholipase A2 (PLA2), from -8.0 to -10.7 kcal/mol for interleukin-1 receptor-associated kinase 4 (IRAK-4), and from -8.6 to -11.6 kcal/mol for the cyclooxygenase 2 (COX2) with Gln-Phe-Tyr model seems to be the most prominent. Results from pharmacophore, drug-likeness and ADMET prediction analyses clearly evidenced the usability of the peptides to be developed as an effective drug, beneficial for COVID-19 treatment.

Identifying the Hot Spot Residues of the SARS-CoV-2 Main Protease Using MM-PBSA and Multiple Force Fields

In this study, we investigated the binding affinities between the main protease of SARS-CoV-2 virus (Mpro) and its various ligands to identify the hot spot residues of the protease. To benchmark the influence of various force fields on hot spot residue identification and binding free energy calculation, we performed MD simulations followed by MM-PBSA analysis with three different force fields: CHARMM36, AMBER99SB, and GROMOS54a7. We performed MD simulations with 100 ns for 11 protein–ligand complexes. From the series of MD simulations and MM-PBSA calculations, it is identified that the MM-PBSA estimations using different force fields are weakly correlated to each other. From a comparison between the force fields, AMBER99SB and GROMOS54a7 results are fairly correlated while CHARMM36 results show weak or almost no correlations with the others. Our results suggest that MM-PBSA analysis results strongly depend on force fields and should be interpreted carefully. Additionally, we identified the hot spot residues of Mpro, which play critical roles in ligand binding through energy decomposition analysis. It is identified that the residues of the S4 subsite of the binding site, N142, M165, and R188, contribute strongly to ligand binding. In addition, the terminal residues, D295, R298, and Q299 are identified to have attractive interactions with ligands via electrostatic and solvation energy. We believe that our findings will help facilitate developing the novel inhibitors of SARS-CoV-2.

Potential Inhibitors Identification of Severe Acute Respiratory Syndrome-Related Coronavirus 2 (SARS-CoV-2) Angiotensin-Converting Enzyme 2 and Main Protease from Anatolian Traditional Plants

Background: The 2019 novel coronavirus disease (COVID-19) has caused a global health catastrophe by affecting the whole human population around the globe. Unfortunately, there is no specific medication or treatment for COVID-19 currently available. Objective: It’s extremely necessary to apply effective drug treatment in order to end the pandemic period and return daily life to normal. In terms of the urgency of treatment, rather than focusing on the discovery of novel compounds, it is critical to explore the effects of existing herbal agents with proven antiviral properties on the virus. Method: Molecular docking studies were carried out with three different methods, Glide extra precision (XP) docking, Induced Fit docking (IFD), and Molecular Mechanics/Generalized Born Surface Area (MM/GBSA), to determine the potential effects of 58 phytochemicals in the content of Rosmarinus officinalis, Thymbra spicata, Satureja thymbra, and Stachys lavandulifolia plants -have antiviral and antibacterial effects- against Main Protease (Mpro) and Angiotensin Converting Enzyme 2 (ACE2) enzymes. Results: 7 compounds stand out among all molecules by showing very high binding affinities. According to our findings, the substances chlorogenic acid, rosmarinic acid, and rosmanol exhibit extremely significant binding affinities for both Mpro and ACE2 enzymes. Furthermore, it was discovered that carnosic acid and alpha-cadinol showed potential anti-Mpro activity, whereas caffeic acid and carvacrol had promising anti-ACE2 activity. Conclusion: Chlorogenic acid, rosmarinic acid, rosmanol, carnosic acid, alpha-cadinol, caffeic acid, and carvacrol compounds have been shown to be powerful anti-SARS-COV-2 agents in docking simulations against Mpro and ACE2 enzymes, as well as ADME investigations.

In Silico Design of Dual-Binding Site Anti-Cholinesterase Phytochemical Heterodimers as Treatment Options for Alzheimer's Disease

The number of patients with neurodegenerative diseases, particularly Alzheimer’s disease (AD), continues to grow yearly. Cholinesterase inhibitors (ChEIs) represent the first-line symptomatic drug treatment for mild-to-moderate AD; however, there is an unmet need to produce ChEIs with improved efficacy and reduced side effects. Herein, phytochemicals with reported anti-acetylcholinesterase (AChE) activity were ranked in silico for their anti-AChE potential. Ligands with a similar or higher binding affinity to AChE than galantamine were then selected for the design of novel dual-binding site heterodimeric drugs. In silico molecular docking of heterodimers with the target enzymes, AChE and butyrylcholinesterase (BuChE), were performed, and anti-cholinesterase binding affinities were compared with donepezil. Drug-likeliness properties and toxicity of the heterodimers were assessed using the SwissADME and ProTox-II webservers. Nine phytochemicals displayed similar or higher binding affinities to AChE than galantamine: sanguinarine > huperzine A > chelerythrine > yohimbine > berberine > berberastine > naringenin > akuammicine > carvone. Eleven heterodimeric ligands were designed with phytochemicals separated by four- or five-carbon alkyl-linkers. All heterodimers were theoretically potent AChE and BuChE dual-binding site inhibitors, with the highest affinity achieved with huperzine-4C-naringenin, which displayed 34% and 26% improved affinity to AChE and BuChE, respectively, then the potent ChEI drug, donepezil. Computational pharmacokinetic and pharmacodynamic screening suggested that phytochemical heterodimers would display useful gastrointestinal absorption and with relatively low predicted toxicity. Collectively, the present study suggests that phytochemicals could be garnered for the provision of novel ChEIs with enhanced drug efficacy and low toxicity.