Structural Basis for Substrate Binding, Catalysis and Inhibition of Breast Cancer Target Mitochondrial Creatine Kinase by Covalent Inhibitor via Cryo-EM.

Mitochondrial creatine kinases are key players in maintaining energy homeostasis in cells by working in conjunction with cytosolic creatine kinases for energy transport from mitochondria to cytoplasm. High levels of MtCK observed in Her2+ breast cancer and inhibition of breast cancer cell growth by substrate analog, cyclocreatine, indicate dependence of cancer cells on the 'energy shuttle' for cell growth and survival. Hence, understanding the key mechanistic features of creatine kinases and their inhibition plays an important role in the development of cancer therapeutics. Herein, we present the mutational and structural investigation on understudied ubiquitous mitochondrial creatine kinase (uMtCK). Our cryo-EM structures and biochemical data on uMtCK showed closure of the loop comprising residue His61 is specific to and relies on creatine binding and the reaction mechanism of phosphoryl transfer depends on electrostatics in the active site. In addition, the previously identified covalent inhibitor CKi showed inhibition in breast cancer BT474 cells, however our biochemical and structural data indicated that CKi is not a potent inhibitor for breast cancer due to strong dependency on the covalent link formation and inability to induce conformational changes upon binding.

Active site remodeling in tumor-relevant IDH1 mutants drives distinct kinetic features and potential resistance mechanisms.

Mutations in human isocitrate dehydrogenase 1 (IDH1) drive tumor formation in a variety of cancers by replacing its conventional activity with a neomorphic activity that generates an oncometabolite. Little is understood of the mechanistic differences among tumor-driving IDH1 mutants. We previously reported that the R132Q mutant unusually preserves conventional activity while catalyzing robust oncometabolite production, allowing an opportunity to compare these reaction mechanisms within a single active site. Here, we employ static and dynamic structural methods and observe that, compared to R132H, the R132Q active site adopts a conformation primed for catalysis with optimized substrate binding and hydride transfer to drive improved conventional and neomorphic activity over R132H. This active site remodeling reveals a possible mechanism of resistance to selective mutant IDH1 therapeutic inhibitors. This work enhances our understanding of fundamental IDH1 mechanisms while pinpointing regions for improving inhibitor selectivity.

Active site remodeling in tumor-relevant IDH1 mutants drives distinct kinetic features and potential resistance mechanisms.

Mutations in human isocitrate dehydrogenase 1 (IDH1) drive tumor formation in a variety of cancers by replacing its conventional activity with a neomorphic activity that generates an oncometabolite. Little is understood of the mechanistic differences among tumor-driving IDH1 mutants. We previously reported that the R132Q mutant uniquely preserves conventional activity while catalyzing robust oncometabolite production, allowing an opportunity to compare these reaction mechanisms within a single active site. Here, we employed static and dynamic structural methods and found that, compared to R132H, the R132Q active site adopted a conformation primed for catalysis with optimized substrate binding and hydride transfer to drive improved conventional and neomorphic activity over R132H. This active site remodeling revealed a possible mechanism of resistance to selective mutant IDH1 therapeutic inhibitors. This work enhances our understanding of fundamental IDH1 mechanisms while pinpointing regions for improving inhibitor selectivity.

Active site remodeling in tumor-relevant IDH1 mutants drives distinct kinetic features and potential resistance mechanisms.

Mutations in human isocitrate dehydrogenase 1 (IDH1) drive tumor formation in a variety of cancers by replacing its conventional activity with a neomorphic activity that generates an oncometabolite. Little is understood of the mechanistic differences among tumor-driving IDH1 mutants. We previously reported that the R132Q mutant uniquely preserves conventional activity while catalyzing robust oncometabolite production, allowing an opportunity to compare these reaction mechanisms within a single active site. Here, we employed static and dynamic structural methods and found that, compared to R132H, the R132Q active site adopted a conformation primed for catalysis with optimized substrate binding and hydride transfer to drive improved conventional and neomorphic activity over R132H. This active site remodeling revealed a possible mechanism of resistance to selective mutant IDH1 therapeutic inhibitors. This work enhances our understanding of fundamental IDH1 mechanisms while pinpointing regions for improving inhibitor selectivity.

TMEM2 is a bona fide hyaluronidase possessing intrinsic catalytic activity.

Transmembrane protein 2 (TMEM2) was originally identified as a membrane-anchored protein of unknown function. We previously demonstrated that TMEM2 can degrade hyaluronan (HA). Furthermore, we showed that induced global knockout of Tmem2 in adult mice results in rapid accumulation of incompletely degraded HA in bodily fluids and organs, supporting the identity of TMEM2 as a cell surface hyaluronidase. In spite of these advances, no direct evidence has been presented to demonstrate the intrinsic hyaluronidase activity of TMEM2. Here, we directly establish the catalytic activity of TMEM2. The ectodomain of TMEM2 (TMEM2ECD) was expressed as a His-tagged soluble protein and purified by affinity and size-exclusion chromatography. Both human and mouse TMEM2ECD robustly degrade fluorescein-labeled HA into 5 to 10 kDa fragments. TMEM2ECD exhibits this HA-degrading activity irrespective of the species of TMEM2 origin and the position of epitope tag insertion. The HA-degrading activity of TMEM2ECD is more potent than that of HYAL2, a hyaluronidase which, like TMEM2, has been implicated in cell surface HA degradation. Finally, we show that TMEM2ECD can degrade not only fluorescein-labeled HA but also native high-molecular weight HA. In addition to these core findings, our study reveals hitherto unrecognized confounding factors, such as the quality of reagents and the choice of assay systems, that could lead to erroneous conclusions regarding the catalytic activity of TMEM2. In conclusion, our results demonstrate that TMEM2 is a legitimate functional hyaluronidase. Our findings also raise cautions regarding the choice of reagents and methods for performing degradation assays for hyaluronidases.

Engineering gain-of-function mutants of a WW domain by dynamics and structural analysis.

Proteins gain optimal fitness such as foldability and function through evolutionary selection. However, classical studies have found that evolutionarily designed protein sequences alone cannot guarantee foldability, or at least not without considering local contacts associated with the initial folding steps. We previously showed that foldability and function can be restored by removing frustration in the folding energy landscape of a model WW domain protein, CC16, which was designed based on Statistical Coupling Analysis (SCA). Substitutions ensuring the formation of five local contacts identified as "on-path" were selected using the closest homolog native folded sequence, N21. Surprisingly, the resulting sequence, CC16-N21, bound to Group I peptides, while N21 did not. Here, we identified single-point mutations that enable N21 to bind a Group I peptide ligand through structure and dynamic-based computational design. Comparison of the docked position of the CC16-N21/ligand complex with the N21 structure showed that residues at positions 9 and 19 are important for peptide binding, whereas the dynamic profiles identified position 10 as allosterically coupled to the binding site and exhibiting different dynamics between N21 and CC16-N21. We found that swapping these positions in N21 with matched residues from CC16-N21 recovers nature-like binding affinity to N21. This study validates the use of dynamic profiles as guiding principles for affecting the binding affinity of small proteins.

Design of novel cyanovirin-N variants by modulation of binding dynamics through distal mutations.

We develop integrated co-evolution and dynamic coupling (ICDC) approach to identify, mutate, and assess distal sites to modulate function. We validate the approach first by analyzing the existing mutational fitness data of TEM-1 ß-lactamase and show that allosteric positions co-evolved and dynamically coupled with the active site significantly modulate function. We further apply ICDC approach to identify positions and their mutations that can modulate binding affinity in a lectin, cyanovirin-N (CV-N), that selectively binds to dimannose, and predict binding energies of its variants through Adaptive BP-Dock. Computational and experimental analyses reveal that binding enhancing mutants identified by ICDC impact the dynamics of the binding pocket, and show that rigidification of the binding residues compensates for the entropic cost of binding. This work suggests a mechanism by which distal mutations modulate function through dynamic allostery and provides a blueprint to identify candidates for mutagenesis in order to optimize protein function.

Development of a Robust High-Throughput Screening Platform for Inhibitors of the Striatal-Enriched Tyrosine Phosphatase (STEP).

Many human diseases are the result of abnormal expression or activation of protein tyrosine kinases (PTKs) and protein tyrosine phosphatases (PTPs). Not surprisingly, more than 30 tyrosine kinase inhibitors (TKIs) are currently in clinical use and provide unique treatment options for many patients. PTPs on the other hand have long been regarded as "undruggable" and only recently have gained increased attention in drug discovery. Striatal-enriched tyrosine phosphatase (STEP) is a neuron-specific PTP that is overactive in Alzheimer's disease (AD) and other neurodegenerative and neuropsychiatric disorders, including Parkinson's disease, schizophrenia, and fragile X syndrome. An emergent model suggests that the increase in STEP activity interferes with synaptic function and contributes to the characteristic cognitive and behavioral deficits present in these diseases. Prior efforts to generate STEP inhibitors with properties that warrant clinical development have largely failed. To identify novel STEP inhibitor scaffolds, we developed a biophysical, label-free high-throughput screening (HTS) platform based on the protein thermal shift (PTS) technology. In contrast to conventional HTS using STEP enzymatic assays, we found the PTS platform highly robust and capable of identifying true hits with confirmed STEP inhibitory activity and selectivity. This new platform promises to greatly advance STEP drug discovery and should be applicable to other PTP targets.

Structural basis for the association of PLEKHA7 with membrane-embedded phosphatidylinositol lipids.

PLEKHA7 (pleckstrin homology domain containing family A member 7) plays key roles in intracellular signaling, cytoskeletal organization, and cell adhesion, and is associated with multiple human cancers. The interactions of its pleckstrin homology (PH) domain with membrane phosphatidyl-inositol-phosphate (PIP) lipids are critical for proper cellular localization and function, but little is known about how PLEKHA7 and other PH domains interact with membrane-embedded PIPs. Here we describe the structural basis for recognition of membrane-bound PIPs by PLEHA7. Using X-ray crystallography, nuclear magnetic resonance, molecular dynamics simulations, and isothermal titration calorimetry, we show that the interaction of PLEKHA7 with PIPs is multivalent, distinct from a discrete one-to-one interaction, and induces PIP clustering. Our findings reveal a central role of the membrane assembly in mediating protein-PIP association and provide a roadmap for understanding how the PH domain contributes to the signaling, adhesion, and nanoclustering functions of PLEKHA7.

Local Interactions That Contribute Minimal Frustration Determine Foldability.

Earlier experiments suggest that the evolutionary information (conservation and coevolution) encoded in protein sequences is necessary and sufficient to specify the fold of a protein family. However, there is no computational work to quantify the effect of such evolutionary information on the folding process. Here we explore the role of early folding steps for sequences designed using coevolution and conservation through a combination of computational and experimental methods. We simulated a repertoire of native and designed WW domain sequences to analyze early local contact formation and found that the N-terminal ß-hairpin turn would not form correctly due to strong non-native local contacts in unfoldable sequences. Through a maximum likelihood approach, we identified five local contacts that play a critical role in folding, suggesting that a small subset of amino acid pairs can be used to solve the "needle in the haystack" problem to design foldable sequences. Thus, using the contact probability of those five local contacts that form during the early stage of folding, we built a classification model that predicts the foldability of a WW sequence with 81% accuracy. This classification model was used to redesign WW domain sequences that could not fold due to frustration and make them foldable by introducing a few mutations that led to the stabilization of these critical local contacts. The experimental analysis shows that a redesigned sequence folds and binds to polyproline peptides with a similar affinity as those observed for native WW domains. Overall, our analysis shows that evolutionary-designed sequences should not only satisfy the folding stability but also ensure a minimally frustrated folding landscape.

Conformational States of the Cytoprotective Protein Bcl-xL.

Bcl-xL is a major inhibitor of apoptosis, a fundamental homeostatic process of programmed cell death that is highly conserved across evolution. Because it plays prominent roles in cancer, Bcl-xL is a major target for anticancer therapy and for studies aimed at understanding its structure and activity. Although Bcl-xL is active primarily at intracellular membranes, most studies have focused on soluble forms of the protein lacking both the membrane-anchoring C-terminal tail and the intrinsically disordered loop, and this has resulted in a fragmented view of the protein's biological activity. Here, we describe the conformation of full-length Bcl-xL. Using NMR spectroscopy, molecular dynamics simulations, and isothermal titration calorimetry, we show how the three structural elements affect the protein's structure, dynamics, and ligand-binding activity in both its soluble and membrane-anchored states. The combined data provide information about the molecular basis for the protein's functionality and a view of its complex molecular mechanisms.

Structurally plastic NEMO and oligomerization prone IKK2 subunits define the behavior of human IKK2:NEMO complexes in solution.

The human IκB Kinase (IKK) is a multisubunit protein complex of two kinases and one scaffolding subunit that controls induction of transcription factor NF-κB activity. IKK behaves as an entity of aberrantly high apparent molecular weight in solution. Recent X-ray crystallographic and cryo-electron microscopy structures of individual catalytic subunits (IKK1/IKKα and IKK2/IKKß) reveal that they are both stably folded dimeric proteins that engage in extensive homo-oligomerization through unique surfaces that are required for activation of their respective catalytic activities. The NEMO/IKKγ subunit is a predominantly coiled coil protein that is required for activation of IKK through the canonical NF-κB signaling pathway. Here we report size-exclusion chromatography, multi-angle light scattering, analytical centrifugation, and thermal denaturation analyses of full-length human recombinant NEMO as well as deletion and disease-linked variants. We observe that NEMO is predominantly a dimer in solution, although by virtue of its modular coiled coil regions NEMO exhibits complicated solution dynamics involving portions that are mutually antagonistic toward homodimerization. This behavior causes NEMO to behave as a significantly larger sized particle in solution. Analyses of NEMO in complex with IKK2 indicate that NEMO preserves this structurally dynamic character within the multisubuit complex and provides the complex-bound IKK2 further propensity toward homo-oligomerization. These observations provide critical information on the structural plasticity of NEMO subunit dimers which helps clarify its role in diseases and in IKK regulation through oligomerization-dependent phosphorylation of catalytic IKK2 subunit dimers.

An acidic residue buried in the dimer interface of isocitrate dehydrogenase 1 (IDH1) helps regulate catalysis and pH sensitivity.

Isocitrate dehydrogenase 1 (IDH1) catalyzes the reversible NADP+-dependent conversion of isocitrate to α-ketoglutarate (αKG) to provide critical cytosolic substrates and drive NADPH-dependent reactions like lipid biosynthesis and glutathione regeneration. In biochemical studies, the forward reaction is studied at neutral pH, while the reverse reaction is typically characterized in more acidic buffers. This led us to question whether IDH1 catalysis is pH-regulated, which would have functional implications under conditions that alter cellular pH, like apoptosis, hypoxia, cancer, and neurodegenerative diseases. Here, we show evidence of catalytic regulation of IDH1 by pH, identifying a trend of increasing kcat values for αKG production upon increasing pH in the buffers we tested. To understand the molecular determinants of IDH1 pH sensitivity, we used the pHinder algorithm to identify buried ionizable residues predicted to have shifted pKa values. Such residues can serve as pH sensors, with changes in protonation states leading to conformational changes that regulate catalysis. We identified an acidic residue buried at the IDH1 dimer interface, D273, with a predicted pKa value upshifted into the physiological range. D273 point mutations had decreased catalytic efficiency and, importantly, loss of pH-regulated catalysis. Based on these findings, we conclude that IDH1 activity is regulated, at least in part, by pH. We show this regulation is mediated by at least one buried acidic residue ∼12 Å from the IDH1 active site. By establishing mechanisms of regulation of this well-conserved enzyme, we highlight catalytic features that may be susceptible to pH changes caused by cell stress and disease.

Water Networks and Correlated Motions in Mutant Isocitrate Dehydrogenase 1 (IDH1) Are Critical for Allosteric Inhibitor Binding and Activity.

Point mutations in human isocitrate dehydrogenase 1 (IDH1) can drive malignancies, including lower-grade gliomas and secondary glioblastomas, chondrosarcomas, and acute myeloid leukemias. These mutations, which usually affect residue R132, ablate the normal activity of catalyzing the NADP+-dependent oxidation of isocitrate to α-ketoglutarate (αKG) while also acquiring a neomorphic activity of reducing αKG to d-2-hydroxyglutarate (D2HG). Mutant IDH1 can be selectively therapeutically targeted due to structural differences that occur in the wild type (WT) versus mutant form of the enzyme, though the full mechanisms of this selectivity are still under investigation. Here we probe the mechanistic features of the neomorphic activity and selective small molecule inhibition through a new lens, employing WaterMap and molecular dynamics simulations. These tools identified a high-energy path of water molecules connecting the inhibitor binding site with the αKG and NADP+ binding sites in mutant IDH1. This water path aligns spatially with the α10 helix from WT IDH1 crystal structures. Mutating residues at the termini of this water path specifically disrupted inhibitor binding and/or D2HG production, revealing additional key residues to consider in optimizing druglike molecules against mutant IDH1. Taken together, our findings from molecular simulations and mutant enzyme kinetic assays provide insight into how disrupting water paths through enzyme active sites can impact not only inhibitor potency but also substrate recognition and activity.

Molecular Inhibitor of QSOX1 Suppresses Tumor Growth In Vivo.

Quiescin sulfhydryl oxidase 1 (QSOX1) is an enzyme overexpressed by many different tumor types. QSOX1 catalyzes the formation of disulfide bonds in proteins. Because short hairpin knockdowns (KD) of QSOX1 have been shown to suppress tumor growth and invasion in vitro and in vivo, we hypothesized that chemical compounds inhibiting QSOX1 enzymatic activity would also suppress tumor growth, invasion, and metastasis. High throughput screening using a QSOX1-based enzymatic assay revealed multiple potential QSOX1 inhibitors. One of the inhibitors, known as "SBI-183," suppresses tumor cell growth in a Matrigel-based spheroid assay and inhibits invasion in a modified Boyden chamber, but does not affect viability of nonmalignant cells. Oral administration of SBI-183 inhibits tumor growth in 2 independent human xenograft mouse models of renal cell carcinoma. We conclude that SBI-183 warrants further exploration as a useful tool for understanding QSOX1 biology and as a potential novel anticancer agent in tumors that overexpress QSOX1.

Calcium sensing by the STIM1 ER-luminal domain.

Stromal interaction molecule 1 (STIM1) monitors ER-luminal Ca2+ levels to maintain cellular Ca2+ balance and to support Ca2+ signalling. The prevailing view has been that STIM1 senses reduced ER Ca2+ through dissociation of bound Ca2+ from a single EF-hand site, which triggers a dramatic loss of secondary structure and dimerization of the STIM1 luminal domain. Here we find that the STIM1 luminal domain has 5-6 Ca2+-binding sites, that binding at these sites is energetically coupled to binding at the EF-hand site, and that Ca2+ dissociation controls a switch to a second structured conformation of the luminal domain rather than protein unfolding. Importantly, the other luminal-domain Ca2+-binding sites interact with the EF-hand site to control physiological activation of STIM1 in cells. These findings fundamentally revise our understanding of physiological Ca2+ sensing by STIM1, and highlight molecular mechanisms that govern the Ca2+ threshold for activation and the steep Ca2+ concentration dependence.

Adapted ATPase domain communication overcomes the cytotoxicity of p97 inhibitors.

The AAA+ ATPase p97 regulates ubiquitin-dependent protein homeostasis and has been pursued as a cancer drug target. The ATP-competitive inhibitor CB-5083 and allosteric inhibitor NMS-873 are the most advanced p97 inhibitors described to date. Previous studies have reported that their cytotoxicity can be readily overcome and involves single p97 mutations in the linker between the D1 and D2 ATPase domains and within D2. We report here that the proline 472 to leucine (P472L) mutation, in the D1-D2 linker and identified in CB-5083-resistant cells, desensitizes p97 to both inhibitor classes. This mutation does not disrupt the distinct D2-binding sites of the inhibitors. Instead, P472L changes ATPase domain communication within the p97 hexamer. P472L enhances cooperative D2 ATP binding and hydrolysis. This mechanism alters the function of the D1-D2 linker in the control of D2 activity involving the ATP-bound state of D1. Although increased D2 activity is sufficient to desensitize the P472L mutant to NMS-873, the mutant's desensitization to CB-5083 also requires D1 ATPase domain function. Our study highlights the remarkable adaptability of p97 ATPase domain communication that enables escape from mechanistically distinct classes of cytotoxic p97 inhibitors.

Oxidized analogs of Di(1H-indol-3-yl)methyl-4-substituted benzenes are NR4A1-dependent UPR inducers with potent and safe anti-cancer activity.

Di(1H-indol-3-yl)(4-trifluoromethylphenyl)methane (DIM-Ph-4-CF3) is an analog of orphan nuclear receptor 4A1 (NR4A1) ligand cytosporone B. We have synthesized several oxidation products of DIM-Ph-4-CF3, focusing on analogs with electron-withdrawing or donating groups at their phenyl ring 4-positions, and examined their anti-cancer activity and mechanism-of-action. Mesylates (DIM-Ph-4-X+ OMs-s) having CF3, CO2Me and Cl groups were more effective inhibitors of cancer cell viability than their precursors. 19F NMR spectroscopy and differential scanning calorimetry strongly indicated interactions of DIM-Ph-4-CF3+ OMs- with the NR4A1 ligand binding domain, and compound-induced apoptosis of prostate cancer cells was dependent on NR4A1. DIM-Ph-4-CF3+ OMs- showed robust inhibition of LNCaP prostate cancer xenografts with no apparent toxicity. In vitro and in vivo, DIM-Ph-4-CF3+ OMs- activated proapoptotic unfolded protein response (UPR) signaling in prostate cancer cells. Independently of DIM-Ph-4-CF3+ OMs-, the bulk of NR4A1 localized to the cytoplasm in various cancer cell lines, suggesting a cytoplasmic mechanism-of-action of DIM-Ph-4-CF3+ OMs- in UPR induction and cell death. In summary, the data suggest that oxidized analogs of DIM-Ph-4-CF3 possess potent and safe anti-cancer activity which is mediated through UPR signaling downstream of NR4A1 binding.

Physiological temperatures reduce dimerization of dengue and Zika virus recombinant envelope proteins.

The spread of dengue (DENV) and Zika virus (ZIKV) is a major public health concern. The primary target of antibodies that neutralize DENV and ZIKV is the envelope (E) glycoprotein, and there is interest in using soluble recombinant E (sRecE) proteins as subunit vaccines. However, the most potent neutralizing antibodies against DENV and ZIKV recognize epitopes on the virion surface that span two or more E proteins. Therefore, to create effective DENV and ZIKV vaccines, presentation of these quaternary epitopes may be necessary. The sRecE proteins from DENV and ZIKV crystallize as native-like dimers, but studies in solution suggest that these dimers are marginally stable. To better understand the challenges associated with creating stable sRecE dimers, we characterized the thermostability of sRecE proteins from ZIKV and three DENV serotypes, DENV2-4. All four proteins irreversibly unfolded at moderate temperatures (46-53 °C). At 23 °C and low micromolar concentrations, DENV2 and ZIKV were primarily dimeric, and DENV3-4 were primarily monomeric, whereas at 37 °C, all four proteins were predominantly monomeric. We further show that the dissociation constant for DENV2 dimerization is very temperature-sensitive, ranging from <1 µm at 25 °C to 50 µm at 41 °C, due to a large exothermic enthalpy of binding of -79 kcal/mol. We also found that quaternary epitope antibody binding to DENV2-4 and ZIKV sRecE is reduced at 37 °C. Our observation of reduced sRecE dimerization at physiological temperature highlights the need for stabilizing the dimer as part of its development as a subunit vaccine.

Diabetes reversal by inhibition of the low-molecular-weight tyrosine phosphatase.

Obesity-associated insulin resistance plays a central role in type 2 diabetes. As such, tyrosine phosphatases that dephosphorylate the insulin receptor (IR) are potential therapeutic targets. The low-molecular-weight protein tyrosine phosphatase (LMPTP) is a proposed IR phosphatase, yet its role in insulin signaling in vivo has not been defined. Here we show that global and liver-specific LMPTP deletion protects mice from high-fat diet-induced diabetes without affecting body weight. To examine the role of the catalytic activity of LMPTP, we developed a small-molecule inhibitor with a novel uncompetitive mechanism, a unique binding site at the opening of the catalytic pocket, and an exquisite selectivity over other phosphatases. This inhibitor is orally bioavailable, and it increases liver IR phosphorylation in vivo and reverses high-fat diet-induced diabetes. Our findings suggest that LMPTP is a key promoter of insulin resistance and that LMPTP inhibitors would be beneficial for treating type 2 diabetes.