Rational Approach toward COVID-19's Main Protease Inhibitors: A Hierarchical Biochemoinformatics Analysis.

This study investigated the potential of selected compounds as inhibitors of SARS-CoV-2 Mpro through pharmacokinetic and toxicological analyses, molecular docking, and molecular dynamics simulations. In silico molecular docking simulations revealed promising ligands with favorable binding affinities for Mpro, ranging from -6.2 to -9.5 kcal/mol. Moreover, molecular dynamics simulations demonstrated the stability of protein-ligand complexes over 200 ns, maintaining protein secondary structures. MM-PBSA analysis revealed favorable interactions between ligands and Mpro, with negative binding energy values. Hydrogen bond formation capacity during molecular dynamics was confirmed, indicating consistent interactions with Mpro catalytic residues. Based on these findings, selected ligands show promise for future studies in developing COVID-19 treatments.

Design and Identification of Inhibitors for the Spike-ACE2 Target of SARS-CoV-2.

When an epidemic started in the Chinese city of Wuhan in December 2019, coronavirus was identified as the cause. Infection by the virus occurs through the interaction of viral S protein with the hosts' angiotensin-converting enzyme 2. By leveraging resources such as the DrugBank database and bioinformatics techniques, ligands with potential activity against the SARS-CoV-2 spike protein were designed and identified in this investigation. The FTMap server and the Molegro software were used to determine the active site of the Spike-ACE2 protein's crystal structure. Virtual screening was performed using a pharmacophore model obtained from antiparasitic drugs, obtaining 2000 molecules from molport®. The ADME/Tox profiles were used to identify the most promising compounds with desirable drug characteristics. The binding affinity investigation was then conducted with selected candidates. A molecular docking study showed five structures with better binding affinity than hydroxychloroquine. Ligand_003 showed a binding affinity of -8.645 kcal·mol-1, which was considered an optimal value for the study. The values presented by ligand_033, ligand_013, ligand_044, and ligand_080 meet the profile of novel drugs. To choose compounds with favorable potential for synthesis, synthetic accessibility studies and similarity analyses were carried out. Molecular dynamics and theoretical IC50 values (ranging from 0.459 to 2.371 µM) demonstrate that these candidates are promising for further tests. Chemical descriptors showed that the candidates had strong molecule stability. Theoretical analyses here show that these molecules have potential as SARS-CoV-2 antivirals and therefore warrant further investigation.

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.

Identification of Potential COX-2 Inhibitors for the Treatment of Inflammatory Diseases Using Molecular Modeling Approaches.

Non-steroidal anti-inflammatory drugs are inhibitors of cyclooxygenase-2 (COX-2) that were developed in order to avoid the side effects of non-selective inhibitors of COX-1. Thus, the present study aims to identify new selective chemical entities for the COX-2 enzyme via molecular modeling approaches. The best pharmacophore model was used to identify compounds within the ZINC database. The molecular properties were determined and selected with Pearson's correlation for the construction of quantitative structure-activity relationship (QSAR) models to predict the biological activities of the compounds obtained with virtual screening. The pharmacokinetic/toxicological profiles of the compounds were determined, as well as the binding modes through molecular docking compared to commercial compounds (rofecoxib and celecoxib). The QSAR analysis showed a fit with R = 0.9617, R2 = 0.9250, standard error of estimate (SEE) = 0.2238, and F = 46.2739, with the tetra-parametric regression model. After the analysis, only three promising inhibitors were selected, Z-964, Z-627, and Z-814, with their predicted pIC50 (-log IC50) values, Z-814 = 7.9484, Z-627 = 9.3458, and Z-964 = 9.5272. All candidates inhibitors complied with Lipinski's rule of five, which predicts a good oral availability and can be used in in vitro and in vivo tests in the zebrafish model in order to confirm the obtained in silico data.

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.

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.

Hierarchical Virtual Screening of Potential New Antibiotics from Polyoxygenated Dibenzofurans against Staphylococcus aureus Strains.

Staphylococcus aureus is a microorganism with high morbidity and mortality due to antibiotic-resistant strains, making the search for new therapeutic options urgent. In this context, computational drug design can facilitate the drug discovery process, optimizing time and resources. In this work, computational methods involving ligand- and structure-based virtual screening were employed to identify potential antibacterial agents against the S. aureus MRSA and VRSA strains. To achieve this goal, tetrahydroxybenzofuran, a promising antibacterial agent according to in vitro tests described in the literature, was adopted as the pivotal molecule and derivative molecules were considered to generate a pharmacophore model, which was used to perform virtual screening on the Pharmit platform. Through this result, twenty-four molecules were selected from the MolPort® database. Using the Tanimoto Index on the BindingDB web server, it was possible to select eighteen molecules with greater structural similarity in relation to commercial antibiotics (methicillin and oxacillin). Predictions of toxicological and pharmacokinetic properties (ADME/Tox) using the eighteen most similar molecules, showed that only three exhibited desired properties (LB255, LB320 and LB415). In the molecular docking study, the promising molecules LB255, LB320 and LB415 showed significant values in both molecular targets. LB320 presented better binding affinity to MRSA (-8.18 kcal/mol) and VRSA (-8.01 kcal/mol) targets. Through PASS web server, the three molecules, specially LB320, showed potential for antibacterial activity. Synthetic accessibility (SA) analysis performed on AMBIT and SwissADME web servers showed that LB255 and LB415 can be considered difficult to synthesize and LB320 is considered easy. In conclusion, the results suggest that these ligands, particularly LB320, may bind strongly to the studied targets and may have appropriate ADME/Tox properties in experimental studies.

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.

Potential inhibitors of the enzyme acetylcholinesterase and juvenile hormone with insecticidal activity: study of the binding mode via docking and molecular dynamics simulations.

Models validation in QSAR, pharmacophore, docking and others can ensure the accuracy and reliability of future predictions in design and selection of molecules with biological activity. In this study, pyriproxyfen was used as a pivot/template to search the database of the Maybridge Database for potential inhibitors of the enzymes acetylcholinesterase and juvenile hormone as well. The initial virtual screening based on the 3D shape resulted in 2000 molecules with Tanimoto index ranging from 0.58 to 0.88. A new reclassification was performed on the overlapping of positive and negative charges, which resulted in 100 molecules with Tanimoto's electrostatic score ranging from 0.627 to 0.87. Using parameters related to absorption, distribution, metabolism and excretion and the pivot molecule, the molecules selected in the previous stage were evaluated regarding these criteria, and 21 were then selected. The pharmacokinetic and toxicological properties were considered and for 12 molecules, the DEREK software not fired any alert of toxicity, which were thus considered satisfactory for prediction of biological activity using the Web server PASS. In the molecular docking with insect acetylcholinesterase, the Maybridge3_002654 molecule had binding affinity of -11.1 kcal/mol, whereas in human acetylcholinesterase, the Maybridge4_001571molecule show in silico affinity of -10.2 kcal/mol, and in the juvenile hormone, the molecule MCULE-8839595892 show in silico affinity value of -11.6 kcal/mol. Subsequent long-trajectory molecular dynamics studies indicated considerable stability of the novel molecules compared to the controls.AbbreviationsQSARquantitative structure-activity relationshipsPASSprediction of activity spectra for substancesCommunicated by Ramaswamy H. Sarma.