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.

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.

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.

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.

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.

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.

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.

Relative impact of 3- and 5-hydroxyl groups of cytosporone B on cancer cell viability.

A novel and the shortest route, thus far, for preparing cytosporone B (Csn-B) is reported. Csn-B and two analogs were used to probe the importance of hydroxyl groups at the 3- and 5-positions of the Csn-B benzene ring in inhibiting the viability of human H460 lung cancer and LNCaP prostate cancer cells, inducing H460 cell apoptosis, and interacting with the NR4A1 (TR3) ligand-binding domain (LBD). These studies indicate that Csn-B and 5-Me-Csn-B, having a phenolic hydroxyl at the 3-position of their aromatic rings, had similar activities in inhibiting cancer cell viability and in inducing apoptosis, whereas 3,5-(Me)2-Csn-B was unable to do so. These results are in agreement with ligand-binding experiments showing that the interaction with the NR4A1 LBD required the presence of the 3-hydroxyl group.

Utrophin ABD binds to F-actin in an open conformation.

Structural analyses of actin binding regions comprising tandem calponin homology domains alone and when bound to F-actin have revealed a number of different conformations with calponin homology domains in 'open' and 'closed' positions. In an attempt to resolve these issues we have examined the properties of the utrophin actin binding domain in open and closed conformations in order to verify the conformation when bound to F-actin. Locking the actin binding domain in a closed conformation using engineered cysteine residues in each calponin homology domain reduced the affinity for F-actin without affecting the stoichiometry furthermore differential scanning calorimetry experiments revealed a reduction in melting temperature on binding to actin. The data suggest the amino-terminal utrophin actin binding domain is in an open conformation in solution and when bound to F-actin.

Targeted drug delivery to tumor vasculature by a carbohydrate mimetic peptide.

Although numerous carbohydrates play significant roles in mammalian cells, carbohydrate-based drug discovery has not been explored due to the technical difficulty of chemically synthesizing complex carbohydrate structures. Previously, we identified a series of carbohydrate mimetic peptides and found that a 7-mer peptide, designated I-peptide, inhibits hematogenous carbohydrate-dependent cancer cell colonization. During analysis of the endothelial surface receptor for I-peptide, we found that I-peptide bound to annexin 1 (Anxa1). Because Anxa1 is a highly specific tumor vasculature surface marker, we hypothesized that an I-peptide-like peptide could target anticancer drugs to the tumor vasculature. This study identifies IFLLWQR peptide, designated IF7, as homing to tumors. When synthetic IF7 peptide was conjugated to fluorescent Alexa 488 (A488) and injected intravenously into tumor-bearing mice, IF7-A488 targeted tumors within minutes. IF7 conjugated to the potent anticancer drug SN-38 and injected intravenously into nude mice carrying human colon HCT116 tumors efficiently suppressed tumor growth at low dosages with no apparent side effects. These results suggest that IF7 serves as an efficient drug delivery vehicle by targeting Anxa1 expressed on the surface of tumor vasculature. Given its extremely specific tumor-targeting activity, IF7 may represent a clinically relevant vehicle for anticancer drugs.

A mechanism of Rap1-induced stabilization of endothelial cell--cell junctions.

Activation of Rap1 small GTPases stabilizes cell--cell junctions, and this activity requires Krev Interaction Trapped gene 1 (KRIT1). Loss of KRIT1 disrupts cardiovascular development and causes autosomal dominant familial cerebral cavernous malformations. Here we report that native KRIT1 protein binds the effector loop of Rap1A but not H-Ras in a GTP-dependent manner, establishing that it is an authentic Rap1-specific effector. By modeling the KRIT1-Rap1 interface we designed a well-folded KRIT1 mutant that exhibited a ~40-fold-reduced affinity for Rap1A and maintained other KRIT1-binding functions. Direct binding of KRIT1 to Rap1 stabilized endothelial cell-cell junctions in vitro and was required for cardiovascular development in vivo. Mechanistically, Rap1 binding released KRIT1 from microtubules, enabling it to locate to cell--cell junctions, where it suppressed Rho kinase signaling and stabilized the junctions. These studies establish that the direct physical interaction of Rap1 with KRIT1 enables the translocation of microtubule-sequestered KRIT1 to junctions, thereby supporting junctional integrity and cardiovascular development.

Vaccinia virus virulence factor N1L is a novel promising target for antiviral therapeutic intervention.

The 14 kDa homodimeric N1L protein is a potent vaccinia and variola (smallpox) virulence factor. It is not essential for viral replication, but it causes a strong attenuation of viral production in culture when deleted. The N1L protein is predicted to contain the BH3-like binding domain characteristic of Bcl-2 family proteins, and it is able to bind the BH3 peptides. Its overexpression has been reported to prevent infected cells from committing apoptosis. Therefore, interfering with the N1L apoptotic blockade may be a legitimate therapeutic strategy affecting the viral growth. By using in silico ligand docking and an array of in vitro assays, we have identified submicromolar (600 nM) N1L antagonists belonging to the family of polyphenols. Their affinity is comparable to that of the BH3 peptides (70-1000 nM). We have also identified the natural polyphenol resveratrol as a moderate N1L inhibitor. Finally, we show that our ligands efficiently inhibit growth of vaccinia virus.