Prediction of ligand binding mode among multiple cross-docking poses by molecular dynamics simulations

Optimization of Novel α7 Nicotinic Acetylcholine Receptor Positive Allosteric Modulators and the Discovery of a Preclinical Development Candidate Molecule (RGH-560)

During optimization of a previously identified lead compound, attempts were made to optimize the reactive indole structural element, the suboptimal metabolic stability, as well as the low kinetic solubility. It was concluded that the indole was important for in vitro activity. With the aim of further improvements, more thorough modifications were also carried out. As a result, a new chemotype (the azetidinespirochromone family) was identified, which proved to be 1 order of magnitude less lipophilic retaining the same high level of in vitro potency as the lead series itself, however, with improved metabolic stability and kinetic solubility. Compound 53 showed the most balanced physicochemical and pharmacological profile with significant in vivo efficacy in the scopolamine-induced amnesia test. Based on these promising results, cognitive enhancement through the positive modulation of α7 nAChRs appears to be a viable approach. Compound 53 was selected to be a preclinical development candidate (as RGH-560).

A Study on the Interactions of Proteinase K with Myricetin and Myricitrin by Multi-Spectroscopy and Molecular Modeling

Myricetin (MYR) and myricitrin (MYT) are well recognized for their nutraceutical value, such as antioxidant, hypoglycemic, and hypotensive effects. In this work, fluorescence spectroscopy and molecular modeling were adopted to investigate the conformational and stability changes of proteinase K (PK) in the presence of MYR and MYT. The experimental results showed that both MYR and MYT could quench fluorescence emission via a static quenching mechanism. Further investigation demonstrated that both hydrogen bonding and van der Waals forces play significant roles in the binding of complexes, which is consistent with the conclusions of molecular modeling. Synchronous fluorescence spectroscopy, Förster resonance energy transfer, and site-tagged competition experiments were performed to prove that the binding of MYR or MYT to PK could alter its micro-environment and conformation. Molecular docking results revealed that either MYR or MYT spontaneously interacted with PK at a single binding site via hydrogen bonding and hydrophobic interactions, which is consistent with the results of spectroscopic measurements. A 30 ns molecular dynamics simulation was conducted for both PK-MYR and PK-MYT complexes. The calculation results showed that no large structural distortions or interaction changes occurred during the entire simulation time span. The average RMSD changes of PK in PK-MYR and PK-MYT were 2.06 and 2.15 Å, respectively, indicating excellent stability of both complexes. The molecular simulation results suggested that both MYR and MYT could interact with PK spontaneously, which is in agreement with spectroscopic results. This agreement between experimental and theoretical results indicates that the method herein could be feasible and worthwhile for protein–ligand complex studies.

In Vitro Assessment on Designing Novel Antibiofilms of Pseudomonas aeruginosa Using a Computational Approach

An anti-biofilm that can inhibit the matrix of biofilm formation is necessary to prevent recurrent and chronic Pseudomonas aeruginosa infection. This study aimed to design compounds with a new mechanism through competitive inhibitory activity against phosphomannomutase/phosphoglucomutase (PMM/PGM), using in vitro assessment and a computational (in silico) approach. The active site of PMM/PGM was assessed through molecular redocking using L-tartaric acid as the native ligand and other small molecules, such as glucaric acid, D-sorbitol, and ascorbic acid. The docking program set the small molecules to the active site, showing a stable complex formation. Analysis of structural similarity, bioavailability, absorption, distribution, metabolism, excretion, and toxicity properties proved the potential application of ligands as an anti-biofilm. In vitro assessment with crystal violet showed that the ligands could reach up to 95.87% inhibition at different concentrations. The nitrocellulose membrane and scanning electron microscopic visualization showed that the untreated P. aeruginosa biofilm was denser than the ligand-treated biofilm.

3,4,3'-Tri-O-methylellagic acid as an anticancer agent: in vitro and in silico studies

We report a natural product compound isolated from Syzygium polycephalum known as 3,4,3'-tri-O-methylellagic acid (T-EA) as a candidate drug for cancer treatment. The characterization of the isolated T-EA compound was carried out using various spectroscopic methods. The in vitro evaluation showcased the inhibition activity of T-EA towards the T47D and HeLa cell lines with EC₅₀ values of 55.35 ± 6.28 μg mL^-1 and 12.57 ± 2.22 μg mL^-1, respectively. Meanwhile, the in silico evaluation aimed to understand the interaction of T-EA with enzymes responsible for cancer regulation at the molecular level by targeting the hindrance of cyclin-dependent kinase 9 (CDK9) and sirtuin 1 (SIRT1) enzymes. T-EA showed a binding free energy towards the SIRT1 protein of ΔG _{bind (MM-GBSA)}: -30.98 ± 0.25 kcal mol^-1 and ΔG _{bind (MM-PBSA)}: -24.07 ± 0.30 kcal mol^-1, while that of CDK9 was ΔG _{bind (MM-GBSA)}: -29.50 ± 0.22 kcal mol^-1 and ΔG _{bind (MM-PBSA)}: -25.87 ± 0.40 kcal mol^-1. The obtained results from this research could be considered as important information on 3,4,3'-tri-O-methylellagic acid as a drug to treat cervical and breast cancers.

A Versatile Class of 1,4,4-Trisubstituted Piperidines Block Coronavirus Replication In Vitro

There is a clear need for novel antiviral concepts to control SARS-CoV-2 infection. Based on the promising anti-coronavirus activity observed for a class of 1,4,4-trisubstituted piperidines, we here conducted a detailed analysis of the structure–activity relationship of these structurally unique inhibitors. Despite the presence of five points of diversity, the synthesis of an extensive series of analogues was readily achieved by Ugi four-component reaction from commercially available reagents. After evaluating 63 analogues against human coronavirus 229E, four of the best molecules were selected and shown to have micromolar activity against SARS-CoV-2. Since the action point was situated post virus entry and lying at the stage of viral polyprotein processing and the start of RNA synthesis, enzymatic assays were performed with CoV proteins involved in these processes. While no inhibition was observed for SARS-CoV-2 nsp12-nsp7-nsp8 polymerase, nsp14 N7-methyltransferase and nsp16/nsp10 2’-O-methyltransferase, nor the nsp3 papain-like protease, the compounds clearly inhibited the nsp5 main protease (M^pro). Although the inhibitory activity was quite modest, the plausibility of binding to the catalytic site of M^pro was established by in silico studies. Therefore, the 1,4,4-trisubstituted piperidines appear to represent a novel class of non-covalent CoV M^pro inhibitors that warrants further optimization and development.

One-Pot Synthesis and in Silico Molecular Docking Studies of Arylselanyl Hydrazides as Potential Antituberculosis Agents

The present study reports a simple two-step method for the synthesis of arylselanyl hydrazide derivatives using hypophosphorous acid and polyethylene glycol (H₃ PO₂ /PEG-400) as an alternative reducing system and hydrazine hydrate (NH₂ NH₂ ⋅xH₂ O/50-60 %). This single-vessel procedure was employed with methyl acrylate 2a and methyl bromoacetate 2b using diaryl diselenides to generate the nucleophile species to produce, respectively, 3-(arylselanyl)propane-hydrazides 4a-e and 2-(arylselanyl)acetohydrazides 5a-e with good yields by accelerating the reduction of -Se-Se- bond, when compared to available methods. The synthesized molecules are structurally similar to the isoniazid (INH). Therefore, we perform in silico molecular docking studies, using the lactoperoxidase enzyme, in order to verify whether the INH Se derivatives could interact in a similar way to INH at the active site of the mammalian enzyme. The in silico results indicated a similar type of interaction of the arylselanyl hydrazide derivatives with that of INH. In view of the similar in silico interaction of the selenium derivatives of INH, the arylselanyl hydrazide derivatives reported here should be tested against Mycobacterium tuberculosis in vitro.

Binding and Activation of Serotonergic G-Protein Coupled Receptors by the Multimodal Antidepressant Vortioxetine

G-protein coupled receptors (GPCRs) are important pharmacological targets. Despite substantial progress, important questions still remain concerning the details of activation: how can a ligand act as an agonist in one receptor but as an antagonist in a homologous receptor, and how can agonists activate a receptor despite lacking polar functional groups able to interact with helix 5 as is the case for the related adrenergic receptors? Studying vortioxetine (VXT), an important multimodal antidepressant drug, may elucidate both questions. Herein, we present a thorough in silico analysis of VXT binding to 5-HT_1A, 5-HT_1B, and 5-HT₇ receptors and compare it with available experimental data. We are able to rationalize the differential mode of action of VXT at different receptors, but also, in the case of the 5-HT_1A receptor, we observe the initial steps of activation that inform about an activation mechanism that does not involve polar interaction with helix 5. The results extend our current understanding of agonist and antagonist action at aminergic GPCRs.

The dolabellane diterpenes as potential inhibitors of the SARS-CoV-2 main protease: molecular insight of the inhibitory mechanism through computational studies

An investigation has been carried out on natural products from dolabellane derivatives to understand their potential in inhibiting the SARS-CoV-2 main protease (3CLpro) using an in silico approach. Inhibition of the 3CLpro enzyme is a promising target in stopping the replication of the SARS-CoV-2 virus through inhibition of the subsite binding pocket. The redocking process aims to determine the 3CLpro active sites. The redocking requirement showed a good pose with an RMSD value of 1.39 Å. The combination of molecular docking and MD simulation shows the results of DD13 as a candidate which had a good binding affinity (kcal mol-1) to inhibit the 3CLpro enzyme activity. Prediction of binding free energy (kcal mol-1) of DD13 using the Molecular Mechanics-Poisson Boltzmann/Generalized Born Surface Area (MM-PB/GBSA) approach shows the results ΔG bind(MM-GBSA): -52.33 ± 0.34 and ΔG bind(MM-PBSA): -43.52 ± 0.42. The key residues responsible for the inhibition mechanism are Hie41, Ser46, Met49, Asn142, Cys145, Hie163, Met165, and Gln189. Additionally, pharmacokinetic prediction recommended that DD13 had promising criteria as a drug candidate. The results demonstrated in this study provide theoretical information to obtain a potential inhibitor against the SARS-CoV-2 main protease.

Docking-generated multiple ligand poses for bootstrapping bioactivity classifying Machine Learning: Repurposing covalent inhibitors for COVID-19-related TMPRSS2 as case study

In the present work we introduce the use of multiple docked poses for bootstrapping machine learning-based QSAR modelling. Ligand-receptor contact fingerprints are implemented as descriptor variables. We implemented this method for the discovery of potential inhibitors of the serine protease enzyme TMPRSS2 involved the infectivity of coronaviruses. Several machine learners were scanned, however, Xgboost, support vector machines (SVM) and random forests (RF) were the best with testing set accuracies reaching 90%. Three potential hits were identified upon using the method to scan known untested FDA approved drugs against TMPRSS2. Subsequent molecular dynamics simulation and covalent docking supported the results of the new computational approach.

In silico drug repositioning on F508del-CFTR: A proof-of-concept study on the AIFA library

Computational drug repositioning is of growing interest to academia and industry, for its ability to rapidly screen a huge number of candidates in silico (exploiting comprehensive drug datasets) together with reduced development cost and time. The potential of drug repositioning has not been fully evaluated yet for cystic fibrosis (CF), a disease mainly caused by deletion of Phe 508 (F508del) of the cystic fibrosis transmembrane conductance regulator (CFTR) protein. F508del-CFTR is thus withheld in the endoplasmic reticulum and rapidly degraded by the ubiquitin/proteasome system. CF is still a fatal disease. Nowadays, it is treatable by some CFTR-rescuing drugs, but new-generation drugs with stronger therapeutic benefits and fewer side effects are still awaited. In this manuscript we report about the results of a pilot computational drug repositioning screening in search of F508del-CFTR-targeted drugs performed on AIFA library by means of a dedicated computational pipeline and surface plasmon resonance binding assay to experimentally validate the computational findings.

Exploration of stilbenoid trimers as potential inhibitors of sirtuin1 enzyme using a molecular docking and molecular dynamics simulation approach

A combination of molecular docking and molecular dynamics simulation (250 ns) has been carried out to study the interaction of stilbenoid trimer compounds with the SIRT1 enzyme as the target protein. SIRT1 expression regulates cellular stress responses that lead to the development of cancer. Redocking showed a good native ligand pose with an RMSD value of 1.40 Å at the receptor active site's coordinates. The molecular docking score uses a grid score functional (kcal mol⁻¹), which shows results of 1NS: 79.56, TS1: −26.83, TS2: −87.77, and TS3: −83.67. The TS2 and TS3 candidates were chosen for further analysis because they had a lower grid score than the native ligand (1NS). Furthermore, prediction of binding free energy (kcal mol⁻¹) using the Quantum Mechanics/generalized Born Surface Area (QM/MM-GBSA) method shows the results of 1NS: −31.52 ± 0.39, TS2: −58.99 ± 0.34, and TS3: −43.38 ± 0.35. These results indicate that the TS2 and TS3 compounds have good potential as inhibitors of the SIRT1 enzyme. Additionally, the amino acid residues were responsible for the inhibition mechanism through hydrogen bond interactions at the molecular level, including ASP22, PHE91, PRO11, ILE165, ASP166, and VAL230. The observations made in this study provide theoretical information for exploring the stilbenoid trimers as anticancer agents by targeting the SIRT1 enzyme.

A combination of molecular docking and molecular dynamics simulation (250 ns) has been carried out to study the interaction of stilbenoid trimer compounds with the SIRT1 enzyme as the target protein.