Identification of FDA approved drugs and nucleoside analogues as potential SARS-CoV-2 A1pp domain inhibitor: An in silico study

Machine learning and classical MD simulation to identify inhibitors against the P37 envelope protein of monkeypox virus

Monkeypox virus (MPXV) outbreak is a serious public health concern that requires international attention. P37 of MPXV plays a pivotal role in DNA replication and acts as one of the promising targets for antiviral drug design. In this study, we intent to screen potential analogs of existing FDA approved drugs of MPXV against P37 using state-of-the-art machine learning and computational biophysical techniques. AlphaFold2 guided all-atoms molecular dynamics simulations optimized P37 structure is used for molecular docking and binding free energy calculations. Similar to members of Phospholipase-D family , the predicted P37 structure also adopts a β-α-β-α-β sandwich fold, harbouring strongly conserved HxKxxxxD motif. The binding pocket comprises of Tyr48, Lys86, His115, Lys117, Ser130, Asn132, Trp280, Asn240, His325, Lys327 and Tyr346 forming strong hydrogen bonds and dense hydrophobic contacts with the screened analogs and is surrounded by positively charged patches. Loops connecting the two domains and C-terminal region exhibit high degree of flexibility. In some structural ensembles, the partial disorderness in the C-terminal region is presumed to be due to its low confidence score, acquired during structure prediction. Transition from loop to β-strands (244-254 aa) in P37-Cidofovir and its analog complexes advocates the need for further investigations. MD simulations support the accuracy of the molecular docking results, indicating the potential of analogs as potent binders of P37. Taken together, our results provide preferable understanding of molecular recognition and dynamics of ligand-bound states of P37, offering opportunities for development of new antivirals against MPXV. However, the need of in vitro and in vivo assays for confirmation of these results still persists.Communicated by Ramaswamy H. Sarma.

Predicting FFAR4 agonists using structure-based machine learning approach based on molecular fingerprints

Free Fatty Acid Receptor 4 (FFAR4), a G-protein-coupled receptor, is responsible for triggering intracellular signaling pathways that regulate various physiological processes. FFAR4 agonists are associated with enhancing insulin release and mitigating the atherogenic, obesogenic, pro-carcinogenic, and pro-diabetogenic effects, normally associated with the free fatty acids bound to FFAR4. In this research, molecular structure-based machine-learning techniques were employed to evaluate compounds as potential agonists for FFAR4. Molecular structures were encoded into bit arrays, serving as molecular fingerprints, which were subsequently analyzed using the Bayesian network algorithm to identify patterns for screening the data. The shortlisted hits obtained via machine learning protocols were further validated by Molecular Docking and via ADME and Toxicity predictions. The shortlisted compounds were then subjected to MD Simulations of the membrane-bound FFAR4-ligand complexes for 100 ns each. Molecular analyses, encompassing binding interactions, RMSD, RMSF, RoG, PCA, and FEL, were conducted to scrutinize the protein-ligand complexes at the inter-atomic level. The analyses revealed significant interactions of the shortlisted compounds with the crucial residues of FFAR4 previously documented. FFAR4 as part of the complexes demonstrated consistent RMSDs, ranging from 3.57 to 3.64, with minimal residue fluctuations 5.27 to 6.03 nm, suggesting stable complexes. The gyration values fluctuated between 22.8 to 23.5 nm, indicating structural compactness and orderliness across the studied systems. Additionally, distinct conformational motions were observed in each complex, with energy contours shifting to broader energy basins throughout the simulation, suggesting thermodynamically stable protein-ligand complexes. The two compounds CHEMBL2012662 and CHEMBL64616 are presented as potential FFAR4 agonists, based on these insights and in-depth analyses. Collectively, these findings advance our comprehension of FFAR4's functions and mechanisms, highlighting these compounds as potential FFAR4 agonists worthy of further exploration as innovative treatments for metabolic and immune-related conditions.

In silico investigation of potential phytoconstituents against ligand- and voltage-gated ion channels as antiepileptic agents

The most promising anticonvulsant phytocompounds were explored in this work using docking, molecular dynamic (MD) simulation, and Molecular Mechanics-Poisson-Boltzmann Surface Area (MM-PBSA) approaches. A total of 70 phytochemicals were screened against α-amino-3-hydroxyl-5-methyl-4-isoxazole propionic acid (AMPA), N-methyl-d-aspartate (NMDA), voltage-gated sodium ion channels (VGSC), and carbonic anhydrase enzyme II (CA II) receptors, and the docking results were compared to the reference drug phenytoin. Amentoflavone displayed the highest affinity for AMPA and VGSC receptors, with docking scores of - 10.4 and - 10.1 kcal/mol, respectively. Oliganthin H-NMDA and epigallocatechin-3-gallate-CA II complexes showed docking scores of - 10.9 and - 6.9 kcal/mol, respectively. All four complexes depicted a high dock score compared to the phenytoin complex at the binding site of the corresponding proteins. The MD simulation investigated the stabilities and favorable conformation of apoproteins and ligand/reference-bound complexes. The results revealed that proteins AMPA, VGSC, and CA II were more efficiently stabilized by lead phytochemicals than phenytoin binding. Additionally, principal component analysis and MM-PBSA results suggested that these lead phytocompounds have good compactness and strong binding free energy. Further, physicochemical and pharmacokinetic studies revealed that these final lead phytochemicals would be suitable for oral intake, have sufficient intestinal permeability, and have the ability to cross the blood-brain barrier (BBB). Comprehensively, this study predicted amentoflavone as the best lead phytochemical out of the 70 anticonvulsant phytocompounds that can be used to treat epilepsy.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-03948-1.

Structural and energetic analysis of NS5 protein inhibition by small molecules in Japanese encephalitis virus using machine learning and steered molecular dynamics approach

One of the viral diseases that affect millions of people around the world, particularly in developing countries, is Japanese encephalitis (JE). In this study, the conserved protein of this virus, that is, non-structural protein 5 (NS5), was used as a target protein for this study, and a compound library of 749 antiviral molecules was screened against NS5. The current study employed machine learning-based virtual screening combined with molecular docking. Here, three hits (24360, 123519051 and 213039) had lower binding energies (< -8 kcal/mol) than the control, S-Adenosyl-L-homocysteine (SAH). All the compounds showed significant H-bond interactions with functional residues, which were also observed by the control. Molecular dynamics simulation, MM/GBSA for binding free energy analysis, principal component analysis and free energy landscape were also performed to study the stability of the complex formation. All three compounds had similar root mean square deviation trends, which were comparable to the control, SAH. Post-MD, the 123519051-receptor complex had the highest number of H-bonds (4 to 5) after the control, out of which three exhibited the highest percentage occupancy (50%, 24% and 79%). Both docking and MD, 123519051 showed an H-bond with the residue Gly111, which was also found for the control-protein complex. 123519051 showed the lowest binding free energy with ΔGbind of -89 kJ/mol. Steered molecular dynamics depicted that 123519051 had the maximum magnitude of dissociation (1436.43 kJ/mol/nm), which was more than the control, validating its stable complex formation. This study concluded that 123519051 is a binder and could inhibit the protein NS5 of JE.Communicated by Ramaswamy H. Sarma.

Computational exploration of microsomal cytochrome P450 3A1 enzyme modulation by phytochemicals of Cichorium intybus L.: Insights into drug metabolism

Drug metabolism plays a crucial role in drug fate, including therapeutic inactivation or activation, as well as the formation of toxic compounds. This underscores the importance of understanding drug metabolism in drug discovery and development. Considering the substantial costs associated with traditional drug development methods, computational approaches have emerged as valuable tools for predicting the metabolic fate of drug candidates. With this in mind, the present study aimed to investigate the potential mechanisms underlying the modulation of microsomal cytochrome P450 3A1 (CYP3A1) enzyme activity by various phytochemicals found in Cichorium intybus L., commonly known as chicory. To achieve this goal, several in silico methods, including molecular docking and molecular dynamics (MD) simulation, were employed to explore computationally the microsomal CYP3A1 enzyme. Schrodinger software was utilized for the molecular docking study, which involved the interaction analysis between CYP3A1 and 28 phytoconstituents of Cichorium intybus. Virtual screening of 28 compounds from chicory led to the identification of the top five ranked compounds. These compounds were evaluated for drug-likeness properties, pharmacokinetic profiles, and predicted binding affinities to CYP3A1. Caffeoylshikimic acid and cichoric acid emerged as promising candidates due to their favorable characteristics, including good oral bioavailability and high binding affinities to CYP3A1. Molecular dynamics simulations were conducted to assess the stability of caffeoylshikimic acid within the CYP3A1 binding pocket. The results demonstrated that caffeoylshikimic acid maintained stable interactions with the enzyme throughout the simulation, suggesting its potential as an effective modulator of CYP3A1 activity. The findings of this study have the potential to provide valuable insights into the complex molecular mechanisms by which Cichorium intybus L. acts on hepatocytes and modulates CYP3A1 enzyme expression or activity. By elucidating the impact of these phytochemicals on drug metabolism, this research contributes to our understanding of how chicory may interact with drugs and influence their efficacy and safety profiles.

Molecular docking and MD simulations reveal protease inhibitors block the catalytic residues in Prp8 intein of Aspergillus fumigatus: a potential target for antimycotics

Resistance to azoles and amphotericin B especially in Aspergillus fumigatus is a growing concern towards the treatment of invasive fungal infection. At this critical juncture, intein splicing would be a productive, and innovative target to establish therapies against resistant strains. Intein splicing is the central event for the activation of host protein, essential for the growth and survival of various microorganisms including A. fumigatus. The splicing process is a four-step protease-like nucleophilic cascade. Thus, we hypothesise that protease inhibitors would successfully halt intein splicing and potentially restrict the growth of the aforementioned pathogen. Using Rosetta Fold and molecular dynamics simulations, we modelled Prp8 intein structure; resembling classic intein fold with horse shoe shaped splicing domain. To fully comprehend the active site of Afu Prp8 intein, C1, T62, H65, H818, N819 from intein sequences and S820, the first C-extein residue are selected. Molecular docking shows that two FDA-approved drugs, i.e. Lufotrelvir and Remdesivir triphosphate efficiently interact with Prp8 intein from the assortment of 212 protease inhibitors. MD simulation portrayed that Prp8 undergoes conformational change upon ligand binding, and inferred the molecular recognition and stability of the docked complexes. Per-residue decomposition analysis confirms the importance of F: block R802, V803, and Q807 binding pocket in intein splicing domain towards recognition of inhibitors, along with active site residues through strong hydrogen bonds and hydrophobic contacts. However, in vitro and in vivo assays are required to confirm the inhibitory action on Prp8 intein splicing; which may pave the way for the development of new antifungals for A. fumigatus.Communicated by Ramaswamy H. Sarma.

Exploration of phytochemical compounds against Marburg virus using QSAR, molecular dynamics, and free energy landscape

Marburg virus disease (MVD) is caused by the Marburg virus, a one-of-a-kind zoonotic RNA virus from the genus Filovirus. Thus, this current study employed AI-based QSAR and molecular docking-based virtual screening for identifying potential binders against the target protein (nucleoprotein (NP)) of the Marburg virus. A total of 2727 phytochemicals were used for screening, out of which the top three compounds (74977521, 90470472, and 11953909) were identified based on their predicted bioactivity (pIC50) and binding score (< - 7.4 kcal/mol). Later, MD simulation in triplicates and trajectory analysis were performed which showed that 11953909 and 74977521 had the most stable and consistent complex formations and had the most significant interactions with the highest number of hydrogen bonds. PCA (principal component analysis) and FEL (free energy landscape) analysis indicated that these compounds had favourable energy states for most of the conformations. The total binding free energy of the compounds using the MM/GBSA technique showed that 11953909 (ΔGTOTAL = - 30.78 kcal/mol) and 74977521 (ΔGTOTAL = - 30 kcal/mol) had the highest binding affinity with the protein. Overall, this in silico pipeline proposed that the phytochemicals 11953909 and 74977521 could be the possible binders of NP. This study aimed to find phytochemicals inhibiting the protein's function and potentially treating MVD.

Bioinformatics insights into CENP-T and CENP-W protein-protein interaction disruptive amino acid substitution in the CENP-T-W complex

Kinetochores are multi-protein assemblies present at the centromere of the human chromosome and play a crucial role in cellular mitosis. The CENP-T and CENP-W chains form a heterodimer, which is an integral part of the inner kinetochore, interacting with the linker DNA on one side and the outer kinetochore on the other. Additionally, the CENP-T-W dimer interacts with other regulatory proteins involved in forming inner kinetochores. The specific roles of different amino acids in the CENP-W at the protein-protein interaction (PPI) interface during the CENP-T-W dimer formation remain incompletely understood. Since cell division goes awry in diseases like cancer, this CENP-T-W partnership is a potential target for new drugs that could restore healthy cell division. We employed molecular docking, binding free energy calculations, and molecular dynamics (MD) simulations to investigate the disruptive effects of amino acids substitutions in the CENP-W chain on CENP-T-W dimer formation. By conducting a molecular docking study and analysing hydrogen bonding interactions, we identified key residues in CENP-W (ASN-46, ARG-53, LEU-83, SER-86, ARG-87, and GLY-88) for further investigation. Through site-directed mutagenesis and subsequent binding free energy calculations, we refined the selection of mutant. We chose four mutants (N46K, R53K, L83K, and R87E) of CENP-W to assess their comparative potential in forming CENP-T-W dimer. Our analysis from 250 ns long revealed that the substitution of LEU83 and ARG53 residues in CENP-W with the LYS significantly disrupts the formation of CENP-T-W dimer. In conclusion, LEU83 and ARG53 play a critical role in CENP-T and CENP-W dimerization which is ultimately required for cellular mitosis. Our findings not only deepen our understanding of cell division but also hint at exciting drug-target possibilities.

An Update on the Current State of SARS-CoV-2 Mac1 Inhibitors

Non-structural protein 3 (nsp3) from all coronaviruses (CoVs) contains a conserved macrodomain, known as Mac1, that has been proposed as a potential therapeutic target for CoVs due to its critical role in viral pathogenesis. Mac1 is an ADP-ribose binding protein and ADP-ribosylhydrolase that promotes replication and blocks IFN responses, though the precise mechanisms it uses to carry out these functions remain unknown. Over the past 3 years following the onset of COVID-19, several groups have used high-throughput screening with multiple assays and chemical modifications to create unique chemical inhibitors of the SARS-CoV-2 Mac1 protein. Here, we summarize the current efforts to identify selective and potent inhibitors of SARS-CoV-2 Mac1.

Hesperidin potentially interacts with the catalytic site of gamma-secretase and modifies notch sensitive genes and cancer stemness marker expression in colon cancer cells and colonosphere

Gamma secretase (GS) produces Notch Intracellular Domain (NICD) by trans-membrane cleavage of notch receptor. The NICD enters the nucleus and activates the notch signaling pathway (NSP) by activating notch-responsive gene transcription. Hyperactivation of NSP is related to cancer aggressiveness, therapy resistance, and poor therapy outcome, and decreased overall disease-free survival in patients. Till date, none of the GS inhibitors (GSI) has been clinically approved due to their toxicity in patients. Thus in the present study, we explored the GS catalytic site binding potential of hesperidin (natural flavone glycoside) and its effect on notch responsive gene expression in HCT-116 cells. Molecular docking, MM-GBSA binding energy calculations, and molecular dynamics (MD) simulation experiments were performed to study the GS catalytic site binding potential of hesperidin. The compound showed better GS catalytic site binding potential at the active site compared to experimentally validated GSI, N-N-(3, 5-Difluorophenacetyl)-L-alanyl-S-phenylglycine t-butyl ester (DAPT) in molecular docking and MM-GBSA experiments. MD simulation results showed that hesperidin forms stable and energetically favorable complex with gamma secretase in comparison to standard inhibitor (DAPT)-GS complex. Further, in vitro experiments showed that hesperidin inhibited cell growth and sphere formation potential in HCT-116 cells. Further, hesperidin treatment altered notch responsive genes (Hes1, Hey1, and E-cad) and cancer stemness/self-renewal markers expression at transcription levels. In conclusion, hesperidin produces toxicity in HCT-116 cells and decreases colonosphere formation by inhibiting transcription of notch signaling pathway target genes and stemness markers.Communicated by Ramaswamy H. Sarma.

Exploiting the potential of natural polyphenols as antivirals against monkeypox envelope protein F13 using machine learning and all-atoms MD simulations

The re-emergence of monkeypox (MPX), in the era of COVID-19 pandemic is a new global menace. Regardless of its leniency, there are chances of MPX expediting severe health deterioration. The role of envelope protein, F13 as a critical component for production of extracellular viral particles makes it a crucial drug target. Polyphenols, exhibiting antiviral properties have been acclaimed as an effective alternative to the traditional treatment methods for management of viral diseases. To facilitate the development of potent MPX specific therapeutics, herein, we have employed state-of-the-art machine learning techniques to predict a highly accurate 3-dimensional structure of F13 as well as identify binding hotspots on the protein surface. Additionally, we have effectuated high-throughput virtual screening methodology on 57 potent natural polyphenols having antiviral activities followed by all-atoms molecular dynamics (MD) simulations, to substantiate the mode of interaction of F13 protein and polyphenol complexes. The structure-based virtual screening based on Glide SP, XP and MM/GBSA scores enables the selection of six potent polyphenols having higher binding affinity towards F13. Non-bonded contact analysis, of pre- and post- MD complexes propound the critical role of Glu143, Asp134, Asn345, Ser321 and Tyr320 residues in polyphenol recognition, which is well supported by per-residue decomposition analysis. Close-observation of the structural ensembles from MD suggests that the binding groove of F13 is mostly hydrophobic in nature. Taken together, this structure-based analysis from our study provides a lead on Myricetin, and Demethoxycurcumin, which may act as potent inhibitors of F13. In conclusion, our study provides new insights into the molecular recognition and dynamics of F13-polyphenol bound states, offering new promises for development of antivirals to combat monkeypox. However, further in vitro and in vivo experiments are necessary to validate these results.

Evaluation of phytoconstituents of Tinospora cordifolia against K417N and N501Y mutant spike glycoprotein and main protease of SARS-CoV-2- an in silico study

Coronavirus disease 2019 (COVID-19) caused appalling conditions over the globe, which is currently faced by the entire human population. One of the primary reasons behind the uncontrollable situation is the lack of specific therapeutics. In such conditions, drug repurposing of available drugs (viz. Chloroquine, Lopinavir, etc.) has been proposed, but various clinical and preclinical investigations indicated the toxicity and adverse side effects of these drugs. This study explores the inhibition potency of phytochemicals from Tinospora cordifolia (Giloy) against SARS CoV-2 drugable targets (spike glycoprotein and Mpro proteins) using molecular docking and MD simulation studies. ADMET, virtual screening, MD simulation, postsimulation analysis (RMSD, RMSF, Rg, SASA, PCA, FES) and MM-PBSA calculations were carried out to predict the inhibition efficacy of the phytochemicals against SARS CoV-2 targets. Tinospora compounds showed better binding affinity than the corresponding reference. Their binding affinity ranges from -9.63 to -5.68 kcal/mole with spike protein and -10.27 to -7.25 kcal/mole with main protease. Further 100 ns exhaustive simulation studies and MM-PBSA calculations supported favorable and stable binding of them. This work identifies Nine Tinospora compounds as potential inhibitors. Among those, 7-desacetoxy-6,7-dehydrogedunin was found to inhibit both spike (7NEG) and Mpro (7MGS and 6LU7) proteins, and Columbin was found to inhibit selected spike targets (7NEG and 7NX7). In all the analyses, these compounds performed well and confirms the stable binding. Hence the identified compounds, advocated as potential inhibitors can be taken for further in vitro and in vivo experimental validation to determine their anti-SARS-CoV-2 potential.Communicated by Ramaswamy H. Sarma.

Multifaceted 3D-QSAR analysis for the identification of pharmacophoric features of biphenyl analogues as aromatase inhibitors

Aromatase, a cytochrome P450 enzyme, is responsible for the conversion of androgens to estrogens, which fuel the multiplication of cancerous cells. Inhibition of estrogen biosynthesis by aromatase inhibitors (AIs) is one of the highly advanced therapeutic approach available for the treatment of estrogen-positive breast cancer. Biphenyl moiety aids lipophilicity to the conjugated scaffold and enhances the accessibility of the ligand to the target. The present study is focused on the investigation of, the mode of binding of biphenyl with aromatase, prediction of ligand-target binding affinities, and pharmacophoric features essential for favorable for aromatase inhibition. A multifaceted 3D-QSAR (SOMFA, Field and Gaussian) along with molecular docking, molecular dynamic simulations and pharmacophore mapping were performed on a series of biphenyl bearing molecules (1-33) with a wide range of aromatase inhibitory activity (0.15-920 nM). Among the generated 3D-QSAR models, the Force field-based 3D-QSAR model (R² = 0.9151) was best as compared to SOMFA and Gaussian Field (R²=0.7706, 0.9074, respectively). However, all the generated 3D-QSAR models were statistically fit, robust enough, and reliable to explain the variation in biological activity in relation to pharmacophoric features of dataset molecules. A four-point pharmacophoric features with three acceptor sites (A), one aromatic ring (R) features, AAAR_1, were obtained with the site and survival score values 0.890 and 4.613, respectively. The generated 3D-QSAR plots in the study insight into the structure-activity relationship of dataset molecules, which may help in the designing of potent biphenyl derivatives as newer inhibitors of aromatase.Communicated by Ramaswamy H. Sarma.

Design, synthesis, and evaluation of 4'-phosphonomethoxy pyrimidine ribonucleosides as potential anti-influenza agents

Influenza viruses belong to the Orthomyxoviridae family and cause acute respiratory distress in humans. The developed drug resistance toward existing drugs and the emergence of viral mutants that can escape vaccines mandate the search for novel antiviral drugs. Herein, the synthesis of epimeric 4'-methyl-4'-phosphonomethoxy [4'-C-Me-4'-C-(O-CH₂ P═O)] pyrimidine ribonucleosides, their phosphonothioate [4'-C-Me-4'-C-(O-CH₂ P═S)] derivatives, and their evaluation against an RNA viral panel are described. Selective formation of the α- l-lyxo epimer, [4'-C-(α)-Me-4'-C-(β)-(O-CH₂ -P(═O)(OEt)₂ )] over the β- d-ribo epimer [4'-C-(β)-Me-4'-C-(α)-(O-CH₂ -P(═O)(OEt)₂ )] was explained by DFT equilibrium geometry optimizations studies. Pyrimidine nucleosides having the [4'-C-(α)-Me-4'-C-(β)-(O-CH₂ -P(═O)(OEt)₂ )] framework showed specific activity against influenza A virus. Significant anti-influenza virus A (H1N1 California/07/2009 isolate) was observed with the 4'-C-(α)-Me-4'-C-(β)-O-CH₂ -P(═O)(OEt)₂ -uridine derivative 1 (EC₅₀  = 4.56 mM, SI₅₀  > 56), 4-ethoxy-2-oxo-1(2H)-pyrimidin-1-yl derivative 3 (EC₅₀  = 5.44 mM, SI₅₀  > 43) and the cytidine derivative 2 (EC₅₀  = 0.81 mM, SI₅₀  > 13), respectively. The corresponding thiophosphonates 4'-C-(α)-Me-4'-C-(β)-(O-CH₂ -P( S)(OEt)₂ ) and thionopyrimidine nucleosides were devoid of any antiviral activity. This study shows that the 4'-C-(α)-Me-4'-(β)-O-CH₂ -P(═O)(OEt)₂ ribonucleoside can be further optimized to provide potent antiviral agents.

Unlocking SGK1 inhibitor potential of bis-[1-N,7-N, pyrazolo tetraethoxyphthalimido{-4-(3,5-Dimethyl-4-(spiro-3-methylpyazolo)-1,7-dihydro-1H-dipyrazolo[3,4-b;4',3'-e]pyridin-8-yl)}]p-disubstituted phenyl compounds: a computational study

SGK1 (Serum and Glucocorticoid Regulated Kinase 1), a serine/threonine kinase that is activated by various stimuli, including serum and glucocorticoids. It controls inflammation, apoptosis, hormone release, neuro-excitability and cell proliferation, all of which play an important role in cancer progression and metastasis. SGK1 was recently proposed as a potential drug target for cancer, diabetes, and neurodegenerative diseases. In this study, molecular docking, physiochemical, toxicological properties and molecular dynamic simulation of the Bis-[1-N,7-N, Pyrazolo tetraethoxyphthalimido{-4-(3,5-Dimethyl-4-(spiro-3-methylpyazolo)-1,7-dihydro-1H-dipyrazolo[3,4-b;4',3'-e]pyridin-8-yl)}]p-disubstituted phenyl compoundsand reference EMD638683 against new SGK1 target protein. Compared to the reference inhibitor EMD638683, we choose the best compounds (series 2-6) based on the binding energy (in the range from -11.0 to -10.6 kcal/mol). With the exception of compounds 2 and 6, none of the compounds posed a risk for AMES toxicity or carcinogenicity due to their toxicological properties. 100 ns MD simulation accompanied by MM/PBSA energy calculations and PCA. According to MD simulation results, the binding of compounds 3, 4 and 5 stabilizes the SGK1 structure and causes febrile conformational changes compared to EMD638683. As a result of this research, the final selected compounds 3, 4 and 5 can be used as scaffolds to develop promising SGK1 inhibitors for the treatment of related diseases such as cancer.Communicated by Ramaswamy H. Sarma.

Identification of the effective α-amylase inhibitors from Dalbergia odorifera: Virtual screening, spectroscopy, molecular docking, and molecular dynamic simulation

Inhibiting the activity of α-amylase has been considered as one efficient way to prevent and treat type 2 diabetes recently. Dalbergia odorifera, a kind of Leguminosae plant, has a positive therapeutic effect on type 2 diabetes, possibly contributing by some constituents that can inhibit the activity of α-amylase. In this study, we found that eriodictyol was one potential constituent through virtual screening. The interaction mode between eriodictyol and α-amylase was elucidated by molecular docking, multi-spectroscopic analysis, and molecular dynamic simulation. The results revealed that eriodictyol quenched the intrinsic fluorescence of α-amylase, and the quenching mode was static quenching. Eriodictyol could spontaneously interact with α-amylase, mostly stabilized and influenced by the hydrophobic interaction, while the binding sites (n) was 1.13 ± 0.07 and binding constant (K_b) was (1.43 ± 0.14) × 10⁵ at 310 K, respectively. In addition, FT-IR and CD had been applied to identify that eriodictyol can trigger the conformational change of α-amylase. Taken together, the results provided some experimental data for developing new α-amylase inhibitors from Dalbergia odorifera, which may further prevent and treat diabetes and diabetes complications.

In Silico Identification of New Anti-SARS-CoV-2 Main Protease (Mpro) Molecules with Pharmacokinetic Properties from Natural Sources Using Molecular Dynamics (MD) Simulations and Hierarchical Virtual Screening

Infectious agents such as SARS-CoV, MERS-CoV, and SARS-CoV-2 have emerged in recent years causing epidemics with high mortality rates. The quick development of novel therapeutic compounds is required in the fight against such pathogenic agents. Unfortunately, the traditional drug development methods are time-consuming and expensive. In this study, computational algorithms were utilized for virtual screening of a library of natural compounds in the ZINC database for their affinity towards SARS-CoV-2 Mpro. Compounds such as cinanserin, nelfinavir, baicalin, baicalein, candesartan cilexetil, chloroquine, dipyridamole, and hydroxychloroquine have the ability to prevent SARS-CoV-2 Mpro from facilitating COVID 19 infection; thus, they treat COVID 19. However, these drugs majorly act to reduce the symptoms of the disease. No anti-viral drug against COVID 19 virus infection has been discovered and approved. Therefore, this study sought to explore natural inhibitors of SARS-CoV-2 Mpro to develop a pharmacophore model for virtual screening of natural compounds in the ZINC database as potential candidates for SARS-CoV-2 Mpro inhibitors and as therapeutic molecules against COVID 19. This study undertook in silico methods to identify the best anti-viral candidates targeting SAR-CoV-2 Mpro from natural sources in the ZINC database. Initially, reported anti-SARS-CoV-2 Mpro molecules were integrated into designing a pharmacophore model utilizing PharmaGist. Later, the pharmacophore model was loaded into ZINCPHARMER and screened against the ZINC database to identify new probable drug candidates. The root means square deviation (RMSD) values of the potential drug candidates informed the selection of some of them, which were docked with SARS-CoV-2 Mpro to comprehend their interactions. From the molecular docking results, the top four candidates (ZINC000254823011, ZINC000072307130, ZINC000013627512, and ZINC000009418994) against SARS-CoV-2 Mpro, with binding energies ranging from –8.2 kcal/mol to –8.6 kcal/mol, were examined for their oral bioavailability and other pharmacokinetic properties. Consequently, ZINC000072307130 emerged as the only orally bioavailable drug candidate with desirable pharmacokinetic properties. This candidate drug was used to perform MD simulations, and the outcomes revealed that ZINC000072307130 formed a stable complex with the viral main protease. Consequently, ZINC000072307130 emerges as a potential anti-SARS-CoV-2 Mpro inhibitor for the production of new COVID 19 drugs.

Withaniasomnifera phytochemicals possess SARS-CoV-2 RdRp and human TMPRSS2 protein binding potential

Abstract

Coronavirus disease-19 (COVID-19) pandemic caused by severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) has infected approximately 26 million people and caused more than 6 million deaths globally. Spike (S)-protein on the outer surface of the virus uses human trans-membrane serine protease-2 (TMPRSS2) to gain entry into the cell. Recent reports indicate that human dipeptidyl peptidase-4 inhibitors (DPP4 or CD26) could also be utilized to check the S-protein mediated viral entry into COVID-19 patients. RNA dependent RNA polymerase (RdRp) is another key virulence protein of SARS-CoV-2 life cycle. The study aimed to identify the potential anti-SARS-CoV-2 inhibitors present in Withania somnifera (Solanaceae) using computer aided drug discovery approach. Molecular docking results showed that flavone glycoside, sugar alcohol, and flavonoid present in W. somnifera showed - 11.69, - 11.61, - 10.1, - 7.71 kcal/mole binding potential against S-protein, CD26, RdRp, and TMPRSS2 proteins. The major standard inhibitors of the targeted proteins (Sitagliptin, VE607, Camostat mesylate, and Remdesivir) showed the - 7.181, - 6.6, - 5.146, and - 7.56 kcal/mole binding potential. Furthermore, the lead phytochemicals and standard inhibitors bound and non-bound RdRp and TMPRSS2 proteins were subjected to molecular dynamics (MD) simulation to study the complex stability and change in protein conformation. The result showed energetically favorable and stable complex formation in terms of RMSD, RMSF, SASA, Rg, and hydrogen bond formation. Drug likeness and physiochemical properties of the test compounds exhibited satisfactory results. Taken together, the present study suggests the presence of potential anti-SARS-CoV-2 phytochemicals in W. somnifera that requires further validation in in vitro and in vivo studies.

Graphical Abstract

Supplementary information

The online version contains supplementary material available at 10.1007/s42535-022-00404-4.

Drug Resistance Mechanism of M46I-Mutation-Induced Saquinavir Resistance in HIV-1 Protease Using Molecular Dynamics Simulation and Binding Energy Calculation

Drug-resistance-associated mutation in essential proteins of the viral life cycle is a major concern in anti-retroviral therapy. M46I, a non-active site mutation in HIV-1 protease has been clinically associated with saquinavir resistance in HIV patients. A 100 ns molecular dynamics (MD) simulation and MM-PBSA calculations were performed to study the molecular mechanism of M46I-mutation-based saquinavir resistance. In order to acquire deeper insight into the drug-resistance mechanism, the flap curling, closed/semi-open/open conformations, and active site compactness were studied. The M46I mutation significantly affects the energetics and conformational stability of HIV-1 protease in terms of RMSD, RMSF, Rg, SASA, and hydrogen formation potential. This mutation significantly decreased van der Waals interaction and binding free energy (∆G) in the M46I–saquinavir complex and induced inward flap curling and a wider opening of the flaps for most of the MD simulation period. The predominant open conformation was reduced, but inward flap curling/active site compactness was increased in the presence of saquinavir in M46I HIV-1 protease. In conclusion, the M46I mutation induced structural dynamics changes that weaken the protease grip on saquinavir without distorting the active site of the protein. The produced information may be utilized for the discovery of inhibitor(s) against drug-resistant HIV-1 protease.

A candidate triple-negative breast cancer vaccine design by targeting clinically relevant cell surface markers: an integrated immuno and bio-informatics approach

Triple-negative breast cancer (TNBC) is an aggressive, metastatic/invasive sub-class of breast cancer (BCa). Cell surface protein-derived multi-epitope vaccine-mediated targeting of TNBC cells could be a better strategy against the disease. Literature-based identified potential cell surface markers for TNBC cells were subjected to expression pattern and survival analysis in BCa patient sample using TCGA database. The cytotoxic and helper T-lymphocytes antigenic epitopes in the test proteins were identified, selected and fused together with the appropriate linkers and an adjuvant, to construct the multi-epitope vaccine (MEV). The immune profile, physiochemical property (PP) and world population coverage of the MEV was studied. Immune simulation, cloning in a suitable vector, molecular docking (against Toll-like receptors, MHC (I/II) molecules), and molecular dynamics simulations of the MEV was performed. Cell surface markers were differentially expressed in TNBC samples and showed poor survival in TNBC patients. Satisfactory PP and WPC (up to 89 and 99%) was observed. MEV significant stable binding with the immune molecules and induced the immune cells in silico. The designed vaccine has capability to elicit immune response which could be utilized to target TNBC alone/combination with other therapy. The experimental studies are required to check the efficacy of the vaccine.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-022-03140-3.

The Conserved Macrodomain Is a Potential Therapeutic Target for Coronaviruses and Alphaviruses

Emerging and re-emerging viral diseases pose continuous public health threats, and effective control requires a combination of non-pharmacologic interventions, treatment with antivirals, and prevention with vaccines. The COVID-19 pandemic has demonstrated that the world was least prepared to provide effective treatments. This lack of preparedness has been due, in large part, to a lack of investment in developing a diverse portfolio of antiviral agents, particularly those ready to combat viruses of pandemic potential. Here, we focus on a drug target called macrodomain that is critical for the replication and pathogenesis of alphaviruses and coronaviruses. Some mutations in alphavirus and coronaviral macrodomains are not tolerated for virus replication. In addition, the coronavirus macrodomain suppresses host interferon responses. Therefore, macrodomain inhibitors have the potential to block virus replication and restore the host’s protective interferon response. Viral macrodomains offer an attractive antiviral target for developing direct acting antivirals because they are highly conserved and have a structurally well-defined (druggable) binding pocket. Given that this target is distinct from the existing RNA polymerase and protease targets, a macrodomain inhibitor may complement current approaches, pre-empt the threat of resistance and offer opportunities to develop combination therapies for combating COVID-19 and future viral threats.

Lack of efficacy of mono-mode of action therapeutics in COVID-19 therapy - How the lack of predictive power of preclinical cell and animal studies leads developments astray

The diverse experiences regarding the failure of tested drugs in the fight against COVID‐19 made it clear that one should at least question the requirement to apply classical preclinical development strategies that demand cell and animal efficacy models to be tested before going into clinical trials. Most animals are not susceptible to infection with SARS‐CoV‐2, and so this led to one‐sided virus replication experiments in cells and the use of animal models that have little in common with the complex pathogenesis of COVID‐19 in humans. Therefore, non‐clinical development strategies were designed to meet regulatory requirements, but they did not truly reflect the situation in the clinic. This has led the search for effective agents astray in many cases. As proof of this statement, we now bring together the results of such required preclinical experiments and compare with the results in clinical trials. Two clear conclusions that can be drawn from the experience to date: The required preclinical models are unsuitable for the development of innovative treatments medical devices in the case of COVID‐19 and mono‐action strategies (e.g. direct antivirals) are of very little or no benefit to patients under randomized,blinded conditions. Our hypothesis is that the complex situation of COVID‐19 may benefit from multi‐mode drugs. Here, the molecular class of aptamers could be a solution.

Identification of Natural Inhibitors Against SARS-CoV-2 Drugable Targets Using Molecular Docking, Molecular Dynamics Simulation, and MM-PBSA Approach

The present study explores the SARS-CoV-2 drugable target inhibition efficacy of phytochemicals from Indian medicinal plants using molecular docking, molecular dynamics (MD) simulation, and MM-PBSA analysis. A total of 130 phytochemicals were screened against SARS-CoV-2 Spike (S)-protein, RNA-dependent RNA polymerase (RdRp), and Main protease (M^pro). Result of molecular docking showed that Isoquercetin potentially binds with the active site/protein binding site of the Spike, RdRP, and Mpro targets with a docking score of -8.22, -6.86, and -9.73 kcal/mole, respectively. Further, MS 3, 7-Hydroxyaloin B, 10-Hydroxyaloin A, showed -9.57, -7.07, -8.57 kcal/mole docking score against Spike, RdRP, and M^pro targets respectively. The MD simulation was performed to study the favorable confirmation and energetically stable complex formation ability of Isoquercetin and 10-Hydroxyaloin A phytochemicals in M^pro-unbound/ligand bound/standard inhibitor bound system. The parameters such as RMSD, RMSF, Rg, SASA, Hydrogen-bond formation, energy landscape, principal component analysis showed that the lead phytochemicals form stable and energetically stabilized complex with the target protein. Further, MM-PBSA analysis was performed to compare the Gibbs free energy of the M^pro-ligand bound and standard inhibitor bound complexes. The analysis revealed that the His-41, Cys145, Met49, and Leu27 amino acid residues were majorly responsible for the lower free energy of the complex. Drug likeness and physiochemical properties of the test compounds showed satisfactory results. Taken together, the study concludes that that the Isoquercetin and 10-Hydroxyaloin A phytochemical possess significant efficacy to bind SARS-Cov-2 M^pro active site. The study necessitates further in vitro and in vivo experimental validation of these lead phytochemicals to assess their anti-SARS-CoV-2 potential.

An overview of potential inhibitors targeting non-structural proteins 3 (PLpro and Mac1) and 5 (3CLpro/Mpro) of SARS-CoV-2

There is an urgent need to develop effective treatments for coronavirus disease 2019 (COVID-19), which is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The rapid spread of SARS-CoV-2 has resulted in a global pandemic that has not only affected the daily lives of individuals but also had a significant impact on the global economy and public health. Although extensive research has been conducted to identify inhibitors targeting SARS-CoV-2, there are still no effective treatment strategies to combat COVID-19. SARS-CoV-2 comprises two important proteolytic enzymes, namely, the papain-like proteinase, located within non-structural protein 3 (nsp3), and nsp5, both of which cleave large replicase polypeptides into multiple fragments that are required for viral replication. Moreover, a domain within nsp3, known as the macrodomain (Mac1), also plays an important role in viral replication. Inhibition of their functions should be able to significantly interfere with the replication cycle of the virus, and therefore these key proteins may serve as potential therapeutic targets. The functions of the above viral targets and their corresponding inhibitors have been summarized in the current review. This review provides comprehensive updates of nsp3 and nsp5 inhibitor development and would help advance the discovery of novel anti-viral therapeutics against SARS-CoV-2.

Potential diabetes overtreatment and risk of adverse events among older adults in Ontario: a population-based study

AIMS/HYPOTHESIS

More than 25% of older adults (age ≥75 years) have diabetes and may be at risk of adverse events related to treatment. The aim of this study was to assess the prevalence of intensive glycaemic control in this group, potential overtreatment among older adults and the impact of overtreatment on the risk of serious events.

METHODS

We conducted a retrospective, population-based cohort study of community-dwelling older adults in Ontario using administrative data. Participants were ≥75 years of age with diagnosed diabetes treated with at least one anti-hyperglycaemic agent between 2014 and 2015. Individuals were categorised as having intensive or conservative glycaemic control (HbA_1c <53 mmol/mol [<7%] or 54-69 mmol/mol [7.1-8.5%], respectively), and as undergoing treatment with high-risk (i.e. insulin, sulfonylureas) or low-risk (other) agents. We measured the composite risk of emergency department visits, hospitalisations, or death within 30 days of reaching intensive glycaemic control with high-risk agents.

RESULTS

Among 108,620 older adults with diagnosed diabetes in Ontario, the mean (± SD) age was 80.6 (±4.5) years, 49.7% were female, and mean (± SD) diabetes duration was 13.7 (±6.3) years. Overall, 61% of individuals were treated to intensive glycaemic control and 21.6% were treated to intensive control using high-risk agents. Using inverse probability treatment weighting with propensity scores, intensive control with high-risk agents was associated with nearly 50% increased risk of the composite outcome compared with conservative glycaemic control with low-risk agents (RR 1.49, 95% CI 1.08, 2.05).

CONCLUSIONS/INTERPRETATION

Our findings underscore the need to re-evaluate glycaemic targets in older adults and to reconsider the use of anti-hyperglycaemic medications that may lead to hypoglycaemia, especially in setting of intensive glycaemic control.