Skip to main content

Genetic and process engineering strategies for enhanced recombinant N-glycoprotein production in bacteria

A Correction to this article was published on 12 November 2021

This article has been updated

Abstract

Background

The production of N-linked glycoproteins in genetically amenable bacterial hosts offers great potential for reduced cost, faster/simpler bioprocesses, greater customisation, and utility for distributed manufacturing of glycoconjugate vaccines and glycoprotein therapeutics. Efforts to optimize production hosts have included heterologous expression of glycosylation enzymes, metabolic engineering, use of alternative secretion pathways, and attenuation of gene expression. However, a major bottleneck to enhance glycosylation efficiency, which limits the utility of the other improvements, is the impact of target protein sequon accessibility during glycosylation.

Results

Here, we explore a series of genetic and process engineering strategies to increase recombinant N-linked glycosylation, mediated by the Campylobacter-derived PglB oligosaccharyltransferase in Escherichia coli. Strategies include increasing membrane residency time of the target protein by modifying the cleavage site of its secretion signal, and modulating protein folding in the periplasm by use of oxygen limitation or strains with compromised oxidoreductase or disulphide-bond isomerase activity. These approaches achieve up to twofold improvement in glycosylation efficiency. Furthermore, we also demonstrate that supplementation with the chemical oxidant cystine enhances the titre of glycoprotein in an oxidoreductase knockout strain by improving total protein production and cell fitness, while at the same time maintaining higher levels of glycosylation efficiency.

Conclusions

In this study, we demonstrate that improved protein glycosylation in the heterologous host could be achieved by mimicking the coordination between protein translocation, folding and glycosylation observed in native host such as Campylobacter jejuni and mammalian cells. Furthermore, it provides insight into strain engineering and bioprocess strategies, to improve glycoprotein yield and titre, and to avoid physiological burden of unfolded protein stress upon cell growth. The process and genetic strategies identified herein will inform further optimisation and scale-up of heterologous recombinant N-glycoprotein production.

Background

Attachment of carbohydrates to proteins is the most abundant post-translational modification [1,2,3]. Protein glycosylation occurs in all Domains of life, and in eukaryotes more than half of all proteins are predicted to be glycosylated [4, 5]. It generally involves the transfer of glycans onto the amide side chain of asparagine (N-linked) or the hydroxyl group of serine or threonine (O-linked) amino acid residues. This extra layer of molecular information affects a variety of protein features, from structural dynamics such as folding and stability, to involvement in complex cellular physiology such as cell interactions and pathogenicity [6,7,8,9].

Eukaryotes and a small number of bacterial species from the epsilon subdivision of the Proteobacteria share a homologous process by which they perform N-linked protein glycosylation [10]. Both require assembly of a glycan on a lipid anchor, known as a lipid-linked oligosaccharide (LLO), followed by transfer of the synthesised glycan en bloc by an oligosaccharyltransferase (OTase) onto the asparagine residue of the acceptor protein. However, there are some key variations in reaction components and pathway locations (Fig. 1A, B) [11,12,13]. In eukaryotes, N-linked glycosylation is initiated on the cytosolic side of the rough endoplasmic reticulum (ER) membrane (Fig. 1B). There, glycosyltransferases (GTs) assemble a conserved eukaryotic heptasaccharide structure on a polyprenol diphosphate moiety known as dolichol. This LLO is then flipped to the luminal face of the ER and a further seven sugars are added by different GTs in the ER lumen before the glycan is transferred to the target protein by the OTase complex. In prokaryotes, the best characterised N-linked glycosylation pathway is pgl (protein glycosylation) in the bacterium Campylobacter jejuni (Fig. 1A) [11, 14, 15]. In this pathway, a Cj heptasaccharide structure is sequentially assembled by specific GTs onto an undecaprenol diphosphate lipid carrier on the cytosolic face of the inner membrane (IM), flipped onto the periplasmic face of the IM, and transferred to the target protein. A conserved enzyme in eukaryotic and prokaryotic N-linked glycosylation is the OTase, which catalyses covalent attachment of the glycan to the acceptor sequon of the target protein [16, 17]. The consensus acceptor sequon in Archaea and Eukaryota is N-X-S/T (X ≠ Proline), while in bacteria, acidic amino acids at the -2 position are required (D/E-X1-N-X2-S/T, X1, 2 ≠ Proline) [18, 19]. Except for some single-celled protists such as Leishmania major and Trypanosoma brucei, eukaryotic OTases are hetero-oligomeric, comprised of seven to nine proteins [10, 20]. The catalytic sub-unit of eukaryotic OTase is STT3 and in mammals two STT3 isoforms (STT3A and B) have distinct roles during N-linked glycosylation. The STT3A-dependent complex associates with the translocation machinery to optimise co-translational glycan transfer to protein, whilst STT3B-dependent glycosylation occurs following translocation [21, 22]. Bacterial N-linked OTases, such as PglB from C. jejuni, consist of a single transmembrane protein homologue of eukaryotic STT3, that can carry out both co- and post-translocational glycosylation [23,24,25,26].

Fig. 1
figure1

Schematic overview of native N-linked protein glycosylation pathway and proposed strategies in this study to improve sequon accessibility of recombinant target protein to PglB during glycosylation in glycoengineered E. coli. N-linked glycosylation in bacteria/C. jejuni (A) compared to early-stage of eukaryotic/mammalian N-linked glycosylation pathway (B). A In C. jejuni, undecaprenol-linked glycan is synthesised by glycosyltransferases in the cytoplasm, flipped to the periplasm by flippase (PglK), and then transferred co or post-translocationally to a target protein by PglB. B In mammalian N-linked glycosylation, dolichol-linked glycan is synthesised both in the cytosol and endoplasmic reticulum (ER) lumen, and glycan is transferred co- or post-translocationally (later not shown) by different STT3 isoforms. C, D Two different strategies are proposed to enhance heterologous protein glycosylation in E. coli. C Approach based on increasing PglB interaction with target protein in the membrane. Increased membrane residence of target protein was achieved by introducing signal peptide mutations with poor processivity or uncleaved signal. D Approach based on increasing sequon accessibility of target proteins to PglB via modulation of target protein folding state during or after translocation. Using a disulphide-bonded model protein, protein folding and maturation in the periplasm are delayed by producing the protein under suboptimal conditions for disulphide formation such as under oxygen-depleted conditions, in the absence of oxidoreductases, or under chemical redox treatment. Solid arrow = increased reaction, dashed arrow = reduced reaction

Benefiting from improved pharmacological quality, about 70% of marketed recombinant human therapeutic proteins contain (mostly N-linked) glycans [27,28,29]. Production in mammalian cells, predominantly in the Chinese Hamster Ovary (CHO) cell, enables the coupling of complex mammalian glycans to glycoproteins at an industrial scale [30, 31]. However, costly manufacturing processes, heterogeneous products, and risk of viral contamination remain a significant challenge for production in mammalian cell lines [32,33,34]. As glycosylation is conserved in all eukaryotic cells, other eukaryotic expression systems such as yeast, plant and insect cells have also been employed to perform this modification [35, 36]. These systems offer the advantage of relatively faster production and lower risk of mammalian virus introduction, if non-mammalian or non-humanised glycans are acceptable within the glycoprotein product [37,38,39]. However, the production of protein containing a novel glycan or purely heterologous glycoprotein often requires alteration to native pathways which can potentially affect cell fitness, as many glycosylation pathways are essential for viability in eukaryotic cells [40,41,42].

Functional transfer of the pgl locus from C. jejuni into E. coli, a bacterial workhorse for recombinant protein production, opened up a new area of bacterial glycoengineering [15]. Glyco-competent E. coli carrying heterologous protein glycosylation machinery is an attractive platform owing to its potential as a non-virulent, rapidly growing host, with low fermentation costs [43,44,45]. From an engineering standpoint, the absence of a native protein glycosylation pathway in E. coli facilitates rational design to introduce orthogonal glycosylation pathways, without interference from endogenous GTs, hence leading to the production of homogenous products [45, 46]. In addition, one could expect a minimised effect of glycoengineering on the fitness of E. coli cells, as the host does not depend on protein glycosylation for cell viability. Bacterial glycoengineering also benefits from the use of the key enzyme Cj PglB OTase, which has been known to have relaxed substrate specificities in term of glycans and target proteins [18, 23, 24, 47]. Transfer of glycans with various size and composition has been demonstrated along with a wide range of protein recipients spanning from bacterial to eukaryotic origin. The only caveat is that the glycans require the presence of an acetamido group in the reducing-end sugar [47], and the target protein needs to contain the sequon in a structurally exposed and flexible region [18, 26]. To date, glyco-competent E. coli has been extensively developed in the manufacture of novel recombinant bacterial vaccines and glycoconjugates [46, 48,49,50,51,52]. Further, promising progress has been made to engineer glyco-competent E. coli to produce authentic mammalian glycans and glycoproteins [53,54,55,56], and to explore and engineer bacterial OTase with a greater substrate specificity [57,58,59].

Nevertheless, a common challenge of recombinant N-linked glycoprotein production in E. coli is inefficient glycosylation [26, 44, 45, 60]. To overcome this challenge some strains were improved by eliminating competing pathways either in glycan biosynthesis or glycan destination [61,62,63,64,65]. Other efforts focused on the genomic integration and expression of pglB and other glycosyltransferases, led to ~ twofold enhancement in glycosylation efficiency [65,66,67]. Glycosylation efficiency was also increased by both metabolic pathway engineering (1.2- to 1.8-fold) and process optimisation (1.3- to 1.5-fold) [33, 43, 68, 69]. Overall, glycosylation efficiency is reported to largely depend upon sequon accessibility, and so it is essential to understand how the OTase and target proteins interact during co- and post-translocational glycosylation [23,24,25,26]. In vitro studies have shown that folded protein can be glycosylated more efficiently by PglB (~ 15-fold) if the sequon is located on an exposed, flexible and unstructured region [26]. Sequon accessibility can also be enhanced in vitro by interfering with protein folding or destabilisation [23]. In contrast, coupling of glycosylation and the Sec-translocation pathway in E. coli improved AcrA and PEB3 glycosylation (two to fourfold) [26]. However, glycosylation efficiency of these recombinant targets remained much lower in E. coli compared to C. jejuni (up to 15-fold). Under certain fermentation conditions previously evaluated, increased glycosylation efficiency occurs in parallel with a decrease in protein production, rather than an increase in total amount of glycoprotein per cell (yield) [48]. Further, these conditions can also lead to growth defects, negatively affecting total glycoprotein per culture volume (titre) [48]. So, in order to develop alternative genetic and process strategies it is important to evaluate glycoprotein yield and titre as well as glycosylation efficiency.

Mammalian STT3A interacts directly with Sec61 translocon complex, thus increasing the local concentration of OTase-target protein substrate and allowing rapid recognition of the sequon in relaxed-unfolded regions of the protein during translocation before the protein diffuses away from the membrane [22]. Moreover, recent studies demonstrated the presence of an oxidoreductase-like protein in the eukaryotic OTase complex—Ost3p/Ost6p in yeast or N33/Tusc3 in human STT3A and B containing OTases. These oxidoreductase-like proteins modulate folding of disulphide-containing proteins in the ER lumen prior to glycosylation leading to increased sequon accessibility [70, 71]. However, to date the involvement of folding modulators has not been reported in the bacterial N-glycoprotein pathway.

Here, using E. coli as a host, we explored several strategies to maximise N-linked glycosylation including modulation of expression, sequon folding state and accessibility of the target protein, in order to increase interaction with PglB (Fig. 1C, D). The target proteins selected varied in their structure, disulphide-bond content and sequon position in order to explore the impact of different protein structural contexts upon the strategies employed. Our initial study sought to explore the influence of modulating target production level upon glycosylation efficiency. We next investigated if increasing membrane residency time of the target protein, thereby increasing local interaction between OTase PglB and the target protein could enhance glycosylation (Fig. 1C). To further explore the effect of protein folding modulation on glycosylation, we used model disulphide-bond containing proteins (Fig. 1D). These were produced under sub-optimal conditions for disulphide formation, such as under oxygen-depleted condition or in the absence of oxidoreductase (ΔdsbB or ΔdsbC), in order to mimic the ER lumen of eukaryotes. As production of recombinant protein under sub-optimal conditions can result in reduced cell viability and total protein production, we explored supplementation with the chemical oxidant cystine to recover disulphide bond-containing protein production yields in the oxidoreductase mutant (ΔdsbB). The various genetic and process conditions presented here demonstrated enhanced glycosylation efficiency, yield and titre, informing potential strategies to improve N-linked glycosylation in E. coli for the production of various disulphide containing and non-disulphide target glycoproteins.

Results

Selection of model proteins to assess glycosylation efficiency

The four model proteins used in this study, N-glycosylation reporter protein NGRP, anti-β-galactosidase single-chain Fv scFv13R4 and R4CM, and bovine pancreatic ribonuclease RNase A were selected based on their diversity of protein structure, disulphide-bond content, and sequon location within the protein (Fig. 2A, B). These model proteins have been used previously to study N-linked glycosylation in E. coli [23, 58, 65, 72]. For simplicity all were modified by the addition of a single bacterial N-linked glycosylation sequon to assess glycosylation efficiency. A C-terminal hexahistidine tag was also added for detection of glycosylated and non-glycosylated forms by immunoblotting (anti-His). NGRP is a truncated version of the native periplasmic C. jejuni Cj0114 protein [73] designed to improve the performance of glycopeptide structural characterisation by mass spectrometry [58]. NGRP contains either a sequon within the C-terminal region and lacks disulphide-bonds (Fig. 2A, B, i). The scFv13R4 protein was modified by addition of a C-terminal DQNAT glycosylation sequon (Fig. 2B, ii). This sequon was repositioned closer to the cysteine residue involved in disulphide-bond formation to generate the R4CM variant (Fig. 2A ii, B iii). RNase A was modified to replace S32 with D converting the nearby native eukaryotic sequon (NLT) to the bacterial consensus sequon (DRNLT) [23]. This sequon is also adjacent to a cysteine residue involved in disulphide-bond formation (Fig. 2A iii, B iv). Varying the location of this sequon has been shown to influence glycosylation [23, 72]. scFv13R4 and R4CM are predicted to have a single disulphide-bond [74] whilst RNase A contains four, three of which are non-consecutive.

Fig. 2
figure2

Structural variation of model glycoproteins. A Structural models of (i) NGRP, (ii) scFv13R4 and scFv13R4CM were generated by Phyre2 based on protein homology prediction (≥ 99% confidence) [113]. (iii) X-ray crystallographic structure of RNase A (PDB code 3WMR). Ribbon model of the proteins was drawn by UCSF Chimera [114]. Position of the sequon variants (D/E-X-N-X-S/T, X ≠ P) within the protein is indicated. Disulphide-bonds are highlighted in yellow. C-terminal sequon (DQNAT) of scFv13R4 is not displayed in the protein model. B Linear representation of proteins; (i) NGRP, (ii) scFv13R4, (iii) scFv13R4CM, and (iv) RNase A. Position of disulphide-bonds (C–C) and glycosylation sites (N) are indicated with amino acid positions

Expression attenuation of the NGRP model target has limited impact upon glycosylation efficiency

It was previously shown that NGRP was produced in E. coli as a soluble periplasmic protein [58]. Initially in this study, we aimed to explore the impact of NGRP expression level upon glycosylation efficiency using inducible expression vector pDEST-ORS (Fig. 3A). In this vector, expression of the target gene is regulated both transcriptionally by the IPTG-inducible tac promoter (Ptac) and translationally by the PPDA-inducible orthogonal riboswitch (ORS) located within the 5ʹ UTR [75, 76]. The SecB-dependent signal peptide of PelB was integrated at the N-terminus of the target protein to direct secretion into the periplasm. The 5ʹ coding sequence context of genes has previously been shown to influence riboswitch conformation and gene expression levels [77]. Here, three variants containing different synonymous codons at the first eight amino acids within the PelB leader sequence were used, termed wild type (WT), PelB 1 and PelB 2-NGRP (Fig. 3A). pDEST-ORS-NGRP was transformed into E. coli Top10F’ expressing the C. jejuni N-glycosylation machinery from pACYCpgl (glyco-competent). Cells were cultivated in triplicate and expression induced over a range of different inducer concentrations (100 μM IPTG and 0, 2, 8, 40, 100, 200 and 400 μM PPDA) (“Methods” section). To analyse glycoprotein production, proteins were collected from periplasmic fractions. Total target protein was quantified by Western blot densitometry in which pre-determined purified protein (between 15 to 100 ng) was used as a standard for quantification (“Methods” section). The quantification data were plotted as the yield mg/g of dry cell weight (DCW) following the normalisation and conversion of measured OD600 into DCW (“Methods” section). Glycosylation efficiency (%) was determined from the proportion of glycosylated to total NGRP (G1/G1 + G0). NGRP expression in glyco-competent E. coli produced a single band indicative of glycosylated NGRP (~ 22 kDa) with slower migration in the gel compared to non-glycosylated NGRP (~ 20 kDa) (Fig. 3B). Previous study has confirmed by mass spectrometry that this altered mobility of NGRP is due to N-linked glycosylation [58].

Fig. 3
figure3

Glycosylation of NGRP in glyco-competent E. coli. A Organisation of pDEST-ORS expression vector used in this study. Target gene (NGRP) was fused with Sec signal peptide PelB (N-terminal) and Hexahistidine-tag (C-terminal). Expression of the target gene was regulated transcriptionally by Ptac via IPTG induction and translationally by orthogonal riboswitch (ORS) via PPDA titration. IPTG = Isopropyl β-d-1-thiogalactopyranoside, PPDA = Pyrimido [4,5-d] pyrimidine-2,4-diamine. Three synonymous nucleotide sequence variants of PelB-NGRP N-terminal codon (PelB 1, PelB 2, and PelB WT or wild-type) are shown. These variants were tested to explore the impact of the 5ʹ codon context upon the riboswitch-dependent regulatory function. B Western blot analysis of periplasmic fractions of glyco-competent (GC) and non-glycocompetent (Non-GC) strains of E. coli. Anti-His antibody was used to detect the NGRP. Arrows indicate non-glycosylated (G0) and glycosylated (G1) NGRP. C, D Quantitative Western blot analysis (Densitometry) of NGRP located in the periplasmic fractions from the two signal peptide variants, C PelB 1 and D PelB 2-NGRP. Proteins were produced with increasing PPDA inducer concentrations (100 μM IPTG and 0, 2, 8, 40, 100, 200, 400 μM PPDA). Total proteins were quantified using a pre-determined purified NGRP standard curve (15 ng to 100 ng) (“Methods”). The data were converted into mg/g of dry cell weight (DCW) based on normalisation and calculation with measured OD600 of the samples. Glycosylated (yellow bar) and non-glycosylated (green bar) protein are as shown (left y-axis). % Glycosylation (% G1/G0 + G1) is indicated (black circle, right y-axis). Data were processed from three biological replicates; error bars represent standard deviation from mean values. Representative Western blots are shown as insets

For all constructs no NGRP production was detected in the absence of inducer confirming tight control of basal expression (Additional file 1: Figure S1A). Production of NGRP containing PelB-WT showed poor titratability upon induction (Additional file 1: Figure S1B–D) but titration of NGRP production was observed with increasing PPDA inducer in PelB 1 and PelB 2-NGRP constructs giving a dynamic range of expression of 6 and 15-fold respectively (expression level at Max/Min inducer concentration) (Fig. 3C, D). Interestingly, NGRP was glycosylated with similar efficiency (PelB 1-NGRP = 45 ± 4%, PelB 2-NGRP = 49 ± 5%) irrespective of the total amount of NGRP produced (PelB 1-NGRP = 3.9 ± 0.9 to 22.7 ± 1.5 mg/g DCW, PelB 2-NGRP = 1.2 ± 0.4 to 17.7 ± 1.5 mg/g DCW) (Fig. 3C, D). Correlation analysis (Pearson) showed no strong or significant relationship between glycosylation efficiency and the protein production level of NGRP (Additional file 1: Figure S2). This appears to indicate that the cell glycosylation capacity i.e., lipid-glycan substrate availability and glycosylation enzyme activity do not limit NGRP glycosylation, at least within the parameters of this experiment. Cell growth was analysed based on the optical density of cell culture at the end of incubation or induction. The result showed that increased NGRP production level has minimal but significant impact on cell growth (Pearson correlation PelB 1-NGRP, r = − 0.7645, P = 0.0453 PelB 2-NGRP, r = − 0.8042, P = 0.0292) (Additional file 1: Figure S3). As the PelB 1 variant of leader sequence gave highest maximum expression and improved dynamic range of gene expression, this variant then was used within the pDEST-ORS construct for the subsequent experiments.

Design of a signal peptide cleavage site mutant to increase membrane residency time of target protein during glycosylation

As the gene titration analysis of NGRP indicated that glycosylation efficiency was independent of the amount of protein produced, we sought to explore if glycosylation efficiency was dependent upon the target protein residency time within the membrane, and therefore whether co-translocational glycosylation could be enhanced. To do this we constructed a system to allow membrane trapping of the target protein through modification of the PelB signal peptide cleavage site. We reasoned that by decreasing the recognition and processivity of the cleavage site by signal peptidase I (SPaseI) (Additional file 1: Figure S4), that the target protein would be trapped in the membrane leading to increased residency time between the OTase and the target protein-Sec complex. In addition to NGRP, we also used scFv13R4 as a target protein, to explore how the strategy would work with structurally distinct proteins.

The SPaseI-dependent cleavage site is located at the C terminal end of the signal peptide after the recognition motif Ala-X-Ala (X any amino acids) (Additional file 1: Figure S4) [78]. Substituting Ala with Thr at -1 from the cutting site has been shown to reduce processing [78, 79]. Another key amino acid is a proline residue around − 4 to − 6, whose structural turn promotes signal peptide alignment with SPaseI, and deletion or replacement is reported to inhibit signal peptide cleavage [79, 80]. In silico analysis (SignalP 4.1) of the resultant modified PelB signal peptide Leu-Thr-Met-Thr (TMT) motif showed a predicted 40–60% decrease in cleavage site processing compared to the wild type Pro-Ala-Met-Ala (wt) motif (Additional file 1: Figure S5A-B and C-D, table for Y-score) [81]. With only residues around the cleavage site modified in the TMT variant, 90% of Sec signal sequence was left unchanged (N, H and C-region) [82], retaining the signal for protein sorting and delivery to the translocon (Additional file 1: Figures S4, S5). PCR site-directed mutagenesis then was performed to generate a TMT mutant of PelB NGRP and scFv13R4 in the pDEST-ORS construct (“Methods” section).

To investigate the effect of the cleavage site mutations upon protein localisation and glycosylation, NGRP and scFv13R4 wt and TMT variants were expressed in glyco-competent (GC) and non glyco-competent (Non-GC, control) E. coli strains (“Methods” section). We hypothesised there would be four possible NGRP or scFv13R4 isoforms generated based on their signal peptide processivity and glycosylation state (Fig. 4A). These predicted forms were processed (a-form), unprocessed (b-form), processed glycosylated (c-form), and unprocessed glycosylated (d-form). Glycosylated and unprocessed protein will migrate slower than non-glycosylated and processed protein due to additional glycan (single, ~ 1.5 kDa) and un-cleaved signal peptide (~ 2 kDa). As expected, the fastest migrating species were observed in the periplasm, (Fig. 4B, C, lanes 1 and 5), and are predicted to be correctly processed (a-form). In contrast, the mutated signal peptide cleavage site clearly impacts upon the release of target protein into the periplasm fraction (Fig. 4B, C, lanes 3 and 7). Further, proteins with the TMT mutation were only observed located in the membrane fraction (Fig. 4B, C, lanes 4 and 8), which are most likely unprocessed (b-form). Glycosylated isoforms were observed as an additional higher molecular weight band in the periplasm or membrane of glyco-competent cell (predicted c or d-form) (Fig. 4B, lanes 5 and 8, and Fig. 4C, lanes 5, 6 and 8). Glycosylation was confirmed using a polyclonal CjNgp antiserum [83] but as it detects both glycosylated and non-glycosylated NGRP (Additional file 1: Figure S6) it could not be used to probe the degree of glycosylation in NGRP. However, the glycosylation status of scFv13R4 could be clearly observed (Fig. 4C). We noticed that the gel migration of some predicted isoforms could not be clearly resolved or with apparent different migration across the different lanes of the gel, (e.g., Fig. 4B, lanes 5, 6 and 8, and Fig. 4C, lanes 5 and 8). To rule out that band separation was due to differences in migration in different gel lanes, we mixed the samples (lane 5 and 8) and ran them together for an extended electrophoresis running time. By this re-running, all 4 isoforms (a–d forms) were detected as expected and the difference in migration rate was evident (Additional file 1: Figure S8). No glycosylated isoforms were observed in the non glyco-competent controls (Fig. 4B, C, lanes 1–4). Although analysis of the membrane fraction from both glyco-competent and noncompetent cells producing scFv13R4 with the wt signal peptide cleavage site produced bands (Fig. 4C, lanes 2 & 6) with similar migration as those in the periplasm fractions (a-form) (Fig. 4C, lanes 1 and 5). However, we suspect this band is correctly processed but insoluble scFv13R4 from the periplasm that has been associated with the membrane during the fractionation process, indeed this observation has been previously reported for this protein [76].

Fig. 4
figure4

Production of NGRP and scFv13R4 isoforms containing signal peptide cleavage site variants. A Schematic representation of four predicted isoforms of NGRP or scFv13R4 (a-d-forms) based on their signal peptide processivity and protein glycosylation. B, C Western blot analysis of membrane (M) and periplasmic (P) expression of wild type (wt) and signal peptide cleavage mutant (TMT) NGRP (B) and scFv13R4 (C) in glyco-competent (GC) and non glyco-competent (Non-GC) E. coli. B, C All NGRP and scFv13R4 isoforms were detected by anti-His antibody. Glycosylated scFv13R4 was detected by CjNgp antibodies. Predicted a-d isoforms within the bands are indicated. The lanes in panel B are from the same blot while white vertical lines indicates a non-adjacent lanes. The protein migration band relative to each other in different lanes in the blot is unchanged. Additional file 1: Figure S7 shows the uncropped Western blot image of B. DG Quantitative Western blot analysis (Densitometry) of membrane (TMT) and periplasmic (wt) localised DE NGRP and FG scFv13R4 produced in glyco-competent E. coli. Proteins were produced under different induction conditions (100 μM IPTG and 0, 8, 40 and 200 μM PPDA). Anti-His antibody was used to detect the proteins (glycosylated and non-glycosylated). Total proteins were quantified using a pre-determined purified NGRP or scFv13R4 standard curve (5 ng to 100 ng) (“Methods”). The data were converted into mg/g of dry cell weight (DCW) based on normalisation and calculation with measured OD600 of the samples. Glycosylated (yellow bar) and non-glycosylated (green bar) protein as shown (left y-axis). % Glycosylation (% G1/G0 + G1) is indicated (black circle, right y-axis). All data (DG) were processed from three biological replicates. Error bars indicate standard deviation from mean values

In agreement with previous findings, the ability of PglB to perform co-translocational glycosylation was supported here by observation of glycosylated pre-protein TMT-NGRP and TMT-scFv13R4 (d-form) in the membrane of glyco-competent cells (Fig. 4B, C, lanes 8) [26]. Whereas wt-NGRP seems to be glycosylated post-translocation or just before the protein release into the periplasm after signal peptide processing, since no glycosylation of the membrane bound wt-NGRP pre-protein (d-form) was detected (Fig. 4B, lane 6). Glycosylation of membrane bound NGRP was only observed for the TMT variant where an additional band in the membrane was shown, suggesting that glycan is mostly acquired during an extended residency time within the membrane (Fig. 4B, lane 8). In subsequent experiments, we performed quantitative analysis of glycosylation efficiency for the target proteins in both the membrane and in the periplasm.

Glycosylation of NGRP and scFv13R4 during extended membrane residency time

Our data demonstrated that target proteins containing the modified cleavage site (TMT) were not released to the periplasm and were retained in the membrane. To evaluate how secretion attenuation contributed to glycosylation, the wt and TMT variants of NGRP and scFv13R4 were produced in glyco-competent E. coli, under different inducer concentrations (100 μM IPTG and 0, 8, 40, 200 μM PPDA) (“Methods” section). Figure 4D–G shows the quantitative Western blot analysis of protein detected in the periplasm and membrane fractions for the wt and TMT variants. We observed that total target protein production increased upon induction, but with only modest dynamic ranges observed for scFv13R4 (wt = 1.6-fold, from 7.5 ± 2.3 to 12.3 ± 3.8 mg/g DCW; TMT = 1.4-fold, from 0.39 ± 0.07 to 0.56 ± 0.17 mg/g DCW) in contrast to NGRP which displayed greater dynamic ranges (wt = 4.3-fold, from 3.3 ± 0.7 to 11.7 ± 0.6 mg/g DCW; TMT = 2.4-fold, from 4.7 ± 0.5 to 10.1 ± 0.9 mg/g DCW). Both membrane (TMT) and periplasmic (wt) fractions of NGRP yielded similar amounts of target protein under the same induction condition (Fig. 4D, E; Additional file 1: Figure S9A). In contrast, the membrane fraction for TMT-scFv13R4 had significantly lower target protein than the periplasmic fraction wt-scFv13R4 (~ 15- to 20-fold reduction) (Fig. 4F, G; Additional file 1: Figure S9B). One possible explanation is that the two exemplar target proteins display different intrinsic stability and proteolytic susceptibility within the membrane, which consequentially reduces the relative protein yield in the membrane compared to that observed in the periplasm fraction. As previously observed (Fig. 3C, D), the proportion of glycosylated NGRP was minimally impacted by the protein production levels, both in the periplasm fraction (wt) and in the membrane fraction (TMT) (average glycosylation efficiency wt = 29 ± 6%, TMT = 31 ± 3%) (Fig. 4D, E). Indeed, correlation analysis (Pearson) showed no significant relationship between protein level and glycosylation efficiency for wt and TMT-NGRP (Additional file 1: Figure S10A and B). Increased NGRP production levels had no negative impact upon cell growth at lower induction condition (0–40 μM PPDA) (OD600 wt ~ 8, TMT ~ 11) (Additional file 1: Figure S11A, B). Only at the highest induction levels (200 μM PPDA), was a reduction in cell growth observed (OD600 wt ~ 5, TMT ~ 3) (Additional file 1: Figure S11A, B). Overall, we did observe cell culture dependent effects, as lower NGRP glycosylation efficiency and protein yields were achieved in this experiment (Fig. 4) compared to the previous result (Fig. 3) (glycosylation efficiency mean 45 ± 4% vs 29 ± 6%, and total yield 23.94 ± 2.21 vs 11.7 ± 0.6 mg/g DCW). The previous experiment was run under microculture cultivation conditions (vs. flask) and shorter induction time (3 h vs 6 h) (see “Methods” section).

In contrast to the limited impact of titration upon glycosylation efficiency for NGRP, increased scFv13R4 production levels led to decreased glycosylation efficiency (wt = 1.6-fold, from 96 to 63%; TMT = 1.4-fold, from 77 to 47%) (Fig. 4F, G). Correlation analysis (Pearson) showed a significant inverse relationship between protein expression level and glycosylation efficiency for wt-scFv13R4 (r = − 0.9920, P = 0.008), while for TMT no significant correlation was observed) (Additional file 1: Figure S10C, D). Cell growth analysis showed that production of both the scFv13R4 variants (wt and TMT) led to reduced cell growth at higher protein production levels (wt OD600 5.4 ± 0.2 to 2.4 ± 0.6; TMT OD600 3.0 ± 0.5 to 2.1 ± 0.2) (Additional file 1: Figure S11C, D). Therefore, in addition to the effect of protein production level, reduced cell fitness could also possibly contribute to the decreasing glycosylation activity in the cell [48].

Interestingly, glycosylation efficiency of the membrane-bound TMT-scFv13R4 was always consistently higher than for the periplasmic wt-scFV13R4 under the same induction conditions (96–63% vs 77–47%, 0.015 ≤ P ≤ 0.08) (Fig. 4, F and G, and Additional file 1: Figure S12A). Whereas the glycosylation efficiency observed for the NGRP variants was similar under the same induction conditions (0–40 μM PPDA, wt 34%–28% vs. TMT 34%–29%). The exception to this observed trend was under the highest induction conditions (200 μM PPDA) which showed lower glycosylation efficiency of wt-NGRP (22%) compared to TMT-NGRP (29%) (Fig. 4D, E, and Additional file 1: Figure S12B). The lower production level of TMT-scFv13R4, when compared to the wt-scFv13R4, is a potential contributing factor to the observed increase glycosylation efficiency. However, analysis of the relationship between protein yield and glycosylation efficiency for wt-scFv13R4 by linear regression (Additional file 1: Figure S13), indicates that at lower protein production levels (≤ 4.15 mg/g DCW) the target protein would be predicted to be fully glycosylated. Therefore, it would be expected that TMT-scFv13R4 would be fully glycosylated at the yields observed in the membrane (from 0.39 ± 0.07 to 0.56 ± 0.17 mg/g DCW). Therefore, protein production levels alone are likely not the only factor that can affect TMT-scFv13R4 glycosylation, and other factors, such as sequon accessibility during membrane retention could limit glycosylation. From this analysis it is clear that protein glycosylation is a combination of glycosylation both during protein translocation through the membrane (co-translocation) and after translocation in the periplasm (post-translocation). So, to achieve maximal glycosylation efficiency optimal conditions for both processes are required. Therefore, the extended membrane residency time of the TMT variant compared to the wt, appears to be a likely contributing factor upon the observed increase in glycosylation.

Increased glycosylation efficiency of membrane-bound scFv13R4-TMT indicated that extended membrane residency time may enhance protein glycosylation compared to the when protein is released into the periplasm (scFv13R4-wt). However, since glycosylation efficiency of NGRP was not affected by enhanced membrane retention, we suspected that the structural context of the sequon within NGRP and scFv13R4 influences their respective abilities to accept glycan from PglB. Unlike NGRP, scFv13R4 contains disulphide-bonds, formation of which are important for correct protein folding and which might be impeded by membrane retention [84, 85]. Indeed, the different effects of membrane retention upon protein folding and stability of NGRP and scFv13R4 might also influence the different protein yields observed between periplasmic (wt) and membrane production (TMT) of both proteins (Fig. 4D–G). Additionally, having more complex secondary and tertiary structure than NGRP, scFv13R4 might require longer time to fold which consequently may extend the unfolded glycosylation-competent state in the membrane [85]. If formation of structural elements (e.g., disulphide-bonds) in the target protein could be impeded to modulate protein folding and glycosylation, we hypothesised that perhaps specific conditions that impede structure formation could also be exploited to enhance protein glycosylation. Thus, we explored the impact of varying both genetic and process (growth) conditions that might lead to attenuation of disulphide-bond formation, to improve production of the target recombinant glycoprotein.

Expressing disulphide-bond containing proteins under oxygen-depleted condition enhances glycosylation efficiency

Intramolecular disulphide-bonds are essential to maintain the folding and stability of many proteins [86]. In E. coli and many Gram-negative bacteria, their formation is catalysed enzymatically by the dedicated disulphide-bond forming (DSB) pathway, primarily under aerobic conditions [87,88,89,90]. The key enzymes in this pathway are the periplasmic thioredoxin-like protein DsbA and the transmembrane recycling factor protein DsbB. DsbA oxidises protein substrates to give rise to a disulphide-bond, whilst reduced DsbA is quickly recycled to its active oxidised form by DsbB. DsbB is then re-oxidised by ubiquinone via the electron transport chain, in which the final electron acceptor is oxygen during aerobic growth [90]. However, during growth in low oxygen or anaerobic conditions, menaquinone mediates the transfer of electrons to alternative acceptors such as nitrate or small organic molecules such as fumarate. This leads to a reduced oxidation rate for both DsbAB relative to the aerobic pathway [87, 88].

Here we varied oxygen availability to modulate disulphide-bond formation and protein folding and monitored consequent protein glycosylation levels. To vary the oxygen transfer to cells in liquid media, cultures were grown in shake flasks with different culture to flask volume ratios (5:50, 10:50, and 25:50 ml) (“Methods” section). As described scFv13R4, scFv13R4CM, and also RNase A were modified with glycosylation sequons and used as the target disulphide-bond containing proteins and NGRP (DQNAT) acted as a non-disulphide bond containing control protein.

A modest but significant increase in glycosylation efficiency was observed when disulphide-bond containing proteins were produced in higher culture to flask volume ratio (5:50 to 25:50 ml; scFv13R4 62 to 79%, P = 0.005; scFv13R4CM 41 to 52%, P = 0.028; RNase A 20 to 34%, P = 0.017) (Fig. 5A–C). A rise in glycosylation efficiency was also examined for the control NGRP when the cells were grown at higher culture to flask volume ratio (5:50 to 10:50 ml, 25% to 40%, P = 0.004) (Fig. 5D). However, in contrast to the trend observed for disulphide bond-containing proteins, NGRP glycosylation efficiency decreased when the cells were grown at the highest culture to flask volume ratio (10:50 to 25:50 ml, 40% to 26%, P = 0.017—Fig. 5D). Each protein was produced at generally similar levels in these different conditions (mean yield scFv13R4 = 4.7 ± 1.0 mg/g DCW, scFv13R4CM = 4.8 ± 0.6 mg/g DCW, RNase A = 2.9 ± 0.6 mg/g DCW, NGRP = 4.0 ± 0.8 mg/g DCW). scFv13R4 protein production saw a small, but not significant, decrease at higher culture ratios (1.3-fold, 5:50 to 25:50), and RNase A also saw a small, but not significant, increase (1.3-fold, 5:50 to 25:50) (Fig. 5A, C). In addition, correlation analysis (Pearson) showed that the change in glycosylation efficiency was not significantly associated with the production level of disulphide-bond proteins during growth in different culture to flask volume ratio (Additional file 1: Figure S14).

Fig. 5
figure5

Effect of culture conditions and oxygen availability upon target protein glycosylation in E. coli. AD Effect of culture to flask volume ratio (oxygen transfer efficiency) upon target protein glycosylation in E. coli. Quantitative Western blot analysis (densitometry) of A scFv13R4, B scFv13R4CM, C RNase A and D non-disulphide control protein NGRP located in periplasm of glyco-competent E. coli in shake flask under three different cultures to flask volume ratio (5:50, 10:50, and 25:50 ml). EH Quantitative Western blot analysis (densitometry) of E scFv13R4, F scFv13R4CM, G RNase A and H NGRP control non-disulphide bond-containing protein, detected in the periplasm of glyco-competent E. coli under different oxygen levels culture (3% and 15% O2). AH Total proteins were quantified using pre-determined purified scFv13R4CM, RNase, or NGRP standard curve (5 ng to 100 ng). The data were converted into mg/g of dry cell weight (DCW) based on normalisation and calculation with measured OD600 of the samples. AH Glycosylated (yellow bar) and non-glycosylated (green bar) protein as shown (left y-axis). % Glycosylation (% G1/G0 + G1) is indicated (black circle, right y-axis). Statistical analysis was conducted by unpaired t-test with Welch’s correction to control sample at lowest culture to flask volume ratio 5:50 (AD) or to control normoxic culture (E–H) (P < 0.05*, < 0.01**, for % glycosylation; P < 0.05, for total protein mg/g DCW). All data were processed from three biological replicates. Error bars indicate standard deviation from mean values

In order to further explore the effect of lower oxygen culture conditions upon protein glycosylation, cells were grown in a hypoxic chamber with 3% (hypoxic) or 15% O2 (normoxic) (“Methods” section). Production under hypoxic conditions also increased glycosylation efficiency of disulphide bond-containing proteins, however only scFv13R4CM showed a significant change (15% to 3% O2; % glycosylation scFv13R4 73% to 82%, scFv13R4CM 51% to 61% P = 0.027, RNase A 26% to 32%, NGRP 30% to 28%), which may also be attributed to the observed greater reduction in scFv13R4CM protein yield (scFv13R4CM from 5.9 ± 0.9 to 2.9 ± 0.5 mg/g DCW, P = 0.005) (Fig. 5E–H). To test if glycosylation of disulphide bond-containing proteins was more sensitive to the change in oxygen levels compared to glycosylation of non-disulphide bond-containing proteins (NGRP), we quantified the fold-change of glycosylation between the proteins expressed in hypoxic (treatment) and normoxic (control) conditions. This fold-change is defined as relative glycosylation efficiency (RGE) was calculated as follows:

$${\text{RGE}} = \left( {\frac{{{\text{Glycosylation}}\;{\text{efficiency}}\;{\text{under}}\;{\text{treatment}}\;\left( {\text{G1/G1 + G0}} \right)}}{{{\text{Glycosylation}}\;{\text{efficiency}}\;{\text{of}}\;{\text{the}}\;{\text{control}}\;\left( {\text{G1/G1 + G0}} \right)}}} \right)$$

Therefore, even if glycosylation of non-disulphide bond-containing protein were to change in a similar direction as the model disulphide bond-containing protein, the degree of change could be compared using the RGE metric. The RGE results showed that disulphide bond-containing proteins experienced a modest increase in glycosylation fold-change during hypoxic treatment (scFv13R4 1.1-fold, scFv13R4CM 1.2-fold, RNase A 1.3-fold) but this was significantly different to the change in glycosylation observed for NGRP under the same conditions (0.92-fold) (P = 2.0 × 10–4, 3.3 × 10–6, 8.3 × 10–5 respectively) (Additional file 1: Figure S15).

While the glycosylation efficiency of disulphide bond-containing proteins increased when produced under oxygen-depleted conditions, cell growth was also dramatically reduced (Table 1). The negative impact upon growth seems independent of target protein and most likely related to cellular oxygen requirements, since both cells expressing disulphide or non-disulphide bond-containing proteins were negatively affected. Volumetric production (titre) of recombinant protein was decreased as a result of reduced cell growth during production in oxygen depleted conditions (Table 1).

Table 1 Comparative cell growth and volumetric productivity (protein titre) during recombinant disulphide-bond protein production in oxygen-depleted conditions

Production and glycosylation of disulphide bond-containing proteins in the absence of oxidoreductase DsbB

Studies have shown that in the absence of DsbA or DsbB, in order to maintain cell viability and disulphide-bond formation either additional oxygen or media supplementation with a strong oxidant is required [91]. Nevertheless, the rate of disulphide-bond formation driven directly by oxygen supplementation in these knockout strains does not fully rescue the reaction catalysed by the DSB pathway [88, 91]. The previous experiments demonstrated that producing disulphide bond-containing proteins under suboptimal conditions for disulphide-bond formation might delay protein folding, leading to enhanced interaction time with OTase and the sequon, and greater glycan transfer. Therefore, we hypothesised that when disulphide bond-containing proteins were produced in glyco-competent E. coli lacking one of DsbAB enzymes, such as in the ΔdsbB strain, glycosylation would be enhanced. To test this, a glyco-competent ΔdsbB strain was transformed with expression vectors encoding the model disulphide bond-containing proteins and control NGRP, and glycosylation analysed following standard shake flask cultivation (“Methods” section).

Densitometric analysis of periplasmic proteins demonstrated increased glycosylation efficiency of scFv13R4 and scFv13R4CM in the ΔdsbB strain compared to the wild-type (wt) strain (scFv13R4 73% to 88%, P = 0.047; scFv13R4CM 50% to 80%, P = 1.1 × 10–4) (Fig. 6A, B). However, expression in the ΔdsbB strain showed a negative effect for glycosylation of RNase A (25% to 10%, P = 0.025) (Fig. 6C). In parallel a significant reduction in protein production (50–90%) was observed for all three disulphide-containing model proteins in ΔdsbB relative to the wt strain (scFv13R4 from 3.5 ± 0.4 to 1.0 ± 0.3 mg/g DCW, P = 1.7 × 10–3; scFv13R4CM from 5.0 ± 0.6 to 2.7 ± 0.2 mg/g DCW, P = 1.4 × 10–2; RNase A from 3.0 ± 0.1 to 0.24 ± 0.05 mg/g DCW, P = 5.9 × 10–7). We suspected that the proteins might be prone to degradation following an extended delay in disulphide-bond formation and protein folding in the absence of the DsbB enzyme [89, 90]. Indeed, consistent with this notion, RNase A production was most severely affected since it has four disulphide-bonds compared to one for scFv13R4 and scFv13R4CM (Fig. 2A ii–iii and B ii–iv). However, as glycosylation efficiency can be impacted by lower amounts of target protein substrate, increased glycosylation efficiency of scFv13R4 and scFv13R4CM cannot be attributed solely to the delay of disulphide-bond formation during expression in the ΔdsbB strain. No significant change in total protein and glycosylation efficiency of control NGRP was observed in the ΔdsbB strain in comparison to the wt strain (wt = 8.5 ± 1.4 mg/g DCW, ΔdsbB = 7.2 ± 1.2 mg/g DCW) (Fig. 6D). The result indicates the specific effect of DsbB activity upon production and glycosylation of the target disulphide bond-containing proteins.

Fig. 6
figure6

Glycosylation of disulphide bond-containing proteins in oxidoreductase mutant (ΔdsbB) of E. coli. Quantitative Western blot analysis (densitometry) of A scFv13R4, B scFv13R4CM, C RNase A and D control non-disulphide bond-containing protein NGRP expressed in periplasmic of glyco-competent E. coli wild-type (wt) or ΔdsbB strain. Total proteins (AD) were quantified using pre-determined purified scFv13R4CM, RNase A, or NGRP standard curve (5 ng to 100 ng). The data were converted into yield mg/g of dry cell weight (DCW) based on normalisation and calculation with measured OD600 of the samples Glycosylated (yellow bar) and non-glycosylated (green bar) protein as shown (left y-axis). % Glycosylation (% G1/G0 + G1) is indicated (black circle, right y-axis). Additional file 1: Figure S16, A and B showed representative Western blot images for scFv13R4CM and RNase A produced in wt and ΔdsbB strains. Representative Western blot image for NGRP and scFv13R4 have been shown before in Figs. 3B, C, and 4B, C. Statistical analysis was conducted by unpaired t-test with Welch’s correction to control sample expressed in wt strain (P < 0.05*, < 0.01**, for % glycosylation; P < 0.05, < 0.01, < 0.0001, for total protein mg/g DCW). All data were processed from three biological replicates. Error bars indicate standard deviation from mean values

Table 2 summarises the comparative growth and volumetric production (recombinant protein titre) in wild-type and ΔdsbB strains producing model glycoproteins. Reduction in cell growth was observed when scFv13R4 and scFv13R4CM were produced in ΔdsbB relative to the wt strain. The reduction in cell growth along with protein production per cell (protein yield mg/g DCW) results in lower scFv13R4 and scFv13R4CM protein titres (Fig. 6A, B; Table 2). In contrast, both wild-type and ΔdsbB strains expressing RNase A showed similar growth (Table 2). Therefore, the observed reduction in RNase A protein titre (Table 2) was predominately impacted by the decrease in protein per cell (yield) (Fig. 6C, D). No difference in the cell growth was observed between wild-type and ΔdsbB expressing non-disulphide bond-containing protein NGRP as well as the protein titre produced between both strains.

Table 2 Comparative cell growth and volumetric production (protein titre) during recombinant disulphide-bond protein production in the wild-type (wt) and ΔdsbB strain

Supplementation with chemical oxidant (cystine) improves production of disulphide bond-containing glycoproteins in the ΔdsbB strain

The above analysis shows that producing disulphide bond-containing proteins in the absence of DsbB reduces recombinant protein production (yield and titre) and cell growth. Thus, we investigated if supplementation with a small molecule oxidant (e.g. cystine) could maintain enhanced glycosylation efficiency but also recover the total target protein production levels and cell growth in the ΔdsbB strain. Moreover, alternative oxidants may provide different rates of modulation for disulphide-bond formation and protein folding compared to DsbAB enzymes. In this way glycosylation efficiency could be improved without heavily sacrificing protein production levels. To perform this experiment, wild-type (wt) and ΔdsbB strains of glyco-competent E. coli were prepared as before, but now 100 μM of cystine was added into the culture medium along with inducers (“Methods” section). Total production levels for all disulphide bond-containing proteins increased in the ΔdsbB strain upon supplementation with cystine (scFv13R4 2.0-fold, from 1.2 ± 0.2 to 2.4 ± 0.1 mg/g DCW, P = 9.6 × 10–4; scFv13R4CM 2.1-fold, from 1.3 ± 0.2 to 2.7 ± 0.6 mg/g DCW, P = 0.047; RNase A 5.0-fold, from 0.78 ± 0.11 to 3.4 ± 0.8 mg/g DCW, P = 1.4 × 10–5) (Fig. 7A–C). In addition, glycosylation efficiency was almost identical with and without cystine treatment even though increased protein yields were observed (average glycosylation efficiency ± cystine treatment; scFv13R4 ΔdsbB = 80 ± 2%, scFv13R4CM ΔdsbB = 74 ± 4%, RNase A ΔdsbB = 12 ± 3%) (Fig. 7A–C). Therefore, resulting in an increase in the absolute amounts of glycoprotein under this treatment. The scFv13R4 produced in ΔdsbB strain was also evaluated for the binding activity against its cognate antigen β-galactosidase (“Methods” section). The result showed increased specific binding activity of scFv13R4 from the ΔdsbB strain compared to that produced in the wild-type strain (ΔdsbB 0 cys = 1.8-fold, ΔdsbB 100 cys = 1.4-fold) (Additional file 1: Figure S17). It has been reported previously that in vitro glycosylation of scFv13R4 with Cj heptasaccharide does not affect protein activity [92]. However, glycosylation enhanced activity of another exemplar single chain Fv (murine 3D5 scFv) when produced in vivo (in glyco-competent E. coli), has been reported [72]. As such, this present analysis is consistent and indicates enhanced activity of the scFv13R4 glycoconjugate when produced via in vivo glycosylation.

Fig. 7
figure7

Impact of cystine supplementation upon glycosylation of recombinant proteins in ΔdsbB strain and glycosylation of RNase A in disulphide-bond isomerase mutant (ΔdsbC). (A-D) Quantitative Western blot analysis (densitometry) of A scFv13R4, B scFv13R4CM, C RNase A and D NGRP non-disulphide control protein produced in the periplasm of glyco-competent E. coli wild-type (wt) or ΔdsbB strain supplemented with or without 100 μM cystine during protein expression. E, F Quantitative Western blot analysis (densitometry) of E RNase A and F control non-disulphide bond-containing protein NGRP expressed in periplasmic of glyco-competent E. coli ΔdsbC. Total proteins (AF) were quantified using pre-determined purified scFv13R4CM, RNase, or NGRP standard curve (5 ng to 100 ng). The data were converted into mg/g of dry cell weight (DCW) based on normalisation and calculation with measured OD600 of the samples. Glycosylated (yellow bar) and non-glycosylated (green bar) protein as shown (left y-axis). % Glycosylation (% G1/G0 + G1) is indicated (black circle, right y-axis). Statistical analysis was conducted by unpaired t-test with Welch’s correction to control sample expressed in wt or ΔdsbB strain without cystine treatment (AD) or to control sample expressed in wt strain (E, F) (P < 0.05*, < 0.01**, for % glycosylation; P < 0.05, < 0.001, for normalised total protein). All data were processed from three biological replicates. Error bars indicate standard deviation from mean values

Cystine function appears to be redundant in the presence of DsbAB, as there was no improvement in the production of disulphide bond-containing proteins in the wild-type strain (Fig. 7A–C). While RNase A production in ΔdsbB strain was considerably improved by cystine supplementation, glycosylation in this strain was still negatively affected relative to the wild-type (glycosylation efficiency wt = 27 ± 5%, ΔdsbB = 12 ± 3%) (Fig. 7C). No difference in NGRP production and glycosylation were observed when produced in ΔdsbB strain with or without cystine (average yield = 5.2 ± 0.5 mg/g DCW, average glycosylation efficiency = 45 ± 4%) (Fig. 7D). The result confirmed the specific effect of cystine treatment only upon disulphide bond-containing proteins.

In addition to improvement in total protein yield and absolute amounts of glycoprotein, cystine treatment also produced a better growth phenotype for the ΔdsbB strain producing disulphide bond-containing proteins, while the growth level was not fully restored to the level observed for the untreated wild-type (Table 3). Improvement in protein production per cell (yield) and cell growth resulted in increasing disulphide-bond protein titre produced in ΔdsbB strain under cystine treatment (Table 3). The wild-type strain expressing disulphide bond-containing proteins grew similarly with or without cystine supplementation (Table 3). Further, cystine addition appeared to have no impact upon the growth of the ΔdsbB strain producing non-disulphide bond-containing protein NGRP and the protein titre was produced similarly (Table 3). These data indicate that stable production of recombinant disulphide bond-containing proteins is responsible for the growth improvement exhibited in ΔdsbB strain, rather than a more global impact upon cell redox state due to the absence of DsbB.

Table 3 Comparative cell growth and volumetric production (protein titre) during recombinant disulphide-bond protein production in the wild-type (wt) and ΔdsbB strain supplemented with small molecule oxidant (cystine)

It is worthy of mention that we also noticed a difference in glycosylation efficiency improvement among disulphide bond-containing proteins during the expression in the ΔdsbB strain compared to the production under oxygen-depleted conditions. The increase in glycosylation efficiency of scFv13R4CM in ΔdsbB strain relative to wt strain was higher (up to 1.6-fold) (Figs. 6B and 7B) than compared to production under low oxygen transfer conditions (25/50 vs. 5/50 flask to culture vol ratio) (1.3-fold) (Fig. 5B). Both conditions gave a similar improvement in glycosylation efficiency for scFv13R4 (1.2- to 1.3-fold) (Figs. 5A, 6A, 7A). While glycosylation efficiency of RNase A was enhanced during low oxygen transfer (1.7-fold) (Fig. 5C) and hypoxic growth/low oxygen level (1.3-fold) (Fig. 5G). It is unclear why RNase A glycosylation was negatively affected during expression in ΔdsbB, even more considering the expression level was also greatly reduced (Figs. 6C and 7C). We speculate that without DsbAB assistance, RNase A forms an intermediate conformation which has reduced interaction with OTase. Different strategies to enhance glycosylation may be required to optimise folding-dependent glycosylation of RNase A in vivo.

Production of RNase A in the absence of disulphide-bond isomerase (ΔdsbC) improves glycosylation efficiency

Cystine supplementation was shown to increase the production of RNase A in a ΔdsbB strain of E. coli (yield and titre) (Fig. 7C). However, glycosylation efficiency remained low compared with production in the wild-type strain. Due to the presence of non-consecutive disulphides, RNase A is more likely to be sensitive to incorrect pairing during its folding in the periplasm [90, 93]. In the periplasm, mispaired disulphide-bonds are recognised by the disulphide-bond isomerase DsbC, which catalyses reshuffling reactions of the mispaired bonds [90, 94]. Previously, Kowarik et al. demonstrated in vitro that rapid oxidation of RNase A generated mixed disulphide isoforms as preferable substrates for bacterial OTase (PglB) [23]. Another in vitro study showed the formation of intermediate disulphide isoforms were dominantly observed within RNase A in the absence of DsbC [93]. Indeed, the activity of the oxidoreductase subunit in some eukaryotic OTases leads to transient mixed disulphide forms of the target proteins, promoting an intermediate folding state, and consequently increasing their window for glycosylation [71, 95].

To investigate if RNase A glycosylation could be improved in the same way in vivo (in bacteria) by attempting to increase intermediate mixed disulphide forms, protein production was performed in a ΔdsbC strain of glyco-competent E. coli. Analysis of the periplasmic fraction showed increasing glycosylation efficiency of RNase A produced in ΔdsbC (29% to 53%, P = 0.001) (Fig. 7E). Total protein was also higher compared with production in the wt strain (1.4-fold, from 3.1 ± 0.3 to 4.3 ± 0.5 mg/g DCW, P = 0.042) (Fig. 7E). NGRP production and glycosylation was similar in both wt and ΔdsbC strain (protein = 4.4 ± 0.6 mg/g DCW, % glycosylation = 48 ± 3%) indicating no effect upon the non-disulphide bond-containing protein (Fig. 7F). Production of scFv13R4 and scFv13R4CM in the ΔdsbC strain resulted in compromised protein production levels (scFv13R4, 1.6-fold, from 4.8 ± 0.4 to 2.9 ± 0.3 mg/g DCW, P = 0.004; scFv13R4CM, 1.2-fold, from 4.8 ± 0.1 to 4.1 ± 0.4 mg/g DCW), and minor decreases in protein glycosylation efficiency (scFv13R4, from 70 to 65%; scFv13R4CM, from 49 to 39%, P = 0.002) (Additional file 1: Figure S18). While the isomerase activity is not necessary for proteins with only single disulphide-bonds or those possessing consecutive disulphide-bonds, it has been reported that DsbC displays a chaperone function to improve stability and solubility of various scFv model proteins [96,97,98].

To confirm that RNase A produced in ΔdsbC strain increases production of intermediate folding form of RNase, we analysed the activity of RNase produced from this strain (“Methods” section). A recent study reported that glycosylation of RNase A produced in glyco-competent E. coli had a similar activity with its non-glycosylated version [99]. Therefore, any decrease in the activity of RNase A could be associated with the increased presence of pre-folded intermediate forms. As we expected, the activity assay showed 3.3-fold reduction of RNase A produced in the periplasm of ΔdsbC compared to the one produced in the wt strain (Additional file 1: Figure S19). Going forward, protein folding modulation strategies could be coordinated (e.g., by titrating the expression or activity of folding modulators) to increase glycoprotein yield, glycosylation and to retain biological activity.

Surprisingly, the ΔdsbC strain expressing RNase A grew better than the wt strain expressing the same protein which also resulted in increased recombinant protein titre (Table 4). Both wild-type and ΔdsbC expressing NGRP showed little or no difference in the cell growth and protein titre (Table 4). We also observed that the expression of scFv13R4 but not the CM variant in the ΔdsbC strain negatively impacted cell growth, and consequently the recombinant protein titre (Additional file 2: Table S6). Together with the result from glycoprotein production in the ΔdsbB strain, we suggest different proteins could be glycosylated more efficiently in vivo based upon modulation of their distinct folding requirements to place them in a specific unfolded or intermediate form for optimal glycosylation.

Table 4 Comparative cell growth and volumetric production (protein titre) during RNase A production in the wild-type (wt) and ΔdsbC strain

Discussion

Our initial study with NGRP glycosylation indicated that for this particular protein, production levels had little impact upon glycosylation efficiency. Availability of lipid base, sugar precursor, and glycan conjugating enzyme could in principle be limiting factors for glycosylation, for example in the scenario where target protein production exceeds cell glycosylation capacity [33, 48, 64, 66, 69]. However, if the enzyme and lipid-glycan substrate are not rate limiting and protein glycosylation efficiency is observed similarly across different expression levels, we propose that protein glycosylation in this context could depend predominantly upon sequon accessibility to the glycosylation enzyme (PglB/oligosaccharyltransferase). In this simple way, if the sequon is not accessible for glycosylation once target protein has folded or diffused away from the membrane, the degree of glycosylation will only be dependent upon the transient unfolded state or membrane residency time during translocation, in which PglB “sees” the sequon and not the level of target protein production.

To explore the effect of membrane residency time upon protein glycosylation, we monitored the glycosylation state of proteins located in the periplasm and inner membrane. Using NGRP and scFv13R4 as the model proteins, in which the target protein was trapped in the membrane by generating a poorly processed signal cleavage site (TMT variant) attached to the target proteins. Improvement in glycosylation efficiency of TMT-scFv13R4, located in the membrane fraction, over the wt-scFv13R4 located in the periplasm fraction, seemed to agree with the hypothesis in this study. By increasing membrane residency of target protein scFv13R4, glycosylation efficiency was improved by up to 25–40%. Following extended time in the membrane, local interaction of the protein substrate and the OTase would likely increase, leading to enhanced glycosylation efficiency. However, the effectiveness of this strategy seems to be protein specific, as some factors such as protein size and structural complexity might also influence protein folding rate during membrane retention, and therefore influence sequon accessibility for glycosylation. For example, NGRP is a small protein with a predicted structure consisting of only a single α-helical chain (d =  ~ 10–12 Å). For this protein, folding could be initiated very early whilst inside the translocon tunnel (d =  ~ 20 Å) (Fig. 2A, i) [100, 101]. Previous studies showed that some proteins with more complex architecture and complex tertiary structural modifications, such as disulphide-bonds, underwent folding suppression within the translocon [84, 85]. Moreover, for many of these proteins, folding would continue for some period after translocation [102,103,104]. Taken together, we assume that simpler proteins like NGRP do not experience the same spatial and structural constraints as scFv13R4 upon folding in the translocon during extended residency time. These constraints suggest that, scFv13R4 is likely to remain unfolded during membrane retention, meaning that the sequon remains available for glycan addition for a greater length of time, resulting in greater interaction with OTase and more efficient glycosylation.

We then further investigated if modulation of target protein folding could enhance protein glycosylation. By monitoring the glycosylation efficiency of model disulphide bond-containing proteins, we found that producing the protein under suboptimal condition for disulphide-bond formation could enhance glycosylation efficiency. In this experiment, the rate of disulphide-bond synthesis was adjusted by varying the oxidants (oxygen, chemical oxidants, oxidoreductases DsbAB and DsbC) through changing culture conditions and E. coli production strain genetic backgrounds. Owing to the structural diversity of model disulphide bond-containing proteins, we also identified specific environmental and biochemical conditions to optimise glycosylation efficiency of different proteins. Sequon location was shown to play a key role is deciding which methodology is most likely to result in improved glycosylation efficiency. If the sequon is located around highly structured regions this can lead to reduced glycosylation efficiency, as shown previously [23, 71, 105, 106]. For example, we show that when located near residues involved in disulphide-bond formation, the sequon is less well glycosylated by PglB. For scFv13R4CM where the sequon is proximal to these regions, expression in the ΔdsbB strain may provide an extension in the time required for disulphide-formation, therefore increasing the sequon exposure for glycan transfer by PglB, leading to ~ 1.6-fold (from ~ 50% to 80%) enhanced glycosylation efficiency (Figs. 6B and 7B). In another example, RNase A glycosylation efficiency was improved ~ 1.9-fold (from ~ 28% to 53%) when the protein was produced in a ΔdsbC strain (Fig. 7E). As a protein with multiple disulphides, RNase A production in the ΔdsbC strain may fold into a specific intermediate form, leading to the sequon becoming more accessible for PglB compared to the one produced in ΔdsbB or under oxygen-depleted conditions. Lastly, scFv13R4 glycosylation efficiency was observed to be less sensitive, and displayed relatively similar levels across all the different disulphide attenuation treatments (1.2- to 1.3-fold or from ~ 65% to 80%). Having a distal location from cysteine residues, the availability of the R4 sequon may be less influenced by disulphide-bond formation and hence a limited effect upon interaction with PglB.

Production of disulphide bond-containing glycoproteins under these suboptimal conditions required for their folding could be problematic in terms of achieving high glycoprotein yields and titres. Indeed, cell growth was dramatically reduced during cultivation in oxygen-depleted conditions, possibly due to general metabolic burden caused by oxygen limitation. On the other hand, fitness issues during protein production in the ΔdsbB strain seemed to be related to protein specific stress conditions, since strains producing the control protein NGRP grew similarly under these conditions. While increased glycosylation efficiency was observed during the production of disulphide-bond containing protein RNase A in the ΔdsbC strain, the produced proteins were shown to be less active. Total production of disulphide bond-containing proteins was reduced 50–90% when expressed in the ΔdsbB strain. However, supplementation with the chemical oxidant cystine, was shown to rescue protein yield and cell viability whilst maintaining glycosylation efficiency of the ΔdsbB strain, therefore increasing protein titre and the absolute amount of glycoprotein produced. This result indicated the potential for using different oxidants or to tune expression levels of the native oxidoreductase to modulate folding-dependent glycosylation of target proteins and to balance with host cell physiology. To further optimise the strategy, improvement in total protein production, glycosylation efficiency, protein activity, and the effect upon the cell growth need to be considered. Therefore, those studies will be informative for the selection of strains, medium, cultivation conditions, and bioprocess design for glycoprotein production. This includes options to separate biomass production and target protein expression to avoid the physiological burden of unfolded protein stress during cell growth and to improve volumetric production (protein titre). However, it is also worthy of mention that for some production strategies, that enhanced glycosylation efficiency can be a main consideration, especially if protein purity is essential and the downstream protein purification step can be a bottleneck [48]. Moreover, if the amount of protein needed for an application is rather low, for example for use as a glycoconjugate vaccine (25–100 μg needed per vaccination), improved glycosylation efficiency may compensate for the reduction in total protein production titres [48].

As the use of different folding attenuation conditions to produce the model proteins in this study showed different pros and cons in terms of protein glycosylation efficiency, yields, and cell growth. The production strategy could also be prioritised dependent on the ultimate application and the requirement for maximising the total amount of glycoprotein (titre) vs. maximising the proportion of glycosylated product (glycosylation efficiency). The highest glycoprotein titre of the scFv model proteins used here tended to be achieved by production in the wild-type strain and under more oxygenated conditions (Fig. 8A, B). On the other hand, for production of scFv13R4 and scFv13R4CM, enhanced glycosylation efficiency was obtained by production in the ΔdsbB strain, in which improved titre could be achieved by cystine supplementation during the cultivation (Fig. 8A, B). Production of RNase A in the ΔdsbC strain resulted in increased glycosylation efficiency and glycoprotein titre, while at cost of reduced protein activity (Fig. 8C, and Additional file 1: Figure S19). Future study could explore how optimisation of folding modulators could benefit glycosylation efficiency, glycoprotein titre, and maintain biological activity of the product. Production of glycosylated NGRP in terms of glycosylation efficiency seemed to give similar results irrespective of production in the wild-type or oxidoreductase strains (Fig. 8D). The optimised condition for production of glycosylated NGRP was flask cultivation at a 10/50 culture to flask volume ratio, which gave improvement both in glycoprotein titre and glycosylation efficiency, regardless the type of producer strains (wt/ΔdsbBdsbC) (Fig. 8D). It is worthy of mention that increased NGRP glycoprotein yield and titre was achieved at higher induction expression (200 μM PPDA at the experiment in the Fig. 6D compared to 40 μM PPDA in any other experiments, “Methods” section) without significantly affecting glycosylation efficiency (Fig. 8D).

Fig. 8
figure8

Summary of glycosylation efficiency (%) and glycoprotein titre (mg/L of glycosylated protein) of A scFv13R4, B scFv13R4CM, C RNase A, and D NGRP during production in glyco-competent E. coli K-12 by use of different cultivation conditions to modulate target protein folding. Glycoprotein titres were converted from the yield (mg/g DCW), cell growth, and total protein titre results given in the Figs. 5, 6, 7, and Tables 1, 2, 3, 4 and Additional file 2: Table S6 (“Methods” section). Colour symbols indicate the experiment or figure/table sources for the data (Fig; 1 = Table 1, 2 = Table 2, 3 = Table 3, 4 = Table 4, 5 = Additional file 2: Table S6). A similar wt and ΔdsbB cultivation were run in four different batches of experiment, which were (i) wt oxygen transfer experiment (10/50, Fig. 5A–D), (ii) wt against oxidoreductase knockout ΔdsbB (Fig. 6A–D), (iii) wt and ΔdsbB 0 cystine treatment (Fig. 7A–D), (iv) wt against oxidoreductase knockout ΔdsbC (Fig. 7E, F, Additional file 1: Figure S18A, B). A variation of inducer concentration for NGRP expression was used in the experiment of Fig. 6D (200 μM instead of 40 μM PPDA, “Methods” section)

The final steps of protein translocation through the secretion machinery and subsequent folding could be viewed as a competing reaction with glycosylation, since the sequon will potentially be less available for glycan transfer upon target protein folding. In eukaryotic cells and the bacterium C. jejuni, glycosylation is coordinated temporally and spatially with protein translocation and folding pathways to increase glycosylation efficiency, and further has been mechanistically demonstrated in higher eukaryotic (mammalian) cells [22, 26]. This is demonstrated by coupling of transmembrane OTase activities with Sec translocation. The subunits of some eukaryotic heteroligomeric Otases possess oxidoreductase activity which interfere with disulphide-bond formation of the target proteins, leading to formation of transient intermediate protein conformations more accessible for glycosylation [71, 95]. However, expression of heterologous glycosylation machinery in E. coli, such as the pgl pathway from C. jejuni, appears to be uncoordinated and less efficient at protein glycosylation than in the native host. One explanation for this reduced glycosylation efficiency of the heterologous host is the potential lack of accessory proteins for modulating the translocation process and protein folding environment thereby affecting target protein sequon accessibility for the OTase [26]. Informed by glycosylation events in the native system, efforts to modulate sequon interaction with PglB could be a potential approach to optimise glycoprotein production in E. coli, the global aim of this present study.

Based on the general principle of folding-dependent glycosylation, the strategy performed here could be similarly adopted by using specific conditions for the production of non-disulphide bond-containing proteins. Recently, DeLisa and co-workers [92] demonstrated the versatility of cell-free glycoprotein production systems for rapid screening of functional Otases and glycosyltransferases from different bacteria and eukaryotes. Due to the open nature and modularity of this platform, screening could be rapidly and simply expanded using collections of chaperones and folding catalysts to improve glycosylation of new target proteins. Incorporating native or engineered chaperone candidates from E. coli or even C. jejuni may lead to identification of accessory proteins promoting more efficient protein glycosylation [107,108,109,110]. Moreover, the rapid characterisation of folding modulator-target protein combinations in vitro would allow early validation before further optimisation in vivo.

Lastly, the study presented here also may also show the potential of increasing membrane residency time of the target protein during Sec translocation to enhance glycosylation efficiency. While this study provides a proof of concept for this strategy, the use of uncleaved signal peptide would provide practical challenges in the downstream process of protein production, such as membrane purification and signal peptide removal. For further work, the use of alternative protein delivery mechanisms to the periplasm via the signal recognition particle (SRP) pathway may provide further opportunities to coordinate translation, translocation and glycosylation [24]. Since translocation occurs co-translationally via the SRP, modifying codon usage of the gene encoding the target protein could be used to influence transport rate [111], and hence glycosylation efficiency.

Conclusion

In this study, we demonstrated that improved glycosylation efficiency in the heterologous host could be achieved by mimicking the coordination between protein translocation, folding and glycosylation observed in native such as Campylobacter jejuni and mammalian hosts. Furthermore, it provides insight into strain engineering and bioprocess strategies, to improve glycoprotein yield, titre and to avoid physiological burden of unfolded protein stress to cell growth. In addition to existing strategies to enhance protein glycosylation in E. coli, the work presented here provides new approaches that can be used in order to develop more robust recombinant glycoengineering platforms.

Methods

Strains, plasmids, and general growth conditions

Cloning strains: E. coli NEB5α (NEB). Expression strain: E. coli Top10F’ (Thermo Fisher); E. coli K-12 wild -type (wt), dsbB::kandsbB) and dsbC::kandsbC) with inactivated dsbB or dsbC due to insertion of a kanamycin-resistance gene were used as indicated (Keio collection) [112]. Plasmids: Four genes encoding target proteins; Cj0114 NGRP, anti-β-galactosidase single-chain Fv scFv13R4 and R4CM, and bovine pancreatic ribonuclease RNase A, all modified with N-terminal PelB signal peptide, C-terminal hexahistidine (His6) tag, and single bacterial glycosylation sequon were synthesised by gBlocks (Integrated DNA Technologies). The genes were subcloned into ampicillin resistance determinant-containing pDEST-ORS expression vector [76]. Chloramphenicol resistance determinant-containing pACYC184 plasmid carrying the entire C. jejuni pgl locus (pACYCpgl) [15] was transformed into the expression strain to allow glycoprotein production (glyco-competent). General growth conditions: E. coli NEB5α strains were routinely grown in LB medium (0.5% yeast extract, 0.5% NaCl, 1.0% Bactotryptone) to enrich constructed plasmids. For glycoprotein production studies, expression strains of E. coli were grown in LB medium as starter cultures and in TB medium (2.7% yeast extract, 4.5% glycerol, 1.3% Bactotryptone) with 0.2% glucose for protein expression. All media were supplemented with antibiotics as required (ampicillin 100 μg/ml, kanamycin 50 μg/ml, and chloramphenicol 25 μg/ml). Protein expression was induced by the addition of IPTG (Isopropyl β-d-1-thiogalactopyranoside) and PPDA (Pyrimido [4,5-d] pyrimidine-2,4-diamine) as indicated.

In silico analysis of protein structure and translocation

Homology modelling of NGRP, scFv13R4, and scFv13R4CM was performed using Phyre2 (protein homology/analogy recognition engine) [113]. Ribbon models were drawn by UCSF Chimera [114]. Signal peptide profile and protein translocation of NGRP/scFv13R4 wt and TMT were predicted using SignalP4.1 [81].

PCR site-directed mutagenesis

Primers used to generate mutant versions of the PelB cleavage site (TMT) are listed in Additional file 2: Table S1 and were supplied by Integrated DNA Technologies. Inserts NGRP or scFv13R4 and linearised plasmid (backbone) pDEST-ORS-PelB TMT mutant were generated from the template pDEST-ORS-PelB wt NGRP or scFv13R4 by PCR: 98 °C for 30 s; 35 cycles of 10 s at 98 °C, 30 s at the primer annealing temperature (gradient), and extension at 72oCfor 30 s per 1 kb; with a final extension step of 72 °C for 10 min. PCR products were treated with DpnI (NEB) and purified using a gel extraction kit (Qiagen). Insert NGRP or scFv13R4 was assembled into backbone pDEST-ORS-PelB TMT following NEB builder Hi-fi mix protocol (NEB) to create pDEST-ORS-NGRP/scFv13R4 TMT. Plasmid products were cloned into E. coli NEB5α and isolated by miniprep according to the manufacturers (Qiagen) protocol. All constructs were confirmed by DNA sequencing.

Growth conditions for glycoprotein production

To prepare inocula for glycoprotein production tests, expression strains of E. coli bearing pACYCpgl (glyco-competent) were transformed with pDEST-ORS constructs specific for each target protein. Freshly plated colony transformants were picked and grown overnight (37 °C, 200 rpm) in LB medium supplemented with antibiotics. The starter culture then was transferred (1% v/v) into TB medium (with added glucose and antibiotics) in a shake flask and grown to OD600 = 0.6 (37 °C, 200 rpm). At this point, cells were induced with IPTG (100 μM) and PPDA (0–400 μM or 0–200 μM for titration assay; 400 μM during membrane residency test of wt and TMT variant; 40 μM for NGRP during expression in oxygen-depleted conditions, in ΔdsbC strain, and in ΔdsbB strain with cystine supplementation; 200 μM for NGRP during expression in ΔdsbB strain; 0 μM for scFv13R4 and scFv13R4CM, and 200 μM for RNase A in all conditions for redox attenuation experiment) and then divided into three for different cultivation strategies. Microculture: 1 ml of cultures in 96-well deep well plate, 37 °C 1000 rpm for 3 h (titration assay) or 6 h (membrane residency test of wt and TMT variant). Shake flask: 10 ml (standard condition) or 5, 10 and 25 ml (treatment conditions, as indicated) of culture in 50 ml shake flask, 37 °C 200 rpm for 4 h or 6 h (titration assay of wt and TMT variant only). 6-well: 5 ml of cultures in 6-well plate, incubated in a Don Whitley VA500 microaerobic cabinet set to 3% or 15% O2, 37 °C 100 rpm for 4 h.

Cell fractionation

Cell fractionation was performed from the final culture of glycoprotein production experiments after the specified time post-induction. Cultures were normalised to the same ODV 1–5 (OD600 x ml = 1–5) to collects periplasmic and membrane fractions. The cultures with specific ODV were centrifuged (17,000g, 4 °C, 1 min) and cell pellets resuspended in 250 µl Buffer 1 (500 mM Sucrose, 5 mM EDTA, 100 mM Tris–acetate pH 8.2), followed by addition of lysozyme (0.16 mg/ml) and Milli-Q water (250 μl). Cells were incubated on ice for 5 min, and 10 µl of 1 M MgSO4 added to stabilise the spheroplast. Periplasm (supernatant) was collected from spheroplast (pellet) by centrifugation (17,000g, 4 °C, 10 min). To produce membrane fractions, spheroplasts were washed with 500 µl Buffer 2 (250 mM sucrose, 10 mM MgSO4, 50 mM Tris–acetate pH 8.2), resuspended in 500 µl Buffer 3 (2.5 mM EDTA, 0.1% sodium deoxycholate, 2U benzonase, 50 mM Tris–acetate pH 8.2), and stored overnight at − 80 °C. The following day, spheroplasts were lysed by freeze-thawing and cytoplasm (supernatant) and membrane fractions (pellet) separated by centrifugation (17,000g, 4 °C, 30 min). Pellets were washed with 500 μl Buffer 3 and resuspended in 500 μl of Buffer 3 + SDS (1%) to obtain the membrane fraction.

Western blot analysis

The amount of glycoprotein produced in periplasmic and membrane fractions was analysed by Western blot. Samples were resuspended in SDS-PAGE loading buffer (ThermoFisher), supplemented with 75 mM dithiothreitol, and boiled for 10 min. Equal amounts of sample were separated by SDS PAGE (Biorad) and Western blotted onto PVDF 0.2 µM pore size membrane using a Turbo Transblot apparatus (Biorad). The membrane was blocked with 5% (w/v) milk powder in PBS (30 min, 50 rpm shaking, room temperature). Membranes were incubated overnight (50 rpm, 4 °C) with mouse monoclonal anti-His antibody (1:3000 dilution in PBS 5% milk, Pierce) and rabbit polyclonal anti-Cj glycan (CjNgp) primary antibody(1:500 dilution in PBS 5% milk) as appropriate. Membranes were washed three times with PBS, followed by incubation with anti-mouse and anti-rabbit secondary antibody (1:30,000 dilution in PBS 5%, Li-Cor) (50 rpm, in the dark at room temperature, 30 min), and then washed three times with PBS. The protein bands were visualised with a Li-Cor Odyssey scanner. Proteins were quantified based on densitometry analysis using Li-Cor Image Studio 5.0. Pre-determined purified protein from each sample (scFv13R4, RNase A, or NGRP) was used as a standard for quantification (5 ng to 100 ng). All data were produced in biological triplicate.

Data processing and statistical analysis

Data were processed and analysed using Microsoft Excel and GraphPad Prism9. Each data point was the average from three biological replicates, and error bars show calculated standard deviation. For Western blot quantification, protein samples were measured from each biological replicate of isolated cell fraction (periplasmic or membrane). Additional file 1: Figure S20 shows a representative Western blots images for blot reproducibility of PelB 2-NGRP (representative from Fig. 3C) expressed in periplasmic fraction from three biological replicates. A standard curve was constructed from pre-determined purified protein (scFv13R4, RNase A, or NGRP: 5–100 ng), and used to calibrate Western blot signal intensity of protein samples. The data was fitted with linear regression, and the linear equation from the standard curve was used to normalise and convert the signal into ng of the protein. Representative standard curve for NGRP and the use of standard curve for protein quantification are shown in Additional file 1: Figure S21. Conversion of protein quantification data into the yield mg/g of dry cell weight (DCW) based on normalisation and calculation from the measured OD600 of samples. Dry cell weight was determined following the previous protocol [76], and the conversion factor of 0.35 mg/ml per OD600 was acquired. Protein titre mg/l was determined from protein quantification data (ng) converted and normalised from OD600 of collected samples and and volume of sample (μl) for Western blot. Statistical analysis was conducted by unpaired t-test with Welch’s correction, and P value < 0.05 was considered statistically significant.

scFv13R4 binding assay

scFv13R4 binding activity against its cognate antigen β-galactosidase was performed following bio-dot protocol [115]. 2 μl of β-galactosidase (0.3 mg/ml) dots were spotted onto a nitrocellulose membrane (Amersham Hydrobond), and the membrane then left to dry for 5–15 min. To block non-specific binding site, the membrane was incubated with 5% (w/v) milk powder in PBS for 1 h (50 rpm, RT). The membrane then was left to dry for 2 h in 40 °C. To prepare binding mix, each 2 μl of protein sample from isolated periplasmic fraction was added onto the dot of β-galactosidase, with 2 μl of protein buffer was added instead as a negative control. The membrane then left to dry again for 15 min. Afterwards, the membrane was probed with primary (mouse monoclonal anti-his) and secondary (anti-mouse) antibody to detect the scFv13R4 following the similar protocol with the Western blot. The protein dots were visualised with the Li-Cor Odyssey scanner and the signal intensity was analysed by densitometry using Image Studio 5.0 Software. Specific binding activity of was measured by normalising the densitometric signal with the amount of protein sample. All data were produced in biological triplicate.

RNase A activity assay

RnaseAlert lab test kit (Invitrogen) was used to assay RNase A activity of the protein sample. The reaction mix consist of 10 μl of RNase A substrate (20× dilution with RNase A buffer solution), 5 μl of protein sample from periplasmic fraction (1000×), and 85 μl of nuclease-free water. As a negative control, 5 μl of protein buffer or 5 μl of protein sample from periplasmic fraction of non RNase A producer strain was added instead into the reaction mix. The reaction mix was transferred into 96-well plate, and then incubated in the Clariostar plate reader for 1 h 37 °C. The RNase A substrate from the kit is a modified RNA nucleotide that emits a green fluorescence when it is degraded by RNase. The fluorescence activity of degraded substrate was recorded overtime during the incubation (RFU at 490/520 nm). RNase A specific activity was quantified by normalising the rate of substrate degradation (RFU/minute) with the amount of protein sample. All data were produced in biological triplicate.

Availability of data and materials

All data generated or analysed in this study are included in this manuscript or in Additional files. Raw data are available from the corresponding author upon reasonable request.

Change history

References

  1. 1.

    Nothaft H, Szymanski CM. Protein glycosylation in bacteria: sweeter than ever. Nat Rev Microbiol. 2010;8:765–78.

    CAS  PubMed  Google Scholar 

  2. 2.

    Aebi M. N-linked protein glycosylation in the ER. Biochim Biophys Acta Mol Cell Res. 2013;1833:2430–7.

    CAS  Google Scholar 

  3. 3.

    Eichler J, Koomey M. Sweet new roles for protein glycosylation in prokaryotes. Trends Microbiol. 2017;25:662–72.

    CAS  PubMed  Google Scholar 

  4. 4.

    Apweiler R, Hermjakob H, Sharon N. On the frequency of protein glycosylation, as deduced from analysis of the SWISS-PROT database. Biochim Biophys Acta Gen Subj. 1999;1473:4–8.

    CAS  Google Scholar 

  5. 5.

    Abu-Qarn M, Eichler J, Sharon N. Not just for Eukarya anymore: protein glycosylation in Bacteria and Archaea. Curr Opin Struct Biol. 2008;18:544–50.

    CAS  PubMed  Google Scholar 

  6. 6.

    Solá RJ, Griebenow K. Glycosylation of therapeutic proteins: an effective strategy to optimize efficacy. BioDrugs. 2011;24:9–21.

    Google Scholar 

  7. 7.

    Xu C, Ng DTW. Glycosylation-directed quality control of protein folding. Nat Rev Mol Cell Biol. 2015;16:742–52.

    CAS  PubMed  Google Scholar 

  8. 8.

    Varki A. Biological roles of glycans. Glycobiology. 2017;27:3–49.

    CAS  PubMed  Google Scholar 

  9. 9.

    Poole J, Day CJ, Itzstein MV, Paton JC, Jennings MP. Glycointeractions in bacterial pathogenesis. Nat Rev Microbiol. 2018;16:440–52.

    CAS  PubMed  Google Scholar 

  10. 10.

    Larkin A, Imperiali B. The expanding horizons of asparagine-linked glycosylation. Biochemistry. 2011;50:4411–26.

    CAS  PubMed  Google Scholar 

  11. 11.

    Szymanski CM, Logan SM, Linton D, Wren BW. Campylobacter—a tale of two protein glycosylation systems. Trends Microbiol. 2003;11:233–8.

    CAS  PubMed  Google Scholar 

  12. 12.

    Schwarz F, Lizak C, Fan YY, Fleurkens S, Kowarik M, Aebi M. Relaxed acceptor site specificity of bacterial oligosaccharyltransferase in vivo. Glycobiology. 2011;21:45–54.

    CAS  PubMed  Google Scholar 

  13. 13.

    Baker JL, Çelik E, DeLisa MP. Expanding the glycoengineering toolbox: the rise of bacterial N-linked protein glycosylation. Trends Biotechnol. 2013;31:313–23.

    CAS  PubMed  Google Scholar 

  14. 14.

    Szymanski CM, Yao R, Ewing CP, Trust TJ, Guerry P. Evidence for a system of general protein glycosylation in Campylobacter jejuni. Mol Microbiol. 1999;32:1022–30.

    CAS  PubMed  Google Scholar 

  15. 15.

    Wacker M, Linton D, Hitchen PG, Nita-Lazar M, Haslam SM, North SJ, Panico M, Morris HR, Dell A, Wren BW, Aebi M. N-linked glycosylation in Campylobacter jejuni and its functional transfer into E. coli. Science. 2002;298:1790–3.

    CAS  PubMed  Google Scholar 

  16. 16.

    Szymanski CM, Wren BW. Protein glycosylation in bacterial mucosal pathogens. Nat Rev Microbiol. 2005;3:225–37.

    CAS  PubMed  Google Scholar 

  17. 17.

    Mohorko E, Glockshuber R, Aebi M. Oligosaccharyltransferase: the central enzyme of N-linked protein glycosylation. J Inherit Metab Dis. 2011;34:869–78.

    CAS  PubMed  Google Scholar 

  18. 18.

    Kowarik M, Young NM, Numao S, Schulz BL, Hug I, Callewaert N, Mills DC, Watson DC, Hernandez M, Kelly JF, Wacker M, Aebi M. Definition of the bacterial N-glycosylation site consensus sequence. EMBO J. 2006;25:1957–66.

    CAS  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Zielinska DF, Gnad F, Schropp K, Wiśniewski JR, Mann M. Mapping N-glycosylation sites across seven evolutionarily distant species reveals a divergent substrate proteome despite a common core machinery. Mol Cell. 2012;46:542–8.

    CAS  PubMed  Google Scholar 

  20. 20.

    Kelleher DJ, Gilmore R. An evolving view of the eukaryotic oligosaccharyltransferase. Glycobiology. 2006;16:47–62.

    Google Scholar 

  21. 21.

    Ruiz-Canada C, Kelleher DJ, Gilmore R. Cotranslational and posttranslational N-glycosylation of polypeptides by distinct mammalian OST isoforms. Cell. 2009;136:272–83.

    CAS  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Shrimal S, Cherepanova NA, Gilmore R. Cotranslational and posttranslocational N-glycosylation of proteins in the endoplasmic reticulum. Semin Cell Dev Biol. 2015;41:71–8.

    CAS  PubMed  Google Scholar 

  23. 23.

    Kowarik M, Numao S, Feldman MF, Schulz BJ, Callewaert N, Kiermaier E, Catrein I, Aebi M. N-linked glycosylation of folded. Science. 2006;314:1148–50.

    CAS  PubMed  Google Scholar 

  24. 24.

    Fisher AC, Haitjema CH, Guarino C, Çelik E, Endicott CE, Reading CA, Merritt JH, Ptak AC, Zhang S, DeLisa MP. Production of secretory and extracellular N-linked glycoproteins in Escherichia coli. Appl Environ Microbiol. 2011;77:871–81.

    CAS  PubMed  Google Scholar 

  25. 25.

    Lizak C, Gerber S, Numao S, Aebi M, Locher KP. X-ray structure of a bacterial oligosaccharyltransferase. Nature. 2011;474:350–5.

    CAS  PubMed  Google Scholar 

  26. 26.

    Silverman JM, Imperiali B. Bacterial N-glycosylation efficiency is dependent on the structural context of target sequons. J Biol Chem. 2016;291:22001–10.

    CAS  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Sethuraman N, Stadheim TA. Challenges in therapeutic glycoprotein production. Curr Opin Biotechnol. 2006;17:341–6.

    CAS  PubMed  Google Scholar 

  28. 28.

    Walsh G. Biopharmaceutical benchmarks 2014. Nat Biotechnol. 2014;32:992–1000.

    CAS  PubMed  Google Scholar 

  29. 29.

    Walsh G. Biopharmaceutical benchmarks 2018. Nat Biotechnol. 2018;36:1136–45.

    CAS  PubMed  Google Scholar 

  30. 30.

    Kim JY, Kim Y-G, Lee GM. CHO cells in biotechnology for production of recombinant proteins: current state and further potential. Appl Microbiol Biotechnol. 2012;93:917–30.

    CAS  PubMed  Google Scholar 

  31. 31.

    Noh SM, Sathyamurthy M, Lee GM. Development of recombinant Chinese hamster ovary cell lines for therapeutic protein production. Curr Opin Chem Eng. 2013;2:391–7.

    Google Scholar 

  32. 32.

    Palmberger D, Wilson IBH, Berger I, Grabherr R, Rendic D. Sweetbac: a new approach for the production of mammalianised glycoproteins in insect cells. PLoS ONE. 2012;7:e34226.

    CAS  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Pandhal J, Woodruff LBA, Jaffe S, Desai P, Ow SY, Noirel J, Gill RT, Wright PC. Inverse metabolic engineering to improve Escherichia coli as an N -glycosylation host. Biotechnol Bioeng. 2013;110:2482–93.

    CAS  PubMed  Google Scholar 

  34. 34.

    Burnett MJB, Burnett AC. Therapeutic recombinant protein production in plants: challenges and opportunities. Plants People Planet. 2020;2:121–32.

    Google Scholar 

  35. 35.

    Schmidt FR. Recombinant expression systems in the pharmaceutical industry. Appl Microbiol Biotechnol. 2004;65:363–72.

    CAS  PubMed  Google Scholar 

  36. 36.

    O’Flaherty R, Bergin A, Flampouri E, Mota LM, Obaidi I, Quigley A, Xie Y, Butler M. Mammalian cell culture for production of recombinant proteins: a review of the critical steps in their biomanufacturing. Biotechnol Adv. 2020;43:107552.

    PubMed  Google Scholar 

  37. 37.

    Vieira Gomes A, Souza Carmo T, Silva Carvalho L, Mendonça Bahia F, Parachin N. Comparison of yeasts as hosts for recombinant protein production. Microorganisms. 2018;6:38.

    PubMed Central  Google Scholar 

  38. 38.

    Schillberg S, Raven N, Spiegel H, Rasche S, Buntru M. Critical analysis of the commercial potential of plants for the production of recombinant proteins. Front Plant Sci. 2019;10:720.

    PubMed  PubMed Central  Google Scholar 

  39. 39.

    Yee CM, Zak AJ, Hill BD, Wen F. The coming age of insect cells for manufacturing and development of protein therapeutics. Ind Eng Chem Res. 2018;57:10061–70.

    CAS  PubMed  PubMed Central  Google Scholar 

  40. 40.

    Lehle L, Strahl S, Tanner W. Protein glycosylation, conserved from yeast to man: a model organism helps elucidate congenital human diseases. Angew Chemie Int Ed. 2006;45:6802–18.

    CAS  Google Scholar 

  41. 41.

    Choi B-K, Bobrowicz P, Davidson RC, Hamilton SR, Kung DH, Li H, Miele RG, Nett JH, Wildt S, Gerngross TU. Use of combinatorial genetic libraries to humanize N-linked glycosylation in the yeast Pichia pastoris. Proc Natl Acad Sci U S A. 2010;107:7113.

    CAS  Google Scholar 

  42. 42.

    De Pourcq K, De Schutter K, Callewaert N. Engineering of glycosylation in yeast and other fungi: current state and perspectives. Appl Microbiol Biotechnol. 2010;87:1617–31.

    PubMed  Google Scholar 

  43. 43.

    Pandhal J, Ow SY, Noirel J, Wright PC. Improving N-glycosylation efficiency in Escherichia coli using shotgun proteomics, metabolic network analysis, and selective reaction monitoring. Biotechnol Bioeng. 2011;108:902–12.

    CAS  PubMed  Google Scholar 

  44. 44.

    Jaffé SRP, Strutton B, Levarski Z, Pandhal J, Wright PC. Escherichia coli as a glycoprotein production host: Recent developments and challenges. Curr Opin Biotechnol. 2014;30:205–10.

    PubMed  Google Scholar 

  45. 45.

    Natarajan A, Jaroentomeechai T, Li M, Glasscock CJ, DeLisa MP. Metabolic engineering of glycoprotein biosynthesis in bacteria. Emerg Top Life Sci. 2018;2:419–32.

    CAS  PubMed  Google Scholar 

  46. 46.

    Keys TG, Aebi M. Engineering protein glycosylation in prokaryotes. Curr Opin Syst Biol. 2017;5:23–31.

    Google Scholar 

  47. 47.

    Wacker M, Feldman MF, Callewaert N, Kowarik M, Clarke BR, Pohl NL, Hernandez M, Vines ED, Valvano MA, Whitfield C, Aebi M. Substrate specificity of bacterial oligosaccharyltransferase suggests a common transfer mechanism for the bacterial and eukaryotic systems. Proc Natl Acad Sci. 2006;103:7088–93.

    CAS  PubMed  PubMed Central  Google Scholar 

  48. 48.

    Ihssen J, Kowarik M, Dilettoso S, Tanner C, Wacker M, Thöny-Meyer L. Production of glycoprotein vaccines in Escherichia coli. Microb Cell Fact. 2010;9:61.

    PubMed  PubMed Central  Google Scholar 

  49. 49.

    Iwashkiw JA, Fentabil MA, Faridmoayer A, Mills DC, Peppler M, Czibener C, Ciocchini AE, Comerci D, J., Ugalde, J. E., Feldman, M. F. . Exploiting the Campylobacter jejuni protein glycosylation system for glycoengineering vaccines and diagnostic tools directed against brucellosis. Microb Cell Fact. 2012;11:13.

    CAS  PubMed  PubMed Central  Google Scholar 

  50. 50.

    Cuccui J, Thomas RM, Moule MG, D’Elia RV, Laws TR, Mills DC, Williamson D, Atkins TP, Prior JL, Wren BW. Exploitation of bacterial N-linked glycosylation to develop a novel recombinant glycoconjugate vaccine against Francisella tularensis. Open Biol. 2013;3:130002–130002.

    PubMed  PubMed Central  Google Scholar 

  51. 51.

    Nothaft H, Davis B, Lock YY, Perez-Munoz ME, Vinogradov E, Walter J, Coros C, Szymanski CM. Engineering the Campylobacter jejuni N-glycan to create an effective chicken vaccine. Sci Rep. 2016;6:26511.

    CAS  PubMed  PubMed Central  Google Scholar 

  52. 52.

    Kay E, Cuccui J, Wren BW. Recent advances in the production of recombinant glycoconjugate vaccines. NPJ Vaccines. 2019;4:16.

    PubMed  PubMed Central  Google Scholar 

  53. 53.

    Schwarz F, Huang W, Li C, Schulz BL, Lizak C, Palumbo A, Numao S, Neri D, Aebi M, Wang LX. A combined method for producing homogeneous glycoproteins with eukaryotic N-glycosylation. Nat Chem Biol. 2010;6:264–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  54. 54.

    Valderrama-Rincon JD, Fisher AC, Merritt JH, Fan Y-Y, Reading CA, Chhiba K, Heiss C, Azadi P, Aebi M, DeLisa MP. An engineered eukaryotic protein glycosylation pathway in Escherichia coli. Nat Chem Biol. 2012;8:434–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  55. 55.

    Srichaisupakit A, Ohashi T, Misaki R, Fujiyama K. Production of initial-stage eukaryotic N-glycan and its protein glycosylation in Escherichia coli. J Biosci Bioeng. 2015;119:399–405.

    CAS  PubMed  Google Scholar 

  56. 56.

    Yates LE, Mills DC, DeLisa MP. Bacterial glycoengineering as a biosynthetic route to customized glycomolecules. Adv Biochem Eng Biotechnol. 2018;175:167–220.

    Google Scholar 

  57. 57.

    Ollis AA, Zhang S, Fisher AC, DeLisa MP. Engineered oligosaccharyltransferases with greatly relaxed acceptor-site specificity. Nat Chem Biol. 2014;10:816–22.

    CAS  PubMed  PubMed Central  Google Scholar 

  58. 58.

    Mills DC, Jervis AJ, Abouelhadid S, Yates LE, Cuccui J, Linton D, Wren BW. Functional analysis of N-linking oligosaccharyl transferase enzymes encoded by deep-sea vent proteobacteria. Glycobiology. 2016;26:398–409.

    CAS  PubMed  Google Scholar 

  59. 59.

    Ihssen J, Haas J, Kowarik M, Wiesli L, Wacker M, Schwede T, Thöny-Meyer L. Increased efficiency of Campylobacter jejuni N-oligosaccharyltransferase PglB by structure-guided engineering. Open Biol. 2015;5:140227.

    PubMed  PubMed Central  Google Scholar 

  60. 60.

    Merritt JH, Ollis AA, Fisher AC, DeLisa MP. Glycans-by-design: engineering bacteria for the biosynthesis of complex glycans and glycoconjugates. Biotechnol Bioeng. 2013;110:1550–64.

    CAS  PubMed  Google Scholar 

  61. 61.

    Feldman MF, Wacker M, Hernandez M, Hitchen PG, Marolda CL, Kowarik M, Morris HR, Dell A, Valvano MA, Aebi M. Engineering N-linked protein glycosylation with diverse O antigen lipopolysaccharide structures in Escherichia coli. Proc Natl Acad Sci USA. 2005;102:3016–21.

    CAS  PubMed  PubMed Central  Google Scholar 

  62. 62.

    Linton D, Dorrell N, Hitchen PG, Amber S, Karlyshev AV, Morris HR, Dell A, Valvano MA, Aebi M, Wren BW. Functional analysis of the Campylobacter jejuni N-linked protein glycosylation pathway. Mol Microbiol. 2005;55:1695–703.

    CAS  PubMed  Google Scholar 

  63. 63.

    Guarino C, Delisa MP. A prokaryote-based cell-free translation system that efficiently synthesizes glycoproteins. Glycobiology. 2012;22:596–601.

    CAS  PubMed  Google Scholar 

  64. 64.

    Wayman JA, Glasscock C, Mansell TJ, DeLisa MP, Varner JD. Improving designer glycan production in Escherichia coli through model-guided metabolic engineering. Metab Eng Commun. 2019;9:e00088.

    PubMed  PubMed Central  Google Scholar 

  65. 65.

    Yates LE, Natarajan A, Li M, Hale ME, Mills DC, DeLisa MP. Glyco-recoded Escherichia coli: recombineering-based genome editing of native polysaccharide biosynthesis gene clusters. Metab Eng. 2019;53:59–68.

    CAS  PubMed  Google Scholar 

  66. 66.

    Pandhal J, Desai P, Walpole C, Doroudi L, Malyshev D, Wright PC. Systematic metabolic engineering for improvement of glycosylation efficiency in Escherichia coli. Biochem Biophys Res Commun. 2012;419:472–6.

    CAS  PubMed  Google Scholar 

  67. 67.

    Strutton B, Jaffé SRP, Pandhal J, Wright PC. Producing a glycosylating Escherichia coli cell factory: the placement of the bacterial oligosaccharyl transferase pglB onto the genome. Biochem Biophys Res Commun. 2018;495:686–92.

    CAS  PubMed  Google Scholar 

  68. 68.

    Ding N, Yang C, Sun S, Han L, Ruan Y, Guo L, Hu X, Zhang J. Increased glycosylation efficiency of recombinant proteins in Escherichia coli by auto-induction. Biochem Biophys Res Commun. 2017;485:138–43.

    CAS  PubMed  Google Scholar 

  69. 69.

    Strutton B, Jaffe SRP, Evans CA, Fowler GJS, Dobson PD, Pandhal J, Wright PC. Engineering pathways in central carbon metabolism help to increase glycan production and improve N-type glycosylation of recombinant proteins in E. coli. Bioengineering. 2019;6:27.

    CAS  PubMed Central  Google Scholar 

  70. 70.

    Yusuf SNHM, Bailey U-M, Tan NY, Jamaluddin MF, Schulz BL. Mixed disulfide formation in vitro between a glycoprotein substrate and yeast oligosaccharyltransferase subunits Ost3p and Ost6p. Biochem Biophys Res Commun. 2013;432:438–43.

    Google Scholar 

  71. 71.

    Mohorko E, Owen RL, Malojčić G, Brozzo MS, Aebi M, Glockshuber R. Structural basis of substrate specificity of human oligosaccharyl transferase subunit N33/Tusc3 and its role in regulating protein N-glycosylation. Structure. 2014;22:590–601.

    CAS  PubMed  Google Scholar 

  72. 72.

    Lizak C, Fan Y, Weber TC, Aebi M. N-Linked glycosylation of antibody fragments in Escherichia coli. Bioconjug Chem. 2011;22:488–96.

    CAS  PubMed  Google Scholar 

  73. 73.

    Scott NE, Parker BL, Connoly AM, Paulech J, Edwards AVG, Crossett B, Falconer L, Kolarich D, Djordjevic SP, Højrup P, Packer NH, Larsen MR, Cordwell SJ. Simultaneous glycan-peptide characterization using hydrophilic interaction chromatography and parallel fragmentation by CID, higher energy collisional dissociation, and electron transfer dissociation MS applied to the N-linked glycoproteome of C. jejuni. Mol Cell Proteom. 2011;10:M000031-MCP201.

    Google Scholar 

  74. 74.

    Martineau P, Betton JM. In vitro folding and thermodynamic stability of an antibody fragment selected in Vivo for high expression levels in Escherichia coli cytoplasm. J Mol Biol. 1999;292:921–9.

    CAS  PubMed  Google Scholar 

  75. 75.

    Dixon N, Duncan JN, Geerlings T, Dunstan MS, McCarthy JEG, Leys D, Micklefield J. Reengineering orthogonally selective riboswitches. Proc Natl Acad Sci USA. 2010;107:2830–5.

    CAS  PubMed  PubMed Central  Google Scholar 

  76. 76.

    Horga LG, Halliwell S, Castiñeiras TS, Wyre C, Matos CFRO, Yovcheva DS, Kent R, Morra R, Williams SG, Smith DC, Dixon N. Tuning recombinant protein expression to match secretion capacity. Microb Cell Fact. 2018;17:199.

    CAS  PubMed  PubMed Central  Google Scholar 

  77. 77.

    Kent R, Halliwell S, Young K, Swainston N, Dixon N. Rationalizing context-dependent performance of dynamic RNA regulatory devices. ACS Synth Biol. 2018;7:1660–8.

    CAS  PubMed  Google Scholar 

  78. 78.

    Fikes JD, Barkocy-Gallagher GA, Klapper DG, Bassford PJ Jr. Maturation of Eschericia coli maltose-binding protein by signal peptidase I in vivo. J Biol Chem. 1990;265:3417–23.

    CAS  PubMed  Google Scholar 

  79. 79.

    Shen LM, Lee J, Cheng S, Jutte H, Kuhn A, Dalbey RE. Use of Site-directed mutagenesis to define the limits of sequence variation tolerated for processing of the M13 procoat protein by the Escherichia coli leader peptidase. Biochemistry. 1991;30:11775–81.

    CAS  PubMed  Google Scholar 

  80. 80.

    Auclair SM, Bhanu MK, Kendall DA. Signal peptidase I: cleaving the way to mature proteins. Protein Sci. 2012;21:13–25.

    CAS  PubMed  Google Scholar 

  81. 81.

    Petersen TN, Brunak S, Von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8:785–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  82. 82.

    Käll L, Krogh A, Sonnhammer ELL. Advantages of combined transmembrane topology and signal peptide prediction-the Phobius web server. Nucleic Acids Res. 2007;35:429–32.

    Google Scholar 

  83. 83.

    Weaver, D. N-linked glycosylation in Campylobacter jejuni and Campylobacter fetus N-glycans as targets for antibody-based detection. Ph.D. thesis, The University of Manchester; 2017.

  84. 84.

    Chen W, Helenius A. Role of ribosome and translocon complex during folding of influenza hemagglutinin in the endoplasmic reticulum of living cells. Mol Biol Cell. 2000;11:765–72.

    CAS  PubMed  PubMed Central  Google Scholar 

  85. 85.

    Ellgaard L, McCaul N, Chatsisvili A, Braakman I. Co- and post-translational protein folding in the ER. Traffic. 2016;17:615–38.

    CAS  PubMed  Google Scholar 

  86. 86.

    Mamathambika BS, Bardwell JC. Disulfide-linked protein folding pathways. Annu Rev Cell Dev Biol. 2008;24:211–35.

    CAS  PubMed  Google Scholar 

  87. 87.

    Bader M, Muse W, Zander T, Bardwell J. Reconstitution of a protein disulfide catalytic system. J Biol Chem. 1998;273:10302–7.

    CAS  PubMed  Google Scholar 

  88. 88.

    Bader M, Muse W, Ballou DP, Gassner C, Bardwell JCA. Oxidative protein folding is driven by the electron transport system. Cell. 1999;98:217–27.

    CAS  PubMed  Google Scholar 

  89. 89.

    Merdanovic M, Clausen T, Kaiser M, Huber R, Ehrmann M. Protein quality control in the bacterial periplasm. Annu Rev Microbiol. 2011;65:149–68.

    CAS  PubMed  Google Scholar 

  90. 90.

    Ke N, Berkmen M. Production of disulfide-bonded proteins in Escherichia coli. Curr Protoc Mol Biol. 2014;108:16.1B.1-16.1B.21.

    Google Scholar 

  91. 91.

    Meehan BM, Landeta C, Boyd D, Beckwith J. The disulfide bond formation pathway is essential for anaerobic growth of Escherichia coli. J Bacteriol. 2017;199:1–9.

    Google Scholar 

  92. 92.

    Jaroentomeechai T, Stark JC, Natarajan A, Glasscock CJ, Yates LE, Hsu KJ, Mrksich M, Jewett MC, DeLisa MP. Single-pot glycoprotein biosynthesis using a cell-free transcription-translation system enriched with glycosylation machinery. Nat Commun. 2018;9:2686.

    PubMed  PubMed Central  Google Scholar 

  93. 93.

    Bader MW, Xie T, Yu CA, Bardwell JCA. Disulfide bonds are generated by quinone reduction. J Biol Chem. 2000;275:26082–8.

    CAS  PubMed  Google Scholar 

  94. 94.

    Shevchik VE, Condemine G, Robert-Baudouy J. Characterization of DsbC, a periplasmic protein of Erwinia chrysanthemi and Eschericia coli with disulfide isomerase activity. EMBO J. 1994;13:2007–12.

    CAS  PubMed  PubMed Central  Google Scholar 

  95. 95.

    Cherepanova NA, Shrimal S, Gilmore R. Oxidoreductase activity is necessary for N-glycosylation of cysteine-proximal acceptor sites in glycoproteins. J Cell Biol. 2014;206:525–39.

    CAS  PubMed  PubMed Central  Google Scholar 

  96. 96.

    Chen J, Song J-L, Zhang S, Wang YC, D-F., Wang, C-C. . Chaperonic activity of DsbC. J Biol Chem. 1999;274:19601–5.

    CAS  PubMed  Google Scholar 

  97. 97.

    Jurado P, Ritz D, Beckwith J, De Lorenzo V, Fernández LA. Production of functional single-chain Fv antibodies in the cytoplasm of Escherichia coli. J Mol Biol. 2002;320:1–10.

    CAS  PubMed  Google Scholar 

  98. 98.

    Zhang Z, Li Z-H, Wang F, Fang M, Yin C-C, Zhou Z-Y, Lin Q, Huang H-L. Overexpression of DsbC and DsbG markedly improves soluble and functional expression of single-chain Fv antibodies in Escherichia coli. Protein Expr Purif. 2002;26:218–28.

    CAS  PubMed  Google Scholar 

  99. 99.

    Li M, Zheng X, Shanker S, Jaroentomeechai T, Moeller TD, Hulbert SW, Kocer I, Byrne J, Cox EC, Fu Q, Zhang S, Labonte JW, Gray JJ, DeLisa MP. Shotgun scanning glycomutagenesis: a simple and efficient strategy for constructing and characterizing neoglycoproteins. Preprint at bioRxiv; 2021.

  100. 100.

    Bhushan S, Gartmann M, Halic M, Armache JP, Jarasch A, Mielke T, Berninghausen O, Wilson DN, Beckmann R. α-Helical nascent polypeptide chains visualized within distinct regions of the ribosomal exit tunnel. Nat Struct Mol Biol. 2010;17:313–7.

    CAS  PubMed  Google Scholar 

  101. 101.

    Tu LW, Deutsch C. A folding zone in the ribosomal exit tunnel for Kv1.3 helix formation. J Mol Biol. 2010;396:1346–60.

    CAS  PubMed  PubMed Central  Google Scholar 

  102. 102.

    Miot M, Betton JM. Protein quality control in the bacterial periplasm. Microb Cell Fact. 2004;65:149–68.

    Google Scholar 

  103. 103.

    Kadokura H, Beckwith J. Detecting folding intermediates of a protein as it passes through the bacterial translocation channel. Cell. 2009;138:1164–73.

    CAS  PubMed  PubMed Central  Google Scholar 

  104. 104.

    Braakman I, Hebert DN. Protein folding in the endoplasmic reticulum. Cold Spring Harb Perspect Biol. 2013;5:a013201.

    PubMed  PubMed Central  Google Scholar 

  105. 105.

    Allen S, Naim HY, Bulleid NJ. Intracellular folding of tissue-type plasminogen activator. Effects of disulfide bond formation on N-linked glycosylation and secretion. J Biol Chem. 1995;270:4797–804.

    CAS  PubMed  Google Scholar 

  106. 106.

    McGinnes LW, Morrison TG. Disulfide bond formation is a determinant of glycosylation site usage in the hemagglutinin-neuraminidase glycoprotein of Newcastle disease virus. J Virol. 1997;71:3083–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  107. 107.

    Kadokura H, Tian H, Zander T, Bardwell JCA, Beckwith J. Snapshots of DsbA in action: detection of proteins in the process of oxidative folding. Science. 2004;303:534–7.

    CAS  PubMed  Google Scholar 

  108. 108.

    Antonoaea R, Fϋrst M, Nishiyama KI, Mϋller M. The periplasmic chaperone PpiD interacts with secretory proteins exiting from the SecYEG translocon. Biochemistry. 2008;47:5649–56.

    CAS  PubMed  Google Scholar 

  109. 109.

    Matern Y, Barion B, Behrens-Kneip S. PpiD is a player in the network of periplasmic chaperones in Escherichia coli. BMC Microbiol. 2010;10:251.

    PubMed  PubMed Central  Google Scholar 

  110. 110.

    Kale A, Phansopa C, Suwannachart C, Craven CJ, Rafferty JB, Kelly DJ. The virulence factor PEB4 (Cj0596) and the periplasmic protein Cj1289 are two structurally related SurA-like chaperones in the human pathogen Campylobacter jejuni. J Biol Chem. 2011;286:21254–65.

    CAS  PubMed  PubMed Central  Google Scholar 

  111. 111.

    Yu CH, Dang Y, Zhou Z, Wu C, Zhao F, Sachs MS, Liu Y. Codon Usage influences the local rate of translation elongation to regulate co-translational protein folding. Mol Cell. 2015;59:744–54.

    CAS  PubMed  PubMed Central  Google Scholar 

  112. 112.

    Baba T, Ara T, Hasegawa M, Takai Y, Okumura Y, Baba M, Datsenko KA, Tomita M, Wanner BL, Mori H. Construction of Escherichia coli K-12 in-frame, single-gene knockout mutants: the Keio collection. Mol Syst Biol. 2006;2(2006):0008.

    PubMed  Google Scholar 

  113. 113.

    Kelly LA, Mezulis S, Yates CM, Wass MN, Sternberg MJE. The Phyre2 web portal for protein modelling, prediction, and analysis. Nat Protoc. 2015;10:845–58.

    Google Scholar 

  114. 114.

    Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, Ferrin TE. UCSF Chimera—a visualization system for exploratory research and analysis. J Comput Chem. 2004;25:1605–12.

    CAS  PubMed  Google Scholar 

  115. 115.

    Morra R, Del Carratore F, Muhamadali H, Horga LG, Halliwell S, Goodacre R, Breitling R, Dixon N. Translation stress positively regulates mscl-dependent excretion of cytoplasmic proteins. MBio. 2018;9:1–18.

    Google Scholar 

Download references

Acknowledgements

We would like to thank Adrian Jervis for help with glycoprotein analysis, SynBioChem for access to the Keio strain collection, and the members of Dixon lab for technical assistance and comments on the manuscript. Concept figures (Figs. 1 and 3A) were created with BioRender illustration tool.

Funding

FP received a PhD scholarship from LPDP Indonesia Endowment Fund for Education (S-903/LPDP.3/2016). ND was supported by a BBSRC David Phillips Fellowship (BB/K014773/1).

Author information

Affiliations

Authors

Contributions

ND and DL designed and coordinated the study. FP planned and performed the experiments. ND and FP analysed the data. FP and ND wrote the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Neil Dixon.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors disclose no conflicts.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Additional figures

Additional file 2: Table S1.

Primers used in PCR site-directed mutagenesis to generate TMT mutant of PelB NGRP and scFv13R4 in pDEST-ORS construct. Table S2. The effect of oxygen transfer conditions (culture to flask volume ratio) to cell growth and protein titre (mg/L) during production of model disulphide-bond proteins in glyco-competent E. coli Top10F’. Table S3. The effect of oxygen level conditions (15% and 3% O2) to cell growth and protein titre (mg/L) during production of model disulphide-bond proteins in glyco-competent E. coli Top10F’. Table S4. The effect of dsbB knockout to cell growth and protein titre (mg/L) during production of model disulphide-bond proteins in glyco-competent E. coli Top10F’. Table S5. The effect of cystine supplementation to cell growth and protein titre (mg/L) during production of model disulphide-bond proteins in glyco-competent E. coli Top10F’ wt and ΔdsbB strains. Table S6. The effect of dsbC knockout to cell growth and protein titre (mg/L) during production of model disulphide-bond proteins in glyco-competent E. coli Top10F’.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Pratama, F., Linton, D. & Dixon, N. Genetic and process engineering strategies for enhanced recombinant N-glycoprotein production in bacteria. Microb Cell Fact 20, 198 (2021). https://doi.org/10.1186/s12934-021-01689-x

Download citation

Keywords

  • Glycoengineering
  • N-glycosylation
  • Bacterial host engineering
  • Protein folding
  • Process optimisation