Streamlining Large Chemical Library Docking with Artificial Intelligence: the PyRMD2Dock Approach.

The present contribution introduces a novel computational protocol called PyRMD2Dock, which combines the Ligand-Based Virtual Screening (LBVS) tool PyRMD with the popular docking software AutoDock-GPU (AD4-GPU) to enhance the throughput of virtual screening campaigns for drug discovery. By implementing PyRMD2Dock, we demonstrate that it is possible to rapidly screen massive chemical databases and identify those with the highest predicted binding affinity to a target protein. Our benchmarking and screening experiments illustrate the predictive power and speed of PyRMD2Dock and highlight its potential to accelerate the discovery of novel drug candidates. Overall, this study showcases the value of combining AI-powered LBVS tools with docking software to enable effective and high-throughput virtual screening of ultralarge molecular databases in drug discovery. PyRMD and the PyRMD2Dock protocol are freely available on GitHub (https://github.com/cosconatilab/PyRMD) as an open-source tool.

Development of novel dipeptide nitriles as inhibitors of rhodesain of Trypanosoma brucei rhodesiense.

In this paper, we developed a new series of dipeptide nitriles that were demonstrated to be reversible rhodesain inhibitors at nanomolar level, with EC50 values against cultured T. b. brucei in the micromolar range. We also proved that our dipeptide nitriles directly bind to the active site of rhodesain acting as competitive inhibitors. Within the most interesting compounds, the dipeptide nitrile 2b showed the highest binding affinity towards rhodesain (Ki = 16 nM) coupled with a good antiparasitic activity (EC50 = 14.1 µM). Moreover, for the dipeptide nitrile 3e, which showed a Ki = 122 nM towards the trypanosomal protease, we obtained the highest antiparasitic activity (EC50 = 8.8 µM). Thus, given the obtained results both compounds could certainly represent new lead compounds for the discovery of new drugs to treat Human African Trypanosomiasis.

FEPrepare: A Web-Based Tool for Automating the Setup of Relative Binding Free Energy Calculations.

Relative binding free energy calculations in drug design are becoming a useful tool in facilitating lead binding affinity optimization in a cost- and time-efficient manner. However, they have been limited by technical challenges such as the manual creation of large numbers of input files to set up, run, and analyze free energy simulations. In this Application Note, we describe FEPrepare, a novel web-based tool, which automates the setup procedure for relative binding FEP calculations for the dual-topology scheme of NAMD, one of the major MD engines, using OPLS-AA force field topology and parameter files. FEPrepare provides the user with all necessary files needed to run a FEP/MD simulation with NAMD. FEPrepare can be accessed and used at https://feprepare.vi-seem.eu/.

PyRMD: A New Fully Automated AI-Powered Ligand-Based Virtual Screening Tool.

Artificial intelligence (AI) algorithms are dramatically redefining the current drug discovery landscape by boosting the efficiency of its various steps. Still, their implementation often requires a certain level of expertise in AI paradigms and coding. This often prevents the use of these powerful methodologies by non-expert users involved in the design of new biologically active compounds. Here, the random matrix discriminant (RMD) algorithm, a high-performance AI method specifically tailored for the identification of new ligands, was implemented in a new fully automated tool, PyRMD. This ligand-based virtual screening tool can be trained using target bioactivity data directly downloaded from the ChEMBL repository without manual intervention. The software automatically splits the available training compounds into active and inactive sets and learns the distinctive chemical features responsible for the compounds' activity/inactivity. PyRMD was designed to easily screen millions of compounds in hours through an automated workflow and intuitive input files, allowing fine tuning of each parameter of the calculation. Additionally, PyRMD features a wealth of benchmark metrics, to accurately probe the model performance, which were used here to gauge its predictive potential and limitations. PyRMD is freely available on GitHub (https://github.com/cosconatilab/PyRMD) as an open-source tool.

Quinolinonyl Non-Diketo Acid Derivatives as Inhibitors of HIV-1 Ribonuclease H and Polymerase Functions of Reverse Transcriptase.

Novel anti-HIV agents are still needed to overcome resistance issues, in particular inhibitors acting against novel viral targets. The ribonuclease H (RNase H) function of the reverse transcriptase (RT) represents a validated and promising target, and no inhibitor has reached the clinical pipeline yet. Here, we present rationally designed non-diketo acid selective RNase H inhibitors (RHIs) based on the quinolinone scaffold starting from former dual integrase (IN)/RNase H quinolinonyl diketo acids. Several derivatives were synthesized and tested against RNase H and viral replication and found active at micromolar concentrations. Docking studies within the RNase H catalytic site, coupled with site-directed mutagenesis, and Mg2+ titration experiments demonstrated that our compounds coordinate the Mg2+ cofactor and interact with amino acids of the RNase H domain that are highly conserved among naïve and treatment-experienced patients. In general, the new inhibitors influenced also the polymerase activity of RT but were selective against RNase H vs the IN enzyme.

Inhibition studies on carbonic anhydrase isoforms I, II, IV and IX with N-arylsubstituted secondary sulfonamides featuring a bicyclic tetrahydroindazole scaffold.

Carbonic Anhydrases (CAs) are pharmaceutically relevant targets for the treatment of several disease conditions. The ubiquitous localization of these enzymes and the high homology shared by the different isoforms represent substantial impediments for the discovery of potential drugs devoid of off-target side effects. As a consequence, substantial efforts are still needed to allow for the full realization of the pharmacological potential of CA modulators. In this contribution, starting from our previous studies, we describe the synthesis of a set of new bicyclic tetrahydroindazoles featuring a secondary sulfonamide. Biological evaluation of the inhibitory activity against the hCA I, II, IV, and IX isoforms allowed drawing a structure-activity relationship profile that was rationalized through theoretical studies. This allowed dissecting the new molecules into the single portions influencing the zinc chelation properties and the selectivity profile thereby offering a new platform for the discovery of new isotype selective CA inhibitors.

Lead Discovery of SARS-CoV-2 Main Protease Inhibitors through Covalent Docking-Based Virtual Screening.

During almost all 2020, coronavirus disease 2019 (COVID-19) pandemic has constituted the major risk for the worldwide health and economy, propelling unprecedented efforts to discover drugs for its prevention and cure. At the end of the year, these efforts have culminated with the approval of vaccines by the American Food and Drug Administration (FDA) and the European Medicines Agency (EMA) giving new hope for the future. On the other hand, clinical data underscore the urgent need for effective drugs to treat COVID-19 patients. In this work, we embarked on a virtual screening campaign against the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Mpro chymotrypsin-like cysteine protease employing our in-house database of peptide and non-peptide ligands characterized by different types of warheads acting as Michael acceptors. To this end, we employed the AutoDock4 docking software customized to predict the formation of a covalent adduct with the target protein. In vitro verification of the inhibition properties of the most promising candidates allowed us to identify two new lead inhibitors that will deserve further optimization. From the computational point of view, this work demonstrates the predictive power of AutoDock4 and suggests its application for the in silico screening of large chemical libraries of potential covalent binders against the SARS-CoV-2 Mpro enzyme.

Structure-Activity Relationship Studies on Oxazolo[3,4-a]pyrazine Derivatives Leading to the Discovery of a Novel Neuropeptide S Receptor Antagonist with Potent In Vivo Activity.

Neuropeptide S modulates important neurobiological functions including locomotion, anxiety, and drug abuse through interaction with its G protein-coupled receptor known as neuropeptide S receptor (NPSR). NPSR antagonists are potentially useful for the treatment of substance abuse disorders against which there is an urgent need for new effective therapeutic approaches. Potent NPSR antagonists in vitro have been discovered which, however, require further optimization of their in vivo pharmacological profile. This work describes a new series of NPSR antagonists of the oxazolo[3,4-a]pyrazine class. The guanidine derivative 16 exhibited nanomolar activity in vitro and 5-fold improved potency in vivo compared to SHA-68, a reference pharmacological tool in this field. Compound 16 can be considered a new tool for research studies on the translational potential of the NPSergic system. An in-depth molecular modeling investigation was also performed to gain new insights into the observed structure-activity relationships and provide an updated model of ligand/NPSR interactions.

Peptidyl Vinyl Ketone Irreversible Inhibitors of Rhodesain: Modifications of the P2 Fragment.

In this paper, we report the design, synthesis and biological investigation of a series of peptidyl vinyl ketones obtained by modifying the P2 fragment of previously reported highly potent inhibitors of rhodesain, the main cysteine protease of Trypanosoma brucei rhodesiense. Investigation of the structure-activity relationship led us to identify new rhodesain inhibitors endowed with an improved selectivity profile (a selectivity index of up to 22 000 towards the target enzyme), and/or an improved antitrypanosomal activity in the sub-micromolar range.

Pyrrolyl Pyrazoles as Non-Diketo Acid Inhibitors of the HIV-1 Ribonuclease H Function of Reverse Transcriptase.

Due to the biological liability of diketo acid (DKA) chain, we transferred this element of our previously reported anti-HIV-1 pyrrolyl derivatives to a non-DKA scaffold, obtaining a series of pyrrolyl-pyrazole carboxylic acids as new RNase H inhibitors. Among the newly synthesized derivatives, oxyphenylpyrrolyl-pyrazoles demonstrated inhibitory activities within the low micromolar/submicromolar range with compound 11b being the most potent. Interestingly, all tested compounds showed up to 2 orders of magnitude of selectivity for RNase H vs integrase. Docking studies within the RNase H catalytic site, coupled with site-directed mutagenesis, showed the key structural features that could confer the ability to establish specific interactions within RNase H. Furthermore, they proved the ability of our compounds to interact with amino acids highly conserved among HIV-1 subspecies isolated among patients carrying drug-resistant variants. In the end, the newly discovered pyrazole carboxylic acid derivatives feature promising serum stability with respect to their corresponding DKAs.

Development of Novel Benzodiazepine-Based Peptidomimetics as Inhibitors of Rhodesain from Trypanosoma brucei rhodesiense.

Starting from the reversible rhodesain inhibitors 1 a-c, which have Ki values towards the target protease in the low-micromolar range, we have designed a series of peptidomimetics, 2 a-g, that contain a benzodiazepine scaffold as a ß-turn mimetic; they are characterized by a specific peptide sequence for the inhibition of rhodesain. Considering that irreversible inhibition is strongly desirable in the case of a parasitic target, a vinyl ester moiety acting as Michael-acceptor was introduced as the warhead; this portion was functionalized in order to evaluate the size of corresponding enzyme pocket that could accommodate this substituent. With this investigation, we identified an irreversible rhodesain inhibitor (i. e., 2 g) with a k2nd value of 90 000 M-1 min-1 that showed antitrypanosomal activity in the low-micromolar range (EC50 =1.25 µM), this may be considered a promising lead compound in the drug-discovery process for treating human African trypanosomiasis (HAT).

Optimization Strategy of Novel Peptide-Based Michael Acceptors for the Treatment of Human African Trypanosomiasis.

This paper describes an optimization strategy of the highly active vinyl ketone 3 which was recognized as a strong inhibitor of rhodesain of Trypanosoma brucei rhodesiense, endowed with a ksecond value of 67 × 106 M-1 min-1 coupled with a high binding affinity (Ki = 38 pM). We now report a new structure-activity relationship study based on structural variations on the P3, P2, and P1' sites which led us to identify two potent lead compounds, i.e., vinyl ketones 4h and 4k. Vinyl ketone 4h showed an impressive potency toward rhodesain (ksecond = 8811 × 105) coupled to a good antiparasitic activity (EC50 = 3.6 µM), while vinyl ketone 4k proved to possess the highest binding affinity toward the trypanosomal protease (Ki = 0.6 pM) and a submicromolar antiparasitic activity (EC50 = 0.67 µM), thus representing new lead compounds in the drug discovery process for the treatment of Human African Trypanosomiasis.

Novel 2-substituted-benzimidazole-6-sulfonamides as carbonic anhydrase inhibitors: synthesis, biological evaluation against isoforms I, II, IX and XII and molecular docking studies.

Inhibition of Carbonic Anhydrases (CAs) has been clinically exploited for many decades for a variety of therapeutic applications. Within a research project aimed at developing novel classes of CA inhibitors (CAIs) with a proper selectivity for certain isoforms, a series of derivatives featuring the 2-substituted-benzimidazole-6-sulfonamide scaffold, conceived as frozen analogs of Schiff bases and secondary amines previously reported in the literature as CAIs, were investigated. Enzyme inhibition assays on physiologically relevant human CA I, II, IX and XII isoforms revealed a number of potent CAIs, showing promising selectivity profiles towards the transmembrane tumor-associated CA IX and XII enzymes. Computational studies were attained to clarify the structural determinants behind the activities and selectivity profiles of the novel inhibitors.

Identification of novel indole derivatives acting as inhibitors of the Keap1-Nrf2 interaction.

Nine indole derivatives (9a-i) were tested as potential inhibitors of the Keap1-Nrf2 interaction. This class of compounds increases the intracellular levels of the transcription factor Nrf2 and the consequent expression of enzymes encoded by genes containing the antioxidant response element (ARE). In the ARE-luciferase reporter assay only 9e-g revealed to be remarkably more active than t-butylhydroxyquinone (t-BHQ), with 9g standing out as the best performing compound. While 9e and 9f are weak acids, 9g is an ampholyte prevailing as a zwitterion in neutral aqueous solutions. The ability of 9e-g to significantly increase levels of Nrf2, NADPH:quinone oxidoreductase 1, and transketolase (TKT) gave further support to the hypothesis that these compounds act as inhibitors of the Keap1-Nrf2 interaction. Docking simulations allowed us to elucidate the nature of the putative interactions between 9g and Keap1.

Discovery of Pyrido[3',2':5,6]thiopyrano[4,3-d]pyrimidine-Based Antiproliferative Multikinase Inhibitors.

Protein kinases dysregulation is extremely common in cancer cells, and the development of new agents able to simultaneously target multiple kinase pathways involved in angiogenesis and tumor growth may offer several advantages in the treatment of cancer. Herein we report the discovery of new pyridothiopyranopyrimidine derivatives (2-4) showing high potencies in VEGFR-2 KDR inhibition as well as antiproliferative effect on a panel of human tumor cell lines. Investigation on the selectivity profile of the representative 2-anilino-substituted compounds 3b, 3i, and 3j revealed a multiplicity of kinase targets that should account for the potent antiproliferative effect produced by these pyridothiopyranopyrimidine derivatives.

Discovery of a Novel Chemotype of Histone Lysine Methyltransferase EHMT1/2 (GLP/G9a) Inhibitors: Rational Design, Synthesis, Biological Evaluation, and Co-crystal Structure.

Since the discovery of compound BIX01294 over 10 years ago, only a very limited number of nonquinazoline inhibitors of H3K9-specific methyltransferases G9a and G9a-like protein (GLP) have been reported. Herein, we report the identification of a novel chemotype for G9a/GLP inhibitors, based on the underinvestigated 2-alkyl-5-amino- and 2-aryl-5-amino-substituted 3 H-benzo[ e][1,4]diazepine scaffold. Our research efforts resulted in the identification 12a (EML741), which not only maintained the high in vitro and cellular potency of its quinazoline counterpart, but also displayed improved inhibitory potency against DNA methyltransferase 1, improved selectivity against other methyltransferases, low cell toxicity, and improved apparent permeability values in both parallel artificial membrane permeability assay (PAMPA) and blood-brain barrier-specific PAMPA, and therefore might potentially be a better candidate for animal studies. Finally, the co-crystal structure of GLP in complex with 12a provides the basis for the further development of benzodiazepine-based G9a/GLP inhibitors.

4-Substituted Benzenesulfonamides Incorporating Bi/Tricyclic Moieties Act as Potent and Isoform-Selective Carbonic Anhydrase II/IX Inhibitors.

As a part of our efforts to expand chemical diversity in the carbonic anhydrases inhibitors (CAIs), three small series of polyheterocyclic compounds (4-6) featuring the primary benzenesulfonamide moiety linked to bi/tricyclic scaffolds were investigated. Highly effective inhibitors against the target tumor-associated hCA IX (low nanomolar/subnanomolar potency levels) showing significant functional selectivity profile toward hCA I, II, and IV isozymes were identified. Molecular docking studies clarified the reasons behind the activity and selectivity of the new compounds.

New insights in the structure-activity relationships of 2-phenylamino-substituted benzothiopyrano[4,3-d]pyrimidines as kinase inhibitors.

Inhibition of angiogenesis via blocking vascular endothelial growth factor receptor (VEGFR) signaling pathway emerged as an established approach in anticancer therapy. So far, many monoclonal antibodies and ATP-competitive small molecule inhibitors have been clinically validated and approved. In this study, structure-activity relationships (SAR) within the 2-phenylamino-substituted benzothiopyrano[4,3-d]pyrimidine class of kinase inhibitors were further refined by the synthesis and biological evaluation of new compounds 1-21 featuring different substitution patterns on the pendant phenyl moiety, combined with H, OCH3, or Cl at 8-position. Most compounds showed a promising human kinase insert domain receptor (KDR) inhibition profile, with IC50 values in the submicromolar/low micromolar range, and promising antiproliferative activity on human umbilical vein endothelial cells (HUVECs) as well as on a panel of three human tumor cell lines. The angio-kinase selectivity profile was assessed for the most promising compound 16 against a set of six human kinases. Finally, computational studies allowed clarifying at molecular level the interaction pattern established by the compounds with KDR, highlighting key stable cation-π interactions, and thus providing the basis for further designing novel inhibitors.

Dual MET and SMO Negative Modulators Overcome Resistance to EGFR Inhibitors in Human Nonsmall Cell Lung Cancer.

Tyrosine kinase inhibitors (TKIs) of the EGF receptor (EGFR) have provided a significant improvement in the disease outcome of nonsmall cell lung cancer (NSCLC). Unfortunately, resistance to these agents frequently occurs, and it is often related to the activation of the Hedgehog (Hh) and MET signaling cascades driving the epithelial-to-mesenchymal transition (EMT). Because the concomitant inhibition of both Hh and MET pathways restores the sensitivity to anti-EGFR drugs, here we aimed at discovering the first compounds that block simultaneously MET and SMO. By using an "in silico drug repurposing" approach and by validating our predictions both in vitro and in vivo, we identified a set of compounds with the desired dual inhibitory activity and enhanced antiproliferative activity on EGFR TKI-resistant NSCLC. The identification of the known MET TKIs, glesatinib and foretinib, as negative modulators of the Hh pathway, widens their application in the context of NSCLC.

Development of Novel Peptide-Based Michael Acceptors Targeting Rhodesain and Falcipain-2 for the Treatment of Neglected Tropical Diseases (NTDs).

This paper describes the development of a class of peptide-based inhibitors as novel antitrypanosomal and antimalarial agents. The inhibitors are based on a characteristic peptide sequence for the inhibition of the cysteine proteases rhodesain of Trypanosoma brucei rhodesiense and falcipain-2 of Plasmodium falciparum. We exploited the reactivity of novel unsaturated electrophilic functions such as vinyl-sulfones, -ketones, -esters, and -nitriles. The Michael acceptors inhibited both rhodesain and falcipain-2, at nanomolar and micromolar levels, respectively. In particular, the vinyl ketone 3b has emerged as a potent rhodesain inhibitor (k2nd = 67 × 106 M-1 min-1), endowed with a picomolar binding affinity (Ki = 38 pM), coupled with a single-digit micromolar activity against Trypanosoma brucei brucei (EC50 = 2.97 µM), thus being considered as a novel lead compound for the discovery of novel effective antitrypanosomal agents.