Structure-based virtual screening discovers potent and selective adenosine A1 receptor antagonists.

Development of subtype-selective leads is essential in drug discovery campaigns targeting G protein-coupled receptors (GPCRs). Herein, a structure-based virtual screening approach to rationally design subtype-selective ligands was applied to the A1 and A2A adenosine receptors (A1R and A2AR). Crystal structures of these closely related subtypes revealed a non-conserved subpocket in the binding sites that could be exploited to identify A1R selective ligands. A library of 4.6 million compounds was screened computationally against both receptors using molecular docking and 20 A1R selective ligands were predicted. Of these, seven antagonized the A1R with micromolar activities and several compounds displayed slight selectivity for this subtype. Twenty-seven analogs of two discovered scaffolds were designed, resulting in antagonists with nanomolar potency and up to 76-fold A1R-selectivity. Our results show the potential of structure-based virtual screening to guide discovery and optimization of subtype-selective ligands, which could facilitate the development of safer drugs.

Macrocycles in Drug Discovery─Learning from the Past for the Future.

We have analyzed FDA-approved macrocyclic drugs, clinical candidates, and the recent literature to understand how macrocycles are used in drug discovery. Current drugs are mainly used in infectious disease and oncology, while oncology is the major indication for the clinical candidates and in the literature Most macrocyclic drugs bind to targets that have difficult to drug binding sites. Natural products have provided 80-90% of the drugs and clinical candidates, whereas macrocycles in ChEMBL have less complex structures. Macrocycles usually reside in the beyond the Rule of 5 chemical space, but 30-40% of the drugs and clinical candidates are orally bioavailable. Simple bi-descriptor models, i.e., HBD ≤ 7 in combination with either MW < 1000 Da or cLogP > 2.5, distinguished orals from parenterals and can be used as filters in design. We propose that recent breakthroughs in conformational analysis and inspiration from natural products will further improve the de novo design of macrocycles.

Conformational Sampling Deciphers the Chameleonic Properties of a VHL-Based Degrader.

Chameleonicity (the capacity of a molecule to adapt its conformations to the environment) may help to identify orally bioavailable drugs in the beyond-Rule-of-5 chemical space. Computational methods to predict the chameleonic behaviour of degraders have not yet been reported and the identification of molecular chameleons still relies on experimental evidence. Therefore, there is a need to tune predictions with experimental data. Here, we employ PROTAC-1 (a passively cell-permeable degrader), for which NMR and physicochemical data prove the chameleonic behaviour, to benchmark the capacity of two conformational sampling algorithms and selection schemes. To characterize the conformational ensembles in both polar and nonpolar environments, we compute three molecular properties proven to be essential for cell permeability: conformer shape (radius of gyration), polarity (3D PSA), and the number of intramolecular hydrogen bonds. Energetic criteria were also considered. Infographics monitored the simultaneous variation of those properties in computed and NMR conformers. Overall, we provide key points for tuning conformational sampling tools to reproduce PROTAC-1 chameleonicity according to NMR evidence. This study is expected to improve the design of PROTAC drugs and the development of computational sustainable strategies to exploit the potential of new modalities in drug discovery.

Autoantibodies to Disease-Related Proteins in Joints as Novel Biomarkers for the Diagnosis of Rheumatoid Arthritis.

OBJECTIVE: This study was undertaken to develop and characterize a multiplex immunoassay for detection of autoantibodies against peptides derived from proteins known to play a role in development of arthritis and that are also expressed in joints. METHODS: We selected peptides from the human counterpart of proteins expressed in the joints, based on mouse models that showed these to be targeted by pathogenic or regulatory antibodies in vivo. Using bead-based flow immunoassays measuring IgG antibodies, we selected triple helical or cyclic peptides, containing the epitopes, to avoid collinear reactivity. We characterized the analytical performance of the immunoassay and then validated it in 3 independent rheumatoid arthritis (RA) cohorts (n = 2,110), Swedish age- and sex-matched healthy controls, and patients with osteoarthritis (OA), patients with psoriatic arthritis (PsA), and patients with systemic lupus erythematosus (SLE). RESULTS: Screening assays showed 5 peptide antigens that discriminated RA patients from healthy controls with 99% specificity (95% confidence interval [CI] 98-100%). In our validation studies, we reproduced the discriminatory capacity of the autoantibodies in 2 other RA cohorts, showing that the autoantibodies had high discriminatory capacity for RA versus OA, PsA, and SLE. The novel biomarkers identified 22.5% (95% CI 19-26%) of early RA patients seronegative for anti-cyclic citrullinated peptide and rheumatoid factor. The usefulness of the biomarkers in identifying seronegative RA patients was confirmed in validation studies using 2 independent cohorts of RA patients and cohorts of patients with OA, PsA, and SLE. CONCLUSION: A multiplex immunoassay with peptides from disease-related proteins in joints was found to be useful for detection of specific autoantibodies in RA serum. Of note, this immunoassay had high discriminatory capacity for early seronegative RA.

Ultralarge Virtual Screening Identifies SARS-CoV-2 Main Protease Inhibitors with Broad-Spectrum Activity against Coronaviruses.

Drugs targeting SARS-CoV-2 could have saved millions of lives during the COVID-19 pandemic, and it is now crucial to develop inhibitors of coronavirus replication in preparation for future outbreaks. We explored two virtual screening strategies to find inhibitors of the SARS-CoV-2 main protease in ultralarge chemical libraries. First, structure-based docking was used to screen a diverse library of 235 million virtual compounds against the active site. One hundred top-ranked compounds were tested in binding and enzymatic assays. Second, a fragment discovered by crystallographic screening was optimized guided by docking of millions of elaborated molecules and experimental testing of 93 compounds. Three inhibitors were identified in the first library screen, and five of the selected fragment elaborations showed inhibitory effects. Crystal structures of target-inhibitor complexes confirmed docking predictions and guided hit-to-lead optimization, resulting in a noncovalent main protease inhibitor with nanomolar affinity, a promising in vitro pharmacokinetic profile, and broad-spectrum antiviral effect in infected cells.

Key interactions in the trimolecular complex consisting of the rheumatoid arthritis-associated DRB1*04:01 molecule, the major glycosylated collagen II peptide and the T-cell receptor.

OBJECTIVES: Rheumatoid arthritis (RA) is an autoimmune disease strongly associated with the major histocompatibility complex (MHC) class II allele DRB1*04:01, which encodes a protein that binds self-peptides for presentation to T cells. This study characterises the autoantigen-presenting function of DRB1*04:01 (HLA-DRA*01:01/HLA-DRB1*04:01) at a molecular level for prototypic T-cell determinants, focusing on a post-translationally modified collagen type II (Col2)-derived peptide. METHODS: The crystal structures of DRB1*04:01 molecules in complex with the peptides HSP70289-306, citrullinated CILP982-996 and galactosylated Col2259-273 were determined on cocrystallisation. T cells specific for Col2259-273 were investigated in peripheral blood mononuclear cells from patients with DRB1*04:01-positive RA by cytofluorometric detection of the activation marker CD154 on peptide stimulation and binding of fluorescent DRB1*0401/Col2259-273 tetramer complexes. The cDNAs encoding the T-cell receptor (TCR) α-chains and ß-chains were cloned from single-cell sorted tetramer-positive T cells and transferred via a lentiviral vector into TCR-deficient Jurkat 76 cells. RESULTS: The crystal structures identified peptide binding to DRB1*04:01 and potential side chain exposure to T cells. The main TCR recognition sites in Col2259-273 were lysine residues that can be galactosylated. RA T-cell responses to DRB1*04:01-presented Col2259-273 were dependent on peptide galactosylation at lysine 264. Dynamic molecular modelling of a functionally characterised Col2259-273-specific TCR complexed with DRB1*04:01/Col2259-273 provided evidence for differential allosteric T-cell recognition of glycosylated lysine 264. CONCLUSIONS: The MHC-peptide-TCR interactions elucidated in our study provide new molecular insights into recognition of a post-translationally modified RA T-cell determinant with a known dominant role in arthritogenic and tolerogenic responses in murine Col2-induced arthritis.

Structure-Guided Design of G-Protein-Coupled Receptor Polypharmacology.

Many diseases are polygenic and can only be treated efficiently with drugs that modulate multiple targets. However, rational design of compounds with multi-target profiles is rarely pursued because it is considered too difficult, in particular if the drug must enter the central nervous system. Here, a structure-based strategy to identify dual-target ligands of G-protein-coupled receptors is presented. We use this approach to design compounds that both antagonize the A2A adenosine receptor and activate the D2 dopamine receptor, which have excellent potential as antiparkinson drugs. Atomic resolution models of the receptors guided generation of a chemical library with compounds designed to occupy orthosteric and secondary binding pockets in both targets. Structure-based virtual screens identified ten compounds, of which three had affinity for both targets. One of these scaffolds was optimized to nanomolar dual-target activity and showed the predicted pharmacodynamic effect in a rat model of Parkinsonism.

Macrocyclic Peptides Uncover a Novel Binding Mode for Reversible Inhibitors of LSD1.

Lysine-specific demethylase 1 (LSD1) is an epigenetic enzyme which regulates the methylation of Lys4 of histone 3 (H3) and is overexpressed in certain cancers. We used structures of H3 substrate analogues bound to LSD1 to design macrocyclic peptide inhibitors of LSD1. A linear, Lys4 to Met-substituted, 11-mer (4) was identified as the shortest peptide distinctly interacting with LSD1. It was evolved into macrocycle 31, which was >40 fold more potent (K i = 2.3 µM) than 4. Linear and macrocyclic peptides exhibited unexpected differences in structure-activity relationships for interactions with LSD1, indicating that they bind LSD1 differently. This was confirmed by the crystal structure of 31 in complex with LSD1-CoREST1, which revealed a novel binding mode at the outer rim of the LSD1 active site and without a direct interaction with FAD. NMR spectroscopy of 31 suggests that macrocyclization restricts its solution ensemble to conformations that include the one in the crystalline complex. Our results provide a solid basis for the design of optimized reversible LSD1 inhibitors.

Docking Finds GPCR Ligands in Dark Chemical Matter.

High-throughput screening has revealed dark chemical matter, a set of drug-like compounds that has never shown bioactivity despite being extensively assayed. If dark molecules are found active at a therapeutic target, their extraordinary selectivity profiles make excellent starting points for drug development. We explored if ligands of therapeutically relevant G-protein-coupled receptors could be discovered by structure-based virtual screening of the dark chemical matter. Molecular docking screens against crystal structures of the A2A adenosine and the D4 dopamine receptors were carried out, and 53 top-ranked molecules were evaluated experimentally. Two ligands of each receptor were discovered, and the most potent had sub-micromolar affinities. Analysis of bioactivity data showed that the ligands lacked activity at hundreds of off-targets, including several that are associated with adverse effects. Our results demonstrate that virtual screening provides an efficient means to mine the dark chemical space, which could contribute to development of drugs with improved safety profiles.

Dynamics Determine Signaling in a Multicomponent System Associated with Rheumatoid Arthritis.

Strategies that target multiple components are usually required for treatment of diseases originating from complex biological systems. The multicomponent system consisting of the DR4 major histocompatibility complex type II molecule, the glycopeptide CII259-273 from type II collagen, and a T-cell receptor is associated with development of rheumatoid arthritis (RA). We introduced non-native amino acids and amide bond isosteres into CII259-273 and investigated the effect on binding to DR4 and the subsequent T-cell response. Molecular dynamics simulations revealed that complexes between DR4 and derivatives of CII259-273 were highly dynamic. Signaling in the overall multicomponent system was found to depend on formation of an appropriate number of dynamic intramolecular hydrogen bonds between DR4 and CII259-273, together with the positioning of the galactose moiety of CII259-273 in the DR4 binding groove. Interestingly, the system tolerated modifications at several positions in CII259-273, indicating opportunities to use analogues to increase our understanding of how rheumatoid arthritis develops and for evaluation as vaccines to treat RA.

Impact of Dynamically Exposed Polarity on Permeability and Solubility of Chameleonic Drugs Beyond the Rule of 5.

Conformational flexibility has been proposed to significantly affect drug properties outside rule-of-5 (Ro5) chemical space. Here, we investigated the influence of dynamically exposed polarity on cell permeability and aqueous solubility for a structurally diverse set of drugs and clinical candidates far beyond the Ro5, all of which populated multiple distinct conformations as revealed by X-ray crystallography. Efflux-inhibited (passive) Caco-2 cell permeability correlated strongly with the compounds' minimum solvent-accessible 3D polar surface areas (PSA), whereas aqueous solubility depended less on the specific 3D conformation. Inspection of the crystal structures highlighted flexibly linked aromatic side chains and dynamically forming intramolecular hydrogen bonds as particularly effective in providing "chameleonic" properties that allow compounds to display both high cell permeability and aqueous solubility. These structural features, in combination with permeability predictions based on the correlation to solvent-accessible 3D PSA, should inspire drug design in the challenging chemical space far beyond the Ro5.

Cutting Edge: Processing of Oxidized Peptides in Macrophages Regulates T Cell Activation and Development of Autoimmune Arthritis.

APCs are known to produce NADPH oxidase (NOX) 2-derived reactive oxygen species; however, whether and how NOX2-mediated oxidation affects redox-sensitive immunogenic peptides remains elusive. In this study, we investigated a major immunogenic peptide in glucose-6-phosphate isomerase (G6PI), a potential autoantigen in rheumatoid arthritis, which can form internal disulfide bonds. Ag presentation assays showed that presentation of this G6PI peptide was more efficient in NOX2-deficient (Ncf1m1J/m1J mutant) mice, compared with wild-type controls. IFN-γ-inducible lysosomal thiol reductase (GILT), which facilitates disulfide bond-containing Ag processing, was found to be upregulated in macrophages from Ncf1 mutant mice. Ncf1 mutant mice exhibited more severe G6PI peptide-induced arthritis, which was accompanied by the increased GILT expression in macrophages and enhanced Ag-specific T cell responses. Our results show that NOX2-dependent processing of the redox-sensitive autoantigens by APCs modify T cell activity and development of autoimmune arthritis.

Macrocyclic drugs and clinical candidates: what can medicinal chemists learn from their properties?

Macrocycles are ideal in efforts to tackle "difficult" targets, but our understanding of what makes them cell permeable and orally bioavailable is limited. Analysis of approximately 100 macrocyclic drugs and clinical candidates revealed that macrocycles are predominantly used for infectious disease and in oncology and that most belong to the macrolide or cyclic peptide class. A significant number (N = 34) of these macrocycles are administered orally, revealing that oral bioavailability can be obtained at molecular weights up to and above 1 kDa and polar surface areas ranging toward 250 Å(2). Moreover, insight from a group of "de novo designed" oral macrocycles in clinical studies and understanding of how cyclosporin A and model cyclic hexapeptides cross cell membranes may unlock wider opportunities in drug discovery. However, the number of oral macrocycles is still low and it remains to be seen if they are outliers or if macrocycles will open up novel oral druggable space.

3-aminopiperidine-based peptide analogues as the first selective noncovalent inhibitors of the bacterial cysteine protease IdeS.

A series of eight peptides corresponding to the amino acid sequence of the hinge region of IgG and 17 newly synthesized peptide analogues containing a piperidine moiety as a replacement of a glycine residue were tested as potential inhibitors of the bacterial IgG degrading enzyme of Streptococcus pyogenes , IdeS. None of the peptides showed any inhibitory activity of IdeS, but several piperidine-based analogues were identified as inhibitors. Two different analysis methods were used: an SDS-PAGE based assay to detect IgG cleavage products and a surface plasmon resonance spectroscopy based assay to quantify the degree of inhibition. To investigate the selectivity of the inhibitors for IdeS, all compounds were screened against two other related cysteine proteases (SpeB and papain). The selectivity results show that larger analogues that are active inhibitors of IdeS are even more potent as inhibitors of papain, whereas smaller analogues that are active inhibitors of IdeS inhibit neither SpeB nor papain. Two compounds were identified that exhibit high selectivity against IdeS and will be used for further studies.

(E)-alkene and ethylene isosteres substantially alter the hydrogen-bonding network in class II MHC A(q)/glycopeptide complexes and affect T-cell recognition.

The structural basis for antigen presentation by class II major histocompatibility complex (MHC) proteins to CD4(+) T-cells is important for understanding and possibly treating autoimmune diseases. In the work described in this paper, (E)-alkene and ethylene amide-bond isosteres were used to investigate the effect of removing hydrogen-bonding possibilities from the CII259-270 glycopeptide, which is bound by the arthritis-associated murine A(q) class II MHC protein. The isostere-modified glycopeptides showed varying and unexpectedly large losses of A(q) binding that could be linked to the dynamics of the system. Molecular dynamics (MD) simulations revealed that the backbone of CII259-270 and the A(q) protein are able to form up to 11 hydrogen bonds, but fewer than this number are present at any one time. Most of the strong hydrogen-bond interactions were formed by the N-terminal part of the glycopeptide, i.e., in the region where the isosteric replacements were made. The structural dynamics also revealed that hydrogen bonds were strongly coupled to each other; the loss of one hydrogen-bond interaction had a profound effect on the entire hydrogen-bonding network. The A(q) binding data revealed that an ethylene isostere glycopeptide unexpectedly bound more strongly to A(q) than the corresponding (E)-alkene, which is in contrast to the trend observed for the other isosteres. Analysis of the MD trajectories revealed that the complex conformation of this ethylene isostere was structurally different and had an altered molecular interaction pattern compared to the other A(q)/glycopeptide complexes. The introduced amide-bond isosteres also affected the interactions of the glycopeptide/A(q) complexes with T-cell receptors. The dynamic variation of the patterns and strengths of the hydrogen-bond interactions in the class II MHC system is of critical importance for the class II MHC/peptide/TCR signaling system.

Oxazole-modified glycopeptides that target arthritis-associated class II MHC A(q) and DR4 proteins.

The glycopeptide CII259-273, a fragment from type II collagen (CII), can induce tolerance in mice susceptible to collagen-induced arthritis (CIA), which is a validated disease model for rheumatoid arthritis (RA). Here, we describe the design and synthesis of a small series of modified CII259-273 glycopeptides with oxazole heterocycles replacing three potentially labile peptide bonds. These glycopeptidomimetics were evaluated for binding to murine CIA-associated A(q) and human RA-associated DR4 class II major histocompatibility complex (MHC) proteins. The oxazole modifications drastically reduced or completely abolished binding to A(q). Two of the glycopeptidomimetics were, however, well tolerated in binding to DR4 and they also induced strong responses by one or two DR4-restricted T-cell hybridomas. This work contributes to the development of an altered glycopeptide for inducing immunological tolerance in CIA, with the long-term goal of developing a therapeutic vaccine for treatment of RA.

Synthesis and biological evaluation of reversible inhibitors of IdeS, a bacterial cysteine protease and virulence determinant.

Analogues of the irreversible protease inhibitors TPCK and TLCK have been synthesized and tested as inhibitors of the bacterial cysteine protease IdeS excreted by Streptococcuspyogenes. Eight compounds were identified as inhibitors of IdeS in an in vitro assay. The most potent compounds contained an aldehyde function, thus acting as efficient reversible inhibitors, nitrile and azide derivatives showed moderate activity.

Application of gel-phase 19F NMR spectroscopy for optimization of solid-phase synthesis of a hydrophobic peptide from the signal sequence of the mucin MUC1.

This paper describes the manual Fmoc/t-Bu solid-phase synthesis of a difficult nine-residue hydrophobic peptide LLLLTVLTV from one of the signal sequences that flank the tandem repeat of the mucin MUC1. Gel-phase 19F NMR spectroscopy was used as a straightforward method for optimization of the solid-phase synthesis. Different approaches were applied for comparative studies. The strategy based on modified solid-phase conditions using DIC/HOAt for coupling, DBU for Fmoc deprotection, and the incorporation of the pseudo proline dipeptide Fmoc-Leu-Thr(psiMe, Me pro)-OH as a backbone-protecting group was found to be superior according to gel-phase 19F NMR spectroscopy. Implementation of the optimized Fmoc protocol enabled an effective synthesis of signal peptide LLLLTVLTV.

Multivalent sialic acid conjugates inhibit adenovirus type 37 from binding to and infecting human corneal epithelial cells.

Adenovirus type 37 is one of the main causative agents of epidemic keratoconjunctivitis. In a series of publications, we have reported that this virus uses sialic acid as a cellular receptor. Here we demonstrate in vitro that on a molar basis, multivalent sialic acid conjugated to human serum albumin prevents adenovirus type 37 from binding to and infecting human corneal epithelial cells 1000-fold more efficiently than monosaccharidic sialic acid. We also demonstrate that the extraordinary inhibitory effect of multivalent sialic acid is due to the ability of this compound to aggregate virions. We conclude that multivalent sialic acid may be a potential new antiviral drug, for use in the treatment of epidemic keratoconjunctivitis caused by the adenoviruses that use sialic acid as cellular receptor.

Formation of lactones from sialylated MUC1 glycopeptides.

The tumor-associated carbohydrate antigens TN, T, sialyl TN and sialyl T are expressed on mucins in several epithelial cancers. This has stimulated studies directed towards development of glycopeptide-based anticancer vaccines. Formation of intramolecular lactones involving sialic acid residues and suitably positioned hydroxyl groups in neighboring saccharide moieties is known to occur for glycolipids such as gangliosides. It has been suggested that these lactones are more immunogenic and tumor-specific than their native counterparts and that they might find use as cancer vaccines. We have now investigated if lactonization also occurs for the sialyl TN and T antigens of mucins. It was found that the model compound sialyl T benzyl glycoside , and the glycopeptide Ala-Pro-Asp-Thr-Arg-Pro-Ala from the tandem repeat of the mucin MUC1, in which Thr stands for the 2,3-sialyl-T antigen, lactonized during treatment with glacial acetic acid. Compound gave the 1''--> 2' lactone as the major product and the corresponding 1''--> 4' lactone as the minor product. For glycopeptide the 1''--> 4' lactone constitued the major product, whereas the 1''--> 2' lactone was the minor one. When lactonized was dissolved in water the 1''--> 4' lactone underwent slow hydrolysis, whereas the 1''--> 2' remained stable even after a 30 days incubation. In contrast the corresponding 2,6-sialyl-TN glycopeptide did not lactonize in glacial acetic acid.