Instant determination of the artemisinin from various Artemisia annua L. extracts by LC-ESI-MS/MS and their in-silico modelling and in vitro antiviral activity studies against SARS-CoV-2.

INTRODUCTION: Numerous efforts in natural product drug development are reported for the treatment of Coronavirus. Based on the literature, among these natural plants Artemisia annua L. shows some promise for the treatment of SARS-CoV-2. OBJECTIVE: The main objective of our study was to determine artemisinin content by liquid chromatography electrospray ionisation tandem mass spectrometry (LC-ESI-MS/MS), to investigate the in vitro biological activity of artemisinin from the A. annua plants grown in Turkey with various extracted methods, to elaborate in silico activity against SARS-CoV-2 using molecular modelling. METHODOLOGY: Twenty-one different extractions were applied. Direct and sequential extractions studies were compared with ultrasonic assisted maceration, Soxhlet, and ultra-rapid determined artemisinin active molecules by LC-ESI-MS/MS methods. The inhibition of spike protein and main protease (3CL) enzyme activity of SARS-CoV-2 virus was assessed by time resolved fluorescence energy transfer (TR-FRET) assay. RESULTS: Artemisinin content in the range 0.062-0.066%. Artemisinin showed significant inhibition of 3CL protease activity but not Spike/ACE-2 binding. The 50% effective concentration (EC50 ) of artemisinin against SARS-CoV-2 Spike pseudovirus was found greater than 50 µM (EC45 ) in HEK293T cell line whereas the cell viability was 94% of the control (P < 0.01). The immunosuppressive effects of artemisinin on TNF-α production on both pseudovirus and lipopolysaccharide (LPS)-induced THP-1 cells were found significant in a dose dependent manner. CONCLUSION: Further studies of these extracts for COVID-19 treatment will shed light to seek alternative treatment options. Moreover, these natural extracts can be used as an additional treatment option with medicines, as well as prophylactic use can be very beneficial for patients.

Screening of Clinically Approved and Investigation Drugs as Potential Inhibitors of SARS-CoV-2: A Combined in silico and in†vitro Study.

In the current study, we used 7922 FDA approved small molecule drugs as well as compounds in clinical investigation from NIH's NPC database in our drug repurposing study. SARS-CoV-2 main protease as well as Spike protein/ACE2 targets were used in virtual screening and top-100 compounds from each docking simulations were considered initially in short molecular dynamics (MD) simulations and their average binding energies were calculated by MM/GBSA method. Promising hit compounds selected based on average MM/GBSA scores were then used in long MD simulations. Based on these numerical calculations following compounds were found as hit inhibitors for the SARS-CoV-2 main protease: Pinokalant, terlakiren, ritonavir, cefotiam, telinavir, rotigaptide, and cefpiramide. In addition, following 3 compounds were identified as inhibitors for Spike/ACE2: Denopamine, bometolol, and rotigaptide. In order to verify the predictions of in silico analyses, 4 compounds (ritonavir, rotigaptide, cefotiam, and cefpiramide) for the main protease and 2 compounds (rotigaptide and denopamine) for the Spike/ACE2 interactions were tested by in vitro experiments. While the concentration-dependent inhibition of the ritonavir, rotigaptide, and cefotiam was observed for the main protease; denopamine was effective at the inhibition of Spike/ACE2 binding.

Development of Small Molecule MEIS Inhibitors that modulate HSC activity.

Meis1, which belongs to TALE-type class of homeobox gene family, appeared as one of the key regulators of hematopoietic stem cell (HSC) self-renewal and a potential therapeutical target. However, small molecule inhibitors of MEIS1 remained unknown. This led us to develop inhibitors of MEIS1 that could modulate HSC activity. To this end, we have established a library of relevant homeobox family inhibitors and developed a high-throughput in silico screening strategy against homeodomain of MEIS proteins using the AutoDock Vina and PaDEL-ADV platform. We have screened over a million druggable small molecules in silico and selected putative MEIS inhibitors (MEISi) with no predicted cytotoxicity or cardiotoxicity. This was followed by in vitro validation of putative MEIS inhibitors using MEIS dependent luciferase reporter assays and analysis in the ex vivo HSC assays. We have shown that small molecules named MEISi-1 and MEISi-2 significantly inhibit MEIS-luciferase reporters in vitro and induce murine (LSKCD34l°w cells) and human (CD34+, CD133+, and ALDHhi cells) HSC self-renewal ex vivo. In addition, inhibition of MEIS proteins results in downregulation of Meis1 and MEIS1 target gene expression including Hif-1α, Hif-2α and HSC quiescence modulators. MEIS inhibitors are effective in vivo as evident by induced HSC content in the murine bone marrow and downregulation of expression of MEIS target genes. These studies warrant identification of first-in-class MEIS inhibitors as potential pharmaceuticals to be utilized in modulation of HSC activity and bone marrow transplantation studies.

Novel tumor necrosis factor-α (TNF-α) inhibitors from small molecule library screening for their therapeutic activity profiles against rheumatoid arthritis using target-driven approaches and binary QSAR models.

Tumor necrosis factor alpha (TNF-α) is a multifunctional cytokine that acts as a central biological mediator for critical immune functions, including inflammation, infection, and antitumor responses. It plays pivotal role in autoimmune diseases like rheumatoid arthritis (RA). The synthetic antibodies etanercept, infliximab, and adalimumab are approved drugs for the treatment of inflammatory diseases bind to TNF-α directly, preventing its association with the tumor necrosis factor receptor (TNFR). These biologics causes serious side effects such as triggering an autoimmune anti-antibody response or the weakening of the body's immune defenses. Therefore, alternative small-molecule based therapies for TNF-α inhibition is a hot topic both in academia and industry. Most of small-molecule inhibitors reported in the literature target TNF-α, indirectly. In this study, combined in silico approaches have been applied to better understand the important direct interactions between TNF-α and small inhibitors. Our effort executed with the extensive literature review to select the compounds that inhibit TNF-α. High-throughput structure-based and ligand-based virtual screening methods are applied to identify TNF-α inhibitors from 3 different small molecule databases (∼256.000 molecules from Otava drug-like green chemical collection, ∼ 500.000 molecules from Otava Tangible database, ∼2.500.000 Enamine small molecule database) and ∼240.000 molecules from ZINC natural products libraries. Moreover, therapeutic activity prediction, as well as pharmacokinetic and toxicity profiles are also investigated using MetaCore/MetaDrug platform which is based on a manually curated database of molecular interactions, molecular pathways, gene-disease associations, chemical metabolism and toxicity information, uses binary QSAR models. Particular therapeutic activity and toxic effect predictions are based on the ChemTree ability to correlate structural descriptors to that property using recursive partitioning algorithm. Molecular Dynamics (MD) simulations were also performed for selected hits to investigate their detailed structural and dynamical analysis beyond docking studies. As a result, at least one hit from each database were identified as novel TNF-α inhibitors after comprehensive virtual screening, multiple docking, e-Pharmacophore modeling (structure-based pharmacophore modeling), MD simulations, and MetaCore/MetaDrug analysis. Identified hits show predicted promising anti-arthritic activity and no toxicity. Communicated by Ramaswamy H. Sarma.

Molecular modeling and in vitro approaches towards cholinesterase inhibitory effect of some natural xanthohumol, naringenin, and acyl phloroglucinol derivatives.

BACKGROUND: Many natural products, particularly phenolic compounds, have been reported to have a strong inhibition against acetylcholinesterase (AChE) and butyrylcholinesterase (BChE), the key enzymes in the pathology of Alzheimer's disease (AD). HYPOTHESIS: Therefore, we hypothesized that some xanthahumol, naringenin, and acyl phloroglucinol derivatives (1-14) isolated from Humulus lupulus L. (hops) may have an inhibitory potential against AChE and BChE. METHODS: Inhibitory potential of compounds 1-14 were tested against AChE and BChE using ELISA microtiter assay. Different molecular docking simulations, including IFD and GOLD protocols, were implemented to verify the interactions between the ligands and the active site amino acids and also their binding energies inside the catalytic crevices of AChE and BChE. ADME/Tox analysis were used to determine pharmacological activities of the compounds. RESULTS: Among them, 3­hydroxy­xanthohumol (IC50 = 51.25 ±â€¯0.88 µM) and xanthohumol (IC50 = 71.34 ±â€¯2.09 µM), displayed a moderate AChE inhibition in comparison to that of the reference (galanthamine, IC50 = 2.52 ±â€¯0.15 µM). In addition to 3­hydroxy­xanthohumol (IC50 = 63.07 ±â€¯3.76 µM) and xanthohumol (IC50 = 32.67 ±â€¯2.82 µM), 8-prenylnaringenin (IC50 = 86.58 ±â€¯3.74 µM) also showed micromolar-range inhibition against BChE (galanthamine, IC50 = 46.58 ±â€¯0.91 µM). Rest of the compounds were found to be either inactive or having inhibition below 50%. Prediction of pharmacokinetic studies suggested that all the ligands revealed acceptable drug-like profiles. Docking simulations demonstrate not only the prediction of ligand binding energies of the compounds inside the catalytic domains of the targets, but also highlight the critical amino acids contributing to stabilizations of the ligands. CONCLUSION: Our findings revealed that xanthohumol in particular could be considered as lead molecule to explore new cholinesterase inhibitors for AD.

Discovery of selective dengue virus inhibitors using combination of molecular fingerprint-based virtual screening protocols, structure-based pharmacophore model development, molecular dynamics simulations and in vitro studies.

Dengue virus is a major issue of tropical and sub-tropical regions. The proliferation of virus results in immense number of deaths each year because of unavailability of on-shelf drugs. This issue necessitates the design of novel anti-Dengue drugs. The protease enzyme pathway is the critical target for drug design due to its significance in the replication, survival and other cellular activities of Dengue virus. Keeping in mind the worsening situation regarding Dengue virus, approximately eighteen million drug-like compounds from the ZINC small molecule database have been screened against Nonstructural Protein 3 (NS3) previously by our group. In this study, in order to investigate the effect of extended time of molecular dynamics (MD) simulations on structural and dynamical profiles of used complexes, simulation run time is increased from 50-ns to 100-ns for the each system. In addition, a well-known Dengue virus inhibitor (MB21) from literature is used as reference structure (positive control) to compare the proposed molecules. Post-processing MD analyses including Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) calculations were conducted to predict binding free energies of inhibitors from derived trajectory frames of MD simulations. Identified compounds are further directed to Quantum-Polarized Ligand Docking (QPLD), molecular fingerprint-based virtual screening of another small molecule database (Otava Drug Like small molecule database), and Structure-based Pharmacophore Modeling (E-Pharmacophore). Finally, cell proliferation and cytotoxicity tests as well as pre- and post-treatment on HUH7 cells infected with DENV2 NGC strain are applied for four identified hit molecules (ZINC36681949, ZINC44921800, ZINC95518765 and ZINC39500661) to check whether these drugs inhibit DENV2 from entry and/or exit pathways. Based on cell-based Dengue quantification assays, there is no effect seen on pre-treatment of cells with these compounds indicating that the early infection processes of virus is not affected. In contrast, the post-treatment of cells with these compounds after Dengue virus infection has resulted in a significant 1 log PFU/ml reduction of the virus infectious titre.

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/cm2 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.

The effects of pollen, propolis, and caffeic acid phenethyl ester on tyrosine hydroxylase activity and total RNA levels in hypertensive rats caused by nitric oxide synthase inhibition: experimental, docking and molecular dynamic studies.

The objective of the present study was to evaluate the effects of propolis, pollen, and caffeic acid phenethyl ester (CAPE) on tyrosine hydroxylase (TH) activity and total RNA levels of Nω-nitro-L-arginine methyl ester (L-NAME) inhibition of nitric oxide synthase in the heart, adrenal medulla, and hypothalamus of hypertensive male Sprague dawley rats. The TH activity in the adrenal medulla, heart, and hypothalamus of the rats was significantly increased in the L-NAME group vs. control (p < 0.05). Treatment with L-NAME led to a significant increase in blood pressure (BP) in the L-NAME group compared to control (p < 0.05). These data suggest that propolis, pollen, and CAPE may mediate diminished TH activity in the heart, adrenal medulla, and hypothalamus in hypertensive rats. The decreased TH activity may be due to the modulation and synthesis of catecholamines and BP effects. In addition, the binding mechanism of CAPE within the catalytic domain of TH was investigated by means of molecular modeling approaches. These data suggest that the amino acid residues, Glu429 and Ser354 of TH may play a pivotal role in the stabilization of CAPE within the active site as evaluated by molecular dynamics (MD) simulations. Gibbs binding free energy (ΔGbinding) of CAPE in complex with TH was also determined by post-processing MD analysis approaches (i.e. Poisson-Boltzmann Surface Area (MM-PBSA) method).

Effects of propolis, caffeic acid phenethyl ester, and pollen on renal injury in hypertensive rat: An experimental and theoretical approach.

The objective of this study was to evaluate the antioxidant effects of propolis, caffeic acid phenethyl ester (CAPE; active compound in propolis), and pollen on biochemical oxidative stress biomarkers in rat kidney tissue inhibited by Nω -nitro-L-arginine methyl ester (L-NAME). The biomarkers evaluated were paraoxonase (PON1), oxidative stress index (OSI), total antioxidant status (TAS), total oxidant status (TOS), asymmetric dimethylarginine (ADMA), and nuclear factor kappa B (NF-κB). TAS levels and PON1 activity were significantly decreased in kidney tissue samples in the L-NAME-treated group (P < 0.05). The levels of TAS and PONI were higher in the L-NAME plus propolis, CAPE, and pollen groups compared with the L-NAME-treated group. TOS, ADMA, and NF-κB levels were significantly increased in the kidney tissue samples of the L-NAME-treated group (P < 0.05). However, these parameters were significantly lower in the L-NAME plus propolis, CAPE, and pollen groups (P < 0.05) compared with rats administered L-NAME alone (P < 0.05). Furthermore, the binding energy of CAPE within catalytic domain of glutathione reductase (GR) enzyme as well as its inhibitory mechanism was determined using molecular modeling approaches. In conclusion, experimental and theoretical data suggested that oxidative alterations occurring in the kidney tissue of chronic hypertensive rats may be prevented via active compound of propolis, CAPE administration.

First universal pharmacophore model for hERG1 K+ channel activators: acthER.

The intra-cavitary drug blockade of hERG1 channel has been extensively studied, both experimentally and theoretically. Structurally diverse ligands inadvertently block the hERG1 K+ channel currents lead to drug induced Long QT Syndrome (LQTS). Accordingly, designing either hERG1 channel openers or current activators, with the potential to target other binding pockets of the channel, has been introduced as a viable approach in modern anti-arrhythmia drug development. However, reports and investigations on the molecular mechanisms underlying activators binding to the hERG1 channel remain sparse and the overall molecular design principles are largely unknown. Most of the hERG1 activators were discovered during mandatory screening for hERG1 blockade. To fill this apparent deficit, the first universal pharmacophore model for hERG1 K+ channel activators was developed using PHASE. 3D structures of 18 hERG1 K+ channel activators and their corresponding measured binding affinity values were used in the development of pharmacophore models. These compounds spanned a range of structurally different chemotypes with moderate variation in binding affinity. A five sites AAHRR (A, hydrogen-bond accepting, H, hydrophobic, R, aromatic) pharmacophore model has shown reasonable high statistical results compared to the other developed more than 1000 hypotheses. This model was used to construct steric and electrostatic contour maps. The predictive power of the model was tested with 3 external test set compounds as true unknowns. Finally, the pharmacophore model was combined with the previously developed receptor-based model of hERG1 K+ channel to develop and screen novel activators. The results are quite striking and it suggests a greater future role for pharmacophore modeling and virtual drug screening simulations in deciphering complex patterns of molecular mechanisms of hERG1 channel openers at the target sites. The developed model is available upon request and it may serve as basis for the synthesis of novel therapeutic hERG1 activators.

Identification of novel serotonin reuptake inhibitors targeting central and allosteric binding sites: A virtual screening and molecular dynamics simulations study.

The serotonin (5-hydroxytryptamine, 5HT) transporter (SERT) is a member of neurotransmitter sodium symporter (NSS) family, which maintains neurotransmitter by reuptaking 5HT into synapses. Decrease in serotonin concentrations in synaptic clefts have been reported to cause psychological and neurological disorders. Therefore, inhibition of SERT is a potent strategy for the treatment of related diseases such as depression. In this study, approximately 260,000 small molecules from an available chemical database have been virtually screened both at central and allosteric binding sites of SERT to identify potent novel candidate SERT inhibitors. A set of docking algorithms were used to predict binding modes and energies of compounds. Screening analyses led three top-ranked hit compounds (160234, Otava ID: 7118020138; 159166, Otava ID: 7117171303; and 69419, Otava ID: 118671819) for central binding site (S1) and one compound (93507, Otava ID: 6248262) for allosteric binding site (S2). These promising compounds are then subjected to long multiple molecular dynamics (MD) simulations to elucidate their structural and dynamical profiles at the binding cavities of SERT. Higher predicted binding affinities of identified compounds were also confirmed with binding free energy calculations (MM/GBSA) in comparison with the reference central and allosteric binding site inhibitors, paroxetine (8PR) and escitalopram (68P), respectively. To the best of our knowledge, the present work is the first structure-based high throughput virtual screening study reported using recently revealed crystal structure of SERT for screening inhibitors from chemical databases on S1 and S2 binding sites. Small molecule library screening study yielded candidate compounds both at central and allosteric binding site of SERT, and further experimentation may pave the way for developing novel strong inhibitors.

Discovering novel carbonic anhydrase type IX (CA IX) inhibitors from seven million compounds using virtual screening and in vitro analysis.

Carbonic anhydrase type IX (CA IX) enzyme is mostly over expressed in different cancer cell lines and tumor tissues. Potent CA IX inhibitors can be effective for adjusting the pH imbalance in tumor cells. In the present work, we represented the successful application of high throughput virtual screening (HTVS) of large dataset from ZINC database included of ∼7 million compounds to discover novel inhibitors of CA IX. HTVS and molecular docking were performed using consequence Glide/standard precision (SP), extra precision (XP) and induced fit docking (IFD) molecular docking protocols. For each compound, docking code calculates a set of low-energy poses and then exhaustively scans the binding pocket of the target with small compounds. Novel CA IX inhibitor candidates were suggested based on molecular modeling studies and a few of them were tested using in vitro analysis. These compounds were determined as good inhibitors against human CA IX target with Ki in the range of 0.85-1.58 µM. In order to predict the pharmaceutical properties of the selected compounds, ADME (absorption, distribution, metabolism and excretion) analysis was also carried out.