Fig 2: Loss of CYB5R3 enzymatic activity after ufmylation.a Left, crystal structure of human CYB5R3 (PDB 1UMK). FAD and NADH domains are colored yellow and salmon pink, respectively, whereas bound FAD is colored green. Right, close-up view of the domain interface that recognizes FAD. Broken lines indicate possible hydrogen bonds. b Top, Successive HS-AFM images of CYB5R3. Scale, 10 nm; Height, 0–4 nm. Middle, Histograms of the distances between two globular lobes from 194 frames of 1 molecule. Representative images of closed and open conformations with height profiles are shown right. Bottom, Height profiles of simulated AFM images from the crystal structure (PDB ID: 1UMK) (left) and the manually arranged structure (right) with the simulated images. c The closed conformation (left) corresponds to the crystal structure (1UMK) while the open conformation (right) was manually prepared based on the HS-AFM data. CYB5R3 and FAD are shown with surface and space-filling models, respectively. d Color comparison of purified free CYB5R3 with ufmylated CYB5R3. Recombinant CYB5R3 and UFM1-conjugated CYB5R3 were purified in large quantities as described in Fig. S6. e In vitro enzymatic activity of recombinant CYB5R3 (n = 3) and UFM1-conjugated CYB5R3 (n = 3) were measured as described in Methods section. Data are means ± s.e. Statistical analysis was performed by Šidák’s multiple comparison test after one-way ANOVA. f Left panel: CYB5R3-/- UFSP2-/- HEK293T expressing indicated constructs were lysed and then subjected to immunoblot analysis with anti-CYB5R3 antibody. Data shown are representative of three separate experiments. Right panel: the reductase activity of microsomal fractions prepared from the aforementioned cells was measured as described in the Methods section. Data are means ± s.e. Statistical analysis was performed by Šidák’s multiple comparison test after one-way ANOVA. g, h The amounts of saturated, monounsaturated, and polyunsaturated phosphatidylcholine and of phosphatidylethanolamine (g) (n = 3) and of cholesterol (h) (n = 4) in CYB5R3-/- UFSP2-/- HEK293T expressing UFM1, UFL1, UFBP1, and CYB5R3 or CYB5R3K214R. Data are means ± s.e. Source data are provided as a Source Data file.

Fig 3: CYB5R3 is a bona fide substrate for UFM1.a Immunoblot analysis. Wild-type or UFSP2-/- HEK293T expressing UFL1 and UFBP1 were fractionated into microsomal (M) and cytoplasmic (C) fractions followed by immunoblot analysis. Data shown are representative of three separate experiments. b An experimental approach to identify ufmylated proteins. c Domain structure of CYB5R3. Membrane anchor, FAD-binding site, and lysine residue for UFM1 conjugation are indicated. d Immunoblot analysis. The indicated constructs were transfected into parental or UFSP2-/- HEK293T. 48 h after the transfection, the cell lysates were pulled down with Ni-NTA agarose under denaturing conditions followed by immunoblot analysis. Data shown are representative of three separate experiments. e Immunoblot analysis. A series of CYB5R3KR mutants together with MYC-UFM1, MYC-UFL1, and UFBP-MYC were expressed in UFSP2-/- HEK293T. The cell lysates were analyzed as described in d. Data shown are representative of three separate experiments. f Alignment of the region around K214 of CYB5R3. g In vitro conjugation assay. Recombinant UFM1, UBA5, and UFC1 from E. coli and FLAG-UFL1 and UFBP1-MYC from UFC1-knockout HEK293T cells were incubated with a microsomal fraction prepared from UFBP1-deficient HEK293T cells in the presence or absence of ATP for 90 min, and the mixture was subjected to SDS-PAGE followed by immunoblot analysis with CYB5R3 antibody. Data shown are representative of three separate experiments. h Immunoblot analysis. The indicated constructs were transfected into UFSP2-/- HEK293T. 48 h after transfection, cell lysates were subjected to immunoblot analysis. Data shown are representative of three separate experiments. Source data are provided as a Source Data file.

Fig 4: Ufmylated CYB5R3 interacts with UFBP1.a CYB5R3-/- UFSP2-/- HEK293T expressing indicated constructs were lysed and then subjected to immunoprecipitation. b–c In vitro pull-down assay. The amount of each protein bound to GST-UFBP1 was estimated by immunoblot (b). The faster-migrating ufmylated CYB5R3?N26 proteins are most likely degradation products due to long-term storage and freeze–thawing. UFM1 and UFL1 binding to GST-UFBP1 and each mutant was estimated by immunoblot (c). d Top, analysis of the secondary structure of UFBP1 (116–214) using PSIPRED/DISOPRED61,62 and COILS63. Bottom, alignment of UFIM of UFBP1. e Crystal structure of UFM1 complexed with the UFIM of UFBP1 (left) or the UFIM of UBA5 (right, PDB 5HKH). The side chains of the residues involved in the UFM1-UFIM interaction are shown with a stick model. Broken lines indicate possible hydrogen bonds. f Structural comparison between UFBP1 and UBA5 UFIMs complexed with UFM1. The figure was prepared by superimposing the structure of the UFIM moiety of UFBP1-UFM1 onto that of UBA5 UFIM-UFM1. g In vitro pull-down assay. The amount of each protein bound to GST-UFBP1 was estimated by immunoblot. h–i Wild-type (h) and UFBP1-/- HEK293T (i) expressing indicated constructs were lysed and then subjected to immunoprecipitation. j HEK293T expressing indicated constructs were lysed and then subjected to immunoblot. Bar graphs show the relative value of the ratio of endogenous ufmylated CYB5R3 to CYB5R3 was 1 and the relative value of the ratio of ufmylated UFBP1 to UFBP1 as 1 when UFL1 and UFBP1 were overexpressed, respectively. Data are means ± s.e. Statistical analysis was performed by two-sided Welch’s t test. k HEK293T expressing indicated constructs were lysed and then subjected to pull-down assay with Ni-NTA agarose under denaturing conditions, followed by immunoblot. Bar graph shows the relative value of the ratio of MYC-UFM1-conjugated CYB5R3-His-FLAG to CYB5R3-His-FLAG in crude lysates as 1 when UFL1 and UFBP1 were overexpressed. Data are means ± s.e. Statistical analysis was performed by Šidák’s multiple comparison test after one-way ANOVA. Data for immunoblots presented in this Figure are representative of three separate experiments. Source data are provided as a Source Data file.

Fig 5: Ufmylation of CYB5R3 on the ER.a Immunoprecipitation assay. HEK293T expressing UFL1 and UFBP1 were fractionated into microsomal (M) and cytoplasmic (C) fractions. The fractions were denatured with TNE with SDS (final 1%), diluted 10-fold with TNE, and then immunoprecipitated with CYB5R3 antibody followed by immunoblot analyses. Data shown are representative of three separate experiments. Bar graphs show the relative value of the ratio of endogenous ufmylated CYB5R3 to CYB5R3 was 1 when neither overexpression nor knockdown of any E3 components was performed. Data are means ± s.e. Statistical analysis was performed by two-sided Welch’s t test. b Immunoprecipitation assay. UFL1 or UFBP1 in UFSP2-/- HEK293T were knocked down, and the cells were fractionated into microsomal (M) and cytoplasmic (C) fractions. The fractions were analyzed as described in a. Data shown are representative of four separate experiments. Bar graphs show the relative value of the ratio of endogenous ufmylated CYB5R3 to CYB5R3 was 1 when neither overexpression nor knockdown of any E3 components was performed. Data are means ± s.e. Statistical analysis was performed by Šidák’s multiple comparison test after one-way ANOVA. c Immunoblot analysis. The indicated constructs were transfected into CYB5R3-/- UFSP2-/- HEK293T. 48 h after transfection, the cell lysates were fractionated into microsomal (M) and cytoplasmic (C) fractions followed by immunoblot analysis. Data shown are representative of three separate experiments. Asterisk indicates a degradative product of CYB5R3-His-3xFLAG. Since the CYB5R3 in the cytoplasmic fraction (indicated by asterisk) is almost as mobile as ?N26, it is considered to be the N-terminal truncated form. d Immunofluorescence analysis. GFP-tagged wild-type (CYB5R3-GFP) or ER-localization-defective CYB5R3 (CYB5R3?N26-GFP) were transfected into CYB5R3-/- HeLa. 48 h after transfection, the cells were immunostained with anti-PDI antibody. Bars: 10 µm. Data shown are representative of three independent experiments. Source data are provided as a Source Data file.