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.

Discovery of an Allosteric Ligand Binding Site in SMYD3 Lysine Methyltransferase.

SMYD3 is a multifunctional epigenetic enzyme with lysine methyltransferase activity and various interaction partners. It is implicated in the pathophysiology of cancers but with an unclear mechanism. To discover tool compounds for clarifying its biochemistry and potential as a therapeutic target, a set of drug-like compounds was screened in a biosensor-based competition assay. Diperodon was identified as an allosteric ligand; its R and S enantiomers were isolated, and their affinities to SMYD3 were determined (KD =42 and 84 µM, respectively). Co-crystallization revealed that both enantiomers bind to a previously unidentified allosteric site in the C-terminal protein binding domain, consistent with its weak inhibitory effect. No competition between diperodon and HSP90 (a known SMYD3 interaction partner) was observed although SMYD3-HSP90 binding was confirmed (KD =13 µM). Diperodon clearly represents a novel starting point for the design of tool compounds interacting with a druggable allosteric site, suitable for the exploration of noncatalytic SMYD3 functions and therapeutics with new mechanisms of action.

BioMAX - the first macromolecular crystallography beamline at MAX IV Laboratory.

BioMAX is the first macromolecular crystallography beamline at the MAX IV Laboratory 3 GeV storage ring, which is the first operational multi-bend achromat storage ring. Due to the low-emittance storage ring, BioMAX has a parallel, high-intensity X-ray beam, even when focused down to 20 µm × 5 µm using the bendable focusing mirrors. The beam is tunable in the energy range 5-25 keV using the in-vacuum undulator and the horizontally deflecting double-crystal monochromator. BioMAX is equipped with an MD3 diffractometer, an ISARA high-capacity sample changer and an EIGER 16M hybrid pixel detector. Data collection at BioMAX is controlled using the newly developed MXCuBE3 graphical user interface, and sample tracking is handled by ISPyB. The computing infrastructure includes data storage and processing both at MAX IV and the Lund University supercomputing center LUNARC. With state-of-the-art instrumentation, a high degree of automation, a user-friendly control system interface and remote operation, BioMAX provides an excellent facility for most macromolecular crystallography experiments. Serial crystallography using either a high-viscosity extruder injector or the MD3 as a fixed-target scanner is already implemented. The serial crystallography activities at MAX IV Laboratory will be further developed at the microfocus beamline MicroMAX, when it comes into operation in 2022. MicroMAX will have a 1 µm × 1 µm beam focus and a flux up to 1015 photons s-1 with main applications in serial crystallography, room-temperature structure determinations and time-resolved experiments.

Crystal Structure of Non-Structural Protein 10 from Severe Acute Respiratory Syndrome Coronavirus-2.

Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2), causing Coronavirus Disease 19 (COVID-19), emerged at the end of 2019 and quickly spread to cause a global pandemic with severe socio-economic consequences. The early sequencing of its RNA genome revealed its high similarity to SARS, likely to have originated from bats. The SARS-CoV-2 non-structural protein 10 (nsp10) displays high sequence similarity with its SARS homologue, which binds to and stimulates the 3'-to-5' exoribonuclease and the 2'-O-methlytransferase activities of nsps 14 and 16, respectively. Here, we report the biophysical characterization and 1.6 Å resolution structure of the unbound form of nsp10 from SARS-CoV-2 and compare it to the structures of its SARS homologue and the complex-bound form with nsp16 from SARS-CoV-2. The crystal structure and solution behaviour of nsp10 will not only form the basis for understanding the role of SARS-CoV-2 nsp10 as a central player of the viral RNA capping apparatus, but will also serve as a basis for the development of inhibitors of nsp10, interfering with crucial functions of the replication-transcription complex and virus replication.

Unveiling the Biochemistry of the Epigenetic Regulator SMYD3.

SET and MYND domain-containing protein 3 (SMYD3) is a lysine methyltransferase that plays a central role in a variety of cancer diseases, exerting its pro-oncogenic activity by methylation of key proteins, of both nuclear and cytoplasmic nature. However, the role of SMYD3 in the initiation and progression of cancer is not yet fully understood and further biochemical characterization is required to support the discovery of therapeutics targeting this enzyme. We have therefore developed robust protocols for production, handling, and crystallization of SMYD3 and biophysical and biochemical assays for clarification of SMYD3 biochemistry and identification of useful lead compounds. Specifically, a time-resolved biosensor assay was developed for kinetic characterization of SMYD3 interactions. Functional differences in SMYD3 interactions with its natural small molecule ligands SAM and SAH were revealed, with SAM forming a very stable complex. A variety of peptides mimicking putative substrates of SMYD3 were explored in order to expose structural features important for recognition. The interaction between SMYD3 and some peptides was influenced by SAM. A nonradioactive SMYD3 activity assay using liquid chromatography-mass spectrometry (LC-MS) analysis explored substrate features of importance also for methylation. Methylation was notable only toward MAP kinase kinase kinase 2 (MAP3K2_K260)-mimicking peptides, although binary and tertiary complexes were detected also with other peptides. The analysis supported a random bi-bi mechanistic model for SMYD3 methyltransferase catalysis. Our work unveiled complexities in SMYD3 biochemistry and resulted in procedures suitable for further studies and identification of novel starting points for design of effective and specific leads for this potential oncology target.

Biophysical analysis of the dynamics of calmodulin interactions with neurogranin and Ca2+ /calmodulin-dependent kinase II.

Calmodulin (CaM) functions depend on interactions with CaM-binding proteins, regulated by Ca2+. Induced structural changes influence the affinity, kinetics, and specificities of the interactions. The dynamics of CaM interactions with neurogranin (Ng) and the CaM-binding region of Ca2+/calmodulin-dependent kinase II (CaMKII290-309 ) have been studied using biophysical methods. These proteins have opposite Ca2+ dependencies for CaM binding. Surface plasmon resonance biosensor analysis confirmed that Ca2+ and CaM interact very rapidly, and with moderate affinity ( KDSPR=3µM). Calmodulin-CaMKII290-309 interactions were only detected in the presence of Ca2+, exhibiting fast kinetics and nanomolar affinity ( KDSPR=7.1nM). The CaM-Ng interaction had higher affinity under Ca2+-depleted ( KDSPR=480nM,k1=3.4×105M-1s-1 and k-1 = 1.6 × 10-1 s-1 ) than Ca2+-saturated conditions ( KDSPR=19µM). The IQ motif of Ng (Ng27-50 ) had similar affinity for CaM as Ng under Ca2+-saturated conditions ( KDSPR=14µM), but no interaction was seen under Ca2+-depleted conditions. Microscale thermophoresis using fluorescently labeled CaM confirmed the surface plasmon resonance results qualitatively, but estimated lower affinities for the Ng ( KDMST=890nM) and CaMKII290-309 ( KDMST=190nM) interactions. Although CaMKII290-309 showed expected interaction characteristics, they may be different for full-length CaMKII. The data for full-length Ng, but not Ng27-50 , agree with the current model on Ng regulation of Ca2+/CaM signaling.

Multiplexed experimental strategies for fragment library screening against challenging drug targets using SPR biosensors.

Surface plasmon resonance (SPR) biosensor methods are ideally suited for fragment-based lead discovery.  However, generally applicable experimental procedures and detailed protocols are lacking, especially for structurally or physico-chemically challenging targets or when tool compounds are not available. Success depends on accounting for the features of both the target and the chemical library, purposely designing screening experiments for identification and validation of hits with desired specificity and mode-of-action, and availability of orthogonal methods capable of confirming fragment hits. The range of targets and libraries amenable to an SPR biosensor-based approach for identifying hits is considerably expanded by adopting multiplexed strategies, using multiple complementary surfaces or experimental conditions. Here we illustrate principles and multiplexed approaches for using flow-based SPR biosensor systems for screening fragment libraries of different sizes (90 and 1056 compounds) against a selection of challenging targets. It shows strategies for the identification of fragments interacting with 1) large and structurally dynamic targets, represented by acetyl choline binding protein (AChBP), a Cys-loop receptor ligand gated ion channel homologue, 2) targets in multi protein complexes, represented by lysine demethylase 1 and a corepressor (LSD1/CoREST), 3) structurally variable or unstable targets, represented by farnesyl pyrophosphate synthase (FPPS), 4) targets containing intrinsically disordered regions, represented by protein tyrosine phosphatase 1B  (PTP1B), and 5) aggregation-prone proteins, represented by an engineered form of human tau  (tau K18M). Practical considerations and procedures accounting for the characteristics of the proteins and libraries, and that increase robustness, sensitivity, throughput and versatility are highlighted. The study shows that the challenges for addressing these types of targets is not identification of potentially useful fragments per se, but establishing methods for their validation and evolution into leads.

Identification of fragments binding to SARS-CoV-2 nsp10 reveals ligand-binding sites in conserved interfaces between nsp10 and nsp14/nsp16.

Since the emergence of SARS-CoV-2 in 2019, Covid-19 has developed into a serious threat to our health, social and economic systems. Although vaccines have been developed in a tour-de-force and are now increasingly available, repurposing of existing drugs has been less successful. There is a clear need to develop new drugs against SARS-CoV-2 that can also be used against future coronavirus infections. Non-structural protein 10 (nsp10) is a conserved stimulator of two enzymes crucial for viral replication, nsp14 and nsp16, exhibiting exoribonuclease and methyltransferase activities. Interfering with RNA proofreading or RNA cap formation represents intervention strategies to inhibit replication. We applied fragment-based screening using nano differential scanning fluorometry and X-ray crystallography to identify ligands targeting SARS-CoV-2 nsp10. We identified four fragments located in two distinct sites: one can be modelled to where it would be located in the nsp14-nsp10 complex interface and the other in the nsp16-nsp10 complex interface. Microscale thermophoresis (MST) experiments were used to quantify fragment affinities for nsp10. Additionally, we showed by MST that the interaction by nsp14 and 10 is weak and thereby that complex formation could be disrupted by small molecules. The fragments will serve as starting points for the development of more potent analogues using fragment growing techniques and structure-based drug design.

Discovery of the 4-aminopiperidine-based compound EM127 for the site-specific covalent inhibition of SMYD3.

Recent findings support the hypothesis that inhibition of SMYD3 methyltransferase may be a therapeutic avenue for some of the deadliest cancer types. Herein, active site-selective covalent SMYD3 inhibitors were designed by introducing an appropriate reactive cysteine trap into reversible first-generation SMYD3 inhibitors. The 4-aminopiperidine derivative EM127 (11C) bearing a 2-chloroethanoyl group as reactive warhead showed selectivity for Cys186, located in the substrate/histone binding pocket. Selectivity towards Cys186 was retained even at high inhibitor/enzyme ratio, as shown by mass spectrometry. The mode of interaction with the SMYD3 substrate/histone binding pocket was revealed by crystallographic studies. In enzymatic assays, 11C showed a stronger SMYD3 inhibitory effect compared to the reference inhibitor EPZ031686. Remarkably, 11C attenuated the proliferation of MDA-MB-231 breast cancer cell line at the same low micromolar range of concentrations that reduced SMYD3 mediated ERK signaling in HCT116 colorectal cancer and MDA-MB-231 breast cancer cells. Furthermore, 11C (5 µM) strongly decreased the steady-state mRNA levels of genes important for tumor biology such as cyclin dependent kinase 2, c-MET, N-cadherin and fibronectin 1, all known to be regulated, at least in part, by SMYD3. Thus, 11C is as a first example of second generation SMYD3 inhibitors; this agent represents a covalent and a site specific SMYD3 binder capable of potent and prolonged attenuation of methyltransferase activity.

FragMAXapp: crystallographic fragment-screening data-analysis and project-management system.

Crystallographic fragment screening (CFS) has become one of the major techniques for screening compounds in the early stages of drug-discovery projects. Following the advances in automation and throughput at modern macromolecular crystallography beamlines, the bottleneck for CFS has shifted from collecting data to organizing and handling the analysis of such projects. The complexity that emerges from the use of multiple methods for processing and refinement and to search for ligands requires an equally sophisticated solution to summarize the output, allowing researchers to focus on the scientific questions instead of on software technicalities. FragMAXapp is the fragment-screening project-management tool designed to handle CFS projects at MAX IV Laboratory. It benefits from the powerful computing infrastructure of large-scale facilities and, as a web application, it is accessible from everywhere.

FragMAX: the fragment-screening platform at the MAX IV Laboratory.

Advances in synchrotron storage rings and beamline automation have pushed data-collection rates to thousands of data sets per week. With this increase in throughput, massive projects such as in-crystal fragment screening have become accessible to a larger number of research groups. The quality of support offered at large-scale facilities allows medicinal chemistry-focused or biochemistry-focused groups to supplement their research with structural biology. Preparing the experiment, analysing multiple data sets and prospecting for interesting complexes of protein and fragments require, for both newcomers and experienced users, efficient management of the project and extensive computational power for data processing and structure refinement. Here, FragMAX, a new complete platform for fragment screening at the BioMAX beamline of the MAX IV Laboratory, is described. The ways in which users are assisted in X-ray-based fragment screenings and in which the fourth-generation storage ring available at the facility is best exploited are also described.

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.

Introduction of Intrinsic Kinetics of Protein-Ligand Interactions and Their Implications for Drug Design.

Structure-kinetic relationship analyses and identification of dominating interactions for optimization of lead compounds should ideally be based on intrinsic rate constants instead of the more easily accessible observed kinetic constants, which also account for binding-linked reactions. The intrinsic rate constants for sulfonamide inhibitors and pharmacologically relevant isoforms of carbonic anhydrase were determined by a novel surface plasmon resonance (SPR) biosensor-based approach, using chemodynamic analysis of binding-linked pH-dependent effects. The observed association rates ( kaobs) were pH-dependent and correlated with the fraction of deprotonated inhibitor and protonated zinc-bound water molecule. The intrinsic association rate constants ( kaintr) were pH independent and higher than kaobs. By contrast, the observed and intrinsic dissociation rate constants were identical and pH-independent, demonstrating that the observed association and dissociation mechanisms are inherently different. A model accounting for the differences between intrinsic and observed rate constants was developed, useful also for other interactions with binding-linked protonation reactions.

Kinetically Selective Inhibitors of Human Carbonic Anhydrase Isozymes I, II, VII, IX, XII, and XIII.

To get a better understanding of the possibility of developing selective carbonic anhydrase (CA) inhibitors, interactions between 17 benzenesulphonamide ligands and 6 human CAs (full-length CA I, II, VII, and XIII and catalytic domains of CA IX and XII) were characterized using surface plasmon resonance and fluorescent-based thermal shift assays. Kinetics revealed that the strongest binders had subnanomolar affinities with low dissociation rates (i.e., kd values around 1 × 10(-3) s(-1)) or were essentially irreversible. Chemodynamic analysis of the interactions highlighted an intrinsic mechanism of the CA-sulphonamide interaction kinetics and showed that slow dissociation rates were mediated by large hydrophobic contacts. The studied inhibitors demonstrated a high cross-reactivity within the protein family. However, according to chemical phylogenetic analysis developed for kinetic data, several ligands were found to be selective against certain CA isozymes, indicating that it should be possible to develop selective CA inhibitors suitable for clinical use.

Proximity-dependent initiation of hybridization chain reaction.

Sensitive detection of protein interactions and post-translational modifications of native proteins is a challenge for research and diagnostic purposes. A method for this, which could be used in point-of-care devices and high-throughput screening, should be reliable, cost effective and robust. To achieve this, here we design a method (proxHCR) that combines the need for proximal binding with hybridization chain reaction (HCR) for signal amplification. When two oligonucleotide hairpins conjugated to antibodies bind in close proximity, they can be activated to reveal an initiator sequence. This starts a chain reaction of hybridization events between a pair of fluorophore-labelled oligonucleotide hairpins, generating a fluorescent product. In conclusion, we show the applicability of the proxHCR method for the detection of protein interactions and posttranslational modifications in microscopy and flow cytometry. As no enzymes are needed, proxHCR may be an inexpensive and robust alternative to proximity ligation assays.