Identification of putative antiviral bioactive compounds derived from family Asteraceae: An in silico approach

This computational study investigates 21 bioactive compounds from the Asteraceae family as potential inhibitors targeting the Spike protein (S protein) of SARS-CoV-2. Employing in silico methods and simulations, particularly CDOCKER and MM-GBSA, the study identifies two standout compounds, pterodontic acid and cichoric acid, demonstrating robust binding affinities (−46.1973 and −39.4265 kcal/mol) against the S protein. Comparative analysis with Favipiravir underscores their potential as promising inhibitors. Remarkably, these bioactives exhibit favorable ADMET properties, suggesting safety and efficacy. Molecular dynamics simulations validate their stability and interactions, signifying their potential as effective SARS-CoV-2 inhibitors.     

In vitro and in silico evaluation of fucosterol from Sargassum horridum as potential human acetylcholinesterase inhibitor

The fucosterol has been reported numerous biological activities. In this study, the activity in vitro of the fucosterol from Sargassum horridum as potential human acetylcholinesterase inhibitor was evaluated. The structural identification was obtained by nuclear magnetic resonance (NMR) spectroscopy and based on experimental data, we combined docking and molecular dynamics simulations coupled to the molecular-mechanics-generalized-born-surface-area approach to evaluating the structural and energetic basis for the molecular recognition of fucosterol and neostigmine at the binding site of acetylcholinesterase (AChE). In addition, the Lineweaver–Burk plot showed the nature of a non-competitive inhibition. The maximum velocity (V_max) and the constant of Michaelis–Menten (K_m) estimated for fucosterol (0.006 µM) were 0.015 1/V_o (ΔA/h and 6.399 1/[ACh] mM^?1, respectively. While, for neostigmine (0.14 µM), the V_max was 0.022 1/V_o (ΔA/h) and K_m of 6.726 1/[ACh] mM^?1, these results showed a more effective inhibition by fucosterol respect to neostigmine. Structural analysis revealed that neostigmine reaches the AChE binding site reported elsewhere, whereas fucosterol can act as a no-competitive and competitive acetylcholinesterase inhibitor, in agree with kinetic enzymatic experiments. Binding free energy calculations revealed that fucosterol reaches the acetylcholinesterase binding site with higher affinity than neostigmine, which is according to experimental results. Whereas the per-residue decomposition free energy analysis let us identify crucial residues involved in the molecular recognition of ligands by AChE. Results corroborate the ability of theoretical methods to provide crucial information at the atomic level about energetic and structural differences in the binding interaction and affinity from fucosterol with AChE.

Communicated by Ramaswamy H. Sarma     

Marine-derived sea urchin compounds as potential anti-cancer drug candidate against colorectal cancer: In silico and in vitro studies

Sea urchin-derived compounds are potential candidates for the development of effective drugs for the treatment of cancer diseases. In this study, 19 compounds derived from sea urchin (Diadema savignyi) were used to treat colorectal cancer using the HCT116 cell line. However, molecular docking, ADME (absorption, distribution, metabolism, and excretion), toxicity, molecular dynamic (MD) simulation, and molecular mechanics generalized Born surface area (MM-GBSA) were used to confirm the ligand–protein interaction. Interactions of Importin-11 receptor with sea urchin compounds reveal that four compounds have higher binding affinities (ranging from -8.6 to -7.1 kcal/mol). In vitro testing revealed that the CID 6432458 compound was effective (docking score of −8.6 kcal/mol) against the HCT116 cell line. The cytotoxicity of HCT116 has been documented, with an IC50 value of 6.885 ± 4. MTT assay, apoptosis analysis, and cell cycle assay were utilized to examine cell death in colorectal cancer. In the MTT experiment, 15 µM and 20 µM dosages were associated with 77% cell death; however, flow cytometry analysis using the IC50 value revealed that the selected chemical induced greater apoptosis in the HCT116 cell line (58.5%). The gene expression data revealed that the apoptotic gene BAX is expressed at a higher level than the BCL-2 gene. The IPO11 gene was downregulated during treatment. In the experiment involving the cell cycle, the S phase for the 30  µM dose showed 75.1% apoptosis, which was greater than the other concentrations used alone. These in silico and in vitro analysis will not only provide new information about Importin-11 receptor and insight into colorectal cancer but will also facilitate the development of natural compounds in a significant and worthwhile manner.     

Design and evaluation of an apta-nano-sensor to detect Acetamiprid in vitro and in silico

Pesticide detection is a main concern of food safety experts. Therefore, it is urgent to design an accurate, rapid, and cheap test. Biosensors that detect pesticide residues could replace current methods, such as HPLC or GC-MC. This research designs a biosensor based on aptamer (Oligonucleotide ss-DNA) in the receptor role, silver nanoparticles (AgNPs) as optical sensors and salt (NaCl) as the aggregative inducer of AgNPs to detect the presence of Acetamiprid. After optimization, .6 μM aptamer and 100 mM salt were employed. The selectivity and sensitivity of the complex were examined by different pesticides and different Acetamiprid concentrations. To simulate in vitro experimental conditions, bioinformatics software was used as in silico analysis. The results showed the detection of Acetamiprid at the .02 ppm (89.8 nM) level in addition to selectivity. Docking outputs introduced two loops as active sites in aptamer and confirmed aptamer–Acetamiprid bonding. Circular dichroism spectroscopy (CD) confirmed upon Acetamiprid binding, aptamer was folded due to stem-loop formation. Stability of the Apt–Acetamiprid complex in a simulated aqueous media was examined by molecular dynamic studies.     

Virtual screening of small molecules databases for discovery of novel PARP-1 inhibitors: combination of in silico and in vitro studies

Poly(ADP-ribose) polymerase-1 (PARP-1) enzyme has critical roles in DNA replication repair and recombination. Thus, PARP-1 inhibitors play an important role in the cancer therapy. In the current study, we have performed combination of in silico and in vitro studies in order to discover novel inhibitors against PARP-1 target. Structure-based virtual screening was carried out for an available small molecules database. A total of 257,951 ligands from Otava database were screened at the binding pocket of PARP-1 using high-throughput virtual screening techniques. Filtered structures based on predicted binding energy results were then used in more sophisticated molecular docking simulations (i.e. Glide/standard precision, Glide/XP, induced fit docking – IFD, and quantum mechanics polarized ligand docking – QPLD). Potential high binding affinity compounds that are predicted by molecular simulations were then tested by in vitro methods. Computationally proposed compounds as PARP-1 inhibitors (Otava Compound Codes: 7111620047 and 7119980926) were confirmed by in vitro studies. In vitro results showed that compounds 7111620047 and 7119980926 have IC₅₀ values of 0.56 and 63 μM against PARP-1 target, respectively. The molecular mechanism analysis, free energy perturbation calculations using long multiple molecular dynamics simulations for the discovered compounds which showed high binding affinity against PARP-1 enzyme, as well as structure-based pharmacophore development (E-pharmacophore) studies were also studied.     

Eriodictyol and naringenin inhibit the formation of AGEs: An in vitro and molecular interaction study

Maillard reaction occurs between the carbonyl group of reducing sugars and the free amino groups of protein, which eventually results in the formation and accumulation of advanced glycation end products (AGEs) irreversibly. Excessive production of AGEs is associated with many diseases, such as Alzheimer disease, neuropathy, retinopathy, and nephropathy. In this study, the effects of eriodictyol and naringenin on the inhibition of AGEs were studied with bovine serum albumin (BSA)–methylglyoxal (MGO) model by spectroscopic techniques and molecular docking methods. The fluorescence spectroscopy results suggested that eriodictyol and naringenin could inhibit the formation of AGEs. Circular dichroism (CD) studies indicated that eriodictyol and naringenin could stabilize the structure of BSA and inhibit the formation of AGEs. The molecular docking results demonstrated that eriodictyol formed two hydrogen bonds with Lys 350 and Leu 480 and the main forces were hydrogen bonding and hydrophobic interactions. However, naringenin interacted with Arg 484 of BSA, and the main force was hydrophobic interaction. It can be concluded that eriodictyol and naringenin can inhibit the formation of AGEs and eriodictyol has stronger inhibitory activity of AGEs than that of naringenin, which is probably due to the additional hydroxyl group in the position C‐3′ of B ring of eriodictyol.     

Flavonoids as tyrosinase inhibitors in in silico and in vitro models: basic framework of SAR using a statistical modelling approach

Flavonoids are widely distributed in plants and constitute the most common polyphenolic phytoconstituents in the human diet. In this study, the in vitro inhibitory activity of 44 different flavonoids (1–44) against mushroom tyrosinase was studied, and an in silico study and type of inhibition for the most active compounds were evaluated too. Tyrosinase inhibitors block melanogenesis and take part in melanin production or distribution leading to pigmentation diseases. The in vitro study showed that quercetin was a competitive inhibitor (IC₅₀=44.38 ± 0.13 µM) and achieved higher antityrosinase activity than the control inhibitor kojic acid. The in silico results highlight the importance of the flavonoid core with a hydroxyl at C7 as a strong contributor of interference with tyrosinase activity. According to the developed statistical model, the activity of molecules depends on hydroxylation at C3 and methylation at C8, C7, and C3 in the benzo-γ-pyrane ring of the flavonoids.     

In silico and in vitro analysis of an Aspergillus niger chitin deacetylase to decipher its subsite sugar preferences

Chitin deacetylases (CDAs) are found in many different organisms ranging from marine bacteria to fungi and insects. These enzymes catalyze the removal of acetyl groups from chitinous substrates generating various chitosans, linear copolymers consisting of N-acetylglucosamine (GlcNAc) and glucosamine. CDAs influence the degree of acetylation of chitosans as well as their pattern of acetylation, a parameter that was recently shown to influence the physicochemical properties and biological activities of chitosans. The binding site of CDAs typically consists of around four subsites, each accommodating a single sugar unit of the substrate. It has been hypothesized that the subsite preferences for GlcNAc or glucosamine units play a crucial role in the acetylation pattern they generate, but so far, this characteristic was largely ignored and still lacks structural data on the involved residues. Here, we determined the crystal structure of an Aspergillus niger CDA. Then, we used molecular dynamics simulations, backed up with a variety of in vitro activity assays using different well-defined polymeric and oligomeric substrates, to study this CDA in detail. We found that Aspergillus niger CDA strongly prefers a GlcNAc sugar unit at its −1 subsite and shows a weak GlcNAc preference at the other noncatalytic subsites, which was apparent both when deacetylating and N-acetylating oligomeric substrates. Overall, our results show that the combination of in vitro and in silico methods used here enables the detailed analysis of CDAs, including their subsite preferences, which could influence their substrate targets and the characteristics of chitosans produced by these species.     

The role of S477N mutation in the molecular behavior of SARS-CoV-2 spike protein: An in-silico perspective

The attachment of SARA-CoV-2 happens between ACE2 and the receptor binding domain (RBD) on the spike protein. Mutations in this domain can affect the binding affinity of the spike protein for ACE2. S477N, one of the most common mutations reported in the recent variants, is located in the RBD. Today's computational approaches in biology, especially during the SARS-CoV-2 pandemic, assist researchers in predicting a protein's behavior in contact with other proteins in more detail. In this study, we investigated the interactions of the S477N-hACE2 in silico to find the impact of this mutation on its binding affinity for ACE2 and immunity responses using dynamics simulation, protein–protein docking, and immunoinformatics methods. Our computational analysis revealed an increased binding affinity of N477 for ACE2. Four new hydrogen and hydrophobic bonds in the mutant RBD-ACE2 were formed (with S19 and Q24 of ACE2), which do not exist in the wild type. Also, the protein spike structure in this mutation was associated with an increase in stabilization and a decrease in its fluctuations at the atomic level. N477 mutation can be considered as the cause of increased escape from the immune system through MHC-II.     

3D-QSAR studies and shape based virtual screening for identification of novel hits to inhibit MbtA in Mycobacterium tuberculosis

Mycobacterium tuberculosis, the pathogen responsible for tuberculosis, uses various strategies to survive in a variety of host lesions. The re-emergence of multi-drug-resistant strains of M. tuberculosis underlines the necessity to discover new molecules. Inhibitors of aryl acid adenylating enzyme, MbtA, involved in siderophore biosynthesis in M. tuberculosis, are being explored as potential anti tubercular agents. In this study, we have used 3D-QSAR models and shape based virtual screening to identify novel MbtA inhibitors. 3D-QSAR studies were carried out on nucleoside bisubstrate derivatives. Both Comparative Molecular Field Analysis (r²?=?.944 and r²_pred?=?.938) and Comparative Molecular Similarity Indices Analysis (r²?=?.892 and r²_pred?=?.842) models, developed using Gasteiger charges with all fields, predicted efficiently. A total of 13 hits were identified as novel prospective inhibitors for MbtA by utilizing an insilico workflow. Out of 13 hits, five top ranked hits were used for further molecular dynamics studies to gain more insights about the stability of the complexes.     

An in silico pharmacological approach toward the discovery of potent inhibitors to combat drug resistance HIV-1 protease variants

Protease inhibitors (PIs) are crucial drugs in highly active antiretroviral therapy for human immunodeficiency virus-1 (HIV-1) infections. However, resistance owing to mutations challenge the long-term efficacy in the medication of HIV-1-infected individuals. Lopinavir (LPV) and darunavir (DRV), two second-generation drugs are the most potent among PIs, hustling the drug resistance when mutations occur in the active and nonactive site of the protease (PR). Herein, we strive for compounds that can stifle the function of wild-type (WT) HIV-1 PR along with four major single mutants (I54M, V82T, I84V, and L90M) instigating resistance to the PIs using in silico approach. Six common compounds are retrieved from six databases using combined pharmacophore-based and structure-based virtual screening methodology. LPV and DRV are docked and the binding free energy is calculated to set the cut-off value for selecting compounds. Further, to gain insight into the stability of the complexes the molecular dynamics simulation (MDS) is carried out, which uncovers two lead molecules namely NCI-524545 and ZINC12866729. Both the lead molecules connect with WT and mutant HIV-1 PRs through strong and stable hydrogen bond interactions when compared with LPV and DRV throughout the trajectory analysis. Interestingly, NCI-524545 and ZINC12866729 exhibit direct interactions with I50/50′ by replacing the conserved water molecule as evidenced by MDS, which indicates the credible potency of these compounds. Hence, we concluded that NCI-524545 and ZINC12866729 have great puissant to restrain the role of drug resistance HIV-1 PR variants, which can also show better activity through in vivo and in vitro conditions.     

The use of flax (Linum usitatissimum) as a model system for studies on organogenesis in vitro: the effect of different carbohydrates

Three carbohydrate sources, maltose, sucrose and cellobiose, were compared for the respective effects on regeneration from internodal sections of flax. It was found that different patterns of morphogenesis were attributable to the type of carbohydrate and its concentration. Mean shoot regeneration was highest on a medium containing 87 mM maltose with a rate of 2.4 shoots per explant and root formation was stimulated to the greatest degree on a medium with 87 mM sucrose. This medium also resulted in the highest tissue dry weight.     

Biochemical changes induced by grape seed extract and low level laser therapy administration during intraoral wound healing in rat liver: an experimental and in silico study

In the present study, the changes that occur in rat liver tissue as a result of the use of grape seed extract (GSE) and low level laser therapy (LLLT) in intraoral wound (IW) healing are analyzed using biochemical parameters. Diode laser application groups received 8 J/cm² dose LLLT once a day for 4 days (810 nm wavelength, continuous mode, 0.25 W, 9 s). As a result of the biological parameter analysis, it was determined that the oxidative damage caused by the IWs and recovery period on 7th and 14th days could be substantially removed with GSE applications that have antioxidant capacity especially in rat liver tissue. In addition, the active compound of grape seed, catechin is studied in the active site of glycogen synthase kinase 3 (GSK3) target using molecular modeling approaches. Post-processing molecular dynamics (MD) results for catechin is compared with a standard GSK3 inhibitor. MD simulations assisted for better understanding of inhibition mechanism and the crucial amino acids contributing in the ligand binding. These results along with a through free energy analysis of ligands using sophisticated simulations methods are quite striking and it suggests a greater future role for simulation in deciphering complex patterns of molecular mechanism in combination with methods for understanding drug-receptor interactions.     

Accommodating receptor flexibility and free energy calculation to reduce false positive binders in the discovery of natural products blockers of SARS-COV-2 spike RBD-ACE2 interface

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-COV-2), which causes coronavirus disease-19 (COVID-19) has caused more than 2 million deaths around the globe. The high transmissibility rate of the disease is related to the strong interaction between the virus spike receptor-binding domain (RBD) and the human angiotensin-converting enzyme 2 (ACE2) as documented in several reports. In this study, using state-of-the-art computational methods, natural products were screened and their molecular mechanism to disrupt spike RBD-ACE2 recognition was evaluated. There is the sensitivity of results to receptor ensemble docking calculations. Binding free energy and MD simulation are important tools to evaluate the thermodynamics of binding stability and the capacity of top hits to disrupt RBD-ACE2 recognition. The free energy profiles provide a slight decrease in binding affinity of the virus-receptor interaction. Three flavonoids parvisoflavone B (3), alpinumisoflavone (5) and norisojamicin (2) were effective in blocking the viral entry by binding strongly at the spike RBD-ACE2 interface with the inhibition constant of 0.56, 0.78 and 0.93 μM, respectively. The same compounds demonstrated similar effect on free ACE2 protein. Compound (2), also demonstrated ability to bind strongly on free spike RBD. Well-tempered metadynamics established that parvisoflavone B (3) works by binding to three sites namely interface α, β and loop thereby inhibiting viral cell entry. Owing to their desirable pharmacokinetic properties, the presented top hit natural products are suggested for further SARS-COV-2 molecular targets and subsequent in vitro and in vivo evaluations.     

Molecular basis of drug resistance in smoothened receptor: An in silico study of protein resistivity and specificity

Smoothened (SMO) antagonist Vismodegib effectively inhibits the Hedgehog pathway in proliferating cancer cells. In early stage of treatment, Vismodegib exhibited promising outcomes to regress the tumors cells, but ultimately relapsed due to the drug resistive mutations in SMO mostly occurring before (primary mutations G497W) or after (acquired mutations D473H/Y) anti-SMO therapy. This study investigates the unprecedented insights of structural and functional mechanism hindering the binding of Vismodegib with sensitive and resistant mutant variants of SMO (SMO^Mut). Along with the basic dynamic understanding of Vismodegib-SMO complexes, network propagation theory based on heat diffusion principles is first time applied here to identify the modules of residues influenced by the individual mutations. The allosteric modulation by GLY497 residue in Vismodegib bound SMO wild-type (SMO^WT) conformation depicts the interconnections of intermediate residues of SMO with the atom of Vismodegib and identify two important motifs (E-X-P-L) and (Q-A-N-V-T-I-G) mediating this allosteric regulation. In this study a novel computational framework based on the heat diffusion principle is also developed, which identify significant residues of allosteric site causing drug resistivity in SMO^Mut. This framework could also be useful for assessing the potential allosteric sites of different other proteins. Moreover, previously reported novel inhibitor “ZINC12368305,” which is proven to make an energetically favorable complex with SMO^WT is chosen as a control sample to assess the impact of receptor mutation on its binding and subsequently identify the important factors that govern binding disparity between Vismodegib and ZINC12368305 bound SMO^WT/Mut conformations.     

Comparative molecular docking and molecular-dynamic simulation of wild-type- and mutant carboxylesterase with BTA-hydrolase for enhanced binding to plastic

According to the literature review, microbial degradation of polyethylene terephthalate by PETases has been detected effective and eco-friendly. However, the number of microorganisms capable of such feats is limited with some undesirable bioprospecting results. BTA-hydrolase has been already reported capable of degrading polyethylene terephthalate. Therefore, mutation by in silico site-directed mutagenesis means to introduce current isomer of PETase for polyethylene terephthalate degradative capability as a better approach to resolve this issue. This study aimed to use in silico site-directed mutagenesis to convert a carboxylesterase from Archaeoglobus fulgidus to BTA-hydrolase from Thermobifida fusca by replacing six amino acids in specific locations. This work was followed by molecular docking analysis with polyethylene terephthalate and polypropylene to compare their interactions. The best-docked enzyme-substrate complex was further subjected to molecular dynamics simulation to gauge the binding quality of the BTA-hydrolase, wild-type and mutant-carboxylesterase with only polyethylene terephthalate as a substrate. Results of molecular docking revealed lowest binding energy for the wild-type carboxylesterase-polypropylene complex (-7.5 kcal/mol). The root-mean-square deviation value was observed stable for BTA-hydrolase. Meanwhile, root-mean-square fluctuation was assessed with higher fluctuation for the mutated residue Lys178. Consequently, the Rg value for BTA-hydrolase-ligand complex (～1.68 nm) was the lowest compared to the mutant and wild-type carboxylesterase. The collective data conveyed that mutations imparted a minimal change in the ability of the mutant carboxylesterase to bind to polyethylene terephthalate.