Identification of Potential Antiviral Inhibitors from Hydroxychloroquine and 1,2,4,5-Tetraoxanes Analogues and Investigation of the Mechanism of Action in SARS-CoV-2.

This study aimed to identify potential inhibitors and investigate the mechanism of action on SARS-CoV-2 ACE2 receptors using a molecular modeling study and theoretical determination of biological activity. Hydroxychloroquine was used as a pivot structure and antimalarial analogues of 1,2,4,5 tetraoxanes were used for the construction and evaluation of pharmacophoric models. The pharmacophore-based virtual screening was performed on the Molport® database (~7.9 million compounds) and obtained 313 structures. Additionally, a pharmacokinetic study was developed, obtaining 174 structures with 99% confidence for human intestinal absorption and penetration into the blood-brain barrier (BBB); posteriorly, a study of toxicological properties was realized. Toxicological predictions showed that the selected molecules do not present a risk of hepatotoxicity, carcinogenicity, mutagenicity, and skin irritation. Only 54 structures were selected for molecular docking studies, and five structures showed binding affinity (ΔG) values satisfactory for ACE2 receptors (PDB 6M0J), in which the molecule MolPort-007-913-111 had the best ΔG value of -8.540 Kcal/mol, followed by MolPort-002-693-933 with ΔG = -8.440 Kcal/mol. Theoretical determination of biological activity was realized for 54 structures, and five molecules showed potential protease inhibitors. Additionally, we investigated the Mpro receptor (6M0K) for the five structures via molecular docking, and we confirmed the possible interaction with the target. In parallel, we selected the TopsHits 9 with antiviral potential that evaluated synthetic accessibility for future synthesis studies and in vivo and in vitro tests.

Natural Products-Based Drug Design against SARS-CoV-2 Mpro 3CLpro.

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has received global attention due to the serious threat it poses to public health. Since the outbreak in December 2019, millions of people have been affected and its rapid global spread has led to an upsurge in the search for treatment. To discover hit compounds that can be used alone or in combination with repositioned drugs, we first analyzed the pharmacokinetic and toxicological properties of natural products from Brazil's semiarid region. After, we analyzed the site prediction and druggability of the SARS-CoV-2 main protease (Mpro), followed by docking and molecular dynamics simulation. The best SARS-CoV-2 Mpro complexes revealed that other sites were accessed, confirming that our approach could be employed as a suitable starting protocol for ligand prioritization, reinforcing the importance of catalytic cysteine-histidine residues and providing new structural data that could increase the antiviral development mainly against SARS-CoV-2. Here, we selected 10 molecules that could be in vitro assayed in response to COVID-19. Two compounds (b01 and b02) suggest a better potential for interaction with SARS-CoV-2 Mpro and could be further studied.

Revisiting the Proposition of Binding Pockets and Bioactive Poses for GSK-3ß Allosteric Modulators Addressed to Neurodegenerative Diseases.

Glycogen synthase kinase-3 beta (GSK-3ß) is an enzyme pertinently linked to neurodegenerative diseases since it is associated with the regulation of key neuropathological features in the central nervous system. Among the different kinds of inhibitors of this kinase, the allosteric ones stand out due to their selective and subtle modulation, lowering the chance of producing side effects. The mechanism of GSK-3ß allosteric modulators may be considered still vague in terms of elucidating a well-defined binding pocket and a bioactive pose for them. In this context, we propose to reinvestigate and reinforce such knowledge by the application of an extensive set of in silico methodologies, such as cavity detection, ligand 3D shape analysis and docking (with robust validation of corresponding protocols), and molecular dynamics. The results here obtained were consensually consistent in furnishing new structural data, in particular by providing a solid bioactive pose of one of the most representative GSK-3ß allosteric modulators. We further applied this to the prospect for new compounds by ligand-based virtual screening and analyzed the potential of the two obtained virtual hits by quantum chemical calculations. All potential hits achieved will be subsequently tested by in vitro assays in order to validate our approaches as well as to unveil novel chemical entities as GSK-3ß allosteric modulators.

Identification of Novel Chemical Entities for Adenosine Receptor Type 2A Using Molecular Modeling Approaches.

Adenosine Receptor Type 2A (A2AAR) plays a role in important processes, such as anti-inflammatory ones. In this way, the present work aimed to search for compounds by pharmacophore-based virtual screening. The pharmacokinetic/toxicological profiles of the compounds, as well as a robust QSAR, predicted the binding modes via molecular docking. Finally, we used molecular dynamics to investigate the stability of interactions from ligand-A2AAR. For the search for A2AAR agonists, the UK-432097 and a set of 20 compounds available in the BindingDB database were studied. These compounds were used to generate pharmacophore models. Molecular properties were used for construction of the QSAR model by multiple linear regression for the prediction of biological activity. The best pharmacophore model was used by searching for commercial compounds in databases and the resulting compounds from the pharmacophore-based virtual screening were applied to the QSAR. Two compounds had promising activity due to their satisfactory pharmacokinetic/toxicological profiles and predictions via QSAR (Diverset 10002403 pEC50 = 7.54407; ZINC04257548 pEC50 = 7.38310). Moreover, they had satisfactory docking and molecular dynamics results compared to those obtained for Regadenoson (Lexiscan®), used as the positive control. These compounds can be used in biological assays (in vitro and in vivo) in order to confirm the potential activity agonist to A2AAR.

Experimental and theoretical study on structure-tautomerism among edaravone, isoxazolone, and their heterocycles derivatives as antioxidants.

Edaravone is a heterocyclic pyrazolone compound. It has pronounced effect against free radicals, however renal and hepatic disorders have been reported. Isoxazolones are considered bioisosteric analogues of pyrazolones and may have comparable properties. Thus, we investigated the structural and electronic influences for edaravone, isoxazolone, and their tautomers on antioxidant process. Structure and tautomerism study among edaravone, isoxazolone and their heterocycles derivatives were related to antioxidant mechanisms by using the hybrid DFT method B3LYP with the basis sets 6-31++G(2d,2p). The C-H tautomer was the most stable and energetically favored among them. Intramolecular N-H-N hydrogen bonds and polar medium were responsible for the low energy differences among all possible tautomers. N-H tautomers in both systems proved to be better antioxidant by SET (single electron transfer), while O-H tautomers were better antioxidant on HAT (homolytic hydrogen atom transfer) mechanism. Theoretical calculation showed that edaravone is more potent than phenylisoxazolone, however, both has similar antioxidant scavenging on experimental DPPH. The carbonyliminic system played a very important role in the antioxidant activity for both studied classes.

In Silico Evaluation of Ibuprofen and Two Benzoylpropionic Acid Derivatives with Potential Anti-Inflammatory Activity.

Inflammation is a complex reaction involving cellular and molecular components and an unspecific response to a specific aggression. The use of scientific and technological innovations as a research tool combining multidisciplinary knowledge in informatics, biotechnology, chemistry and biology are essential for optimizing time and reducing costs in the drug design. Thus, the integration of these in silico techniques makes it possible to search for new anti-inflammatory drugs with better pharmacokinetic and toxicological profiles compared to commercially used drugs. This in silico study evaluated the anti-inflammatory potential of two benzoylpropionic acid derivatives (MBPA and DHBPA) using molecular docking and their thermodynamic profiles by molecular dynamics, in addition to predicting oral bioavailability, bioactivity and toxicity. In accordance to our predictions the derivatives proposed here had the potential capacity for COX-2 inhibition in the human and mice enzyme, due to containing similar interactions with the control compound (ibuprofen). Ibuprofen showed toxic predictions of hepatotoxicity (in human, mouse and rat; toxicophoric group 2-arylacetic or 3-arylpropionic acid) and irritation of the gastrointestinal tract (in human, mouse and rat; toxicophoric group alpha-substituted propionic acid or ester) confirming the literature data, as well as the efficiency of the DEREK 10.0.2 program. Moreover, the proposed compounds are predicted to have a good oral bioavailability profile and low toxicity (LD50 < 700 mg/kg) and safety when compared to the commercial compound. Therefore, future studies are necessary to confirm the anti-inflammatory potential of these compounds.

Toward of Safer Phenylbutazone Derivatives by Exploration of Toxicity Mechanism.

A drug design for safer phenylbutazone was been explored by reactivity and docking studies involving single electron transfer mechanism, as well as toxicological predictions. Several approaches about its structural properties were performed through quantum chemistry calculations at the B3LYP level of theory, together with the 6-31+G(d,p) basis sets. Molecular orbital and ionization potential were associated to electron donation capacity. The spin densities contribution showed a preferential hydroxylation at the para-positions of phenyl ring when compared to other positions. In addition, on electron abstractions the aromatic hydroxylation has more impact than alkyl hydroxylation. Docking studies indicate that six structures 1, 7, 8 and 13⁻15 have potential for inhibiting human as well as murine COX-2, due to regions showing similar intermolecular interactions to the observed for the control compounds (indomethacin and refecoxib). Toxicity can be related to aromatic hydroxylation. In accordance to our calculations, the derivatives here proposed are potentially more active as well safer than phenylbutazone and only structures 8 and 13⁻15 were the most promising. Such results can explain the biological properties of phenylbutazone and support the design of potentially safer candidates.

Time-Dependent Density Functional Theory Analysis of Triphenylamine-Functionalized Graphene Doped with Transition Metals for Photocatalytic Hydrogen Production.

The electronic structures and optical properties of triphenylamine-functionalized graphene (G-TPA) doped with transition metals, using water as a solvent, were theoretically investigated to verify the efficiency of photocatalytic hydrogen production with the use of transition metals. This study was performed by Density Functional Theory and Time-dependent Density Functional Theory through Gaussian 09W software, adopting the B3LYP functional for all structures. The 6-31g(d) basis set was used for H, C and N atoms, and the LANL2DZ basis set for transition metals using the Effective Core Potentials method. Two approaches were adopted: (1) using single metallic dopants (Ni, Pd, Fe, Os and Pt) and (2) using combinations of Ni with the other dopants (NiPd, NiPt, NiFe and NiOs). The DOS spectra reveal an increase of accessible states in the valence shell, in addition to a gap decrease for all dopants. This doping also increases the absorption in the visible region of solar radiation where sunlight is most intense (400 nm to 700 nm), with additional absorption peaks. The results lead us to propose the G-TPA structures doped with Ni, Pd, Pt, NiPt or NiPd to be novel catalysts for the conversion of solar energy for photocatalytic hydrogen production, since they improve the absorption of solar energy in the range of interest for solar radiation; and act as reaction centers, reducing the required overpotential for hydrogen production from water.

The antioxidant properties of salicylate derivatives: A possible new mechanism of anti-inflammatory activity.

The synthesis and antioxidant evaluation by DPPH scavenging of a series of salicylic acid derivatives is described. Gentisic acid and its ester, amide, and amino analogs possess more radical scavenging capacity than salicylic acid and other salicylate derivatives. This property can possibly provide an additional pathway for anti-inflammatory activity through either single electron or hydrogen atom transfer, leading to a new strategy for the design of anti-inflammatory agents.

Flagella interacting with a carbon nanowire with the variation of time and initial temperature.

The system proposed consists of a flagellum relaxing around a static carbon nanowire to mimics behavior of a natural flagellum moving with damped harmonic motion along a wire under van der Waals and electrostatic forces. This flagellum is composed of a C20 nanosphere with different sizes of his tail formed by hydrocarbons. The thermodynamic properties such as molar entropy variation, as well as molar heat dissipation, efficiency and speed were obtained to evaluate which system is most stable by using the variable temperature. This system has a number of carbon atoms ranging from 103-110, with a maximum of 300 ps for each simulation. We had simulated molar entropy variation, energies and efficiency changing with time and initial temperature. The results indicate that among the systems studied, the flagellum with five carbon atoms achieved greater stability and better results in this search.

Evaluation and theoretical study on the anti-inflammatory mechanism of 1-nitro-2-phenylethane.

In this study, 1-nitro-2-phenylethane was evaluated with respect to its effects in edema models of acute inflammation induced with carrageenan, dextran, and croton oil. 1-Nitro-2-phenylethane produced inhibition of rat paw edema induced by carrageenan and dextran at the doses of 25 and 50 mg/kg. The same doses caused an inhibition of croton oil-induced ear edema in mice. Our results suggest that 1-nitro-2-phenylethane has anti-inflammatory activity, probably of peripheral origin, acting in the synthesis and/or release of inflammatory mediators. A conformational study of 1-nitro-2-phenylethane was carried out using density functional theory calculations, showing three different groups of conformers corresponding to energy minimum geometries. The stereoelectronic repulsions are responsible for conformational preferences and the one most stable conformer. The prostaglandin endoperoxide synthase mechanism is related more to electrophilic than nucleophilic properties.

Understanding the molecular aspects of tetrahydrocannabinol and cannabidiol as antioxidants.

An antioxidant mechanism of tetrahydrocannabinol (THC) and cannabidiol (CBD) were compared with a simplified model of α-tocopherol, butylhydroxytoluene and hydroxytoluene in order to understand the antioxidant nature of THC and CBD molecules using DFT. The following electronic properties were evaluated: frontier orbitals nature, ionization potential, O-H bond dissociation energy (BDEOH), stabilization energy, and spin density distribution. An important factor that shows an influence in the antioxidant property of THC is the electron abstraction at the phenol position. Our data indicate that the decrease of the HOMO values and the highest ionization potential values are related to phenol, ether, and alkyl moieties. On the other hand, BDEOH in molecules with the cyclohexenyl group at ortho position of phenol are formed from lower energies than the molecules with an ether group at the meta position. In the light of our results, the properties calculated here predict that THC has a sightly higher antioxidant potential than CBD.

Antioxidant capacity of simplified oxygen heterocycles and proposed derivatives by theoretical calculations.

CONTEXT: Some structural properties can be involved in the antioxidant capacity of several polyphenol derivatives, among them their simplified structures. This study examines the contribution of simplified structure for the antioxidant capacity of some natural and synthetic antioxidants. The resonance structures were related to the π-type electron system of carbon-carbon double bonds between both phenyl rings. Trans-resveratrol, phenyl-benzofuran, phenyl-indenone, and benzylidene-benzofuranone are the best basic antioxidant templates among the simplified derivatives studied here. Additionally, the stilbene moiety was found on the molecules with the best antioxidant capacity. Furthermore, our investigation suggests that these compounds can be used as antioxidant scaffold for designing and developing of new promising derivatives. METHODS: To investigate the structure-antioxidant capacity for sixteen simplified natural and proposed derivatives we have employed density functional theory and used Gaussian 09. Our DFT calculations were performed using the B3LYP functional and the 6-31+G(d,p) basis set. All electron transfer mechanisms were investigated by using values of HOMO, ionization potential, energy affinity, stabilization energies, and spin density distributions.

Analgesic Activity of 5-Acetamido-2-Hydroxy Benzoic Acid Derivatives and an In-Vivo and In-Silico Analysis of Their Target Interactions.

The design, synthesis, and evaluation of novel non-steroidal anti-inflammatory drugs (NSAIDs) with better activity and lower side effects are big challenges today. In this work, two 5-acetamido-2-hydroxy benzoic acid derivatives were proposed, increasing the alkyl position (methyl) in an acetamide moiety, and synthesized, and their structural elucidation was performed using 1H NMR and 13C NMR. The changes in methyl in larger groups such as phenyl and benzyl aim to increase their selectivity over cyclooxygenase 2 (COX-2). These 5-acetamido-2-hydroxy benzoic acid derivatives were prepared using classic methods of acylation reactions with anhydride or acyl chloride. Pharmacokinetics and toxicological properties were predicted using computational tools, and their binding affinity (kcal/mol) with COX-2 receptors (Mus musculus and Homo sapiens) was analyzed using docking studies (PDB ID 4PH9, 5KIR, 1PXX and 5F1A). An in-silico study showed that 5-acetamido-2-hydroxy benzoic acid derivates have a better bioavailability and binding affinity with the COX-2 receptor, and in-vivo anti-nociceptive activity was investigated by means of a writhing test induced by acetic acid and a hot plate. PS3, at doses of 20 and 50 mg/kg, reduced painful activity by 74% and 75%, respectively, when compared to the control group (20 mg/kg). Regarding the anti-nociceptive activity, the benzyl showed reductions in painful activity when compared to acetaminophen and 5-acetamido-2-hydroxy benzoic acid. However, the proposed derivatives are potentially more active than 5-acetamido-2-hydroxy benzoic acid and they support the design of novel and safer derivative candidates. Consequently, more studies need to be conducted to evaluate the different pharmacological actions, the toxicity of possible metabolites that can be generated, and their potential use in inflammation and pain therapy.

Density functional theory (DFT) study of edaravone derivatives as antioxidants.

Quantum chemical calculations at the B3LYP/6-31G* level of theory were employed for the structure-activity relationship and prediction of the antioxidant activity of edaravone and structurally related derivatives using energy (E), ionization potential (IP), bond dissociation energy (BDE), and stabilization energies (ΔE(iso)). Spin density calculations were also performed for the proposed antioxidant activity mechanism. The electron abstraction is related to electron-donating groups (EDG) at position 3, decreasing the IP when compared to substitution at position 4. The hydrogen abstraction is related to electron-withdrawing groups (EDG) at position 4, decreasing the BDE(CH) when compared to other substitutions, resulting in a better antioxidant activity. The unpaired electron formed by the hydrogen abstraction from the C-H group of the pyrazole ring is localized at 2, 4, and 6 positions. The highest scavenging activity prediction is related to the lowest contribution at the carbon atom. The likely mechanism is related to hydrogen transfer. It was found that antioxidant activity depends on the presence of EDG at the C(2) and C(4) positions and there is a correlation between IP and BDE. Our results identified three different classes of new derivatives more potent than edaravone.

Hierarchical Virtual Screening Based on Rocaglamide Derivatives to Discover New Potential Anti-Skin Cancer Agents.

Skin Cancer (SC) is among the most common type of cancers worldwide. The search for SC therapeutics using molecular modeling strategies as well as considering natural plant-derived products seems to be a promising strategy. The phytochemical Rocaglamide A (Roc-A) and its derivatives rise as an interesting set of reference compounds due to their in vitro cytotoxic activity with SC cell lines. In view of this, we performed a hierarchical virtual screening study considering Roc-A and its derivatives, with the aim to find new chemical entities with potential activity against SC. For this, we selected 15 molecules (Roc-A and 14 derivatives) and initially used them in docking studies to predict their interactions with Checkpoint kinase 1 (Chk1) as a target for SC. This allowed us to compile and use them as a training set to build robust pharmacophore models, validated by Pearson's correlation (p) values and hierarchical cluster analysis (HCA), subsequentially submitted to prospective virtual screening using the Molport® database. Outputted compounds were then selected considering their similarities to Roc-A, followed by analyses of predicted toxicity and pharmacokinetic properties as well as of consensus molecular docking using three software. 10 promising compounds were selected and analyzed in terms of their properties and structural features and, also, considering their previous reports in literature. In this way, the 10 promising virtual hits found in this work may represent potential anti-SC agents and further investigations concerning their biological tests shall be conducted.

Structure of dihydrochalcones and related derivatives and their scavenging and antioxidant activity against oxygen and nitrogen radical species.

Quantum mechanical calculations at B3LYP/6-31G** level of theory were employed to obtain energy (E), ionization potential (IP), bond dissociation enthalpy (O-H BDE) and stabilization energies (DE(iso)) in order to infer the scavenging activity of dihydrochalcones (DHC) and structurally related compounds. Spin density calculations were also performed for the proposed antioxidant activity mechanism of 2,4,6-trihydroxyacetophenone (2,4,6-THA). The unpaired electron formed by the hydrogen abstraction from the phenolic hydroxyl group of 2,4,6-THA is localized on the phenolic oxygen at 2, 6, and 4 positions, the C3 and C6 carbon atoms at ortho positions, and the C5 carbon atom at para position. The lowest phenolic oxygen contribution corresponded to the highest scavenging activity value. It was found that antioxidant activity depends on the presence of a hydroxyl at the C2 and C4 positions and that there is a correlation between IP and O-H BDE and peroxynitrite scavenging activity and lipid peroxidation. These results identified the pharmacophore group for DHC.

A comparative theoretical mechanism on simplified flavonoid derivatives and isoxazolone analogous as Michael system inhibitor.

Flavonoids are a big class of natural product and have a wide range of biological activities. Some of these applications depend on its antioxidant capacity. Nevertheless, another mechanism can be involved by means of alkylation reaction on α,ß-unsaturated carbonyl system. This study aimed to evaluate the antioxidant capacity and the chemical reactivity among simplified flavonoid derivatives and isoxazolone analogous as Michael system by using B3LYP functional and 6-311 g(d,p) basis set. Frontier molecular orbital, ionization potential (IP), spin density contributions, and Fukui index explain the antioxidant capacity and reactivity index on isoxazolone and its related derivatives. The best contribution at ß-alkene moiety is related to better reactivity of α,ß-unsaturated carbonyl group. A decrease in antioxidant capacity is related to an increase in the chemical reactivity index. The frontier molecular orbitals show that aurone is more reactive than isoxazolone. In accordance with Fukui index, isoxazolone can be better inhibitor as Michael system when compared to flavonoid derivatives. Graphical abstract.

Essential features for antioxidant capacity of ascorbic acid (vitamin C).

Vitamin C or ascorbic acid is an indispensable micronutrient for human health found principally on citrus species such as lemon and orange fruits and vegetables. It was involved in the production of proteins such as collagen. Its biochemical mechanism is related to its antioxidant capacity; however, its function at the cellular level is still unclear. Several theoretical studies about antioxidant and redox mechanisms for ascorbic acid were suggested; however, no derivative was proposed. Thereby, an electronic study of antioxidant capacity for ascorbic acid derivatives was performed using theoretical chemistry at the DFT/ B3LYP/6-311 + + (2d,2p) level of theory. Simplified derivatives show that enol hydroxyls are more important than any other functional group. The vicinal enolic hydroxyl on ß position is more important for antioxidant capacity of ascorbic than hydroxyl on α position. According to our molecular modifications, the keto-alkene compound showed the best values when compared to ascorbic acid in some molecular characteristics. No lactone derivatives have superior application potential as antioxidant when compared with ascorbic acid. Several structures are possible to be proposed and were related to spin density contributions and the increase of chemical stability. New promising structural derivatives related to ascorbic acid can be developed in the future.

Molecular modeling approaches of selective adenosine receptor type 2A agonists as potential anti-inflammatory drugs.

Adenosine A2A receptor (A2AR) is the predominant receptor in immune cells, where its activation triggers cAMP-mediated immunosuppressive signaling and the underlying inhibition of T cells activation and T cells-induced effects mediated by cAMP-dependent kinase proteins mechanisms. In this study, were used ADME/Tox, molecular docking and molecular dynamics simulations to investigate selective adenosine A2AR agonists as potential anti-inflammatory drugs. As a result, we obtained two promising compounds (A and B) that have satisfactory pharmacokinetic and toxicological properties and were able to interact with important residues of the A2AR binding cavity and during the molecular dynamics simulations were able to keep the enzyme complexed.Communicated by Ramaswamy H. Sarma.