Crystal structures of human and mouse ketohexokinase provide a structural basis for species- and isoform-selective inhibitor design.

A molecular understanding of the proteins involved in fructose metabolism is essential for controlling the current spread of fructose-related obesity, diabetes and related adverse metabolic states in Western populations. Fructose catabolism starts with the phosphorylation of D-fructose to fructose 1-phosphate by ketohexokinase (KHK). KHK exists in two alternatively spliced isoforms: the hepatic and intestinal isoform KHK-C and the peripheral isoform KHK-A. Here, the structure of apo murine KHK (mKHK), which differs from structures of human KHK in overall conformation, is reported. An isoform-selective ligand, which offers a 50-fold higher potency on mKHK and human KHK-A compared with KHK-C, is further characterized. In mKHK, large-scale conformational changes are observed upon ligand binding. The structures suggest a combined strategy for the design of species- and isoform-selective KHK inhibitors.

Molecular basis of TASL recruitment by the peptide/histidine transporter 1, PHT1.

PHT1 is a histidine /oligopeptide transporter with an essential role in Toll-like receptor innate immune responses. It can act as a receptor by recruiting the adaptor protein TASL which leads to type I interferon production via IRF5. Persistent stimulation of this signalling pathway is known to be involved in the pathogenesis of systemic lupus erythematosus (SLE). Understanding how PHT1 recruits TASL at the molecular level, is therefore clinically important for the development of therapeutics against SLE and other autoimmune diseases. Here we present the Cryo-EM structure of PHT1 stabilized in the outward-open conformation. By combining biochemical and structural modeling techniques we propose a model of the PHT1-TASL complex, in which the first 16 N-terminal TASL residues fold into a helical structure that bind in the central cavity of the inward-open conformation of PHT1. This work provides critical insights into the molecular basis of PHT1/TASL mediated type I interferon production.

Dichlorophenylpyridine-Based Molecules Inhibit Furin through an Induced-Fit Mechanism.

Inhibitors of the proprotein convertase furin might serve as broad-spectrum antiviral therapeutics. High cellular potency and antiviral activity against acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have been reported for (3,5-dichlorophenyl)pyridine-derived furin inhibitors. Here we characterized the binding mechanism of this inhibitor class using structural, biophysical, and biochemical methods. We established a MALDI-TOF-MS-based furin activity assay, determined IC50 values, and solved X-ray structures of (3,5-dichlorophenyl)pyridine-derived compounds in complex with furin. The inhibitors induced a substantial conformational rearrangement of the active-site cleft by exposing a central buried tryptophan residue. These changes formed an extended hydrophobic surface patch where the 3,5-dichlorophenyl moiety of the inhibitors was inserted into a newly formed binding pocket. Consistent with these structural rearrangements, we observed slow off-rate binding kinetics and strong structural stabilization in surface plasmon resonance and differential scanning fluorimetry experiments, respectively. The discovered furin conformation offers new opportunities for structure-based drug discovery.

Pathological polyQ expansion does not alter the conformation of the Huntingtin-HAP40 complex.

The abnormal amplification of a CAG repeat in the gene coding for huntingtin (HTT) leads to Huntington's disease (HD). At the protein level, this translates into the expansion of a polyglutamine (polyQ) stretch located at the HTT N terminus, which renders HTT aggregation prone by unknown mechanisms. Here we investigated the effects of polyQ expansion on HTT in a complex with its stabilizing interaction partner huntingtin-associated protein 40 (HAP40). Surprisingly, our comprehensive biophysical, crosslinking mass spectrometry and cryo-EM experiments revealed no major differences in the conformation of HTT-HAP40 complexes of various polyQ length, including 17QHTT-HAP40 (wild type), 46QHTT-HAP40 (typical polyQ length in HD patients), and 128QHTT-HAP40 (extreme polyQ length). Thus, HTT polyQ expansion does not alter the global conformation of HTT when associated with HAP40.

Covalent inhibitor reactivity prediction by the electrophilicity index-in and out of scope.

Drug discovery is an expensive and time-consuming process. To make this process more efficient quantum chemistry methods can be employed. The electrophilicity index is one property that can be calculated by quantum chemistry methods, and if calculated correctly gives insight into the reactivity of covalent inhibitors. Herein we present the usage of the electrophilicity index on three common warheads, i.e., acrylamides, 2-chloroacetamides, and propargylamides. We thoroughly examine the properties of the electrophilicity index, show which pitfalls should be avoided, and what the requirements to successfully apply the electrophilicity index are.

BIreactive: A Machine-Learning Model to Estimate Covalent Warhead Reactivity.

In the past decade, the pharmaceutical industry has paid closer attention to covalent drugs. Differently from standard noncovalent drugs, these compounds can exhibit peculiar properties, such as higher potency or longer duration of target inhibition with a potentially lower dosage. These properties are mainly driven by the reactive functional group present in the compound, the so-called warhead that forms a covalent bond with a specific nucleophilic amino-acid on the target. In this work, we report the possibility to combine ab initio activation energies with machine-learning to estimate covalent compound intrinsic reactivity. The idea behind this approach is to have a precise estimation of the transition state barriers, and thus of the compound reactivity, but with the speed of a machine-learning algorithm. We call this method "BIreactive". Here, we demonstrate this approach on acrylamides and 2-chloroacetamides, two warhead classes that possess different reaction mechanisms. In combination with our recently implemented truncation algorithm, we also demonstrate the possibility to use BIreactive not only for fragments but also for lead-like molecules. The generic nature of this approach allows also the extension to several other warheads. The combination of these factors makes BIreactive a valuable tool for the covalent drug discovery process in a pharmaceutical context.

A Fast Ab Initio Predictor Tool for Covalent Reactivity Estimation of Acrylamides.

Thanks to their unique mode of action, covalent drugs represent an exceptional opportunity for drug design. After binding to a biologically relevant target system, covalent compounds form a reversible or irreversible covalent bond with a nucleophilic amino acid. Due to the inherently large binding energy of a covalent bond, covalent binders exhibit higher potencies and thus allow potentially lower drug dosages. However, a proper balancing of compound reactivity is key for the design of covalent binders, to achieve high levels of target inhibition while minimizing promiscuous covalent binding to nontarget proteins. In this work, we demonstrated the possibility to apply the electrophilicity index concept to estimate covalent compound reactivity. We tested this approach on acrylamides, one of the most prominent classes of covalent warheads. Our study clearly demonstrated that, for compounds with molecular weight (MW) below 250 Da, the electrophilicity index can be directly used to estimate compound reactivity. On the other hand, for leadlike molecules (MW > 250 Da) we implemented a new truncation algorithm that has to be applied before reactivity calculations. This algorithm can ensure the localization of HOMO/LUMO orbitals on the compound warhead and thus a correct estimation of its reactivity. Our results also indicate that caution should be used when employing the electrophilicity index to estimate the reactivity of nonterminal acrylamides. The nonparametric nature of this method and its reasonable computational cost make it a suitable tool to support covalent drug design.

Crystal Structure of CC Chemokine Receptor 2A in Complex with an Orthosteric Antagonist Provides Insights for the Design of Selective Antagonists.

We determined two crystal structures of the chemokine receptor CCR2A in complex with the orthosteric antagonist MK-0812. Full-length CCR2A, stabilized by rubredoxin and a series of five mutations were resolved at 3.3 Å. An N- and C-terminally truncated CCR2A construct was crystallized in an alternate crystal form, which yielded a 2.7 Å resolution structure using serial synchrotron crystallography. Our structures provide a clear structural explanation for the observed key role of residue E2917.39 in high-affinity binding of several orthosteric CCR2 antagonists. By combining all the structural information collected, we generated models of co-structures for the structurally diverse pyrimidine amide class of CCR2 antagonists. Even though the representative Ex15 overlays well with MK-0812, it also interacts with the non-conserved H1213.33, resulting in a significant selectivity over CCR5. Insights derived from this work will facilitate drug discovery efforts directed toward highly selective CCR2 antagonists with potentially superior efficacy.

The cryo-electron microscopy structure of huntingtin.

Huntingtin (HTT) is a large (348 kDa) protein that is essential for embryonic development and is involved in diverse cellular activities such as vesicular transport, endocytosis, autophagy and the regulation of transcription. Although an integrative understanding of the biological functions of HTT is lacking, the large number of identified HTT interactors suggests that it serves as a protein-protein interaction hub. Furthermore, Huntington's disease is caused by a mutation in the HTT gene, resulting in a pathogenic expansion of a polyglutamine repeat at the amino terminus of HTT. However, only limited structural information regarding HTT is currently available. Here we use cryo-electron microscopy to determine the structure of full-length human HTT in a complex with HTT-associated protein 40 (HAP40; encoded by three F8A genes in humans) to an overall resolution of 4 Å. HTT is largely α-helical and consists of three major domains. The amino- and carboxy-terminal domains contain multiple HEAT (huntingtin, elongation factor 3, protein phosphatase 2A and lipid kinase TOR) repeats arranged in a solenoid fashion. These domains are connected by a smaller bridge domain containing different types of tandem repeats. HAP40 is also largely α-helical and has a tetratricopeptide repeat-like organization. HAP40 binds in a cleft and contacts the three HTT domains by hydrophobic and electrostatic interactions, thereby stabilizing the conformation of HTT. These data rationalize previous biochemical results and pave the way for improved understanding of the diverse cellular functions of HTT.

Structure of the Human Protein Kinase ZAK in Complex with Vemurafenib.

The mixed lineage kinase ZAK is a key regulator of the MAPK pathway mediating cell survival and inflammatory response. ZAK is targeted by several clinically approved kinase inhibitors, and inhibition of ZAK has been reported to protect from doxorubicin-induced cardiomyopathy. On the other hand, unintended targeting of ZAK has been linked to severe adverse effects such as the development of cutaneous squamous cell carcinoma. Therefore, both specific inhibitors of ZAK, as well as anticancer drugs lacking off-target activity against ZAK, may provide therapeutic benefit. Here, we report the first crystal structure of ZAK in complex with the B-RAF inhibitor vemurafenib. The cocrystal structure displayed a number of ZAK-specific features including a highly distorted P loop conformation enabling rational inhibitor design. Positional scanning peptide library analysis revealed a unique substrate specificity of the ZAK kinase including unprecedented preferences for histidine residues at positions -1 and +2 relative to the phosphoacceptor site. In addition, we screened a library of clinical kinase inhibitors identifying several inhibitors that potently inhibit ZAK, demonstrating that this kinase is commonly mistargeted by currently used anticancer drugs.

Mode of action of nintedanib in the treatment of idiopathic pulmonary fibrosis.

Idiopathic pulmonary fibrosis (IPF) is a progressive and ultimately fatal disease characterised by fibrosis of the lung parenchyma and loss of lung function. Although the pathogenic pathways involved in IPF have not been fully elucidated, IPF is believed to be caused by repetitive alveolar epithelial cell injury and dysregulated repair, in which there is uncontrolled proliferation of lung fibroblasts and differentiation of fibroblasts into myofibroblasts, which excessively deposit extracellular matrix (ECM) proteins in the interstitial space. A number of profibrotic mediators including platelet-derived growth factor (PDGF), fibroblast growth factor (FGF) and transforming growth factor-ß are believed to play important roles in the pathogenesis of IPF. Nintedanib is a potent small molecule inhibitor of the receptor tyrosine kinases PDGF receptor, FGF receptor and vascular endothelial growth factor receptor. Data from in vitro studies have shown that nintedanib interferes with processes active in fibrosis such as fibroblast proliferation, migration and differentiation, and the secretion of ECM. In addition, nintedanib has shown consistent anti-fibrotic and anti-inflammatory activity in animal models of lung fibrosis. These data provide a strong rationale for the clinical efficacy of nintedanib in patients with IPF, which has recently been demonstrated in phase III clinical trials.

Fast native-SAD phasing for routine macromolecular structure determination.

We describe a data collection method that uses a single crystal to solve X-ray structures by native SAD (single-wavelength anomalous diffraction). We solved the structures of 11 real-life examples, including a human membrane protein, a protein-DNA complex and a 266-kDa multiprotein-ligand complex, using this method. The data collection strategy is suitable for routine structure determination and can be implemented at most macromolecular crystallography synchrotron beamlines.

Crystallizing Membrane Proteins in the Lipidic Mesophase. Experience with Human Prostaglandin E2 Synthase 1 and an Evolving Strategy.

The lipidic mesophase or in meso method for crystallizing membrane proteins has several high profile targets to its credit and is growing in popularity. Despite its success, the method is in its infancy as far as rational crystallogenesis is concerned. Consequently, significant time, effort, and resources are still required to generate structure-grade crystals, especially with a new target type. Therefore, a need exists for crystallogenesis protocols that are effective with a broad range of membrane protein types. Recently, a strategy for crystallizing a prokaryotic α-helical membrane protein, diacylglycerol kinase (DgkA), by the in meso method was reported (Cryst. Growth. Des.2013, 14, 2846-2857). Here, we describe its application to the human α-helical microsomal prostaglandin E2 synthase 1 (mPGES1). While the DgkA strategy proved useful, significant modifications were needed to generate structure-quality crystals of this important therapeutic target. These included protein engineering, using an additive phospholipid in the hosting mesophase, performing multiple rounds of salt screening, and carrying out trials at 4 °C in the presence of a tight binding ligand. The crystallization strategy detailed here should prove useful for generating structures of other integral membrane proteins by the in meso method.

Crystal structure of glucokinase regulatory protein.

Glucokinase (GK) plays a major role in the regulation of blood glucose homeostasis in both the liver and the pancreas. In the liver, GK is controlled by the GK regulatory protein (GKRP). GKRP in turn is activated by fructose 6-phosphate (F6P) and inactivated by fructose 1-phosphate (F1P). Disrupting the GK-GKRP complex increases the activity of GK in the cytosol and is considered an attractive concept for the regulation of blood glucose. We have determined the crystal structure of GKRP in its inactive F1P-bound form. The binding site for F1P is located deeply buried at a domain interface, and H-D exchange experiments confirmed that F1P and F6P compete for this site. The structure of the inactive GKRP-F1P complex provides a starting point for understanding the mechanism of fructose phosphate-dependent GK regulation at an atomic level.

A systematic approach to increase the efficiency of membrane protein production in cell-free expression systems.

High amounts of membrane protein samples are needed for structural or functional analysis and a first bottleneck is often to obtain sufficient production efficiencies. The reduced complexity of protein production in cell-free expression systems results in a frequent correlation of efficiency problems with the essential transcription/translation process. We present a systematic tag variation strategy for the rapid improvement of cell-free expression efficiencies of membrane proteins based on the optimization of translation initiation. A small number of rationally designed short expression tags is attached via overlap PCR to the 5-prime end of the target protein coding sequence. The generated pool of DNA templates is analyzed in a cell-free expression screen and the most efficient template is selected for further preparative scale protein production. The expression tags can be minimized to only a few codons and no further impact on the coding sequence is required. The complete process takes only few days and the synthesized PCR fragments can be used directly as templates for preparative scale cell-free reactions. The strategy is exemplified with the production of a set of G-protein coupled receptors and yield improvements of up to 32-fold were obtained. All proteins were finally synthesized in amounts sufficient for further quality optimization and initial crystallization screens.

Structure of a nanobody-stabilized active state of the ß(2) adrenoceptor.

G protein coupled receptors (GPCRs) exhibit a spectrum of functional behaviours in response to natural and synthetic ligands. Recent crystal structures provide insights into inactive states of several GPCRs. Efforts to obtain an agonist-bound active-state GPCR structure have proven difficult due to the inherent instability of this state in the absence of a G protein. We generated a camelid antibody fragment (nanobody) to the human ß(2) adrenergic receptor (ß(2)AR) that exhibits G protein-like behaviour, and obtained an agonist-bound, active-state crystal structure of the receptor-nanobody complex. Comparison with the inactive ß(2)AR structure reveals subtle changes in the binding pocket; however, these small changes are associated with an 11 Å outward movement of the cytoplasmic end of transmembrane segment 6, and rearrangements of transmembrane segments 5 and 7 that are remarkably similar to those observed in opsin, an active form of rhodopsin. This structure provides insights into the process of agonist binding and activation.

The Implication of the First Agonist Bound Activated GPCR X-ray Structure on GPCR in Silico Modeling.

The very recently published first X-ray structure of the ß2 adrenergic receptor in its active state hosting a small molecule (PDB ID: 3P0G) reveals a lot of information about the G-protein-coupled receptor (GPCR) activation process from a structural point of view. When compared to the inactive state crystal structure of ß2, large differences are seen in the GPCR helical structure at the cytoplasmatic side, whereas very subtle changes occur at the ligand binding site. The observation that there are hardly any differences in the binding site of agonists and inverse agonists implies that in silico predictions of the efficacy of ligands will be very hard. This is illustrated by the example of an already published binding mode of a ß2 agonist, which has been modeled into the inactive state X-ray structure of the ß2 receptor. When comparing the modeled structure to the new activated X-ray structure, quantitative agreement of the binding mode is found, implying that the subtle changes between agonist binding to the activated state and inverse agonist binding to the inactive state can currently not be captured by standard in silico modeling methods.

Molecular recognition of the protein phosphatase 1 glycogen targeting subunit by glycogen phosphorylase.

Disrupting the interaction between glycogen phosphorylase and the glycogen targeting subunit (G(L)) of protein phosphatase 1 is emerging as a novel target for the treatment of type 2 diabetes. To elucidate the molecular basis of binding, we have determined the crystal structure of liver phosphorylase bound to a G(L)-derived peptide. The structure reveals the C terminus of G(L) binding in a hydrophobically collapsed conformation to the allosteric regulator-binding site at the phosphorylase dimer interface. G(L) mimics interactions that are otherwise employed by the activator AMP. Functional studies show that G(L) binds tighter than AMP and confirm that the C-terminal Tyr-Tyr motif is the major determinant for G(L) binding potency. Our study validates the G(L)-phosphorylase interface as a novel target for small molecule interaction.

C3 exoenzymes, novel insights into structure and action of Rho-ADP-ribosylating toxins.

The family of C3-like exoenzymes comprises seven bacterial ADP-ribosyltransferases of different origin. The common hallmark of these exoenzymes is the selective N-ADP-ribosylation of the low molecular mass GTP-binding proteins RhoA, B, and C and inhibition of signal pathways controlled by Rho GTPases. Therefore, C3-like exoenzymes were applied as pharmacological tools for analyses of cellular functions of Rho protein in numerous studies. Recent structural and functional analyses of C3-like exoenzymes provide detailed information on the molecular mechanisms and functional consequences of ADP-ribosylation catalyzed by these toxins. More recently additional non-enzymatic actions of C3-like ADP-ribosyltransferases have been identified showing that C3 transferases from Clostridium botulinum and Clostridium limosum form a GDI-like complex with the Ras-like low molecular mass GTPase Ral without ADP-ribosylation. These results add novel information on the molecular mode of action(s) of C3-like exoenzymes and are discussed in this review.

Crystal structure of the C3bot-RalA complex reveals a novel type of action of a bacterial exoenzyme.

C3 exoenzymes from bacterial pathogens ADP-ribosylate and inactivate low-molecular-mass GTPases of the Rho subfamily. Ral, a Ras subfamily GTPase, binds the C3 exoenzymes from Clostridium botulinum and C. limosum with high affinity without being a substrate for ADP ribosylation. In the complex, the ADP-ribosyltransferase activity of C3 is blocked, while binding of NAD and NAD-glycohydrolase activity remain. Here we report the crystal structure of C3 from C. botulinum in a complex with GDP-bound RalA at 1.8 A resolution. C3 binds RalA with a helix-loop-helix motif that is adjacent to the active site. A quaternary complex with NAD suggests a mode for ADP-ribosyltransferase inhibition. Interaction of C3 with RalA occurs at a unique interface formed by the switch-II region, helix alpha3 and the P loop of the GTPase. C3-binding stabilizes the GDP-bound conformation of RalA and blocks nucleotide release. Our data indicate that C. botulinum exoenzyme C3 is a single-domain toxin with bifunctional properties targeting Rho GTPases by ADP ribosylation and Ral by a guanine nucleotide dissociation inhibitor-like effect, which blocks nucleotide exchange.