Supervised Molecular Dynamics (SuMD) Insights into the mechanism of action of SARS-CoV-2 main protease inhibitor PF-07321332

The chemical structure of PF-07321332, the first orally available Covid-19 clinical candidate, has recently been revealed by Pfizer. No information has been provided about the interaction pattern between PF-07321332 and its biomolecular counterpart, the SARS-CoV-2 main protease (M^pro). In the present work, we exploited Supervised Molecular Dynamics (SuMD) simulations to elucidate the key features that characterise the interaction between this drug candidate and the protease, emphasising similarities and differences with other structurally related inhibitors such as Boceprevir and PF-07304814. The structural insights provided by SuMD will hopefully be able to inspire the rational discovery of other potent and selective protease inhibitors.

Targeting the coronavirus SARS-CoV-2: computational insights into the mechanism of action of the protease inhibitors lopinavir, ritonavir and nelfinavir

Coronavirus SARS-CoV-2 is a recently discovered single-stranded RNA betacoronavirus, responsible for a severe respiratory disease known as coronavirus disease 2019, which is rapidly spreading. Chinese health authorities, as a response to the lack of an effective therapeutic strategy, started to investigate the use of lopinavir and ritonavir, previously optimized for the treatment and prevention of HIV/AIDS viral infection. Despite the clinical use of these two drugs, no information regarding their possible mechanism of action at the molecular level is still known for SARS-CoV-2. Very recently, the crystallographic structure of the SARS-CoV-2 main protease (Mpro), also known as C30 Endopeptidase, was published. Starting from this essential structural information, in the present work we have exploited supervised molecular dynamics, an emerging computational technique that allows investigating at an atomic level the recognition process of a ligand from its unbound to the final bound state. In this research, we provided molecular insight on the whole recognition pathway of Lopinavir, Ritonavir, and Nelfinavir, three potential C30 Endopeptidase inhibitors, with the last one taken into consideration due to the promising in-vitro activity shown against the structurally related SARS-CoV protease.

Re-Exploring the Ability of Common Docking Programs to Correctly Reproduce the Binding Modes of Non-Covalent Inhibitors of SARS-CoV-2 Protease Mpro

In the latest few decades, molecular docking has imposed itself as one of the most used approaches for computational drug discovery. Several docking benchmarks have been published, comparing the performance of different algorithms in respect to a molecular target of interest, usually evaluating their ability in reproducing the experimental data, which, in most cases, comes from X-ray structures. In this study, we elucidated the variation of the performance of three docking algorithms, namely GOLD, Glide, and PLANTS, in replicating the coordinates of the crystallographic ligands of SARS-CoV-2 main protease (Mpro). Through the comparison of the data coming from docking experiments and the values derived from the calculation of the solvent exposure of the crystallographic ligands, we highlighted the importance of this last variable for docking performance. Indeed, we underlined how an increase in the percentage of the ligand surface exposed to the solvent in a crystallographic complex makes it harder for the docking algorithms to reproduce its conformation. We further validated our hypothesis through molecular dynamics simulations, showing that the less stable protein–ligand complexes (in terms of root-mean-square deviation and root-mean-square fluctuation) tend to be derived from the cases in which the solvent exposure of the ligand in the starting system is higher.

On the Value of Using 3D Shape and Electrostatic Similarities in Deep Generative Methods

Multiparameter optimization, the heart of drug design, is still an open challenge. Thus, improved methods for automated compound design with multiple controlled properties are desired. Here, we present a significant extension to our previously described fragment-based reinforcement learning method (DeepFMPO) for the generation of novel molecules with optimal properties. As before, the generative process outputs optimized molecules similar to the input structures, now with the improved feature of replacing parts of these molecules with fragments of similar three-dimensional (3D) shape and electrostatics. We developed and benchmarked a new python package, ESP-Sim, for the comparison of the electrostatic potential and the molecular shape, allowing the calculation of high-quality partial charges (e.g., RESP with B3LYP/6-31G**) obtained using the quantum chemistry program Psi4. By performing comparisons of 3D fragments, we can simulate 3D properties while overcoming the notoriously difficult step of accurately describing bioactive conformations. The new improved generative (DeepFMPO v3D) method is demonstrated with a scaffold-hopping exercise identifying CDK2 bioisosteres. The code is open-source and freely available.

Deciphering the Molecular Recognition Mechanism of Multidrug Resistance Staphylococcus aureus NorA Efflux Pump Using a Supervised Molecular Dynamics Approach

The use and misuse of antibiotics has resulted in critical conditions for drug-resistant bacteria emergency, accelerating the development of antimicrobial resistance (AMR). In this context, the co-administration of an antibiotic with a compound able to restore sufficient antibacterial activity may be a successful strategy. In particular, the identification of efflux pump inhibitors (EPIs) holds promise for new antibiotic resistance breakers (ARBs). Indeed, bacterial efflux pumps have a key role in AMR development; for instance, NorA efflux pump contributes to Staphylococcus aureus (S. aureus) resistance against fluoroquinolone antibiotics (e.g., ciprofloxacin) by promoting their active extrusion from the cells. Even though NorA efflux pump is known to be a potential target for EPIs development, the absence of structural information about this protein and the little knowledge available on its mechanism of action have strongly hampered rational drug discovery efforts in this area. In the present work, we investigated at the molecular level the substrate recognition pathway of NorA through a Supervised Molecular Dynamics (SuMD) approach, using a NorA homology model. Specific amino acids were identified as playing a key role in the efflux pump-mediated extrusion of its substrate, paving the way for a deeper understanding of both the mechanisms of action and the inhibition of such efflux pumps.

Can We Still Trust Docking Results? An Extension of the Applicability of DockBench on PDBbind Database

The number of entries in the Protein Data Bank (PDB) has doubled in the last decade, and it has increased tenfold in the last twenty years. The availability of an ever-growing number of structures is having a huge impact on the Structure-Based Drug Discovery (SBDD), allowing investigation of new targets and giving the possibility to have multiple structures of the same macromolecule in a complex with different ligands. Such a large resource often implies the choice of the most suitable complex for molecular docking calculation, and this task is complicated by the plethora of possible posing and scoring function algorithms available, which may influence the quality of the outcomes. Here, we report a large benchmark performed on the PDBbind database containing more than four thousand entries and seventeen popular docking protocols. We found that, even in protein families wherein docking protocols generally showed acceptable results, certain ligand-protein complexes are poorly reproduced in the self-docking procedure. Such a trend in certain protein families is more pronounced, and this underlines the importance in identification of a suitable protein–ligand conformation coupled to a well-performing docking protocol.

Scaffold Repurposing of in-House Chemical Library toward the Identification of New Casein Kinase 1 δ Inhibitors

Recent studies have highlighted the key role of Casein kinase 1 δ (CK1δ) in the development of several neurodegenerative pathologies, such as Alzheimer's disease (AD), Parkinson's disease (PD), and amyotrophic lateral sclerosis (ALS). So far, CK1δ inhibitors are noncovalent ATP competitive ligands and no drugs are currently available for this molecular target, hence the interest in developing new CK1δ inhibitors. The study aims to identify new inhibitors able to bind the enzyme; by a dual approach in silico/in vitro, the virtual screening has been performed on an in-house chemical library, which was previously designed and synthesized for other targets. The work can, therefore, be seen in the scaffold repurposing logic. The proposed strategy has led to the identification of two hits, having a novel scaffold in the landscape of CK1δ inhibitors and with an activity in the micromolar range.

Revisiting the Allosteric Regulation of Sodium Cation on the Binding of Adenosine at the Human A2A Adenosine Receptor: Insights from Supervised Molecular Dynamics (SuMD) Simulations

One of the most intriguing findings highlighted from G protein-coupled receptor (GPCR) crystallography is the presence, in many members of class A, of a partially hydrated sodium ion in the middle of the seven transmembrane helices (7TM) bundle. In particular, the human adenosine A2A receptor (A2A AR) is the first GPCR in which a monovalent sodium ion was crystallized in a distal site from the canonical orthosteric one, corroborating, from a structural point of view, its role as a negative allosteric modulator. However, the molecular mechanism by which the sodium ion influences the recognition of the A2A AR agonists is not yet fully understood. In this study, the supervised molecular dynamics (SuMD) technique was exploited to analyse the sodium ion recognition mechanism and how its presence influences the binding of the endogenous agonist adenosine. Due to a higher degree of flexibility of the receptor extracellular (EC) vestibule, we propose the sodium-bound A2A AR as less efficient in stabilizing the adenosine during the different steps of binding.

A Computational Workflow for the Identification of Novel Fragments Acting as Inhibitors of the Activity of Protein Kinase CK1δ

Fragment-Based Drug Discovery (FBDD) has become, in recent years, a consolidated approach in the drug discovery process, leading to several drug candidates under investigation in clinical trials and some approved drugs. Among these successful applications of the FBDD approach, kinases represent a class of targets where this strategy has demonstrated its real potential with the approved kinase inhibitor Vemurafenib. In the Kinase family, protein kinase CK1 isoform δ (CK1δ) has become a promising target in the treatment of different neurodegenerative diseases such as Alzheimer’s disease, Parkinson’s disease, and amyotrophic lateral sclerosis. In the present work, we set up and applied a computational workflow for the identification of putative fragment binders in large virtual databases. To validate the method, the selected compounds were tested in vitro to assess the CK1δ inhibition.

Inspecting the Mechanism of Fragment Hits Binding on SARS-CoV-2 Mpro by Using Supervised Molecular Dynamics (SuMD) Simulations

Computational approaches supporting the early characterization of fragment molecular recognition mechanism represent a valuable complement to more expansive and low‐throughput experimental techniques. In this retrospective study, we have investigated the geometric accuracy with which high‐throughput supervised molecular dynamics simulations (HT‐SuMD) can anticipate the experimental bound state for a set of 23 fragments targeting the SARS‐CoV‐2 main protease. Despite the encouraging results herein reported, in line with those previously described for other MD‐based posing approaches, a high number of incorrect binding modes still complicate HT‐SuMD routine application. To overcome this limitation, fragment pose stability has been investigated and integrated as part of our in‐silico pipeline, allowing us to prioritize only the more reliable predictions.

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.

Computational strategies to identify new drug candidates against neuroinflammation

The even more increasing application of computational approaches in these last decades has deeply modified the process of discovery and commercialization of new therapeutic entities. This is especially true in the field of neuroinflammation, in which both the peculiar anatomical localization and the presence of the blood-brain barrier makeit mandatory to finely tune the candidates' physicochemical properties from the early stages of the discovery pipeline. The aim of this review is therefore to provide a general overview to the readers about the topic of neuroinflammation, together with the most common computational strategies that can be exploited to discover and design small molecules controlling neuroinflammation, especially those based on the knowledge of the three-dimensional structure of the biological targets of therapeutic interest. The techniques used to describe the molecular recognition mechanisms, such as molecular docking and molecular dynamics, will therefore be eviscerated, highlighting their advantages and their limitations. Finally, we report several case studies in which computational methods have been applied in drug discovery on neuroinflammation, focusing on the last decade's research.

Are two riboses better than one? The case of the recognition and activation of adenosine receptors

Traditionally, molecular recognition between the orthosteric site of adenosine receptors and their endogenous ligand occurs with a 1:1 stoichiometry. Inspired by previous mechanistic insights derived from Supervised Molecular Dynamics (SuMD) simulations, which suggested an alternative 2:1 binding stoichiometry, we synthesized BRA1, a bis-ribosyl adenosine derivative, tested its ability to bind to and activate the members of the adenosine receptors family, and rationalized its activity through molecular modeling.

Ribose and Non-Ribose A2A Adenosine Receptor Agonists: Do They Share the Same Receptor Recognition Mechanism?

Adenosine receptors have been a promising class of targets for the development of new therapies for several diseases. In recent years, a renewed interest in this field has risen, thanks to the implementation of a novel class of agonists that lack the ribose moiety, once considered essential for the agonistic profile. Recently, an X-ray crystal structure of the A_2A adenosine receptor has been solved, providing insights about the receptor activation from this novel class of agonists. Starting from this structural information, we have performed supervised molecular dynamics (SuMD) simulations to investigate the binding pathway of a non-nucleoside adenosine receptor agonist as well as one of three classic agonists. Furthermore, we analyzed the possible role of water molecules in receptor activation.

Shedding Light on the Molecular Recognition of Sub-Kilodalton Macrocyclic Peptides on Thrombin by Supervised Molecular Dynamics

Macrocycles are attractive structures for drug development due to their favorable structural features, potential in binding to targets with flat featureless surfaces, and their ability to disrupt protein–protein interactions. Moreover, large novel highly diverse libraries of low-molecular-weight macrocycles with therapeutically favorable characteristics have been recently established. Considering the mentioned facts, having a validated, fast, and accurate computational protocol for studying the molecular recognition and binding mode of this interesting new class of macrocyclic peptides deemed to be helpful as well as insightful in the quest of accelerating drug discovery. To that end, the ability of the in-house supervised molecular dynamics protocol called SuMD in the reproduction of the X-ray crystallography final binding state of a macrocyclic non-canonical tetrapeptide—from a novel library of 8,988 sub-kilodalton macrocyclic peptides—in the thrombin active site was successfully validated. A comparable binding mode with the minimum root-mean-square deviation (RMSD) of 1.4 Å at simulation time point 71.6 ns was achieved. This method validation study extended the application domain of the SuMD sampling method for computationally cheap, fast but accurate, and insightful macrocycle–protein molecular recognition studies.

[1,2,4]Triazolo[1,5-c]pyrimidines as Tools to Investigate A3 Adenosine Receptors in Cancer Cell Lines

The A3 adenosine receptor is an interesting target whose role in cancer is controversial. In this work, a structural investigation at the 2-position of the [1,2,4]triazolo[1,5-c]pyrimidine nucleus was performed, finding new potent and selective A3 adenosine receptor antagonists such as the ethyl 2-(4-methoxyphenyl)-5-(methylamino)-[1,2,4]triazolo[1,5-c]pyrimidine-8-carboxylate (20, DZ123) that showed a Ki value of 0.47 nM and an exceptional selectivity profile over the other adenosine receptor subtypes. Computational studies were performed to rationalize the affinity and the selectivity profile of the tested compounds at the A3 adenosine receptor and the A1 and A2A adenosine receptors. Compound 20 was tested on both A3 adenosine receptor positive cell lines (CHO-A3 AR transfected, THP1 and HCT16) and on A3 negative cancer cell lines, showing no effect in the latter and a pro-proliferative effect at a low concentration in the former. These interesting results pave the way to further investigation on both the mechanism involved and potential therapeutic applications.

Towards Multitargeted Ligands as Pain Therapeutics: Dual Ligands of the Cavα2δ-1 Subunit of Voltage-Gated Calcium Channel and the μ-Opioid Receptor

The synthesis and pharmacological activity of a new series of dual ligands combining activities towards the α2δ-1 subunit of voltage-gated calcium channels (Cavα2δ-1) and the μ-opioid receptor (MOR) as novel pain therapeutics are reported. A careful exploration of the pharmacophores related to both targets, which in principle had few common characteristics, led to the design of novel compounds exhibiting both activities. The construction of the dual ligands started from published Cavα2δ-1 ligands, onto which MOR ligand pharmacophoric elements were added. This exercise led to new amino-acidic substances with good affinities on both targets as well as good metabolic and physicochemical profiles and low potential for drug-drug interactions. A representative compound, (2S,4S)-4-(4-chloro-3-(((cis)-4-(dimethylamino)-4-phenylcyclohexyl)methyl)-5-fluorophenoxy)pyrrolidine-2-carboxylic acid, displayed promising analgesic activities in several in vivo pain models as well as a reduced side-effect profile in relation to morphine.

Species differences in opsonization and phagocyte recognition of preclinical poly-2-alkyl-2-oxazoline-coated nanoparticles

Poly(ethylene glycol) (PEG) is widely used in nanomedicine design, but emerging PEG immunogenicity in the general population is of therapeutic concern. As alternative, polyoxazolines are gaining popularity, since "polyoxazolinated" nanoparticles show long-circulating properties comparable to PEGylated nanoparticles in mice. Here, we show species differences in opsonization and differential uptake by monocytes and macrophages of nanoparticles coated with either poly-2-methyl-2-oxazoline or poly-2-ethyl-2-oxazoline. These nanoparticles evade murine opsonization process and phagocytic uptake but porcine ficolin 2 (FCN2), through its S2 binding site, recognizes polyoxazolines, and mediates nanoparticle uptake exclusively by porcine monocytes. In human sera, FCN opsonization is isoform-dependent showing inter-individual variability but both FCN2 and complement opsonization promote nanoparticle uptake by human monocytes. However, nanoparticle uptake by human and porcine macrophages is complement-dependent. These findings advance mechanistic understanding of species differences in innate immune recognition of nanomaterials' molecular patterns, and applicable to the selection and chemical design of polymers for engineering of the next generation of stealth nanoparticles.