Molecular Mechanics Generalized Born Surface Area Research Articles

Computer-aided drug design approaches applied toscreen natural product's structural analogs targeting arginase in Leishmania spp.

Introduction: Leishmaniasis is a disease with high mortality rates and approximately 1.5 million new cases each year. Despite the new approaches and advances to fight the disease, there are no effective therapies. Methods: Hence, this study aims to screen for natural products' structural analogs as new drug candidates against leishmaniasis. We applied Computer-aided drug design (CADD) approaches, such as virtual screening, molecular docking, molecular dynamics simulation, molecular mechanics-generalized Born surface area (MM-GBSA) binding free estimation, and free energy perturbation (FEP) aiming to select structural analogs from natural products that have shown anti-leishmanial and anti-arginase activities and that could bind selectively against the Leishmania arginase enzyme. Results: The compounds 2H-1-benzopyran, 3,4-dihydro-2-(2-methylphenyl)-(9CI), echioidinin, and malvidin showed good results against arginase targets from three parasite species and negative results for potential toxicities. The echioidinin and malvidin ligands generated interactions in the active center at pH 2.0 conditions by MM-GBSA and FEP methods. Conclusions: This work suggests the potential anti-leishmanial activity of the compounds and thus can be further in vitro and in vivo experimentally validated.

The inhibitory potential of chemical constituents of Ficus carica targeting interleukin-6 (IL-6) mediated inflammation.

Inflammation is an innate reaction of the body of an individual when subjected to the noxious factors repeatedly. Pharmacological approaches focused at disrupting cytokine signaling networks have become significant therapeutic alternatives for the treatment of inflammatory illnesses, cancer and autoimmune disorders. High levels of inflammatory mediators, particularly interleukin IL-1, IL-6, IL-18, IL-12, and tumor necrosis factor alpha leads to a cytokine storm in the body. Among all the released cytokines in a patient suffering from inflammatory disorder, IL-6 mediator has a pivotal role in this inflammatory cascade which progresses to a cytokine storm. Therefore, the blockage of the IL-6 inflammatory mediator could be a promising treatment option for the patients with hyper inflammatory conditions. The phytochemicals could provide the new lead compounds against the IL-6 mediator. Ficus carica has been the ideal plant of research and investigation due to its commercial, economic and medical importance. The anti-inflammatory properties of F. carica were further investigated by in silico and in vivo approaches. The docking scores of Cyanidin-3,5-diglucoside, Kaempferol-7-O-rutinoside, Cyanidin-3-rhamnoglucoside, and Rutin are -9.231, -8.921, -8.840, and -8.335 Kcal/mole respectively. The free energy of binding and stability of the docked complexes of these top four phytochemicals with the IL-6 were further analyzed by Molecular Mechanics-Generalized Born Surface Area and Molecular Dynamic simulations, respectively. The in vivo anti-inflammatory carrageenan-induced rat paw edema model was used for the validation of in silico results. The maximum percentage paw edema inhibition with petroleum ether and ethyl acetate was 70.32% and 45.05%, respectively. The in vivo anti-inflammatory activity confirms the anti-inflammatory potential of F. carica. Therefore, it is predicted that Cyanidin-3,5-diglucoside, Kaempferol-7-O-rutinoside, Cyanidin-3-rhamnoglucoside, and Rutin have the potential to inhibit the IL-6 mediator which will aid in mitigating the cytokine storm in patients with acute inflammations.

Computational drug repurposing of Akt-1 allosteric inhibitors for non-small cell lung cancer

Non-small cell lung carcinomas (NSCLC) are the predominant form of lung malignancy and the reason for the highest number of cancer-related deaths. Widespread deregulation of Akt, a serine/threonine kinase, has been reported in NSCLC. Allosteric Akt inhibitors bind in the space separating the Pleckstrin homology (PH) and catalytic domains, typically with tryptophan residue (Trp-80). This could decrease the regulatory site phosphorylation by stabilizing the PH-in conformation. Hence, in this study, a computational investigation was undertaken to identify allosteric Akt-1 inhibitors from FDA-approved drugs. The molecules were docked at standard precision (SP) and extra-precision (XP), followed by Prime molecular mechanics—generalized Born surface area (MM-GBSA), and molecular dynamics (MD) simulations on selected hits. Post XP-docking, fourteen best hits were identified from a library of 2115 optimized FDA-approved compounds, demonstrating several beneficial interactions such as pi–pi stacking, pi-cation, direct, and water-bridged hydrogen bonds with the crucial residues (Trp-80 and Tyr-272) and several amino acid residues in the allosteric ligand-binding pocket of Akt-1. Subsequent MD simulations to verify the stability of chosen drugs to the Akt-1 allosteric site showed valganciclovir, dasatinib, indacaterol, and novobiocin to have high stability. Further, predictions for possible biological interactions were performed using computational tools such as ProTox-II, CLC-Pred, and PASSOnline. The shortlisted drugs open a new class of allosteric Akt-1 inhibitors for the therapy of NSCLC.

Drug repositioning to discover novel ornithine decarboxylase inhibitors against visceral leishmaniasis.

Visceral leishmaniasis (VL) is caused by Leishmania donovani (Ld), and most cases occur in Brazil, East Africa, and India. The treatment for VL is limited and has many adverse effects. The development of safer and more efficacious drugs is urgently needed. Drug repurposing is one of the best processes to repurpose existing drugs. Ornithine decarboxylase (ODC) is an important target against L. donovani in the polyamine biosynthesis pathway. In this study, we have modeled the 3D structure of ODC and performed high-throughput virtual screening of 8630 ZINC database ligands against Leishmania donovani ornithine decarboxylase (Ld ODC), selecting 45 ligands based on their high binding score. It is further validated through molecular docking simulation and the selection of the top two lead molecules (ceftaroline fosamil and rimegepant) for Molecular Dynamics (MD) simulation, Density functional theory (DFT), and molecular mechanics generalized born surface area (MMGBSA) analysis. The results showed that the binding affinities of ceftaroline fosamil, and rimegepant are, respectively, -10.719 and 10.159 kcal/mol. The docking complexes of the two lead compounds, ceftaroline fosamil, and rimegepant, with the target ODC, were found stable during molecular dynamics simulations. Furthermore, the analysis of MMGBSA revealed that these compounds had a high binding free energy. The DFT analysis showed that the top lead molecules were more reactive than the standard drug (pentamidine). In-silico findings demonstrated that ceftaroline fosamil, and rimegepant might be recognized as potent antagonists against ODC for the treatment of VL.

Imidazole and Biphenyl Derivatives as Anti-cancer Agents for Glioma Therapeutics: Computational Drug Repurposing Strategy.

Targeting mutated isocitrate dehydrogenase 1 (mIDH1) is one of the key therapeutic strategies for the treatment of glioma. Few inhibitors, such as ivosidenib and vorasidenib, have been identified as selective inhibitors of mIDH1. However, dose-dependent toxicity and limited brain penetration of the blood-brain barrier remain the major limitations of the treatment procedures using these inhibitors. In the present study, computational drug repurposing strategies were employed to identify potent mIDH1- specific inhibitors from the 11,808 small molecules listed in the DrugBank repository. Tanimoto coefficient (Tc) calculations were initially used to retrieve compounds with structurally similar scaffolds to ivosidenib. The resultant compounds were then subjected to molecular docking to discriminate the binders from the non-binders. The binding affinities and pharmacokinetic properties of the screened compounds were examined using prime Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) and QikProp algorithm, respectively. The conformational stability of these molecules was validated using 100 ns molecular dynamics simulation. Together, these processes led to the identification of three-hit molecules, namely DB12001, DB08026, and DB03346, as potential inhibitors of the mIDH1 protein. Of note, the binding free energy calculations and MD simulation studies emphasized the greater binding affinity and structural stability of the hit compounds towards the mIDH1 protein. The collective evidence from our study indicates the activity of DB12001 against recurrent glioblastoma, which, in turn, highlights the accuracy of our adapted strategy. Hence, we hypothesize that the identified lead molecules could be translated for the development of mIDH1 inhibitors in the near future.

Pharmacophore Modeling Guided by Conformational Dynamics Reveals Potent Anticancer Agents

Targeting the interaction between tumor suppressor p53 and murine double minute 2(MDM2) has been an attractive therapeutic strategy of recent cancer research. There are a few number of MDM2-targeted anticancer drug molecules undergoing clinical trials, yet none of them have been approved so far. In this study, a new approach is employed in which dynamics of MDM2 obtained by elastic network models are used as a guide in the generation of the ligand-based pharmacophore model prior to virtual screening. Hit molecules exhibiting high affinity to MDM2 were captured and tested by rigid and induced-fit molecular docking. The knowledge of the binding mechanism was used while creating the induced-fit docking criteria. Application of Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) method provided an accurate prediction of the binding free energy values. Two leading hit molecules which have shown better docking scores, binding free energy values and drug-like molecular properties were identified. These hits exhibited extra intermolecular interactions with MDM2, indicating a stable complex formation and hence would be further tested in vitro. Finally, the combined computational strategy employed in this study can be a promising tool in drug design for the discovery of potential new hits.

Integrated molecular and quantum mechanical approach to identify novel potent natural bioactive compound against 2′-O-methyltransferase (nsp16) of SARS-CoV-2

With the advent of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) outbreak, efforts are still in progress to find out a functional cure for the infection. Among the various protein targets, nsp16 capping protein is one of the vital targets for drug development as it protects the virus against the host cell nucleases and evading innate immunity. The nsp16 protein forms a heterodimer with a co-factor nsp10 and triggers 2′-O-methyltransferase activity which catalyzes the conversion of S-adenosyl methionine into S-adenosyl homocysteine. The free methyl group is transferred to the 2′-O position on ribose sugar at the 5′ end of mRNA to form the cap-1 structure which is essential for replication of the virus and evading the innate immunity of the host. In this study, we identify a potential lead natural bioactive compound against nsp16 protein by systematic cheminformatic analysis of more than 144k natural compounds. Virtual screening, molecular docking interactions, ADMET profiling, molecular dynamics (MD) simulations, molecular mechanics-generalized born surface area (MM-GBSA), free energy analysis and density functional theory analysis were used to discover the potential lead compound. Our investigation revealed that ZINC8952607 (methyl-[(6-methyl-2,3,4,9-tetrahydro-1H-carbazol-1-yl)aminomethyl]BLAHone) has the greatest binding affinity and best pharmacokinetic parameters due to presence of carbazol and BLAHone (biaryl moiety). Further, time-dependent MD simulation analysis substantiates the stability and rigidness of nsp16 protein even after interaction with the lead compound. We believe that the compound ZINC8952607 might establish as a novel natural drug candidate against CoVID-19 infection. Communicated by Ramaswamy H. Sarma

Structural identification of novel pyrimidine derivatives as epidermal growth factor receptor inhibitors using 3D QSAR, molecular docking, and MMGBSA analysis: a rational approach in anticancer drug design

Non-small cell lung cancer (NSCLC) has evolved into the deadliest in the present scenario. The progression of NSCLC is mainly due to the dysregulation of the tyrosine kinase family's epidermal growth factor receptor (EGFR). Thus, EGFR has been widely studied as a major target in the treatment of NSCLC, but the lack of selectivity and drug resistance limit the use of existing therapeutic agents. Considering the urgent necessity for the advanced development of EGFR inhibitors, we have implemented a three-dimensional structure-activity relationship (3D QSAR), molecular docking, and MMGBSA studies on a series of pyrimidine derivatives. In the 3D QSAR, the comparative molecular field analysis model (CoMFA) was obtained with a correlation coefficient (r<sup>2</sup>) = 0.698, cross-validated correlation coefficient (q<sup>2</sup>) = 0.541, and predictive r<sup>2</sup> (r<sup>2</sup><sub>pred</sub>) = 0.509. The comparative molecular similarity indices analysis (CoMSIA) model also displayed similar results with r<sup>2 </sup>= 0.72, q<sup>2 </sup>= 0.586, and r<sup>2</sup><sub>pred</sub>= 0.495. The statistical parameters fulfill the acceptability criteria of the models. Docking studies revealed the binding interactions of the pyrimidine derivatives with double mutant EGFR<sup>L858R/T790M</sup>. Docking scores of the top two selected compounds 29 and 34 were 92.99 and 92.13, respectively. Analyzing 3D QSAR contour plots and docking results reviewed some important structural attributes of EGFR <sup>L858R/T790M</sup> selective inhibitors, which directed the designing of some new molecules. The designed compounds showed good predictive activity and exhibited higher binding interactions with EGFRL858R/T790M than the reference ligand gefitinib. Moreover, to evaluate the binding of selected top hits from docking and designed compounds, MMGBSA (Molecular Mechanics-Generalized Born Surface Area) analysis was performed, which revealed that the designed compound (N7) showed a good binding affinity with EGFR<sup>L858R/T790M</sup> (dG = -68.59 kcal/mol) as compared to other compounds. Further, in silico ADME predictions revealed the drug-likeness of the designed compounds. Thus, this work will guide researchers in future developments of pyrimidine derivatives as EGFR inhibitors.

A deep learning and docking simulation-based virtual screening strategy enables the rapid identification of HIF-1α pathway activators from a marine natural product database

Artificial Intelligence is hailed as a cutting-edge technology for accelerating drug discovery efforts, and our goal was to validate its potential in predicting pharmacological inhibitors of EGLN1 using a deep learning-based architecture, one of its subsidiaries. Egl nine homolog 1 (EGLN1) inhibition prevents poly ubiquitination-mediated proteosomal destruction HIF-1α. The pharmacological interventions aimed at stabilizing HIF-1α have the potential to be a promising treatment option for a range of human diseases, including ischemic stroke. To unveil a novel EGLN1 inhibitor from marine natural products, a custom-based virtual screening was carried out using a Deep Convolutional Neural Network (DCNN) architecture, docking, and molecular dynamics simulation. The custom DCNN model was optimized and further employed to screen marine natural products from the CMNPD database. The docking was performed as a secondary strategy for screened hits. Molecular dynamics (MD) and molecular mechanics/generalized Born surface area (MM-GBSA) were used to analyze inhibitor binding and identify key interactions. The findings support the claim that deep learning-based virtual screening is a rapid, reliable and accurate method of identifying highly contributing drug candidates (EGLN1 inhibitors). This study demonstrates that deep learning architecture can significantly accelerate drug discovery and development, and provides a solid foundation for using (Z)-2-ethylhex-2-enedioic acid [(Z)-2-ethylhex-2-enedioic acid] as a potential EGLN1 inhibitor for treating various health complications. Communicated by Ramaswamy H. Sarma

Triazole based Schiff bases and their oxovanadium(IV) complexes: Synthesis, characterization, antibacterial assay, and computational assessments

The synthesis and characterization of two new Schiff base ligands containing 1,2,4-triazole moieties and their oxovanadium(IV) complexes have been reported. The ligands and their complexes were studied by ultraviolet–visible (UV–Vis), Fourier transform infrared (FTIR), proton nuclear magnetic resonance (1H NMR), electron paramagnetic resonance (EPR), X-ray diffraction (XRD), conductivity measurement, cyclic voltammetry (CV), and elemental analyses. The molar conductance of oxovanadium(IV) complexes were found to be relatively low, depicting their non-electrolytic nature. The XRD patterns reveal the size of particles to be 47.53 nm and 26.28 nm for the two complexes in the monoclinic crystal system. The molecular structures, geometrical parameters, chemical reactivity, stability, and frontier molecular orbital pictures were determined by density functional theory (DFT) calculations. The theoretical vibrational frequencies and EPR g-factors (1.98) were found to correlate well with the experimental values. A distorted square pyramidal geometry with C2 symmetry of the complexes has been proposed from experimental and theoretical results in a synergistic manner. The antimicrobial sensitivity of the ligands and their metal complexes assayed in vitro against four bacterial pathogens viz. Staphylococcus aureus, Escherichia coli, Pseudomonas aeruginosa, and Salmonella Typhi showed that the oxovanadium(IV) complexes are slightly stronger antibacterial agents than their corresponding Schiff base precursors. The binding affinities obtained from the molecular docking calculations with the receptor proteins of bacterial strains (2EUG, 3UWZ, 4GVF, and 4JVD) showed that the Schiff bases and their oxovanadium(IV) complexes have considerable capacity inferring activeness for effective inhibition. The molecular dynamics simulation of a protein-ligand (4JVD-HL2) complex with the best binding affinity of −12.8 kcal/mol for 100 ns showed acceptable stability of the docked pose and binding free energy of −15.17 ± 2.29 kcal/mol from molecular mechanics-generalized Born surface area (MM-GBSA) calculations indicated spontaneity of the reaction. The outcome of the research shows the complementary role of computational methods in material characterization and provides an interesting avenue to pursue for exploring new triazole based Schiff's bases and their vanadium compounds for better properties.

Evaluation of 6-OxP-CD, an Oxime-based cyclodextrin as a viable medical countermeasure against nerve agent poisoning: Experimental and molecular dynamic simulation studies on its inclusion complexes with cyclosarin, soman and VX.

The ability of the cyclodextrin-oxime construct 6-OxP-CD to bind and degrade the nerve agents Cyclosarin (GF), Soman (GD) and S-[2-[Di(propan-2-yl)amino]ethyl] O-ethyl methylphosphonothioate (VX) has been studied using 31P-nuclear magnetic resonance (NMR) under physiological conditions. While 6-OxP-CD was found to degrade GF instantaneously under these conditions, it was found to form an inclusion complex with GD and significantly improve its degradation (t1/2 ~ 2 hrs) relative over background (t1/2 ~ 22 hrs). Consequently, effective formation of the 6-OxP-CD:GD inclusion complex results in the immediate neutralization of GD and thus preventing it from inhibiting its biological target. In contrast, NMR experiments did not find evidence for an inclusion complex between 6-OxP-CD and VX, and the agent's degradation profile was identical to that of background degradation (t1/2 ~ 24 hrs). As a complement to this experimental work, molecular dynamics (MD) simulations coupled with Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) calculations have been applied to the study of inclusion complexes between 6-OxP-CD and the three nerve agents. These studies provide data that informs the understanding of the different degradative interactions exhibited by 6-OxP-CD with each nerve agent as it is introduced in the CD cavity in two different orientations (up and down). For its complex with GF, it was found that the oxime in 6-OxP-CD lies in very close proximity (PGF⋯OOxime ~ 4-5 Å) to the phosphorus center of GF in the 'downGF' orientation for most of the simulation accurately describing the ability of 6-OxP-CD to degrade this nerve agent rapidly and efficiently. Further computational studies involving the center of masses (COMs) for both components (GF and 6-OxP-CD) also provided some insight on the nature of this inclusion complex. Distances between the COMs (ΔCOM) lie closer in space in the 'downGF' orientation than in the 'upGF' orientation; a correlation that seems to hold true not only for GF but also for its congener, GD. In the case of GD, calculations for the 'downGD' orientation showed that the oxime functional group in 6-OxP-CD although lying in close proximity (PGD⋯OOxime ~ 4-5 Å) to the phosphorus center of the nerve agent for most of the simulation, adopts another stable conformation that increase this distance to ~ 12-14 Å, thus explaining the ability of 6-OxP-CD to bind and degrade GD but with less efficiency as observed experimentally (t1/2 ~ 4 hr. vs. immediate). Lastly, studies on the VX:6-OxP-CD system demonstrated that VX does not form a stable inclusion complex with the oxime-bearing cyclodextrin and as such does not interact in a way that is conducive to an accelerated degradation scenario. Collectively, these studies serve as a basic platform from which the development of new cyclodextrin scaffolds based on 6-OxP-CD can be designed in the development of medical countermeasures against these highly toxic chemical warfare agents.

Design and Synthesis of 3-(β-d-Glucopyranosyl)-4-amino/4-guanidino Pyrazole Derivatives and Analysis of Their Glycogen Phosphorylase Inhibitory Potential.

Glycogen phosphorylase (GP) is a key regulator of glucose levels and, with that, an important target for the discovery of novel treatments against type 2 diabetes. β-d-Glucopyranosyl derivatives have provided some of the most potent GP inhibitors discovered to date. In this regard, C-β-d-glucopyranosyl azole type inhibitors proved to be particularly effective, with 2- and 4-β-d-glucopyranosyl imidazoles among the most potent designed to date. His377 backbone C=O hydrogen bonding and ion-ion interactions of the protonated imidazole with Asp283 from the 280s loop, stabilizing the inactive state, were proposed as crucial to the observed potencies. Towards further exploring these features, 4-amino-3-(β-d-glucopyranosyl)-5-phenyl-1H-pyrazole (3) and 3-(β-d-glucopyranosyl)-4-guanidino-5-phenyl-1H-pyrazole (4) were designed and synthesized with the potential to exploit similar interactions. Binding assay experiments against rabbit muscle GPb revealed 3 as a moderate inhibitor (IC50 = 565 µM), but 4 displayed no inhibition at 625 µM concentration. Towards understanding the observed inhibitions, docking and post-docking molecular mechanics-generalized Born surface area (MM-GBSA) binding free energy calculations were performed, together with Monte Carlo and density functional theory (DFT) calculations on the free unbound ligands. The computations revealed that while 3 was predicted to hydrogen bond with His377 C=O in its favoured tautomeric state, the interactions with Asp283 were not direct and there were no ion-ion interactions; for 4, the most stable tautomer did not have the His377 backbone C=O interaction and while ion-ion interactions and direct hydrogen bonding with Asp283 were predicted, the conformational strain and entropy loss of the ligand in the bound state was significant. The importance of consideration of tautomeric states and ligand strain for glucose analogues in the confined space of the catalytic site with the 280s loop in the closed position was highlighted.

Targeting Olokizumab-Interleukin 6 interaction interface to discover novel IL-6 inhibitors

The IL-6/IL-6R or IL-6/GP130 protein-protein interactions play a significant role in controlling the development of chronic inflammatory diseases, such as rheumatoid arthritis, Castleman disease, psoriasis, and, most recently, COVID-19. Modulating or antagonizing protein-protein interactions of IL6 binding to its receptors by oral drugs promises similar efficacy to biological therapy in patients, namely monoclonal antibodies. In this study, we used a crystal structure of the Fab part of olokizumab in a complex with IL-6 (PDB ID: 4CNI) to uncover starting points for small molecule IL-6 antagonist discovery. Firstly, a structure‑based pharmacophore model of the protein active site cavity was generated to identify possible candidates, followed by virtual screening with a significant database Drugbank. After the docking protocol validation, a virtual screening by molecular docking was carried out and a total of 11 top hits were reported. Detailed analysis of the best scoring molecules was performed with ADME/T analysis and molecular dynamics simulation. Furthermore, the Molecular Mechanics-Generalized Born Surface Area (MM/GBSA) technique has been utilized to evaluate the free binding energy. Based on the finding, one newly obtained compound in this study, namely DB15187, may serve as a lead compound for the discovery of IL-6 inhibitors.Communicated by Ramaswamy H. Sarma

Computational evaluation of bioactive compounds from Viscum album (mistletoe) as inhibitors of p63 for pancreatic cancer treatment

Pancreatic ductal adenocarcinoma is an aggressive malignancy usually detectable at the advanced stage, with a 5-year survival rate of less than 8%. It has been reported that a gene called tumor-protein 63 (TP63) is expressed in an aggressive form of pancreatic cancer with a squamous signature. Thus, inhibiting the activity of p63 can be a means of treating and managing PDA. Different studies have shown that plant constituents are rich and can be a promising source for discovering drug candidates. The extract from mistletoe (Viscum album) is known to contain anticancer compounds; however, the specific molecular mechanism of the bioactive compounds is unknown. This study examines the pancreatic cancer therapeutic potential of the bioactive compounds in the flavonoid and phenolic acid constituents of mistletoe by adopting structural bioinformatics and advanced theoretical chemistry techniques via molecular docking, molecular dynamics simulation, molecular mechanics/generalized Born surface area (MM/GBSA) calculations, pharmacokinetic analysis, and density functional theory analysis. The six best compounds from the flavonoid constituent with the highest binding affinity ranging from −6.8 kcal/mol to −6.7 kcal/mol were selected with the control gemcitabine (−5.5 kcal/mol) for further computational analysis after molecular docking. Furthermore, MM/GBSA calculation showed the highest binding energy for the selected docked compounds, which validates their inhibitory potential. Hence, the molecular dynamics simulation, post-simulation analysis, pharmacokinetics model, and DFT results showed that mistletoe compounds are reliable due to their stable interaction with the target protein and drug-likeness properties. Communicated by Ramaswamy H. Sarma

Structure-based computational screening of 470 natural quercetin derivatives for identification of SARS-CoV-2 Mpro inhibitor.

Coronavirus disease 2019 (COVID-19) is a global pandemic infecting the respiratory system through a notorious virus known as the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Due to viral mutations and the risk of drug resistance, it is crucial to identify new molecules having potential prophylactic or therapeutic effect against SARS-CoV-2 infection. In the present study, we aimed to identify a potential inhibitor of SARS-CoV-2 through virtual screening of a compound library of 470 quercetin derivatives by targeting the main protease-Mpro (PDB ID: 6LU7). The study was carried out with computational techniques such as molecular docking simulation studies (MDSS), molecular dynamics (MD) simulations, and molecular mechanics generalized Born surface area (MMGBSA) techniques. Among the natural derivatives, compound 382 (PubChem CID 65604) showed the best binding affinity to Mpro (-11.1 kcal/mol). Compound 382 interacted with LYS5, TYR126, GLN127, LYS137, ASP289, PHE291, ARG131, SER139, GLU288, and GLU290 of the Mpro protein. The SARS-CoV-2 Mpro-382 complex showed acceptable stability during the 100 ns MD simulations. The SARS-CoV-2 Mpro-382 complex also showed an MM-GBSA binding free energy value of -54.0 kcal/mol. The binding affinity, stability, and free energy results for 382 and Mpro were better than those of the native ligand and the standard inhibitors ledipasvir and cobicistat. The conclusion of our study was that compound 382 has the potential to inhibit SARS-Cov-2 Mpro. However, further investigations such as in-vitro assays are recommended to confirm its in-silico potency.

Structure-based small inhibitors search combined with molecular dynamics driven energies for human programmed cell death-1 (PD-1) protein

Human immune system is specialized in distinguishing normal cells from foreign particles mainly through proteins expressed on immune cells called ‘checkpoints’. Immune checkpoints work as a switch to activate and deactivate immune responses. T cells express one of the immune checkpoint, human programmed cell death-1 (PD-1), which normally operates as an off-switch function to protect the normal cell from T-cell attack. Binding of PD-1 to its ligand, the programmed cell death ligand (PD-L1/2) expressed on myeloid/cancer cells, induce downstream inhibitory signals, leading to tumor immune evasion. Targeting PD-1 or PD-L1 can boost the immune response against cancer cells. To design novel small molecule inhibitors for the PD-1, in silico structure-based screening on pharmacophoric points and molecular docking were performed. Based on the docking score and significant binding interaction with the crucial residues of PD-1 (Thr59, Glu61, Ser62, Glu84, Arg86 and Ala132), compounds were selected from the ZINC20 database, and their dynamic behavior and conformational stability were examined through molecular dynamic simulations. Besides, the Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) method was used to calculate the binding strength of each selected inhibitor complexed with PD-1. The binding energy calculations revealed that these selected inhibitors show a considerable affinity for PD-1. The selected novel inhibitors exhibit excellent drug-like and pharmacokinetic properties (absorption, distribution, metabolism, excretion and toxicity). In conclusion, the identified novel compounds (ZINC1443480030, ZINC1002854123, ZINC988238128, ZINC1481242350, ZINC1001739421, ZINC1220816434 and ZINC1167786692) from the current study can be validated in-vitro as potential PD-1 inhibitors and for discovery of novel drugs against PD-1 in the future. Communicated by Ramaswamy H. Sarma

In vitro and in silico evaluations of actinomycin X2and actinomycin D as potent anti-tuberculosis agents.

Multidrug-resistant tuberculosis (MDR-TB) is one of the world's most devastating contagious diseases and is caused by the MDR-Mycobacterium tuberculosis (MDR-Mtb) bacteria. It is therefore essential to identify novel anti-TB drug candidates and target proteins to treat MDR-TB. Here, in vitro and in silico studies were used to investigate the anti-TB potential of two newly sourced actinomycins, actinomycin-X2 (act-X2) and actinomycin-D (act-D), from the Streptomyces smyrnaeus strain UKAQ_23 (isolated from the Jubail industrial city of Saudi Arabia). The anti-TB activity of the isolated actinomycins was assessed in vitro using the Mtb H37Ra, Mycobacterium bovis (BCG), and Mtb H37Rv bacterial strains, using the Microplate Alamar Blue Assay (MABA) method. In silico molecular docking studies were conducted using sixteen anti-TB drug target proteins using the AutoDock Vina 1.1.2 tool. The molecular dynamics (MD) simulations for both actinomycins were then performed with the most suitable target proteins, using the GROningen MAchine For Chemical Simulations (GROMACS) simulation software (GROMACS 2020.4), with the Chemistry at HARvard Macromolecular Mechanics 36m (CHARMM36m) forcefield for proteins and the CHARMM General Force Field (CGenFF) for ligands. In vitro results for the Mtb H37Ra, BCG, and Mtb H37Rv strains showed that act-X2 had minimum inhibitory concentration (MIC) values of 1.56 ±0.0, 1.56 ±0.0, and 2.64 ±0.07 µg/mL and act-D had MIC values of 1.56 ±0.0, 1.56 ±0.0, and 1.80 ±0.24 µg/mL respectively. The in silico molecular docking results showed that protein kinase PknB was the preferred target for both actinomycins, while KasA and pantothenate synthetase were the least preferred targets for act-X2and act-D respectively. The molecular dynamics (MD) results demonstrated that act-X2 and act-D remained stable inside the binding region of PknB throughout the simulation period. The MM/GBSA (Molecular Mechanics/Generalized Born Surface Area) binding energy calculations showed that act-X2 was more potent than act-D. In conclusion, our results suggest that both actinomycins X2 and D are highly potent anti-TB drug candidates. We show that act-X2is better able to antagonistically interact with the protein kinase PknB target than act-D, and thus has more potential as a new anti-TB drug candidate.

Computational identification of drug-like marine natural products as potential RNA polymerase inhibitors against Nipah virus

Nipah virus (NiV) has been an alarming threat to human populations in southern Asia for more than a decade. It is one of the most deadly viruses in the Mononegavirales order. Despite its high mortality rate and virulence, no chemotherapeutic agent or vaccine is publicly available. Hence, this work was conducted to computationally screen marine natural products database for drug-like potential inhibitors for the viral RNA-dependent RNA polymerase (RdRp). The structural model was subjected to molecular dynamics (MD) simulation to obtain the native ensemble of the protein. The CMNPDB dataset of marine natural products was filtered to retain only compounds following Lipinski's five rules. The molecules were energy minimized and docked into different conformers of the RdRp using AutoDock Vina. The best 35 molecules were rescored by GNINA, a deep learning-based docking software. The resulting nine compounds were evaluated for their pharmacokinetic profiles and medicinal chemistry properties. The best five compounds were subjected to MD simulation for 100 ns, followed by binding free energy estimation via Molecular Mechanics/ Generalized Born Surface Area (MM/GBSA) calculations. The results showed remarkable behavior of five hits as inferred by stable binding pose and orientation to block the exit channel of RNA synthesis products in the RdRp cavity. These hits are promising starting materials for in vitro validation and structural modifications to enhance the pharmacokinetic and medicinal chemistry properties for developing antiviral lead compounds.

Computational and experimental studies of salvianolic acid A targets 3C protease to inhibit enterovirus 71 infection.

Hand, foot, and mouth disease (HFMD) is a common childhood infectious disease caused by enterovirus (EV) infection. EV71 is one of the major pathogens causing hand, foot, and mouth disease and is more likely to cause exacerbation and death than other enteroviruses. Although a monovalent vaccine for EV71 has been developed, there are no clinically available anti-EV71 specific drugs. Here, we performed virtual screening and biological experiments based on the traditional Chinese medicine monomer library. We identified a traditional Chinese medicine monomer, Salvianolic acid A (SA), a polyphenolic compound isolated from Salvia miltiorrhiza. Salvianolic acid A inhibits EV71 virus infection in a concentration-dependent manner, and its antiviral activity is higher than that of other reported natural polyphenols and has a high biosafety. Furthermore, molecular dynamics simulations showed that salvianolic acid A can anchor to E71, a member of the enzyme catalytic triad, and cause H40 to move away from the catalytic center. Meanwhile, molecular mechanics generalized born surface area (MMGBSA) and steered molecular dynamics (SMD) results showed that the P1 group of SA was most easily unbound to the S1 pocket of 3Cpro, which provided theoretical support to further improve the affinity of salvianolic acid A with 3Cpro. These findings suggest that salvianolic acid A is a novel EV71 3Cpro inhibitor with excellent antiviral activity and is a promising candidate for clinical studies.

South African (501Y.V2) and the United Kingdom (B.1.1.7) SARS-CoV-2 Spike (S) Protein Variants Demonstrate a Higher Binding Affinity to ACE2.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) affects the lower respiratory tract by binding to angiotensin-converting enzyme 2 (ACE2) via its S-protein. Recent emerging SARS-CoV-2 variants from the United Kingdom (B.1.1.7) and South Africa (501Y.V2) are spreading worldwide at an alarming rate. The new variants have manifested amino acid substitution K417N, E484K, and N501Y on the RBD domain that binds to ACE2. These mutations may influence the binding of the S-protein to ACE2 and affect viral entry into the host cell. In this study, we modelled the amino acid substitutions on the S-protein and utilised the HADDOCK server to assess the S-protein RBD domain binding with ACE2. Additionally, we calculated the binding affinity of ACE2 to S-protein WT, B.1.1.7 and 501Y.V2 variants using Molecular Mechanics-Generalized Born Surface Area (MM/GBSA). We demonstrated that the S-protein of both variants possesses a higher binding affinity to ACE2 than WT, with the South African 501Y.V2 being more infective than the B.1.1.7 that originated in the United Kingdom. The South African 501Y.V2 variant presents three amino acid substitutions that changed the H-bonding network, resulting in a higher affinity to ACE2, indicating that the 501Y.V2 strain is more infective than the B.1.1.7 strain.