A comprehensive study of SARS-CoV-2 main protease (Mpro) inhibitor-resistant mutants selected in a VSV-based system.

Nirmatrelvir was the first protease inhibitor specifically developed against the SARS-CoV-2 main protease (3CLpro/Mpro) and licensed for clinical use. As SARS-CoV-2 continues to spread, variants resistant to nirmatrelvir and other currently available treatments are likely to arise. This study aimed to identify and characterize mutations that confer resistance to nirmatrelvir. To safely generate Mpro resistance mutations, we passaged a previously developed, chimeric vesicular stomatitis virus (VSV-Mpro) with increasing, yet suboptimal concentrations of nirmatrelvir. Using Wuhan-1 and Omicron Mpro variants, we selected a large set of mutants. Some mutations are frequently present in GISAID, suggesting their relevance in SARS-CoV-2. The resistance phenotype of a subset of mutations was characterized against clinically available protease inhibitors (nirmatrelvir and ensitrelvir) with cell-based, biochemical and SARS-CoV-2 replicon assays. Moreover, we showed the putative molecular mechanism of resistance based on in silico molecular modelling. These findings have implications on the development of future generation Mpro inhibitors, will help to understand SARS-CoV-2 protease inhibitor resistance mechanisms and show the relevance of specific mutations, thereby informing treatment decisions.

CG7630 is the Drosophila melanogaster homolog of the cytochrome c oxidase subunit COX7B.

The mitochondrial respiratory chain (MRC) is composed of four multiheteromeric enzyme complexes. According to the endosymbiotic origin of mitochondria, eukaryotic MRC derives from ancestral proteobacterial respiratory structures consisting of a minimal set of complexes formed by a few subunits associated with redox prosthetic groups. These enzymes, which are the "core" redox centers of respiration, acquired additional subunits, and increased their complexity throughout evolution. Cytochrome c oxidase (COX), the terminal component of MRC, has a highly interspecific heterogeneous composition. Mammalian COX consists of 14 different polypeptides, of which COX7B is considered the evolutionarily youngest subunit. We applied proteomic, biochemical, and genetic approaches to investigate the COX composition in the invertebrate model Drosophila melanogaster. We identified and characterized a novel subunit which is widely different in amino acid sequence, but similar in secondary and tertiary structures to COX7B, and provided evidence that this object is in fact replacing the latter subunit in virtually all protostome invertebrates. These results demonstrate that although individual structures may differ the composition of COX is functionally conserved between vertebrate and invertebrate species.

FAM/USP9x, a deubiquitinating enzyme essential for TGFbeta signaling, controls Smad4 monoubiquitination.

The assembly of the Smad complex is critical for TGFbeta signaling, yet the mechanisms that inactivate or empower nuclear Smad complexes are less understood. By means of siRNA screen we identified FAM (USP9x), a deubiquitinase acting as essential and evolutionarily conserved component in TGFbeta and bone morphogenetic protein signaling. Smad4 is monoubiquitinated in lysine 519 in vivo, a modification that inhibits Smad4 by impeding association with phospho-Smad2. FAM reverts this negative modification, re-empowering Smad4 function. FAM opposes the activity of Ectodermin/Tif1gamma (Ecto), a nuclear factor for which we now clarify a prominent role as Smad4 monoubiquitin ligase. Our study points to Smad4 monoubiquitination and deubiquitination as a way for cells to set their TGFbeta responsiveness: loss of FAM disables Smad4-dependent responses in several model systems, with Ecto being epistatic to FAM. This defines a regulative ubiquitination step controlling Smads that is parallel to those impinging on R-Smad phosphorylation.

7-Amino-[1,2,4]triazolo[1,5-a][1,3,5]triazines as CK1δ inhibitors: Exploring substitutions at the 2 and 5-positions.

CK1δ is a serine-threonine kinase involved in several pathological conditions including neuroinflammation and neurodegenerative disorders like Alzheimer's disease, Parkinson's disease, and Amyotrophic Lateral Sclerosis. Specifically, it seems that an inhibition of CK1δ could have a neuroprotective effect in these conditions. Here, a series of [1,2,4]triazolo[1,5-a][1,3,5]triazines were developed as ATP-competitive CK1δ inhibitors. Both positions 2 and 5 have been explored leading to a total of ten compounds exhibiting IC50s comprised between 29.1 µM and 2.08 µM. Three of the four most potent compounds (IC50 < 3 µM) bear a thiophene ring at the 2 position. All compounds have been submitted to computational studies that identified the chain composed of at least 2 atoms (e.g., nitrogen and carbon atoms) at the 5 position as crucial to determine a key bidentate hydrogen bond with Leu85 of CK1δ. Most potent compounds have been tested in vitro, resulting passively permeable to the blood-brain barrier and, safe and slight neuroprotective on a neuronal cell model. These results encourage to further structural optimize the series to obtain more potent CK1δ inhibitors as possible neuroprotective agents to be tested on models of the above-mentioned neurodegenerative diseases.

TumFlow: An AI Model for Predicting New Anticancer Molecules.

Melanoma is the fifth most common cancer in the United States. Conventional drug discovery methods are inherently time-consuming and costly, which imposes significant limitations. However, the advent of Artificial Intelligence (AI) has opened up new possibilities for simulating and evaluating numerous drug candidates, thereby mitigating the requisite time and resources. In this context, normalizing flow models by employing machine learning techniques to create new molecular structures holds promise for accelerating the discovery of effective anticancer therapies. This manuscript introduces TumFlow, a novel AI model designed to generate new molecular entities with potential therapeutic value in cancer treatment. It has been trained on the NCI-60 dataset, encompassing thousands of molecules tested across 60 tumour cell lines, with an emphasis on the melanoma SK-MEL-28 cell line. The model successfully generated new molecules with predicted improved efficacy in inhibiting tumour growth while being synthetically feasible. This represents a significant advancement over conventional generative models, which often produce molecules that are challenging or impossible to synthesize. Furthermore, TumFlow has also been utilized to optimize molecules known for their efficacy in clinical melanoma treatments. This led to the creation of novel molecules with a predicted enhanced likelihood of effectiveness against melanoma, currently undocumented on PubChem.

Fighting Antimicrobial Resistance: Insights on How the Staphylococcus aureus NorA Efflux Pump Recognizes 2-Phenylquinoline Inhibitors by Supervised Molecular Dynamics (SuMD) and Molecular Docking Simulations.

The superbug Staphylococcus aureus (S. aureus) exhibits several resistance mechanisms, including efflux pumps, that strongly contribute to antimicrobial resistance. In particular, the NorA efflux pump activity is associated with S. aureus resistance to fluoroquinolone antibiotics (e.g., ciprofloxacin) by promoting their active extrusion from cells. Thus, since efflux pump inhibitors (EPIs) are able to increase antibiotic concentrations in bacteria as well as restore their susceptibility to these agents, they represent a promising strategy to counteract bacterial resistance. Additionally, the very recent release of two NorA efflux pump cryo-electron microscopy (cryo-EM) structures in complex with synthetic antigen-binding fragments (Fabs) represents a real breakthrough in the study of S. aureus antibiotic resistance. In this scenario, supervised molecular dynamics (SuMD) and molecular docking experiments were combined to investigate for the first time the molecular mechanisms driving the interaction between NorA and efflux pump inhibitors (EPIs), with the ultimate goal of elucidating how the NorA efflux pump recognizes its inhibitors. The findings provide insights into the dynamic NorA-EPI intermolecular interactions and lay the groundwork for future drug discovery efforts aimed at the identification of novel molecules to fight antimicrobial resistance.

Structural Manipulations of Marine Natural Products Inspire a New Library of 3-Amino-1,2,4-Triazine PDK Inhibitors Endowed with Antitumor Activity in Pancreatic Ductal Adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is one of the main aggressive types of cancer, characterized by late prognosis and drug resistance. Among the main factors sustaining PDAC progression, the alteration of cell metabolism has emerged to have a key role in PDAC cell proliferation, invasion, and resistance to standard chemotherapeutic agents. Taking into account all these factors and the urgency in evaluating novel options to treat PDAC, in the present work we reported the synthesis of a new series of indolyl-7-azaindolyl triazine compounds inspired by marine bis-indolyl alkaloids. We first assessed the ability of the new triazine compounds to inhibit the enzymatic activity of pyruvate dehydrogenase kinases (PDKs). The results showed that most of derivatives totally inhibit PDK1 and PDK4. Molecular docking analysis was executed to predict the possible binding mode of these derivatives using ligand-based homology modeling technique. Evaluation of the capability of new triazines to inhibit the cell growth in 2D and 3D KRAS-wild-type (BxPC-3) and KRAS-mutant (PSN-1) PDAC cell line, was carried out. The results showed the capacity of the new derivatives to reduce cell growth with a major selectivity against KRAS-mutant PDAC PSN-1 on both cell models. These data demonstrated that the new triazine derivatives target PDK1 enzymatic activity and exhibit cytotoxic effects on 2D and 3D PDAC cell models, thus encouraging further structure manipulation for analogs development against PDAC.

Lessons Learnt from COVID-19: Computational Strategies for Facing Present and Future Pandemics.

Since its outbreak in December 2019, the COVID-19 pandemic has caused the death of more than 6.5 million people around the world. The high transmissibility of its causative agent, the SARS-CoV-2 virus, coupled with its potentially lethal outcome, provoked a profound global economic and social crisis. The urgency of finding suitable pharmacological tools to tame the pandemic shed light on the ever-increasing importance of computer simulations in rationalizing and speeding up the design of new drugs, further stressing the need for developing quick and reliable methods to identify novel active molecules and characterize their mechanism of action. In the present work, we aim at providing the reader with a general overview of the COVID-19 pandemic, discussing the hallmarks in its management, from the initial attempts at drug repurposing to the commercialization of Paxlovid, the first orally available COVID-19 drug. Furthermore, we analyze and discuss the role of computer-aided drug discovery (CADD) techniques, especially those that fall in the structure-based drug design (SBDD) category, in facing present and future pandemics, by showcasing several successful examples of drug discovery campaigns where commonly used methods such as docking and molecular dynamics have been employed in the rational design of effective therapeutic entities against COVID-19.

Targeting the I7L Protease: A Rational Design for Anti-Monkeypox Drugs?

The latest monkeypox virus outbreak in 2022 showcased the potential threat of this viral zoonosis to public health. The lack of specific treatments against this infection and the success of viral protease inhibitors-based treatments against HIV, Hepatitis C, and SARS-CoV-2, brought the monkeypox virus I7L protease under the spotlight as a potential target for the development of specific and compelling drugs against this emerging disease. In the present work, the structure of the monkeypox virus I7L protease was modeled and thoroughly characterized through a dedicated computational study. Furthermore, structural information gathered in the first part of the study was exploited to virtually screen the DrugBank database, consisting of drugs approved by the Food and Drug Administration (FDA) and clinical-stage drug candidates, in search for readily repurposable compounds with similar binding features as TTP-6171, the only non-covalent I7L protease inhibitor reported in the literature. The virtual screening resulted in the identification of 14 potential inhibitors of the monkeypox I7L protease. Finally, based on data collected within the present work, some considerations on developing allosteric modulators of the I7L protease are reported.

Thermal Titration Molecular Dynamics (TTMD): Not Your Usual Post-Docking Refinement.

Molecular docking is one of the most widely used computational approaches in the field of rational drug design, thanks to its favorable balance between the rapidity of execution and the accuracy of provided results. Although very efficient in exploring the conformational degrees of freedom available to the ligand, docking programs can sometimes suffer from inaccurate scoring and ranking of generated poses. To address this issue, several post-docking filters and refinement protocols have been proposed throughout the years, including pharmacophore models and molecular dynamics simulations. In this work, we present the first application of Thermal Titration Molecular Dynamics (TTMD), a recently developed method for the qualitative estimation of protein-ligand unbinding kinetics, to the refinement of docking results. TTMD evaluates the conservation of the native binding mode throughout a series of molecular dynamics simulations performed at progressively increasing temperatures with a scoring function based on protein-ligand interaction fingerprints. The protocol was successfully applied to retrieve the native-like binding pose among a set of decoy poses of drug-like ligands generated on four different pharmaceutically relevant biological targets, including casein kinase 1δ, casein kinase 2, pyruvate dehydrogenase kinase 2, and SARS-CoV-2 main protease.

Molecular Modeling Unveils the Effective Interaction of B-RAF Inhibitors with Rare B-RAF Insertion Variants.

The Food and Drug Administration (FDA) has approved MAPK inhibitors as a treatment for melanoma patients carrying a mutation in codon V600 of the BRAF gene exclusively. However, BRAF mutations outside the V600 codon may occur in a small percentage of melanomas. Although these rare variants may cause B-RAF activation, their predictive response to B-RAF inhibitor treatments is still poorly understood. We exploited an integrated approach for mutation detection, tumor evolution tracking, and assessment of response to treatment in a metastatic melanoma patient carrying the rare p.T599dup B-RAF mutation. He was addressed to Dabrafenib/Trametinib targeted therapy, showing an initial dramatic response. In parallel, in-silico ligand-based homology modeling was set up and performed on this and an additional B-RAF rare variant (p.A598_T599insV) to unveil and justify the success of the B-RAF inhibitory activity of Dabrafenib, showing that it could adeptly bind both these variants in a similar manner to how it binds and inhibits the V600E mutant. These findings open up the possibility of broadening the spectrum of BRAF inhibitor-sensitive mutations beyond mutations at codon V600, suggesting that B-RAF V600 WT melanomas should undergo more specific investigations before ruling out the possibility of targeted therapy.

Discovery of the 3-Amino-1,2,4-triazine-Based Library as Selective PDK1 Inhibitors with Therapeutic Potential in Highly Aggressive Pancreatic Ductal Adenocarcinoma.

Pyruvate dehydrogenase kinases (PDKs) are serine/threonine kinases, that are directly involved in altered cancer cell metabolism, resulting in cancer aggressiveness and resistance. Dichloroacetic acid (DCA) is the first PDK inhibitor that has entered phase II clinical; however, several side effects associated with weak anticancer activity and excessive drug dose (100 mg/kg) have led to its limitation in clinical application. Building upon a molecular hybridization approach, a small library of 3-amino-1,2,4-triazine derivatives has been designed, synthesized, and characterized for their PDK inhibitory activity using in silico, in vitro, and in vivo assays. Biochemical screenings showed that all synthesized compounds are potent and subtype-selective inhibitors of PDK. Accordingly, molecular modeling studies revealed that a lot of ligands can be properly placed inside the ATP-binding site of PDK1. Interestingly, 2D and 3D cell studies revealed their ability to induce cancer cell death at low micromolar doses, being extremely effective against human pancreatic KRAS mutated cancer cells. Cellular mechanistic studies confirm their ability to hamper the PDK/PDH axis, thus leading to metabolic/redox cellular impairment, and to ultimately trigger apoptotic cancer cell death. Remarkably, preliminary in vivo studies performed on a highly aggressive and metastatic Kras-mutant solid tumor model confirm the ability of the most representative compound 5i to target the PDH/PDK axis in vivo and highlighted its equal efficacy and better tolerability profile with respect to those elicited by the reference FDA approved drugs, cisplatin and gemcitabine. Collectively, the data highlights the promising anticancer potential of these novel PDK-targeting derivatives toward obtaining clinical candidates for combatting highly aggressive KRAS-mutant pancreatic ductal adenocarcinomas.

Correction: Bassani, D.; Moro, S. Past, Present, and Future Perspectives on Computer-Aided Drug Design Methodologies. Molecules 2023, 28, 3906.

In the published publication [...].

Past, Present, and Future Perspectives on Computer-Aided Drug Design Methodologies.

The application of computational approaches in drug discovery has been consolidated in the last decades. These families of techniques are usually grouped under the common name of "computer-aided drug design" (CADD), and they now constitute one of the pillars in the pharmaceutical discovery pipelines in many academic and industrial environments. Their implementation has been demonstrated to tremendously improve the speed of the early discovery steps, allowing for the proficient and rational choice of proper compounds for a desired therapeutic need among the extreme vastness of the drug-like chemical space. Moreover, the application of CADD approaches allows the rationalization of biochemical and interactive processes of pharmaceutical interest at the molecular level. Because of this, computational tools are now extensively used also in the field of rational 3D design and optimization of chemical entities starting from the structural information of the targets, which can be experimentally resolved or can also be obtained with other computer-based techniques. In this work, we revised the state-of-the-art computer-aided drug design methods, focusing on their application in different scenarios of pharmaceutical and biological interest, not only highlighting their great potential and their benefits, but also discussing their actual limitations and eventual weaknesses. This work can be considered a brief overview of computational methods for drug discovery.

Drugging the Undruggable Trypanosoma brucei Monothiol Glutaredoxin 1.

Trypanosoma brucei is a species of kinetoplastid causing sleeping sickness in humans and nagana in cows and horses. One of the peculiarities of this species of parasites is represented by their redox metabolism. One of the proteins involved in this redox machinery is the monothiol glutaredoxin 1 (1CGrx1) which is characterized by a unique disordered N-terminal extension exclusively conserved in trypanosomatids and other organisms. This region modulates the binding profile of the glutathione/trypanothione binding site, one of the functional regions of 1CGrx1. No endogenous ligands are known to bind this protein which does not present well-shaped binding sites, making it target particularly challenging to target. With the aim of targeting this peculiar system, we carried out two different screenings: (i) a fragment-based lead discovery campaign directed to the N-terminal as well as to the canonical binding site of 1CGrx1; (ii) a structure-based virtual screening directed to the 1CGrx1 canonical binding site. Here we report a small molecule that binds at the glutathione binding site in which the binding mode of the molecule was deeply investigated by Nuclear Magnetic Resonance (NMR). This compound represents an important step in the attempt to develop a novel strategy to interfere with the peculiar Trypanosoma Brucei redox system, making it possible to shed light on the perturbation of this biochemical machinery and eventually to novel therapeutic possibilities.

Qualitative Estimation of Protein-Ligand Complex Stability through Thermal Titration Molecular Dynamics Simulations.

The prediction of ligand efficacy has long been linked to thermodynamic properties such as the equilibrium dissociation constant, which considers both the association and the dissociation rates of a defined protein-ligand complex. In the last 15 years, there has been a paradigm shift, with an increased interest in the determination of kinetic properties such as the drug-target residence time since they better correlate with ligand efficacy compared to other parameters. In this article, we present thermal titration molecular dynamics (TTMD), an alternative computational method that combines a series of molecular dynamics simulations performed at progressively increasing temperatures with a scoring function based on protein-ligand interaction fingerprints for the qualitative estimation of protein-ligand-binding stability. The protocol has been applied to four different pharmaceutically relevant test cases, including protein kinase CK1δ, protein kinase CK2, pyruvate dehydrogenase kinase 2, and SARS-CoV-2 main protease, on a variety of ligands with different sizes, structures, and experimentally determined affinity values. In all four cases, TTMD was successfully able to distinguish between high-affinity compounds (low nanomolar range) and low-affinity ones (micromolar), proving to be a useful screening tool for the prioritization of compounds in a drug discovery campaign.

From the Wuhan-Hu-1 strain to the XD and XE variants: is targeting the SARS-CoV-2 spike protein still a pharmaceutically relevant option against COVID-19?

Since the outbreak of the COVID-19 pandemic in December 2019, the SARS-CoV-2 genome has undergone several mutations. The emergence of such variants has resulted in multiple pandemic waves, contributing to sustaining to date the number of infections, hospitalisations, and deaths despite the swift development of vaccines, since most of these mutations are concentrated on the Spike protein, a viral surface glycoprotein that is the main target for most vaccines. A milestone in the fight against the COVID-19 pandemic has been represented by the development of Paxlovid, the first orally available drug against COVID-19, which acts on the Main Protease (Mpro). In this article, we analyse the structural features of both the Spike protein and the Mpro of the recently reported SARS-CoV-2 variant XE, as well the closely related XD and XF ones, discussing their impact on the efficacy of existing treatments against COVID-19 and on the development of future ones.

Bat coronaviruses related to SARS-CoV-2: what about their 3CL proteases (MPro)?

Despite a huge effort by the scientific community to determine the animal reservoir of SARS-CoV-2, which led to the identification of several SARS-CoV-2-related viruses both in bats and in pangolins, the origin of SARS-CoV-2 is still not clear. Recently, Temmam et al. reported the discovery of bat coronaviruses with a high degree of genome similarity with SARS-CoV-2, especially concerning the RBDs of the S protein, which mediates the capability of such viruses to enter and therefore infect human cells through a hACE2-dependent pathway. These viruses, especially the one named BANAL-236, showed a higher affinity for the hACE2 compared to the original strain of SARS-CoV-2. In the present work, we analyse the similarities and differences between the 3CL protease (main protease, Mpro) of these newly reported viruses and SARS-CoV-2, discussing their relevance relative to the efficacy of existing therapeutic approaches against COVID-19, particularly concerning the recently approved orally available Paxlovid, and the development of future ones.

The Multifaceted Role of GPCRs in Amyotrophic Lateral Sclerosis: A New Therapeutic Perspective?

Amyotrophic lateral sclerosis (ALS) is a degenerating disease involving the motor neurons, which causes a progressive loss of movement ability, usually leading to death within 2 to 5 years from the diagnosis. Much effort has been put into research for an effective therapy for its eradication, but still, no cure is available. The only two drugs approved for this pathology, Riluzole and Edaravone, are onlyable to slow down the inevitable disease progression. As assessed in the literature, drug targets such as protein kinases have already been extensively examined as potential drug targets for ALS, with some molecules already in clinical trials. Here, we focus on the involvement of another very important and studied class of biological entities, G protein-coupled receptors (GPCRs), in the onset and progression of ALS. This workaimsto give an overview of what has been already discovered on the topic, providing useful information and insights that can be used by scientists all around the world who are putting efforts into the fight against this very important neurodegenerating disease.

The Effect of C-Phycocyanin on Microglia Activation Is Mediated by Toll-like Receptor 4.

The blue-green alga Spirulina platensis is rich in phycocyanins, that exhibit a wide range of pharmacological actions. C-phycocyanin (C-PC), in particular, possesses hepatoprotective, nephroprotective, antioxidant, and anticancer effects. Furthermore, several studies have reported both anti- and proinflammatory properties of this pigment. However, the precise mechanism(s) of action of C-PC in these processes remain largely unknown. Therefore, here we explored the C-PC effect in in vitro microglia activation. The effect of C-PC on the expression and release of IL-1ß and TNF-α and the activation of NF-κB was examined in primary microglia by real-time PCR, ELISA, and immunofluorescence. Treatment with C-PC up-regulated the expression and release of IL-1ß and TNF-α. C-PC also promoted the nuclear translocation of the NF-κB transcription factor. Then, to elucidate the molecular mechanisms for the immunoregulatory function of C-PC, we focused on investigating the role of Toll-like receptor 4 (TLR4). Accordingly, several TLR4 inhibitors have been used. Curcumin, ciprofloxacin, L48H37, and CLI-095 that suppresses specifically TLR4 signaling, blocked IL-1ß and TNF-α. Overall, these results indicate the immunomodulatory effect of C-PC in microglia cultures and show for the first time that the molecular mechanism implicated in this effect may involve TLR4 activation.