Search

Detailed genomic and epigenomic analyses of MECOM (the MDS1 and EVI1 complex locus) have revealed that inversion or translocation of chromosome 3 drives inv(3)/t(3;3) myeloid leukemias via structural rearrangement of an enhancer that upregulates transcription of EVI1. Here, we identify a novel, previously unannotated oncogenic RNA-splicing derived isoform of EVI1 that is frequently present in inv(3)/t(3;3) acute myeloid leukemia (AML) and directly contributes to leukemic transformation. This EVI1 isoform is generated by oncogenic mutations in the core RNA splicing factor SF3B1, which is mutated in >30% of inv(3)/t(3;3) myeloid neoplasm patients and thereby represents the single most commonly cooccurring genomic alteration in inv(3)/t(3;3) patients. SF3B1 mutations are statistically uniquely enriched in inv(3)/t(3;3) myeloid neoplasm patients and patient-derived cell lines compared with other forms of AML and promote mis-splicing of EVI1 generating an in-frame insertion of 6 amino acids at the 3′ end of the second zinc finger domain of EVI1. Expression of this EVI1 splice variant enhanced the self-renewal of hematopoietic stem cells, and introduction of mutant SF3B1 in mice bearing the humanized inv(3)(q21q26) allele resulted in generation of this novel EVI1 isoform in mice and hastened leukemogenesis in vivo. The mutant SF3B1 spliceosome depends upon an exonic splicing enhancer within EVI1 exon 13 to promote usage of a cryptic branch point and aberrant 3′ splice site within intron 12 resulting in the generation of this isoform. These data provide a mechanistic basis for the frequent cooccurrence of SF3B1 mutations as well as new insights into the pathogenesis of myeloid leukemias harboring inv(3)/t(3;3).

There has been accumulating evidence for the clinical benefit of chemoradiation therapy (CRT), whereas mechanisms in CRT-recurrent clones derived from the primary tumor are still elusive. Herein, we identified an aberrant BUB1B/BUBR1 expression in CRT-recurrent clones in bladder cancer (BC) by comprehensive proteomic analysis. CRT-recurrent BC cells exhibited a cell-cycle-independent upregulation of BUB1B/BUBR1 expression rendering an enhanced DNA repair activity in response to DNA double-strand breaks (DSBs). With DNA repair analyses employing the CRISPR/cas9 system, we revealed that cells with aberrant BUB1B/BUBR1 expression dominantly exploit mutagenic nonhomologous end joining (NHEJ). We further found that phosphorylated ATM interacts with BUB1B/BUBR1 after ionizing radiation (IR) treatment, and the resistance to DSBs by increased BUB1B/BUBR1 depends on the functional ATM. In vivo, tumor growth of CRT-resistant T24R cells was abrogated by ATM inhibition using AZD0156. A dataset analysis identified FOXM1 as a putative BUB1B/BUBR1-targeting transcription factor causing its increased expression. These data collectively suggest a redundant role of BUB1B/BUBR1 underlying mutagenic NHEJ in an ATM-dependent manner, aside from the canonical activity of BUB1B/BUBR1 on the G2/M checkpoint, and offer novel clues to overcome CRT resistance.

Ubiquitin-like with PHD and RING finger domain-containing protein 1 (UHRF1)-dependent DNA methylation is essential for maintaining cell fate during cell proliferation. Developmental pluripotency-associated 3 (DPPA3) is an intrinsically disordered protein that specifically interacts with UHRF1 and promotes passive DNA demethylation by inhibiting UHRF1 chromatin localization. However, the molecular basis of how DPPA3 interacts with and inhibits UHRF1 remains unclear. We aimed to determine the structure of the mouse UHRF1 plant homeodomain (PHD) complexed with DPPA3 using nuclear magnetic resonance. Induced α-helices in DPPA3 upon binding of UHRF1 PHD contribute to stable complex formation with multifaceted interactions, unlike canonical ligand proteins of the PHD domain. Mutations in the binding interface and unfolding of the DPPA3 helical structure inhibited binding to UHRF1 and its chromatin localization. Our results provide structural insights into the mechanism and specificity underlying the inhibition of UHRF1 by DPPA3.

SeviL is a recently isolated lectin found to bind to the linear saccharides of the ganglioside GM1b (Neu5Ac\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document}α(2-3)Gal\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta$$\end{document}β(1-3)GalNAc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta$$\end{document}β(1-4)Gal\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta$$\end{document}β(1-4)Glc) and its precursor, asialo-GM1 (Gal\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta$$\end{document}β(1-3)GalNAc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta$$\end{document}β(1-4)Gal\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta$$\end{document}β(1-4)Glc). The crystal structures of recombinant SeviL have been determined in the presence and absence of ligand. The protein belongs to the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta$$\end{document}β-trefoil family, but shows only weak sequence similarity to known structures. SeviL forms a dimer in solution, with one binding site per subunit, close to the subunit interface. Molecular details of glycan recognition by SeviL in solution were analysed by ligand- and protein-based NMR techniques as well as ligand binding assays. SeviL shows no interaction with GM1 due to steric hindrance with the sialic acid branch that is absent from GM1b. This unusual specificity makes SeviL of great interest for the detection and control of certain cancer cells, and cells of the immune system, that display asialo-GM1.

Heme oxygenase (HO) converts heme to carbon monoxide, biliverdin, and free iron, products that are essential in cellular redox signaling and iron recycling. In higher plants, HO is also involved in the biosynthesis of photoreceptor pigment precursors. Despite many common enzymatic reactions, the amino acid sequence identity between plant-type and other HOs is exceptionally low (∼19.5%), and amino acids that are catalytically important in mammalian HO are not conserved in plant-type HOs. Structural characterization of plant-type HO is limited to spectroscopic characterization by electron spin resonance, and it remains unclear how the structure of plant-type HO differs from that of other HOs. Here, we have solved the crystal structure of Glycine max (soybean) HO-1 (GmHO-1) at a resolution of 1.06 Å and carried out the isothermal titration calorimetry measurements and NMR spectroscopic studies of its interaction with ferredoxin, the plant-specific electron donor. The high-resolution X-ray structure of GmHO-1 reveals several novel structural components: an additional irregularly structured region, a new water tunnel from the active site to the surface, and a hydrogen-bonding network unique to plant-type HOs. Structurally important features in other HOs, such as His ligation to the bound heme, are conserved in GmHO-1. Based on combined data from X-ray crystallography, isothermal titration calorimetry, and NMR measurements, we propose the evolutionary fine-tuning of plant-type HOs for ferredoxin dependency in order to allow adaptation to dynamic pH changes on the stroma side of the thylakoid membrane in chloroplast without losing enzymatic activity under conditions of fluctuating light.

RBM10, is an RNA binding protein that is important for development by regulating the expression of multiple genes. RBM10 is on the X chromosome, and nonsense and frameshift RBM10 variants cause TARP syndrome in males. In a 4-year-old male, we identified a novel maternally inherited missense RBM10 variant in the RRM2 RNA binding domain, c.965C>T, p.Pro322Leu. His clinical features included intellectual disability, developmental delay, growth restriction, hypotonia, and craniofacial malformations. These features were much milder than those described in previously reported cases of TARP syndrome. By in vitro assays, we found that the mutant p.Pro322Leu RBM10 protein retained its specific RNA binding capacity, while gaining a low-affinity nonspecific RNA binding. It was normally localized to the nucleus, but its expression level was significantly reduced with a significantly short half-life. These results indicated that the p.Pro322Leu missense variant causes a developmental disorder in humans through a unique loss-of-function mechanism.

A convenient analytical system for protein-ligand interactions under crude conditions was developed using native mass spectrometry (MS). As a model protein, Escherichia coli (E. coli) dihydrofolate reductase (DHFR) with and without a histidine tag was used for the study. First, overexpressed DHFR with a His-tag was roughly purified with a Ni-sepharose resin and subjected to native mass spectrometry with or without incubation with an inhibitor, Methotrexate (MTX). Even only with the minimum cleanup by the Ni-sepharose resin, intact ions of DHFR-nicotinamide adenine dinucleotide phosphate (NADPH) and DHFR-NADPH-ligand complexes were successfully observed. By optimizing the preparation procedures of the crude sample for native MS, e.g., avoiding sonication for cell lysis, we successfully observed intact ions of the specific DHFR-NADPH-MTX ternary complex starting with cultivation of E. coli in ≤ 25 mL medium. When the crude DHFR sample was mixed with two, four, or eight candidate compounds, only ions of the specific protein-ligand complex were observed. This indicates that the present system can be used as a rapid and convenient method for the rough determination of binding of specific ligands to the target protein without the time-consuming purification of protein samples. Moreover, it is important to rapidly determine specific interactions with target proteins under conditions similar to those in "real" biological systems. Graphical abstract.

Summary Histone lysine acylations play an important role in the regulation of gene transcription in chromatin. Unlike histone acetyl-lysine, molecular recognition of a recently identified crotonyl-lysine mark is much less understood. Here, we report that the YEATS domain of AF9 preferentially binds crotonyl-lysine over acetyl-lysine in histone H3. Nuclear magnetic resonance structural analysis reveals that crotonyl-lysine of histone H3 lysine 18 is engulfed deep in an aromatic cage of the YEATS domain where the carbonyl oxygen of crotonyl-lysine forms a hydrogen bond with the backbone amide of protein residue Tyr78. The crotonyl-lysine, through its unique electron-rich double-bond side chain, engages π−π aromatic stacking and extended hydrophobic/aromatic interactions with the YEATS domain compared with acetyl-lysine. Our mutational analysis confirmed key protein residues Phe59 and Tyr78 for crotonyl-lysine recognition. Importantly, our findings present a new structural mechanism of protein-protein interactions mediated by histone lysine crotonylation, and show how the cells interpret acyl-lysine marks in different biological contexts.

General transcription factor II E (TFIIE) contains an acid-rich region (residues 378–393) in its α -subunit, comprising 13 acidic and two hydrophobic (Phe387 and Val390) residues. Upon binding to the p62 subunit of TFIIH, the acidic region adopts an extended string-like structure on the basic groove of the pleckstrin homology domain (PHD) of p62, and inserts Phe387 and Val390 into two shallow pockets in the groove. Here, we have examined the dynamics of this interaction by NMR and molecular dynamics (MD) simulations. Although alanine substitution of Phe387 and/or Val390 greatly reduced binding to PHD, the binding mode of the mutants was similar to that of the wild-type, as judged by the chemical-shift changes of the PHD. NMR relaxation dispersion profiles of the interaction exhibited large amplitudes for residues in the C-terminal half-string in the acidic region (Phe387, Glu388, Val390, Ala391, and Asp392), indicating a two-site binding mode: one corresponding to the final complex structure, and one to an off-pathway minor complex. To probe the off-pathway complex structure, an atomically detailed free-energy landscape of the binding mode was computed by all-atom multicanonical MD. The most thermodynamically stable cluster corresponded to the final complex structure. One of the next stable clusters was the off-pathway structure cluster, showing the reversed orientation of the C-terminal half-string on the PHD groove, as compared with the final structure. MD calculations elucidated that the C-terminal half-acidic-string forms encounter complexes mainly around the positive groove region with nearly two different orientations of the string, parallel and antiparallel to the final structure. Interestingly, the most encountered complexes exhibit a parallel-like orientation, suggesting that the string has a tendency to bind around the groove in the proper orientation with the aid of Phe387 and/or Val390 to proceed smoothly to the final complex structure.

Protein-ligand interactions have been commonly studied through static structures of the protein-ligand complex. Recently, however, there has been increasing interest in investigating the dynamics of protein-ligand interactions both for fundamental understanding of the underlying mechanisms and for drug development. NMR is a versatile and powerful tool, especially because it provides site-specific quantitative information. NMR has widely been used to determine the dissociation constant (KD), in particular, for relatively weak interactions. The simplest NMR method is a chemical-shift titration experiment, in which the chemical-shift changes of a protein in response to ligand titration are measured. There are other quantitative NMR methods, but they mostly apply only to interactions in the fast-exchange regime. These methods derive the dissociation constant from population-averaged NMR quantities of the free and bound states of a protein or ligand. In contrast, the recent advent of new relaxation-based experiments, including R2 relaxation dispersion and ZZ-exchange, has enabled us to obtain kinetic information on protein-ligand interactions in the intermediate- and slow-exchange regimes. Based on R2 dispersion or ZZ-exchange, methods that can determine the association rate, kon, dissociation rate, koff, and KD have been developed. In these approaches, R2 dispersion or ZZ-exchange curves are measured for multiple samples with different protein and/or ligand concentration ratios, and the relaxation data are fitted to theoretical kinetic models. It is critical to choose an appropriate kinetic model, such as the two- or three-state exchange model, to derive the correct kinetic information. The R2 dispersion and ZZ-exchange methods are suitable for the analysis of protein-ligand interactions with a micromolar or sub-micromolar dissociation constant but not for very weak interactions, which are typical in very fast exchange. This contrasts with the NMR methods that are used to analyze population-averaged NMR quantities. Essentially, to apply NMR successfully, both the type of experiment and equation to fit the data must be carefully and specifically chosen for the protein-ligand interaction under analysis. In this review, we first explain the exchange regimes and kinetic models of protein-ligand interactions, and then describe the NMR methods that quantitatively analyze these specific interactions.

The Bromodomains and Extra-Terminal domain (BET) proteins direct gene transcription in chromatin, and represent new drug targets for cancer and inflammation. Here we report that the ET domain of the BET protein BRD4 recognizes an amphipathic protein sequence motif through establishing a two-strand antiparallel β-sheet anchored on a hydrophobic cleft of the three-helix bundle. This structural mechanism likely explains BRD4 interactions with numerous cellular and viral proteins such as Kaposi’s sarcoma-associated herpesvirus (KSHV) latency–associated nuclear antigen (LANA), and NSD3 whose interaction with BRD4 via this ET domain mechanism is essential for acute myeloid leukemia (AML) maintenance.

Diffusion ordered spectroscopy (DOSY) is used to determine the translational diffusion coefficients of molecules in solution. However, DOSY is highly susceptible to spurious spectral peaks resulting from thermal convection occurring in the NMR tube. Thermal convection therefore must be suppressed for accurate estimation of translational diffusion coefficients. In this study, we developed a new method to effectively suppress thermal convection using glass capillaries. A total of 6 to 18 capillaries (0.8-mm outer diameter) were inserted into a regular 5-mm NMR tube. The capillaries had minimal effect on magnetic field homogeneity and enabled us to obtain clean DOSY spectra of a mixture of small organic compounds. Moreover, the capillaries did not affect chemical shifts or signal intensities in two-dimensional heteronuclear single quantum coherence spectra. Capillaries are a simple and inexpensive means of suppressing thermal convection and thus can be used in a wide variety of DOSY experiments. Copyright © 2016 John Wiley & Sons, Ltd.

Plugging into the pump Photosynthetic organisms use light to fix carbon dioxide in a process that requires both chemical reducing equivalents and adenosine triphosphate (ATP). Balancing the ratio of these inputs is accomplished by a short circuit in electron flow through photosynthetic complex I, a proton pump that contributes to ATP production but does not increase net reducing equivalents in the cell. Schuller et al. solved a cryo–electron microscopy structure of photosynthetic complex I (see the Perspective by Brandt) and went on to reconstitute electron transfer using the electron carrier protein ferredoxin. Science , this issue p. 257 ; see also p. 230

Nuclear magnetic resonance relaxation dispersion (rd) experiments provide kinetics and thermodynamics information of molecules undergoing conformational exchange. Rd experiments often use a Carr-Purcell-Meiboom-Gill (CPMG) pulse train equally separated by a spin-state selective inversion element (U-element). Even with measurement parameters carefully set, however, parts of 1H&ndash, 15N correlations sometimes exhibit large artifacts that may hamper the subsequent analyses. We analyzed such artifacts with a combination of NMR measurements and simulation. We found that particularly the lowest CPMG frequency (&nu, cpmg) can also introduce large artifacts into amide 1H&ndash, 15N and aromatic 1H&ndash, 13C correlations whose 15N/13C resonances are very close to the carrier frequencies. The simulation showed that the off-resonance effects and miscalibration of the CPMG &pi, pulses generate artifact maxima at resonance offsets of even and odd multiples of &nu, cpmg, respectively. We demonstrate that a method once introduced into the rd experiments for molecules having residual dipolar coupling significantly reduces artifacts. In the method the 15N/13C &pi, pulse phase in the U-element is chosen between x and y. We show that the correctly adjusted sequence is tolerant to miscalibration of the CPMG &pi, pulse power as large as &plusmn, 10% for most amide 15N and aromatic 13C resonances of proteins.

Protein dynamics plays important roles in many biological events, such as ligand binding and enzyme reactions. NMR is mostly used for investigating such protein dynamics in a site-specific manner. Recently, NMR has been actively applied to large proteins and intrinsically disordered proteins, which are attractive research targets. However, signal overlap, which is often observed for such proteins, hampers accurate analysis of NMR data. In this study, we have developed a new methodology called relaxation dispersion difference that can extract conformational exchange parameters from overlapped NMR signals measured using relaxation dispersion spectroscopy. In relaxation dispersion measurements, the signal intensities of fluctuating residues vary according to the Carr-Purcell-Meiboon-Gill pulsing interval, whereas those of non-fluctuating residues are constant. Therefore, subtraction of each relaxation dispersion spectrum from that with the highest signal intensities, measured at the shortest pulsing interval, leaves only the signals of the fluctuating residues. This is the principle of the relaxation dispersion difference method. This new method enabled us to extract exchange parameters from overlapped signals of heme oxygenase-1, which is a relatively large protein. The results indicate that the structural flexibility of a kink in the heme-binding site is important for efficient heme binding. Relaxation dispersion difference requires neither selectively labeled samples nor modification of pulse programs; thus it will have wide applications in protein dynamics analysis.

In the folding of β-lactoglobulin (βLG), a predominantly β-sheet protein, a transient intermediate possessing an excess amount of non-native α-helix is formed within a few milliseconds. To characterize the early folding dynamics of βLG in terms of secondary structure content and compactness, we performed submillisecond-resolved circular dichroism (CD) and small-angle X-ray scattering (SAXS) measurements. Time-resolved CD after rapid dilution of urea showed non-native α-helix formation within 200μs. Time-resolved SAXS showed that the radius of gyration (R(g)) of the intermediate at 300 μs was 23.3±0.7 Å, indicating a considerable collapse from the unfolded state having R(g) of 35.1±7.1 Å. Further compaction to R(g) of 21.2±0.3 Å occurred with a time constant of 28±11 ms. Pair distribution functions showed that the intermediate at 300 μs comprises a single collapsed domain with a small fluctuating domain, which becomes more compact after the second collapse. Kinetic measurements in the presence of 2,2,2-trifluoroethanol showed that the intermediate at several milliseconds possessed an increased amount of α-helix but similar R(g) of 23.0±0.8 Å, suggesting similarity of the shape of the intermediate in different solvents. Consequently, the initial collapse occurs globally to a compact state with a small fluctuating domain irrespective of the non-native α-helical contents. The second collapse of the fluctuating domain occurs in accordance with the reported stabilization of the non-native helix around strand A. The non-native helix around strand A might facilitate the formation of long-range contacts required for the folding of βLG.

Upon HIV-1 infection, a reservoir of latently infected resting T cells prevents the eradication of the virus from patients. To achieve complete depletion, the existing virus-suppressing antiretroviral therapy must be combined with drugs that reactivate the dormant viruses. We previously described a novel chemical scaffold compound, MMQO (8-methoxy-6-methylquinolin-4-ol), that is able to reactivate viral transcription in several models of HIV latency, including J-Lat cells, through an unknown mechanism. MMQO potentiates the activity of known latency-reversing agents (LRAs) or “shock” drugs, such as protein kinase C (PKC) agonists or histone deacetylase (HDAC) inhibitors. Here, we demonstrate that MMQO activates HIV-1 independently of the Tat transactivator. Gene expression microarrays in Jurkat cells indicated that MMQO treatment results in robust immunosuppression, diminishes expression of c-Myc, and causes the dysregulation of acetylation-sensitive genes. These hallmarks indicated that MMQO mimics acetylated lysines of core histones and might function as a bromodomain and extraterminal domain protein family inhibitor (BETi). MMQO functionally mimics the effects of JQ1, a well-known BETi. We confirmed that MMQO interacts with the BET family protein BRD4. Utilizing MMQO and JQ1, we demonstrate how the inhibition of BRD4 targets a subset of latently integrated barcoded proviruses distinct from those targeted by HDAC inhibitors or PKC pathway agonists. Thus, the quinoline-based compound MMQO represents a new class of BET bromodomain inhibitors that, due to its minimalistic structure, holds promise for further optimization for increased affinity and specificity for distinct bromodomain family members and could potentially be of use against a variety of diseases, including HIV infection., This work was supported by funding from the Spanish Ministry of Economy and Competitiveness (MINECO); the European Regional Development Fund (grant BFU2014-52237-P); Fundación para la Investigación y Prevención del SIDA en España (FIPSE 360946/10; to A.J.); ERC Synergy Grant 609989; Centro de Excelencia Severo Ochoa 2013-2017, SEV-2012-0208; and Plan Nacional BFU2012-37168 (to G.J.F.). The work was also supported in part by the research fund of the First Hospital of Jilin University (Changchun, China), grants from the National Institutes of Health (to M.-M.Z.) and the European Research Council under the European Union’s Seventh Framework Programme (FP/2007-2013)/ERC STG 337116 Trxn-PURGE, and Dutch AIDS Fonds grants 2014021 and an ErasmusMC research grant (to T.M.).

The roles of non-native α-helices frequently observed in the initial folding stage of β-sheet proteins have been examined for many years. We herein investigated the residue-level structures of several mutants of bovine β-lactoglobulin (βLG) in quenched-flow pH-pulse labeling experiments. βLG assumes a collapsed intermediate with a non-native α-helical structure (I0) in the early stage of folding, although its native form is predominantly composed of β-structures. The protection profile in I0 of pseudo-wild type (WT*) βLG was found to deviate from the pattern of the “average area buried upon folding” (AABUF). In particular, the level of protection at the region of strand A, at which non-native α-helices form in the I0 state, was significantly low compared to AABUF. G17E, the mutant with an increased helical propensity, showed a similar protection pattern. In contrast, the protection pattern for I0 of E44L, the mutant with an increased β-sheet propensity, was distinct from that of WT* and resembled the AABU...

Jumonji domain-containing protein 6 (JMJD6) is a member of the Jumonji C family of Fe(II) and 2-oxoglutarate (2OG) dependent oxygenases. It possesses unique bi-functional oxygenase activities, acting as both an arginine demethylase and a lysyl-hydroxylase. JMJD6 has been reported to be over-expressed in oral, breast, lung, and colon cancers and plays important roles in regulation of transcription through interactions with transcription regulator BRD4, histones, U2AF65, Luc7L3, and SRSF11. Here, we report a structural mechanism revealed by NMR of JMJD6 recognition by the extraterminal (ET) domain of BRD4 in that a JMJD6 peptide (Lys84-Asn96) adapts an α-helix when bound to the ET domain. This intermolecular recognition is established through JMJD6 interactions with the conserved hydrophobic core of the ET domain, and reinforced by electrostatic interactions of JMJD6 with residues in the inter-helical α1-α2 loop of the ET domain. Notably, this mode of ligand recognition is different from that of ET domain recognition of NSD3, LANA of herpesvirus, and integrase of MLV, which involves formation of an intermolecular amphipathic two- or three- strand antiparallel β sheet. Furthermore, we demonstrate that the association between the BRD4 ET domain and JMJD6 likely requires a protein conformational change induced by single-stranded RNA binding.

Chemical shift perturbations (CSPs) in NMR spectra provide useful information about the interaction of a protein with its ligands. However, in a multiple-ligand-binding system, determining quantitative parameters such as a dissociation constant (K(d) ) is difficult. Here, we used a method we named CS-PCA, a principal component analysis (PCA) of chemical shift (CS) data, to analyze the interaction between bovine β-lactoglobulin (βLG) and 1-anilinonaphthalene-8-sulfonate (ANS), which is a multiple-ligand-binding system. The CSP on the binding of ANS involved contributions from two distinct binding sites. PCA of the titration data successfully separated the CSP pattern into contributions from each site. Docking simulations based on the separated CSP patterns provided the structures of βLG-ANS complexes for each binding site. In addition, we determined the K(d) values as 3.42 × 10⁻⁴ M² and 2.51 × 10⁻³ M for Sites 1 and 2, respectively. In contrast, it was difficult to obtain reliable K(d) values for respective sites from the isothermal titration calorimetry experiments. Two ANS molecules were found to bind at Site 1 simultaneously, suggesting that the binding occurs cooperatively with a partial unfolding of the βLG structure. On the other hand, the binding of ANS to Site 2 was a simple attachment without a significant conformational change. From the present results, CS-PCA was confirmed to provide not only the positions and the K(d) values of binding sites but also information about the binding mechanism. Thus, it is anticipated to be a general method to investigate protein-ligand interactions.

Amyloid fibrils are proteinous aggregates associated with various diseases, including Alzheimer's disease, type II diabetes, and dialysis-related amyloidosis. It is generally thought that, during the progression of these diseases, a precursor peptide or protein assumes a partially denatured structure, which interacts with the fibril seed to change into the final amyloid form. β2-Microglobulin (β2m), associated with dialysis-related amyloidosis, is known to form amyloid fibrils at low pH via a partially structured state. However, the molecular mechanism by which the conformation of β2m changes from the precursor to the final fibril structure is poorly understood. We performed various NMR experiments to characterize acid-denatured β2m. The analysis of the transverse relaxation rates revealed that acid-denatured β2m undergoes a structural exchange with an extensively unfolded form. The results of transferred cross-saturation experiments indicated that residues with a residual structure in the acid-denatured state are associated with the interaction with the fibril seed. Our experimental data suggest the partially structured state to be "activated" to become extensively unfolded, in which state the hydrophobic residues are exposed and associate with the seed. Our results provide general information about the extension of amyloid fibrils.

β-lactoglobulin (LG) contains nine β-strands (strands A–I) and one α-helix. Strands A–H form a β-barrel. At neutral pH, equine LG (ELG) is monomeric, whereas bovine LG (BLG) is dimeric, and the I-strands of its two subunits form an intermolecular β-sheet. We previously constructed a chimeric ELG in which the sequence of the I-strand was replaced with that of BLG. This chimera did not dimerize. For this study, we constructed the new chimera we call Gyuba (which means cow and horse in Japanese). The amino acid sequence of Gyuba includes the sequences of the BLG secondary structures and those of the ELG loops. The crystal structure of Gyuba is very similar to that of BLG and indicates that Gyuba dimerizes via the intermolecular β-sheet formed by the two I-strands. Thus, the entire arrangement of the secondary structural elements is important for LG dimer formation.

Folding experiments have suggested that some proteins have kinetic intermediates with a non-native structure. A simple G ̅o model does not explain such non-native intermediates. Therefore, the folding energy landscape of proteins with non-native intermediates should have characteristic properties. To identify such properties, we investigated the folding of bovine β-lactoglobulin (βLG). This protein has an intermediate with a non-native α-helical structure, although its native form is predominantly composed of β-structure. In this study, we prepared mutants whose α-helical and β-sheet propensities are modified and observed their folding using a stopped-flow circular dichroism apparatus. One interesting finding was that E44L, whose β-sheet propensity was increased, showed a folding intermediate with an amount of β-structure similar to that of the wild type, though its folding took longer. Thus, the intermediate seems to be a trapped intermediate. The high α-helical propensity of the wild-type sequence likely causes the folding pathway to circumvent such time-consuming intermediates. We propose that the role of the non-native intermediate is to control the pathway at the beginning of the folding reaction.

Structural changes of barnase during folding were investigated using time-resolved small-angle X-ray scattering (SAXS). The folding of barnase involves a burst-phase intermediate, sometimes designated as the denatured state under physiological conditions, D phys , and a second hidden intermediate. Equilibrium SAXS measurements showed that the radius of gyration ( R g ) of the guanidine unfolded state (U) is 26.9 ± 0.7 A, which remains largely constant over a wide denaturant concentration range. Time-resolved SAXS measurements showed that the R g value extrapolated from kinetic R g data to time zero, R g,0 , is 24.3 ± 0.1 A, which is smaller than that of U but which is expanded from that of folding intermediates of other proteins with similar chain lengths (19 A). After the burst-phase change, a single-exponential reduction in R g 2 was observed, which corresponds to the formation of the native state for the major component containing the native trans proline isomer. We estimated R g of the minor component of D phys containing the non-native cis proline isomer (D phys,cis ) to be 25.7 ± 0.6 A. Moreover, R g of the major component of D phys containing the native proline isomer (D phys,tra ) was estimated as 23.9 ± 0.2 A based on R g,0 . Consequently, both components of the burst-phase intermediate of barnase (D phys,tra and D phys,cis ) are still largely expanded. It was inferred that D phys possesses the N-terminal helix and the center of the β-sheet formed independently and that the formation of the remainder of the protein occurs in the slower phase.

Bovine beta-lactoglobulin (beta LG) has been one of the most extensively studied proteins in the history of protein science mainly because its abundance in cow's milk makes it readily available to researchers. However, compared to other textbook proteins, progress in the study of beta LG has been slow because of obstacles such as a low reversibility from denaturation linked with thiol-disulfide exchange or monomer-dimer equilibrium preventing a detailed NMR analysis. Recently, the expression of various types of recombinant beta LGs combined with heteronuclear NMR analysis has significantly improved understanding of the physico-chemical properties of beta LG. In this review, we address several topics including pH-dependent structural dynamics, ligand binding, and the complex folding mechanism with non-native intermediates. These unique properties might be brought about by conformational frustration of the beta LG structure, partly attributed to the relatively large molecular size of beta LG. We expect studies with beta LG to continue to reveal various important findings, difficult to obtain with small globular proteins, leading to a more comprehensive understanding of the conformation, dynamics and folding of proteins.

The BET proteins are major transcriptional regulators and have emerged as new drug targets, but their functional distinction has remained elusive. In this study, we report that the BET family members Brd2 and Brd4 exert distinct genomic functions at genes whose transcription they co-regulate during mouse T helper 17 (Th17) cell differentiation. Brd2 is associated with the chromatin insulator CTCF and the cohesin complex to support cis-regulatory enhancer assembly for gene transcriptional activation. In this context, Brd2 binds the transcription factor Stat3 in an acetylation-sensitive manner and facilitates Stat3 recruitment to active enhancers occupied with transcription factors Irf4 and Batf. In parallel, Brd4 temporally controls RNA polymerase II (Pol II) processivity during transcription elongation through cyclin T1 and Cdk9 recruitment and Pol II Ser2 phosphorylation. Collectively, our study uncovers both separate and interdependent Brd2 and Brd4 functions in potentiating the genetic program required for Th17 cell development and adaptive immunity.

Bovine beta-lactoglobulin (betaLG) binds a variety of hydrophobic ligands, though precisely how is not clear. To understand the structural basis of this promiscuous binding, we studied the interaction of betaLG with palmitic acid (PA) using heteronuclear NMR spectroscopy. The titration was monitored using tryptophan fluorescence and a HSQC spectrum confirmed a 1:1 stoichiometry for the PA-betaLG complex. Upon the binding of PA, signal disappearances and large changes in chemical shifts were observed for the residues located at the entrance and bottom of the cavity, respectively. This observation indicates that the lower region makes a rigid connection with PA whereas the entrance is more flexible. The result is in contrast to the binding of PA to intestinal fatty acid-binding protein, another member of the calycin superfamily, in which structural consolidation occurs upon ligand binding. On the other hand, the ability of betaLG to accommodate various hydrophobic ligands resembles that of GroEL, in which a large hydrophobic cavity and flexible binding site confer the ability to bind various hydrophobic substrates. Considering these observations, it is suggested that, in addition to the presence of the hydrophobic cavity, the plasticity of the entrance region makes possible the binding of hydrophobic ligands of various shapes. Thus, in contrast to the specific binding seen for many enzymes, betaLG provides an example of binding with low specificity but high affinity, which may play an important role in protein-ligand and protein-protein networks.

Diffusion ordered spectroscopy (DOSY) is used to determine the translational diffusion coefficients of molecules in solution. However, DOSY is highly susceptible to spurious spectral peaks resulting from thermal convection occurring in the NMR tube. Thermal convection therefore must be suppressed for accurate estimation of translational diffusion coefficients. In this study, we developed a new method to effectively suppress thermal convection using glass capillaries. A total of 6 to 18 capillaries (0.8-mm outer diameter) were inserted into a regular 5-mm NMR tube. The capillaries had minimal effect on magnetic field homogeneity and enabled us to obtain clean DOSY spectra of a mixture of small organic compounds. Moreover, the capillaries did not affect chemical shifts or signal intensities in two-dimensional heteronuclear single quantum coherence spectra. Capillaries are a simple and inexpensive means of suppressing thermal convection and thus can be used in a wide variety of DOSY experiments. Copyright © 2016 John WileySons, Ltd.

NBR1 (neighbor of BRCA1 gene 1) is a protein commonly found in ubiquitin-positive inclusions in neurodegenerative diseases. Due to its high architectural similarity to the well studied autophagy receptor protein p62/SQSTM1, NBR1 has been thought to analogously bind to ubiquitin-marked autophagic substrates via its C-terminal ubiquitin-associated (UBA) domain and deliver them to autophagosomes for degradation. Unexpectedly, we find that NBR1 differs from p62 in its UBA structure and accordingly in its interaction with ubiquitin. Structural differences are observed on helix α-3, which is tilted farther from helix α-2 and extended by approximately one turn in NBR1. This results not only in inhibition of a p62-type self-dimerization of NBR1 UBA but also in a significantly higher affinity for monoubiquitin as compared with p62 UBA. Importantly, the NBR1 UBA-ubiquitin complex structure shows that the negative charge of the side chain in front of the conserved MGF motif in the UBA plays an integral role in the recognition of ubiquitin. In addition, NMR and isothermal titration calorimetry experiments show that NBR1 UBA binds to each monomeric unit of polyubiquitin with similar affinity and by the same surface used for binding to monoubiquitin. This indicates that NBR1 lacks polyubiquitin linkage-type specificity, in good agreement with the nonspecific linkages observed in intracellular ubiquitin-positive inclusions. Consequently, our results demonstrate that the structural differences between NBR1 UBA and p62 UBA result in a much higher affinity of NBR1 for ubiquitin, which in turn suggests that NBR1 may form intracellular inclusions with ubiquitylated autophagic substrates more efficiently than p62.

Relaxation dispersion spectroscopy is one of the most widely used techniques for the analysis of protein dynamics. To obtain a detailed understanding of the protein function from the view point of dynamics, it is essential to fit relaxation dispersion data accurately. The grid search method is commonly used for relaxation dispersion curve fits, but it does not always find the global minimum that provides the best-fit parameter set. Also, the fitting quality does not always improve with increase of the grid size although the computational time becomes longer. This is because relaxation dispersion curve fitting suffers from a local minimum problem, which is a general problem in non-linear least squares curve fitting. Therefore, in order to fit relaxation dispersion data rapidly and accurately, we developed a new fitting program called GLOVE that minimizes global and local parameters alternately, and incorporates a Monte-Carlo minimization method that enables fitting parameters to pass through local minima with low computational cost. GLOVE also implements a random search method, which sets up initial parameter values randomly within user-defined ranges. We demonstrate here that the combined use of the three methods can find the global minimum more rapidly and more accurately than grid search alone.

Chemical shift perturbations (CSPs) in NMR spectra provide useful information about the interaction of a protein with its ligands. However, in a multiple-ligand-binding system, determining quantitative parameters such as a dissociation constant (K(d) ) is difficult. Here, we used a method we named CS-PCA, a principal component analysis (PCA) of chemical shift (CS) data, to analyze the interaction between bovine β-lactoglobulin (βLG) and 1-anilinonaphthalene-8-sulfonate (ANS), which is a multiple-ligand-binding system. The CSP on the binding of ANS involved contributions from two distinct binding sites. PCA of the titration data successfully separated the CSP pattern into contributions from each site. Docking simulations based on the separated CSP patterns provided the structures of βLG-ANS complexes for each binding site. In addition, we determined the K(d) values as 3.42 × 10⁻⁴ M² and 2.51 × 10⁻³ M for Sites 1 and 2, respectively. In contrast, it was difficult to obtain reliable K(d) values for respective sites from the isothermal titration calorimetry experiments. Two ANS molecules were found to bind at Site 1 simultaneously, suggesting that the binding occurs cooperatively with a partial unfolding of the βLG structure. On the other hand, the binding of ANS to Site 2 was a simple attachment without a significant conformational change. From the present results, CS-PCA was confirmed to provide not only the positions and the K(d) values of binding sites but also information about the binding mechanism. Thus, it is anticipated to be a general method to investigate protein-ligand interactions.

NMR spectroscopy is uniquely suited to study protein dynamics over a wide range of time scales at atomic resolution. However, existing NMR relaxation methods require highly serial, lengthy data collection, ultimately limiting their application to short-lived samples, such as proteins in living cells. In recent years, the utility of nonuniform sampling (NUS) NMR methodologies has been increasingly recognized, but their application has been rare in relaxation measurements where highly accurate spectral quantification is demanded. Recently, Matsuki et al. developed a new NUS-processing method, SIFT (Spectroscopy by Integration of Frequency and Time domain information), which is highly robust and faithful in reproducing signals. In this work, we demonstrate the gains that are possible with more aggressive use of frequency domain information than was employed previously. This improvement is crucial for SIFT to be used in accelerating relaxation measurements while preserving full analytical accuracy. By taking the KIX domain of mouse CREB-binding protein (CBP) as an example, we demonstrate that this quantitative NUS processing method enables total 10-fold expedition of the R(2) relaxation dispersion measurements. The advanced SIFT processing should be equally useful for other NMR relaxation measurements.

β-lactoglobulin (LG) contains nine β-strands (strands A–I) and one α-helix. Strands A–H form a β-barrel. At neutral pH, equine LG (ELG) is monomeric, whereas bovine LG (BLG) is dimeric, and the I-strands of its two subunits form an intermolecular β-sheet. We previously constructed a chimeric ELG in which the sequence of the I-strand was replaced with that of BLG. This chimera did not dimerize. For this study, we constructed the new chimera we call Gyuba (which means cow and horse in Japanese). The amino acid sequence of Gyuba includes the sequences of the BLG secondary structures and those of the ELG loops. The crystal structure of Gyuba is very similar to that of BLG and indicates that Gyuba dimerizes via the intermolecular β-sheet formed by the two I-strands. Thus, the entire arrangement of the secondary structural elements is important for LG dimer formation.

Amyloid fibril elongation in denatured proteins involves cycles of coupled binding and misfolding. To gain insights into possible kinetic intermediates, we performed hydrogen/deuterium exchange of amide protons during fibril elongation with β(2)-microglobulin (β(2)-m) at pD=2.5, under which conditions β(2)-m is acid denatured. To study the conformational change in monomeric β(2)-m monitored by NMR spectroscopy, we used (15)N-labeled monomers and nonlabeled seeds. Pulse-labeling hydrogen/deuterium exchange with a quenched-flow apparatus indicated that the rate-limiting intermediate at pD=2.5 is not protected from the exchange, even disrupting a hydrophobic cluster present in the acid-denatured β(2)-m. Significant protection was acquired upon transition to the fibrils. In view of the suggestion that the rate-limiting intermediates are bound to the lateral surface of seed fibrils, weak interactions with a largely unfolded conformation might be useful for their dynamic sliding to the growing ends. The results support a new model of fibril elongation with intermediates bound to the lateral surface of seeds.

Amyloid fibrils elongate seed dependently, with preformed fibrils providing a template for propagation of amyloidogenic conformation. Most seeding experiments use relatively few seed fibrils in comparison with monomers, resembling steady-state enzyme kinetics. Pre-steady-state kinetics should also be useful for characterizing the elongation process. With beta(2)-microglobulin (beta(2)-m), a protein responsible for dialysis-related amyloidosis, we measured the pre-steady-state kinetics of fibril elongation at pH 2.5, conditions under which the monomer is largely unfolded. beta(2)-m has Trp residues at positions 60 and 95. We used three single Trp mutants and fluorescence spectroscopy to study structural change upon fibril elongation. To focus on conformational change in monomers, we prepared seeds with a mutant without a Trp residue. At a fixed concentration of monomeric beta(2)-m, the apparent rate of fibril elongation increased with an increase in the concentration of seeds and then saturated, suggesting the accumulation of a rate-limiting intermediate. Importantly, saturation occurred at a seed/monomer ratio of around 10, as expressed by the concentration of the monomer. Because the number of monomers constituting the seed fibrils is much larger than 10, the results suggest that the elongation process is limited by "non-active-site binding." Spectral analysis indicated that, upon this non-active-site binding, both Trp60 and Trp95 are exposed to the solvent, and then only Trp60 is buried upon transition to the fibrils. We propose a new model of fibril elongation in which non-active-site binding plays a major role.