Allosteric regulation and crystallographic fragment screening of SARS-CoV-2 NSP15 endoribonuclease

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of coronavirus disease 2019 (COVID-19). The NSP15 endoribonuclease enzyme, known as NendoU, is highly conserved and plays a critical role in the ability of the virus to evade the immune system. NendoU is a promising target for the development of new antiviral drugs. However, the complexity of the enzyme's structure and kinetics, along with the broad range of recognition sequences and lack of structural complexes, hampers the development of inhibitors. Here, we performed enzymatic characterization of NendoU in its monomeric and hexameric form, showing that hexamers are allosteric enzymes with a positive cooperative index, and with no influence of manganese on enzymatic activity. Through combining cryo-electron microscopy at different pHs, X-ray crystallography and biochemical and structural analysis, we showed that NendoU can shift between open and closed forms, which probably correspond to active and inactive states, respectively. We also explored the possibility of NendoU assembling into larger supramolecular structures and proposed a mechanism for allosteric regulation. In addition, we conducted a large fragment screening campaign against NendoU and identified several new allosteric sites that could be targeted for the development of new inhibitors. Overall, our findings provide insights into the complex structure and function of NendoU and offer new opportunities for the development of inhibitors.

Of problems and opportunities-How to treat and how to not treat crystallographic fragment screening data

In their recent commentary in Protein Science, Jaskolski et al. analyzed three randomly picked diffraction data sets from fragment-screening group depositions from the PDB and, based on that, they claimed that such data are principally problematic. We demonstrate here that if such data are treated properly, none of the proclaimed criticisms persist.

Novel Starting Points for Human Glycolate Oxidase Inhibitors, Revealed by Crystallography-Based Fragment Screening

Primary hyperoxaluria type I (PH1) is caused by AGXT gene mutations that decrease the functional activity of alanine:glyoxylate aminotransferase. A build-up of the enzyme’s substrate, glyoxylate, results in excessive deposition of calcium oxalate crystals in the renal tract, leading to debilitating renal failure. Oxidation of glycolate by glycolate oxidase (or hydroxy acid oxidase 1, HAO1) is a major cellular source of glyoxylate, and siRNA studies have shown phenotypic rescue of PH1 by the knockdown of HAO1, representing a promising inhibitor target. Here, we report the discovery and optimization of six low-molecular-weight fragments, identified by crystallography-based fragment screening, that bind to two different sites on the HAO1 structure: at the active site and an allosteric pocket above the active site. The active site fragments expand known scaffolds for substrate-mimetic inhibitors to include more chemically attractive molecules. The allosteric fragments represent the first report of non-orthosteric inhibition of any hydroxy acid oxidase and hold significant promise for improving inhibitor selectivity. The fragment hits were verified to bind and inhibit HAO1 in solution by fluorescence-based activity assay and surface plasmon resonance. Further optimization cycle by crystallography and biophysical assays have generated two hit compounds of micromolar (44 and 158 µM) potency that do not compete with the substrate and provide attractive starting points for the development of potent and selective HAO1 inhibitors.

Experiences From Developing Software for Large X-Ray Crystallography-Driven Protein-Ligand Studies

The throughput of macromolecular X-ray crystallography experiments has surged over the last decade. This remarkable gain in efficiency has been facilitated by increases in the availability of high-intensity X-ray beams, (ultra)fast detectors and high degrees of automation. These developments have in turn spurred the development of several dedicated centers for crystal-based fragment screening which enable the preparation and collection of hundreds of single-crystal diffraction datasets per day. Crystal structures of target proteins in complex with small-molecule ligands are of immense importance for structure-based drug design (SBDD) and their rapid turnover is a prerequisite for accelerated development cycles. While the experimental part of the process is well defined and has by now been established at several synchrotron sites, it is noticeable that software and algorithmic aspects have received far less attention, as well as the implications of new methodologies on established paradigms for structure determination, analysis, and visualization. We will review three key areas of development of large-scale protein-ligand studies. First, we will look into new software developments for batch data processing, followed by a discussion of the methodological changes in the analysis, modeling, refinement and deposition of structures for SBDD, and the changes in mindset that these new methods require, both on the side of depositors and users of macromolecular models. Finally, we will highlight key new developments for the presentation and analysis of the collections of structures that these experiments produce, and provide an outlook for future developments.

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 Novel BRD4 Ligand Scaffolds by Automated Navigation of the Fragment Chemical Space

Fragment-based drug discovery (FBDD) is a very effective hit identification method. However, the evolution of fragment hits into suitable leads remains challenging and largely artisanal. Fragment evolution is often scaffold-centric, meaning that its outcome depends crucially on the chemical structure of the starting fragment. Considering that fragment screening libraries cover only a small proportion of the corresponding chemical space, hits should be seen as probes highlighting privileged areas of the chemical space rather than actual starting points. We have developed an automated computational pipeline to mine the chemical space around any specific fragment hit, rapidly finding analogues that share a common interaction motif but are structurally novel and diverse. On a prospective application on the bromodomain-containing protein 4 (BRD4), starting from a known fragment, the platform yields active molecules with nonobvious scaffold changes. The procedure is fast and inexpensive and has the potential to uncover many hidden opportunities in FBDD.

SARS-CoV-2 infects the human kidney and drives fibrosis in kidney organoids

Kidney failure is frequently observed during and after COVID-19, but it remains elusive whether this is a direct effect of the virus. Here, we report that SARS-CoV-2 directly infects kidney cells and is associated with increased tubule-interstitial kidney fibrosis in patient autopsy samples. To study direct effects of the virus on the kidney independent of systemic effects of COVID-19, we infected human-induced pluripotent stem-cell-derived kidney organoids with SARS-CoV-2. Single-cell RNA sequencing indicated injury and dedifferentiation of infected cells with activation of profibrotic signaling pathways. Importantly, SARS-CoV-2 infection also led to increased collagen 1 protein expression in organoids. A SARS-CoV-2 protease inhibitor was able to ameliorate the infection of kidney cells by SARS-CoV-2. Our results suggest that SARS-CoV-2 can directly infect kidney cells and induce cell injury with subsequent fibrosis. These data could explain both acute kidney injury in COVID-19 patients and the development of chronic kidney disease in long COVID.

•

COVID-19 patients present tubulo-interstitial kidney fibrosis compared with controls

•

SARS-CoV-2 infection stimulates profibrotic signaling in human kidney organoids

•

SARS-CoV-2 infection can be inhibited by a protease blocker in human kidney organoids

Achieving Efficient Fragment Screening at XChem Facility at Diamond Light Source

In fragment-based drug discovery, hundreds or often thousands of compounds smaller than ~300 Da are tested against the protein of interest to identify chemical entities that can be developed into potent drug candidates. Since the compounds are small, interactions are weak, and the screening method must therefore be highly sensitive; moreover, structural information tends to be crucial for elaborating these hits into lead-like compounds. Therefore, protein crystallography has always been a gold-standard technique, yet historically too challenging to find widespread use as a primary screen. Initial XChem experiments were demonstrated in 2014 and then trialed with academic and industrial collaborators to validate the process. Since then, a large research effort and significant beamtime have streamlined sample preparation, developed a fragment library with rapid follow-up possibilities, automated and improved the capability of I04-1 beamline for unattended data collection, and implemented new tools for data management, analysis and hit identification. XChem is now a facility for large-scale crystallographic fragment screening, supporting the entire crystals-to-deposition process, and accessible to academic and industrial users worldwide. The peer-reviewed academic user program has been actively developed since 2016, to accommodate projects from as broad a scientific scope as possible, including well-validated as well as exploratory projects. Academic access is allocated through biannual calls for peer-reviewed proposals, and proprietary work is arranged by Diamond's Industrial Liaison group. This workflow has already been routinely applied to over a hundred targets from diverse therapeutic areas, and effectively identifies weak binders (1%-30% hit rate), which both serve as high-quality starting points for compound design and provide extensive structural information on binding sites. The resilience of the process was demonstrated by continued screening of SARS-CoV-2 targets during the COVID-19 pandemic, including a 3-week turn-around for the main protease.

Fragment Screening Reveals Starting Points for Rational Design of Galactokinase 1 Inhibitors to Treat Classic Galactosemia

Classic galactosemia is caused by loss-of-function mutations in galactose-1-phosphate uridylyltransferase (GALT) that lead to toxic accumulation of its substrate, galactose-1-phosphate. One proposed therapy is to inhibit the biosynthesis of galactose-1-phosphate, catalyzed by galactokinase 1 (GALK1). Existing inhibitors of human GALK1 (hGALK1) are primarily ATP-competitive with limited clinical utility to date. Here, we determined crystal structures of hGALK1 bound with reported ATP-competitive inhibitors of the spiro-benzoxazole series, to reveal their binding mode in the active site. Spurred by the need for additional chemotypes of hGALK1 inhibitors, desirably targeting a nonorthosteric site, we also performed crystallography-based screening by soaking hundreds of hGALK1 crystals, already containing active site ligands, with fragments from a custom library. Two fragments were found to bind close to the ATP binding site, and a further eight were found in a hotspot distal from the active site, highlighting the strength of this method in identifying previously uncharacterized allosteric sites. To generate inhibitors of improved potency and selectivity targeting the newly identified binding hotspot, new compounds were designed by merging overlapping fragments. This yielded two micromolar inhibitors of hGALK1 that were not competitive with respect to either substrate (ATP or galactose) and demonstrated good selectivity over hGALK1 homologues, galactokinase 2 and mevalonate kinase. Our findings are therefore the first to demonstrate inhibition of hGALK1 from an allosteric site, with potential for further development of potent and selective inhibitors to provide novel therapeutics for classic galactosemia.

Fragment binding to the Nsp3 macrodomain of SARS-CoV-2 identified through crystallographic screening and computational docking

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) macrodomain within the nonstructural protein 3 counteracts host-mediated antiviral adenosine diphosphate-ribosylation signaling. This enzyme is a promising antiviral target because catalytic mutations render viruses nonpathogenic. Here, we report a massive crystallographic screening and computational docking effort, identifying new chemical matter primarily targeting the active site of the macrodomain. Crystallographic screening of 2533 diverse fragments resulted in 214 unique macrodomain-binders. An additional 60 molecules were selected from docking more than 20 million fragments, of which 20 were crystallographically confirmed. X-ray data collection to ultra-high resolution and at physiological temperature enabled assessment of the conformational heterogeneity around the active site. Several fragment hits were confirmed by solution binding using three biophysical techniques (differential scanning fluorimetry, homogeneous time-resolved fluorescence, and isothermal titration calorimetry). The 234 fragment structures explore a wide range of chemotypes and provide starting points for development of potent SARS-CoV-2 macrodomain inhibitors.

Structure and activation mechanism of the human liver-type glutaminase GLS2

Cancer cells exhibit an altered metabolic phenotype, consuming higher levels of the amino acid glutamine. This metabolic reprogramming depends on increased mitochondrial glutaminase activity to convert glutamine to glutamate, an essential precursor for bioenergetic and biosynthetic processes in cells. Mammals encode the kidney-type (GLS) and liver-type (GLS2) glutaminase isozymes. GLS is overexpressed in cancer and associated with enhanced malignancy. On the other hand, GLS2 is either a tumor suppressor or an oncogene, depending on the tumor type. The GLS structure and activation mechanism are well known, while the structural determinants for GLS2 activation remain elusive. Here, we describe the structure of the human glutaminase domain of GLS2, followed by the functional characterization of the residues critical for its activity. Increasing concentrations of GLS2 lead to tetramer stabilization, a process enhanced by phosphate. In GLS2, the so-called “lid loop” is in a rigid open conformation, which may be related to its higher affinity for phosphate and lower affinity for glutamine; hence, it has lower glutaminase activity than GLS. The lower affinity of GLS2 for glutamine is also related to its less electropositive catalytic site than GLS, as indicated by a Thr225Lys substitution within the catalytic site decreasing the GLS2 glutamine concentration corresponding to half-maximal velocity (K_0.5). Finally, we show that the Lys253Ala substitution (corresponding to the Lys320Ala in the GLS “activation” loop, formerly known as the “gating” loop) renders a highly active protein in stable tetrameric form. We conclude that the “activation” loop, a known target for GLS inhibition, may also be a drug target for GLS2.

•

A positive cooperative mechanism of activation is demonstrated for the liver-type glutaminase.

•

The novel crystal structure for the glutaminase domain of human GLS2 is presented.

•

Site-directed mutagenesis suggest the activation loop as a hotspot for inhibitor design.

The low-cost Shifter microscope stage transforms the speed and robustness of protein crystal harvesting

Despite the tremendous success of X-ray cryo-crystallography in recent decades, the transfer of crystals from the drops in which they are grown to diffractometer sample mounts remains a manual process in almost all laboratories. Here, the Shifter, a motorized, interactive microscope stage that transforms the entire crystal-mounting workflow from a rate-limiting manual activity to a controllable, high-throughput semi-automated process, is described. By combining the visual acuity and fine motor skills of humans with targeted hardware and software automation, it was possible to transform the speed and robustness of crystal mounting. Control software, triggered by the operator, manoeuvres crystallization plates beneath a clear protective cover, allowing the complete removal of film seals and thereby eliminating the tedium of repetitive seal cutting. The software, either upon request or working from an imported list, controls motors to position crystal drops under a hole in the cover for human mounting at a microscope. The software automatically captures experimental annotations for uploading to the user's data repository, removing the need for manual documentation. The Shifter facilitates mounting rates of 100-240 crystals per hour in a more controlled process than manual mounting, which greatly extends the lifetime of the drops and thus allows a dramatic increase in the number of crystals retrievable from any given drop without loss of X-ray diffraction quality. In 2015, the first in a series of three Shifter devices was deployed as part of the XChem fragment-screening facility at Diamond Light Source, where they have since facilitated the mounting of over 120 000 crystals. The Shifter was engineered to have a simple design, providing a device that could be readily commercialized and widely adopted owing to its low cost. The versatile hardware design allows use beyond fragment screening and protein crystallography.

Deliberately Losing Control of C-H Activation Processes in the Design of Small-Molecule-Fragment Arrays Targeting Peroxisomal Metabolism

Combined photochemical arylation, "nuisance effect" (SN Ar) reaction sequences have been employed in the design of small arrays for immediate deployment in medium-throughput X-ray protein-ligand structure determination. Reactions were deliberately allowed to run "out of control" in terms of selectivity; for example the ortho-arylation of 2-phenylpyridine gave five products resulting from mono- and bisarylations combined with SN Ar processes. As a result, a number of crystallographic hits against NUDT7, a key peroxisomal CoA ester hydrolase, have been identified.

Discovery of allosteric binding sites by crystallographic fragment screening

Understanding allosteric regulation of proteins is fundamental to our study of protein structure and function. Moreover, allosteric binding pockets have become a major target of drug discovery efforts in recent years. However, even though the function of almost every protein can be influenced by allostery, it remains a challenge to discover, rationalise and validate putative allosteric binding pockets. This review examines how the discovery and analysis of putative allosteric binding sites have been influenced by the availability of centralised facilities for crystallographic fragment screening, along with newly developed computational methods for modelling low occupancy features. We discuss the experimental parameters required for success, and how new methods could influence the field in the future. Finally, we reflect on the general problem of how to translate these findings into actual ligand development programs.

Fragment Binding to the Nsp3 Macrodomain of SARS-CoV-2 Identified Through Crystallographic Screening and Computational Docking

The SARS-CoV-2 macrodomain (Mac1) within the non-structural protein 3 (Nsp3) counteracts host-mediated antiviral ADP-ribosylation signalling. This enzyme is a promising antiviral target because catalytic mutations render viruses non-pathogenic. Here, we report a massive crystallographic screening and computational docking effort, identifying new chemical matter primarily targeting the active site of the macrodomain. Crystallographic screening of diverse fragment libraries resulted in 214 unique macrodomain-binding fragments, out of 2,683 screened. An additional 60 molecules were selected from docking over 20 million fragments, of which 20 were crystallographically confirmed. X-ray data collection to ultra-high resolution and at physiological temperature enabled assessment of the conformational heterogeneity around the active site. Several crystallographic and docking fragment hits were validated for solution binding using three biophysical techniques (DSF, HTRF, ITC). Overall, the 234 fragment structures presented explore a wide range of chemotypes and provide starting points for development of potent SARS-CoV-2 macrodomain inhibitors.

Crystallographic and electrophilic fragment screening of the SARS-CoV-2 main protease

COVID-19, caused by SARS-CoV-2, lacks effective therapeutics. Additionally, no antiviral drugs or vaccines were developed against the closely related coronavirus, SARS-CoV-1 or MERS-CoV, despite previous zoonotic outbreaks. To identify starting points for such therapeutics, we performed a large-scale screen of electrophile and non-covalent fragments through a combined mass spectrometry and X-ray approach against the SARS-CoV-2 main protease, one of two cysteine viral proteases essential for viral replication. Our crystallographic screen identified 71 hits that span the entire active site, as well as 3 hits at the dimer interface. These structures reveal routes to rapidly develop more potent inhibitors through merging of covalent and non-covalent fragment hits; one series of low-reactivity, tractable covalent fragments were progressed to discover improved binders. These combined hits offer unprecedented structural and reactivity information for on-going structure-based drug design against SARS-CoV-2 main protease.

Demonstration of the utility of DOS-derived fragment libraries for rapid hit derivatisation in a multidirectional fashion

Organic synthesis underpins the evolution of weak fragment hits into potent lead compounds. Deficiencies within current screening collections often result in the requirement of significant synthetic investment to enable multidirectional fragment growth, limiting the efficiency of the hit evolution process. Diversity-oriented synthesis (DOS)-derived fragment libraries are constructed in an efficient and modular fashion and thus are well-suited to address this challenge. To demonstrate the effective nature of such libraries within fragment-based drug discovery, we herein describe the screening of a 40-member DOS library against three functionally distinct biological targets using X-Ray crystallography. Firstly, we demonstrate the importance for diversity in aiding hit identification with four fragment binders resulting from these efforts. Moreover, we also exemplify the ability to readily access a library of analogues from cheap commercially available materials, which ultimately enabled the exploration of a minimum of four synthetic vectors from each molecule. In total, 10-14 analogues of each hit were rapidly accessed in three to six synthetic steps. Thus, we showcase how DOS-derived fragment libraries enable efficient hit derivatisation and can be utilised to remove the synthetic limitations encountered in early stage fragment-based drug discovery.

Aspartate/asparagine-β-hydroxylase crystal structures reveal an unexpected epidermal growth factor-like domain substrate disulfide pattern

AspH is an endoplasmic reticulum (ER) membrane-anchored 2-oxoglutarate oxygenase whose C-terminal oxygenase and tetratricopeptide repeat (TPR) domains present in the ER lumen. AspH catalyses hydroxylation of asparaginyl- and aspartyl-residues in epidermal growth factor-like domains (EGFDs). Here we report crystal structures of human AspH, with and without substrate, that reveal substantial conformational changes of the oxygenase and TPR domains during substrate binding. Fe(II)-binding by AspH is unusual, employing only two Fe(II)-binding ligands (His679/His725). Most EGFD structures adopt an established fold with a conserved Cys1–3, 2–4, 5–6 disulfide bonding pattern; an unexpected Cys3–4 disulfide bonding pattern is observed in AspH-EGFD substrate complexes, the catalytic relevance of which is supported by studies involving stable cyclic peptide substrate analogues and by effects of Ca(II) ions on activity. The results have implications for EGFD disulfide pattern processing in the ER and will enable medicinal chemistry efforts targeting human 2OG oxygenases.

AspH catalyses hydroxylation of asparagine and aspartate residues in epidermal growth factor-like domains (EGFDs). Here, the authors present crystal structures of AspH with and without substrates and show that AspH uses EFGD substrates with a non-canonical disulfide pattern.

Rapid Covalent-Probe Discovery by Electrophile-Fragment Screening

Covalent probes can display unmatched potency, selectivity, and duration of action; however, their discovery is challenging. In principle, fragments that can irreversibly bind their target can overcome the low affinity that limits reversible fragment screening, but such electrophilic fragments were considered nonselective and were rarely screened. We hypothesized that mild electrophiles might overcome the selectivity challenge and constructed a library of 993 mildly electrophilic fragments. We characterized this library by a new high-throughput thiol-reactivity assay and screened them against 10 cysteine-containing proteins. Highly reactive and promiscuous fragments were rare and could be easily eliminated. In contrast, we found hits for most targets. Combining our approach with high-throughput crystallography allowed rapid progression to potent and selective probes for two enzymes, the deubiquitinase OTUB2 and the pyrophosphatase NUDT7. No inhibitors were previously known for either. This study highlights the potential of electrophile-fragment screening as a practical and efficient tool for covalent-ligand discovery.

A Parsimonious Mechanism of Sugar Dehydration by Human GDP-Mannose-4,6-dehydratase

Biosynthesis of 6-deoxy sugars, including l-fucose, involves a mechanistically complex, enzymatic 4,6-dehydration of hexose nucleotide precursors as the first committed step. Here, we determined pre- and postcatalytic complex structures of the human GDP-mannose 4,6-dehydratase at atomic resolution. These structures together with results of molecular dynamics simulation and biochemical characterization of wildtype and mutant enzymes reveal elusive mechanistic details of water elimination from GDP-mannose C5″ and C6″, coupled to NADP-mediated hydride transfer from C4″ to C6″. We show that concerted acid–base catalysis from only two active-site groups, Tyr₁₇₉ and Glu₁₅₇, promotes a syn 1,4-elimination from an enol (not an enolate) intermediate. We also show that the overall multistep catalytic reaction involves the fewest position changes of enzyme and substrate groups and that it proceeds under conserved exploitation of the basic (minimal) catalytic machinery of short-chain dehydrogenase/reductases.

Partial-occupancy binders identified by the Pan-Dataset Density Analysis method offer new chemical opportunities and reveal cryptic binding sites

Crystallographic fragment screening uses low molecular weight compounds to probe the protein surface and although individual protein-fragment interactions are high quality, fragments commonly bind at low occupancy, historically making identification difficult. However, our new Pan-Dataset Density Analysis method readily identifies binders missed by conventional analysis: for fragment screening data of lysine-specific demethylase 4D (KDM4D), the hit rate increased from 0.9% to 10.6%. Previously unidentified fragments reveal multiple binding sites and demonstrate: the versatility of crystallographic fragment screening; that surprisingly large conformational changes are possible in crystals; and that low crystallographic occupancy does not by itself reflect a protein-ligand complex's significance.

Gentle, fast and effective crystal soaking by acoustic dispensing

A high-throughput method is described for crystal soaking using acoustic droplet ejection, and its effectiveness is demonstrated.

The steady expansion in the capacity of modern beamlines for high-throughput data collection, enabled by increasing X-ray brightness, capacity of robotics and detector speeds, has pushed the bottleneck upstream towards sample preparation. Even in ligand-binding studies using crystal soaking, the experiment best able to exploit beamline capacity, a primary limitation is the need for gentle and nontrivial soaking regimens such as stepwise concentration increases, even for robust and well characterized crystals. Here, the use of acoustic droplet ejection for the soaking of protein crystals with small molecules is described, and it is shown that it is both gentle on crystals and allows very high throughput, with 1000 unique soaks easily performed in under 10 min. In addition to having very low compound consumption (tens of nanolitres per sample), the positional precision of acoustic droplet ejection enables the targeted placement of the compound/solvent away from crystals and towards drop edges, allowing gradual diffusion of solvent across the drop. This ensures both an improvement in the reproducibility of X-ray diffraction and increased solvent tolerance of the crystals, thus enabling higher effective compound-soaking concentrations. The technique is detailed here with examples from the protein target JMJD2D, a histone lysine demethylase with roles in cancer and the focus of active structure-based drug-design efforts.

Proper modelling of ligand binding requires an ensemble of bound and unbound states

Although noncovalent binding by small molecules cannot be assumed a priori to be stoichiometric in the crystal lattice, occupancy refinement of ligands is often avoided by convention. Occupancies tend to be set to unity, requiring the occupancy error to be modelled by the B factors, and residual weak density around the ligand is necessarily attributed to `disorder'. Where occupancy refinement is performed, the complementary, superposed unbound state is rarely modelled. Here, it is shown that superior accuracy is achieved by modelling the ligand as partially occupied and superposed on a ligand-free `ground-state' model. Explicit incorporation of this model of the crystal, obtained from a reference data set, allows constrained occupancy refinement with minimal fear of overfitting. Better representation of the crystal also leads to more meaningful refined atomic parameters such as the B factor, allowing more insight into dynamics in the crystal. An outline of an approach for algorithmically generating ensemble models of crystals is presented, assuming that data sets representing the ground state are available. The applicability of various electron-density metrics to the validation of the resulting models is assessed, and it is concluded that ensemble models consistently score better than the corresponding single-state models. Furthermore, it appears that ignoring the superposed ground state becomes the dominant source of model error, locally, once the overall model is accurate enough; modelling the local ground state properly is then more meaningful than correcting all remaining model errors globally, especially for low-occupancy ligands. Implications for the simultaneous refinement of B factors and occupancies, and for future evaluation of the limits of the approach, in particular its behaviour at lower data resolution, are discussed.

The XChemExplorer graphical workflow tool for routine or large-scale protein-ligand structure determination

XChemExplorer (XCE) is a data-management and workflow tool to support large-scale simultaneous analysis of protein-ligand complexes during structure-based ligand discovery (SBLD). The user interfaces of established crystallographic software packages such as CCP4 [Winn et al. (2011), Acta Cryst. D67, 235-242] or PHENIX [Adams et al. (2010), Acta Cryst. D66, 213-221] have entrenched the paradigm that a `project' is concerned with solving one structure. This does not hold for SBLD, where many almost identical structures need to be solved and analysed quickly in one batch of work. Functionality to track progress and annotate structures is essential. XCE provides an intuitive graphical user interface which guides the user from data processing, initial map calculation, ligand identification and refinement up until data dissemination. It provides multiple entry points depending on the need of each project, enables batch processing of multiple data sets and records metadata, progress and annotations in an SQLite database. XCE is freely available and works on any Linux and Mac OS X system, and the only dependency is to have the latest version of CCP4 installed. The design and usage of this tool are described here, and its usefulness is demonstrated in the context of fragment-screening campaigns at the Diamond Light Source. It is routinely used to analyse projects comprising 1000 data sets or more, and therefore scales well to even very large ligand-design projects.

Outcome of the First wwPDB/CCDC/D3R Ligand Validation Workshop

Crystallographic studies of ligands bound to biological macromolecules (proteins and nucleic acids) represent an important source of information concerning drug-target interactions, providing atomic level insights into the physical chemistry of complex formation between macromolecules and ligands. Of the more than 115,000 entries extant in the Protein Data Bank (PDB) archive, ∼75% include at least one non-polymeric ligand. Ligand geometrical and stereochemical quality, the suitability of ligand models for in silico drug discovery and design, and the goodness-of-fit of ligand models to electron-density maps vary widely across the archive. We describe the proceedings and conclusions from the first Worldwide PDB/Cambridge Crystallographic Data Center/Drug Design Data Resource (wwPDB/CCDC/D3R) Ligand Validation Workshop held at the Research Collaboratory for Structural Bioinformatics at Rutgers University on July 30-31, 2015. Experts in protein crystallography from academe and industry came together with non-profit and for-profit software providers for crystallography and with experts in computational chemistry and data archiving to discuss and make recommendations on best practices, as framed by a series of questions central to structural studies of macromolecule-ligand complexes. What data concerning bound ligands should be archived in the PDB? How should the ligands be best represented? How should structural models of macromolecule-ligand complexes be validated? What supplementary information should accompany publications of structural studies of biological macromolecules? Consensus recommendations on best practices developed in response to each of these questions are provided, together with some details regarding implementation. Important issues addressed but not resolved at the workshop are also enumerated.

A generic protocol for protein crystal dehydration using the HC1b humidity controller

Dehydration may change the crystal lattice and affect the mosaicity, resolution and quality of X-ray diffraction data. A dehydrating environment can be generated around a crystal in several ways with various degrees of precision and complexity. This study uses a high-precision crystal humidifier/dehumidifier to provide an airstream of known relative humidity in which the crystals are mounted: a precise yet hassle-free approach to altering crystal hydration. A protocol is introduced to assess the impact of crystal dehydration systematically applied to nine experimental crystal systems. In one case, that of glucose isomerase, dehydration triggering a change of space group from I222 to P21212 was observed. This observation is supported by an extended study of the behaviour of the glucose isomerase crystal structure during crystal dehydration.

An overview of heavy-atom derivatization of protein crystals

Heavy-atom derivatization is one of the oldest techniques for obtaining phase information for protein crystals and, although it is no longer the first choice, it remains a useful technique for obtaining phases for unknown structures and for low-resolution data sets. It is also valuable for confirming the chain trace in low-resolution electron-density maps. This overview provides a summary of the technique and is aimed at first-time users of the method. It includes guidelines on when to use it, which heavy atoms are most likely to work, how to prepare heavy-atom solutions, how to derivatize crystals and how to determine whether a crystal is in fact a derivative.

A poised fragment library enables rapid synthetic expansion yielding the first reported inhibitors of PHIP(2), an atypical bromodomain

Research into the chemical biology of bromodomains has been driven by the development of acetyl-lysine mimetics. The ligands are typically anchored by binding to a highly conserved asparagine residue. Atypical bromodomains, for which the asparagine is mutated, have thus far proven elusive targets, including PHIP(2) whose parent protein, PHIP, has been linked to disease progression in diabetes and cancers. The PHIP(2) binding site contains a threonine in place of asparagine, and solution screening have yielded no convincing hits. We have overcome this hurdle by combining the sensitivity of X-ray crystallography, used as the primary fragment screen, with a strategy for rapid follow-up synthesis using a chemically-poised fragment library, which allows hits to be readily modified by parallel chemistry both peripherally and in the core. Our approach yielded the first reported hit compounds of PHIP(2) with measurable IC50 values by an AlphaScreen competition assay. The follow-up libraries of four poised fragment hits improved potency into the sub-mM range while showing good ligand efficiency and detailed structural data.

8-Substituted Pyrido[3,4-d]pyrimidin-4(3H)-one Derivatives As Potent, Cell Permeable, KDM4 (JMJD2) and KDM5 (JARID1) Histone Lysine Demethylase Inhibitors

We report the discovery of N-substituted 4-(pyridin-2-yl)thiazole-2-amine derivatives and their subsequent optimization, guided by structure-based design, to give 8-(1H-pyrazol-3-yl)pyrido[3,4-d]pyrimidin-4(3H)-ones, a series of potent JmjC histone N-methyl lysine demethylase (KDM) inhibitors which bind to Fe(II) in the active site. Substitution from C4 of the pyrazole moiety allows access to the histone peptide substrate binding site; incorporation of a conformationally constrained 4-phenylpiperidine linker gives derivatives such as 54j and 54k which demonstrate equipotent activity versus the KDM4 (JMJD2) and KDM5 (JARID1) subfamily demethylases, selectivity over representative exemplars of the KDM2, KDM3, and KDM6 subfamilies, cellular permeability in the Caco-2 assay, and, for 54k, inhibition of H3K9Me3 and H3K4Me3 demethylation in a cell-based assay.

Docking and Linking of Fragments To Discover Jumonji Histone Demethylase Inhibitors

Development of tool molecules that inhibit Jumonji demethylases allows for the investigation of cancer-associated transcription. While scaffolds such as 2,4-pyridinedicarboxylic acid (2,4-PDCA) are potent inhibitors, they exhibit limited selectivity. To discover new inhibitors for the KDM4 demethylases, enzymes overexpressed in several cancers, we docked a library of 600,000 fragments into the high-resolution structure of KDM4A. Among the most interesting chemotypes were the 5-aminosalicylates, which docked in two distinct but overlapping orientations. Docking poses informed the design of covalently linked fragment compounds, which were further derivatized. This combined approach improved affinity by ∼ 3 log-orders to yield compound 35 (Ki = 43 nM). Several hybrid inhibitors were selective for KDM4C over the related enzymes FIH, KDM2A, and KDM6B while lacking selectivity against the KDM3 and KDM5 subfamilies. Cocrystal structures corroborated the docking predictions. This study extends the use of structure-based docking from fragment discovery to fragment linking optimization, yielding novel KDM4 inhibitors.

Structural basis of glycogen branching enzyme deficiency and pharmacologic rescue by rational peptide design

Glycogen branching enzyme 1 (GBE1) plays an essential role in glycogen biosynthesis by generating α-1,6-glucosidic branches from α-1,4-linked glucose chains, to increase solubility of the glycogen polymer. Mutations in the GBE1 gene lead to the heterogeneous early-onset glycogen storage disorder type IV (GSDIV) or the late-onset adult polyglucosan body disease (APBD). To better understand this essential enzyme, we crystallized human GBE1 in the apo form, and in complex with a tetra- or hepta-saccharide. The GBE1 structure reveals a conserved amylase core that houses the active centre for the branching reaction and harbours almost all GSDIV and APBD mutations. A non-catalytic binding cleft, proximal to the site of the common APBD mutation p.Y329S, was found to bind the tetra- and hepta-saccharides and may represent a higher-affinity site employed to anchor the complex glycogen substrate for the branching reaction. Expression of recombinant GBE1-p.Y329S resulted in drastically reduced protein yield and solubility compared with wild type, suggesting this disease allele causes protein misfolding and may be amenable to small molecule stabilization. To explore this, we generated a structural model of GBE1-p.Y329S and designed peptides ab initio to stabilize the mutation. As proof-of-principle, we evaluated treatment of one tetra-peptide, Leu-Thr-Lys-Glu, in APBD patient cells. We demonstrate intracellular transport of this peptide, its binding and stabilization of GBE1-p.Y329S, and 2-fold increased mutant enzymatic activity compared with untreated patient cells. Together, our data provide the rationale and starting point for the screening of small molecule chaperones, which could become novel therapies for this disease.

Optimisation of a triazolopyridine based histone demethylase inhibitor yields a potent and selective KDM2A (FBXL11) inhibitor

A potent inhibitor of the JmjC histone lysine demethylase KDM2A (compound 35, pIC50 7.2) with excellent selectivity over representatives from other KDM subfamilies has been developed; the discovery that a triazolopyridine compound binds to the active site of JmjC KDMs was followed by optimisation of the triazole substituent for KDM2A inhibition and selectivity.

Crystal structures of 2OG oxygenases involved in ribosomal protein hydroxylation

Post-translational modifications play diverse biological functions. Hydroxylation of collagen proteins has long been a recognised post-translational modification in eukaryotes. In the case of collagen, hydroxylation of prolyl residues, by 2-oxoglutarate and iron dependent enzymes (2OG oxygenases), in collagen proteins allows for the stabilisation of the collagen triple helix structure through conformational restraint and through the addition of a hydrogen bond donor. Additionally, hydroxylation of lysine side chains of collagen is required for cross-linking collagen (and possibly other proteins) in the extra-cellular matrix. Post-translational prolyl hydroxylation also plays a pivotal role in transcriptional regulation of the hypoxic response, as catalyzed by the hypoxia inducible factor / HIF prolyl hydroxylases (PHDs or EGLN enzymes). Recently, ribosomal protein hydroxylation catalyzed by 2OG- and Fe(II)-dependent oxygenases has been found to be a highly conserved post-translational modification in eukaryotes and prokaryotes (Ge et al and Loenarz et al). We present several crystal structures of 2OG oxygenases involved in ribosomal protein hydroxylation.

When protein crystal dehydration helps

A typical protein crystal contains 30-60% solvent. For a naked crystal, this solvent is distributed between solvent shells, where water and solvent molecules make specific interactions with the crystalline protein, and solvent channels filled with disordered solvent molecules. This internal solvent map of the crystal can be modified by placing the crystal in a dehydrating environment. This may in turn induce changes to the crystal lattice and affect mosaicity, resolution and quality of diffraction data. A dehydrating environment can be generated around a crystal in several ways with various degrees of precision and complexity. In this study we have used the HC1 device (Maatel) to mount crystals an air stream of known relative humidity – a precise yet hassle-free approach to altering crystal hydration. We set out to analyse a range of different crystals to establish usable protocols that will allow one to explore to crystal hydration space, either by preparing samples before synchrotron beamtime or by undertaking the experiments during beamtime. Our results, considered in the light of the literature surrounding crystal dehydration, provide guidance for when dehydration can help diffraction.

Virtual Screening of Histone Lysine Demethylase(JMJD2) identifies new inhibitors

The JmjC domain-containing proteins are hydroxylases that confer posttranslational modifications on histone tails, by removing methylation marks on methylated lysine residues. This serves to either promote or repress gene transcription. The JMJD2A-D family members include the enzyme Jumonji domain 2C (JMJD2C), which specifically demethylates di- and trimethylated histone H3 at Lys 9 or Lys 36.[1] Dysregulation of JMJD2C has been implicated in prostate, colonic, and breast cancer as the demethylase can modify the expression levels of oncogenes.[2] The goal of the present study was to identify potent and selective small-molecule inhibitors of JMJD2C, to be used as chemical biology tools to further investigate the role of JMJD2C in cell proliferation and survival. Using high-resolution crystal structures of the JMJD2 subfamily members as templates, we have performed a small molecule virtual docking screen. From the ~3 million molecules that were docked, this experiment identified 21 compounds as possible leads. These compounds were tested against JMJD2C in enzymatic assays and here we report an overall hit rate of 76%, with 8 compounds demonstrating an IC50 of 176μM to 1.18μM. A molecule containing a salicylate core was selected as a candidate for optimization and thus far we have completed several rounds of iterative target-specific compound docking, hybrid molecule design, compound synthesis and in vitro characterization. Notably, our method demonstrated a substantial increase in potency when we linked two docked fragments together and further derivatized this new scaffold, through which we have successfully derived a 65nM inhibitor of JMJD2C. A compound representing the inhibitor scaffold has been co-crystallized with JMJD2A to a resolution of 2.4 Å. In the crystal structure each asymmetric unit contains two JMJD2A monomers, each bound to a single inhibitor molecule. This complex-structure superposes well with the docked pose for the hybrid series of compounds. We are now focusing our efforts on identifying an inhibitor that is selective for the JMJD2 family over other JmjC domain-containing proteins.

I04-1, a future 'fragment screening' beamline facility at Diamond Light Source

I04-1 is one of the six macromolecular crystallography (MX) beamlines at Diamond Light Source (DLS), the third generation synchrotron light source in the UK. It was built and delivered in 2010 as a stable and reliable fixed-wavelength MX station. It is currently preparing to release its user programme for exploiting fragment screening using X-ray crystallography in structural medicinal chemistry projects. For this purpose, the beamline has been going through several upgrades in order to achieve unattended high-throughput ligand crystallography. The new developments are aiming at improving the flux, stability and reliability of the beamline and its auto-alignment. In parallel, a peripheral laboratory is being set up to provide a facility for medium throughput compound soaking. Jointly with the Structural Genomics Consortium (SGC), a semi-automatic crystal soaking and harvesting scheme, which will provide hundreds of MX samples per day, is being tested at DLS. The beamline can currently process 400 crystals per day. However, the recent upgrades and automation should further improve that throughput. In this presentation, we will summarise the current specifications of the beamline and its new features, the development of a peripheral laboratory for compounds soaking and underline the remaining work.

5-Carboxy-8-hydroxyquinoline is a Broad Spectrum 2-Oxoglutarate Oxygenase Inhibitor which Causes Iron Translocation

2-Oxoglutarate and iron dependent oxygenases are therapeutic targets for human diseases. Using a representative 2OG oxygenase panel, we compare the inhibitory activities of 5-carboxy-8-hydroxyquinoline (IOX1) and 4-carboxy-8-hydroxyquinoline (4C8HQ) with that of two other commonly used 2OG oxygenase inhibitors, N-oxalylglycine (NOG) and 2,4-pyridinedicarboxylic acid (2,4-PDCA). The results reveal that IOX1 has a broad spectrum of activity, as demonstrated by the inhibition of transcription factor hydroxylases, representatives of all 2OG dependent histone demethylase subfamilies, nucleic acid demethylases and γ-butyrobetaine hydroxylase. Cellular assays show that, unlike NOG and 2,4-PDCA, IOX1 is active against both cytosolic and nuclear 2OG oxygenases without ester derivatisation. Unexpectedly, crystallographic studies on these oxygenases demonstrate that IOX1, but not 4C8HQ, can cause translocation of the active site metal, revealing a rare example of protein ligand-induced metal movement.

Crystal structures of malonyl-coenzyme A decarboxylase provide insights into its catalytic mechanism and disease-causing mutations

Malonyl-coenzyme A decarboxylase (MCD) is found from bacteria to humans, has important roles in regulating fatty acid metabolism and food intake, and is an attractive target for drug discovery. We report here four crystal structures of MCD from human, Rhodopseudomonas palustris, Agrobacterium vitis, and Cupriavidus metallidurans at up to 2.3 Å resolution. The MCD monomer contains an N-terminal helical domain involved in oligomerization and a C-terminal catalytic domain. The four structures exhibit substantial differences in the organization of the helical domains and, consequently, the oligomeric states and intersubunit interfaces. Unexpectedly, the MCD catalytic domain is structurally homologous to those of the GCN5-related N-acetyltransferase superfamily, especially the curacin A polyketide synthase catalytic module, with a conserved His-Ser/Thr dyad important for catalysis. Our structures, along with mutagenesis and kinetic studies, provide a molecular basis for understanding pathogenic mutations and catalysis, as well as a template for structure-based drug design.

•

Structures of human and bacterial MCDs were determined at up to 2.3 Å resolution

•

Distinct tetrameric and dimeric MCD oligomerizations were observed

•

Unexpected homology to the GNAT superfamily gives insights into catalytic mechanism

•

The structures provide the molecular basis for the disease-causing mutations in MCD

Squeezing the most from every crystal: the fine details of data collection

Modern synchrotron beamlines offer instrumentation of unprecedented quality, which in turn encourages increasingly marginal experiments, and for these, as much as ever, the ultimate success of data collection depends on the experience, but especially the care, of the experimenter. A representative set of difficult cases has been encountered at the Structural Genomics Consortium, a worldwide structural genomics initiative of which the Oxford site currently deposits three novel human structures per month. Achieving this target relies heavily on frequent visits to the Diamond Light Source, and the variety of crystal systems still demand customized data collection, diligent checks and careful planning of each experiment. Here, an overview is presented of the techniques and procedures that have been refined over the years and that are considered synchrotron best practice.

Structural basis for Cul3 protein assembly with the BTB-Kelch family of E3 ubiquitin ligases

Cullin-RING ligases are multisubunit E3 ubiquitin ligases that recruit substrate-specific adaptors to catalyze protein ubiquitylation. Cul3-based Cullin-RING ligases are uniquely associated with BTB adaptors that incorporate homodimerization, Cul3 assembly, and substrate recognition into a single multidomain protein, of which the best known are BTB-BACK-Kelch domain proteins, including KEAP1. Cul3 assembly requires a BTB protein “3-box” motif, analogous to the F-box and SOCS box motifs of other Cullin-based E3s. To define the molecular basis for this assembly and the overall architecture of the E3, we determined the crystal structures of the BTB-BACK domains of KLHL11 both alone and in complex with Cul3, along with the Kelch domain structures of KLHL2 (Mayven), KLHL7, KLHL12, and KBTBD5. We show that Cul3 interaction is dependent on a unique N-terminal extension sequence that packs against the 3-box in a hydrophobic groove centrally located between the BTB and BACK domains. Deletion of this N-terminal region results in a 30-fold loss in affinity. The presented data offer a model for the quaternary assembly of this E3 class that supports the bivalent capture of Nrf2 and reveals potential new sites for E3 inhibitor design.

Plant growth regulator daminozide is a selective inhibitor of human KDM2/7 histone demethylases

The JmjC oxygenases catalyze the N-demethylation of N(ε)-methyl lysine residues in histones and are current therapeutic targets. A set of human 2-oxoglutarate analogues were screened using a unified assay platform for JmjC demethylases and related oxygenases. Results led to the finding that daminozide (N-(dimethylamino)succinamic acid, 160 Da), a plant growth regulator, selectively inhibits the KDM2/7 JmjC subfamily. Kinetic and crystallographic studies reveal that daminozide chelates the active site metal via its hydrazide carbonyl and dimethylamino groups.

Molecular adaptation of the DegQ protease to exert protein quality control in the bacterial cell envelope

To react to distinct stress situations and to prevent the accumulation of misfolded proteins, all cells employ a number of proteases and chaperones, which together set up an efficient protein quality control system. The functionality of proteins in the cell envelope of Escherichia coli is monitored by the HtrA proteases DegS, DegP, and DegQ. In contrast with DegP and DegS, the structure and function of DegQ has not been addressed in detail. Here, we show that substrate binding triggers the conversion of the resting DegQ hexamer into catalytically active 12- and 24-mers. Interestingly, substrate-induced oligomer reassembly and protease activation depends on the first PDZ domain but not on the second. Therefore, the regulatory mechanism originally identified in DegP should be a common feature of HtrA proteases, most of which encompass only a single PDZ domain. Using a DegQ mutant lacking the second PDZ domain, we determined the high resolution crystal structure of a dodecameric HtrA complex. The nearly identical domain orientation of protease and PDZ domains within 12- and 24-meric HtrA complexes reveals a conserved PDZ1 → L3 → LD/L1/L2 signaling cascade, in which loop L3 senses the repositioned PDZ1 domain of higher order, substrate-engaged particles and activates protease function. Furthermore, our in vitro and in vivo data imply a pH-related function of DegQ in the bacterial cell envelope.

Assessment of radiation damage behaviour in a large collection of empirically optimized datasets highlights the importance of unmeasured complicating effects

The radiation damage behaviour in 43 datasets of 34 different proteins collected over a year was examined, in order to gauge the reliability of decay metrics in practical situations, and to assess how these datasets, optimized only empirically for decay, would have benefited from the precise and automatic prediction of decay now possible with the programs RADDOSE [Murray, Garman & Ravelli (2004). J. Appl. Cryst. 37, 513-522] and BEST [Bourenkov & Popov (2010). Acta Cryst. D66, 409-419]. The results indicate that in routine practice the diffraction experiment is not yet characterized well enough to support such precise predictions, as these depend fundamentally on three interrelated variables which cannot yet be determined robustly and practically: the flux density distribution of the beam; the exact crystal volume; the sensitivity of the crystal to dose. The former two are not satisfactorily approximated from typical beamline information such as nominal beam size and transmission, or two-dimensional images of the beam and crystal; the discrepancies are particularly marked when using microfocus beams (<20 µm). Empirically monitoring decay with the dataset scaling B factor (Bourenkov & Popov, 2010) appears more robust but is complicated by anisotropic and/or low-resolution diffraction. These observations serve to delineate the challenges, scientific and logistic, that remain to be addressed if tools for managing radiation damage in practical data collection are to be conveniently robust enough to be useful in real time.