Skip to main content

Exploring the use of leucine zippers for the generation of a new class of inclusion bodies for pharma and biotechnological applications

Abstract

Background

Inclusion bodies (IBs) are biologically active protein aggregates forming natural nanoparticles with a high stability and a slow-release behavior. Because of their nature, IBs have been explored to be used as biocatalysts, in tissue engineering, and also for human and animal therapies. To improve the production and biological efficiency of this nanomaterial, a wide range of aggregation tags have been evaluated. However, so far, the presence in the IBs of bacterial impurities such as lipids and other proteins coexisting with the recombinant product has been poorly studied. These impurities could strongly limit the potential of IB applications, being necessary to control the composition of these bacterial nanoparticles. Thus, we have explored the use of leucine zippers as alternative tags to promote not only aggregation but also the generation of a new type of IB-like protein nanoparticles with improved physicochemical properties.

Results

Three different protein constructs, named GFP, J-GFP-F and J/F-GFP were engineered. J-GFP-F corresponded to a GFP flanked by two leucine zippers (Jun and Fos); J/F-GFP was formed coexpressing a GFP fused to Jun leucine zipper (J-GFP) and a GFP fused to a Fos leucine zipper (F-GFP); and, finally, GFP was used as a control without any tag. All of them were expressed in Escherichia coli and formed IBs, where the aggregation tendency was especially high for J/F-GFP. Moreover, those IBs formed by J-GFP-F and J/F-GFP constructs were smaller, rougher, and more amorphous than GFP ones, increasing surface/mass ratio and, therefore, surface for protein release. Although the lipid and carbohydrate content were not reduced with the addition of leucine zippers, interesting differences were observed in the protein specific activity and conformation with the addition of Jun and Fos. Moreover, J-GFP-F and J/F-GFP nanoparticles were purer than GFP IBs in terms of protein content.

Conclusions

This study proved that the use of leucine zippers strategy allows the formation of IBs with an increased aggregation ratio and protein purity, as we observed with the J/F-GFP approach, and the formation of IBs with a higher specific activity, in the case of J-GFP-F IBs. Thus, overall, the use of leucine zippers seems to be a good system for the production of IBs with more promising characteristics useful for pharma or biotech applications.

Background

Inclusion bodies (IBs) are protein nanoparticles ranging from 50 to 800 nm formed during the production of recombinant proteins and often localized in the cytoplasmic space of bacterial cells [1], although their formation in the periplasm has also been described [2]. In the past, IBs have been regarded as an undesired byproduct of recombinant protein production processes [3]. However, over the last 15 years it has been broadly proven that proteins forming IBs are biologically active [4,5,6,7]. The presence of native proteins forming such aggregates has prompted to assess their potential as biomaterials for a wide range of applications, including biocatalysis, tissue engineering, and antimicrobial and cancer therapies [7,8,9,10]. They offer important advantages over their soluble counterparts such as high stability [11], slow-release behavior [12, 13], and production through cost-effective processes [14]. Aiming to improve the production and biological efficiency of this nanomaterial, the increase of the aggregation tendency of proteins of interest is a key aspect and has been evaluated through the use of different aggregation tags, being VP1 [6], GFIL8 [15] and ELK16 [16] three representative examples. Besides, the scale-up of IB production and purification protocols have also been improved during the last decade [14, 17]. It has been described that the recombinant protein forming this biomaterial coexists with other proteins such as chaperones, and also with lipids [18, 19]. However, the presence of these impurities has been poorly studied, and consequently there is room to optimize IB composition. This is an important challenge to be solved in terms of IB applicability, since for specific applications it is important to control the exact composition of these bacterial nanoparticles. Thus, exploring ways to control the IB formation process emerges as a central strategy to improve IB purity and better control their physicochemical properties. In the present study, we have used an alternative approach for the generation of a new type of IB-like protein nanoparticles using leucine zippers (LZ) as aggregation-seeding domains with the aim to drive protein aggregation and improve IB properties for industrial applications. In marked contrast to aggregation tags used till present to increase protein aggregation, LZ are protein–protein interaction domains consisting of amphipathic α helices that dimerize in parallel, either as homodimers or heterodimers, to form a coiled-coil [20,21,22]. LZ dimerization motifs have already been explored as protein–protein interaction drivers both in recombinant mammalian [23] and bacterial [24, 25] cells. Thus, LZ specific properties make these peptides promising players to improve, control, and stabilize IB quality, thus boosting their potential for pharma and biotech industrial applications. In this study, we have explored the aggregation profile of a GFP reporter protein fused to Jun and Fos LZ [23] at different positions. The purity and formation of GFP IBs have been studied to determine the possible role of LZ in obtaining purer and better controlled IBs.

Results

Construct design and modelling

In this study three different protein constructs, named GFP, J-GFP-F and J/F-GFP, have been engineered to evaluate the effect of Jun and Fos LZ on the protein aggregation process (Fig. 1a). J-GFP-F is a single fusion protein consisting on a GFP flanked by Jun and Fos at N- and C-terminus, respectively, whereas J/F-GFP is constituted by two proteins (a GFP with Jun LZ at N-terminus (J-GFP) and a GFP with Fos at N-terminus (F-GFP)), simultaneously coexpressed (Fig. 1). GFP protein has been used as a control without any LZ tag. The three-dimensional structural arrangements of the constructs were built by iterative threading, taking advantage of the I-TASSER webserver [26]. The generated models were visually inspected and possible three-dimensional arrangements of the J/F units with respect to the GFP were selected as representative conformational arrangements (Fig. 1b). In all constructs, the J/F subunits tended to assume a helical structure, as expected, but it was only possible to generate highly ordered starting domains for J/F-GFP. All models underwent 250 ns of atomistic molecular dynamics simulations in full water solvent (Additional file 1: Figure S1) showing the J/F elongated domains to be structurally unstable leading to partially disordered arrangements similar to those initially modelled for J-GFP-F. Among the generated models for J-GFP-F, no highly ordered structure was present, and the pool differed by the local arrangement of the J/F subdomains (as exemplified by two structures in Fig. 1b).

Fig. 1
figure 1

Recombinant construct. a Diagram of the protein-based constructs and their possible aggregation process as IBs depending on the presence or not of LZ. b Representative three-dimensional models of the J-GFP-F (two models: a and b) and J/F-GFP constructs (one model each), as generated by iterative threading. Construct domains are color coded as follow: GFP (green), Jun (blue), Fos (red)

Protein production and aggregation

Cell growth was determined under the overexpression of the three different constructs used in this study (Fig. 2). Interestingly, a reduction in the bacterial growth was observed when J-GFP-F and J/F-GFP were produced, especially at 5 h (Fig. 2). Moreover, those proteins that impaired the growth (J-GFP-F and J/F-GFP) were produced at lower levels than their GFP counterpart (Table 1). Analyzing the protein yields at different production time points, it can be observed that there is a time-dependent production of both GFP and J/F-GFP, while for J-GFP-F production values kept similar at different times post-induction (Table 1).

Fig. 2
figure 2

Optical density values of recombinant bacteria cultures in LB medium after protein expression induction. Circles represent E. coli pET22b/GFP (GFP), triangles E. coli pET22b/Jun-GFP-Fos (J-GFP-F), and squares represent E. coli pETDuet-1-Jun-GFP/Fos-GFP (J/F-GFP). Different letters depict significant differences between the growth curves (p ≤ 0.0001)

Table 1 Total protein yield (μg) of each recombinant construct at different times post-induction (1, 3, and 5 h)

Despite of the reduction in protein levels and in cell viability for those proteins carrying Jun and Fos (Fig. 2 and Table 1), all constructs aggregated (Fig. 3). Specifically, the analysis of the aggregation ratio showed differences among the constructs used, being the aggregation ratio higher for J/F-GFP than for GFP and J-GFP-F (Fig. 3 and Additional file 1: Table S1).

Fig. 3
figure 3

Protein aggregation ratio (%) for each construct over time. Black, dark grey, and light grey bars represent GFP, J-GFP-F, and J/F-GFP, respectively. Significant differences for construct (p ≤ 0.05) and for time (p ≤ 0.1; Additional file 1: Table S1)

Characteristics of purified IBs

Aiming to explore the specific nanoarchitectural characteristics of the protein nanoparticles formed using the three different tag combinations, IBs of each construct were purified. In all samples, qualitative and quantitative approaches with high resolution electron microscopy imaging of the ultrastructural morphometric of IBs showed a high number of well-formed nanoparticles with round shape and nanoscale size in the three constructs (Fig. 4). However, slight and important shape and size differences between GFP IBs and those nanoparticles formed with Jun or Fos LZ were detected. First of all, control GFP IBs showed a very homogeneous round shape and smooth surface, suggesting a highly compact structure (Fig. 4). In contrast, Jun-Fos based IBs showed a high variability in shape, significantly more amorphous than GFP ones and with a high percentage of IBs showing rough and porous surface. About the size, GFP nanoparticles showed the highest size with a diameter of 400–500 nm and a large mean area (Fig. 4), whereas Jun-GFP-Fos and Jun-GFP/Fos-GFP IBs were significantly smaller with a lower area and a diameter around 250 nm in both cases (Fig. 4).

Fig. 4
figure 4

Representative FESEM images of the isolated IBs for each construct: GFP IBs, J-GFP-F IBs and J/F-GFP IBs. Bars size: 200 nm. IB mean area (nm2), mean diameter (nm) and roundness (%) was calculated for each construct IB (****p ≤ 0.0001)

The analysis of the protein quality in terms of protein activity of GFP, J-GFP-F, and J/F-GFP IBs indicated that the addition of Jun and Fos improved the functional protein content when flanking the GFP, whereas no differences were observed for J/F-GFP compared to GFP (Fig. 5a). This indicates that the GFP forming J-GFP-F nanoparticles has better quality and, in consequence, higher fluorescence per μg of protein than the two other constructs used. Interestingly, the comparison of the specific fluorescence of GFP forming the IBs with the soluble version of this protein indicated that although the activity of the soluble form of all the proteins was higher than when forming IBs, the difference between soluble and IB-forming proteins was especially greater in the case of GFP (Fig. 5b).

Fig. 5
figure 5

Specific fluorescence of the protein constructs. IB specific fluorescence for the three constructs (a). Specific fluorescence ratio of the insoluble fraction compared to the soluble fraction specific fluorescence (b). Black, dark grey, and light grey bars represent GFP, J-GFP-F, and J/F-GFP, respectively. * shows statistically significant differences compared to the GFP construct (p ≤ 0.05), and different letters also show statistically significant differences (p ≤ 0.05)

The conformational properties of the proteins embedded in the IBs were investigated by Fourier transform infrared (FTIR) spectroscopy [27]. The spectra were collected before and after hydrogen/deuterium (H/D) exchange to allow a better assignment of the Amide I band components to the protein secondary structures [28, 29]. GFP and J/F-GFP IBs displayed comparable IR response, while J-GFP-F IBs showed distinct peak positions and relative intensities of the spectral components assigned to β-sheets (Fig. 6 and Additional file 1: Figure S3). In particular, the absorption spectra (both before and after H/D exchange) of J-GFP-F IBs were more similar to those of the soluble GFP (Additional file 1: Figure S3) and the main β-sheet peak showed a higher downshift after H/D exchange compared to GFP and J/F-GFP IBs (Fig. 6). In the deuterated J-GFP-F, this component occurred at around 1620 cm−1, a peak position near to that observed for the main β-sheet band of soluble GFP, which indeed was observed around 1622 cm−1 (Fig. 6) [30, 31]. Overall, these data suggested that the Jun and Fos motifs modulated the conformational features of the formed IBs.

Fig. 6
figure 6

Second derivatives of the FTIR absorption spectra of GFP, J-GFP-F and J/F-GFP IBs and of soluble GFP. Samples were measured in form of protein films, obtained by solvent evaporation, and after re-hydration by D2O in order to allow H/D exchange. Selected peaks and the typical spectral region of the different protein secondary structures after H/D exchange are indicated. α, α-helices; β, β-sheets; βt, β-turns; RC, random coils

Purity of IBs

Interestingly, J-GFP-F and J/F-GFP nanoparticles have a degree of protein purity more than twice than that observed for parental GFP nanoparticles (Fig. 7). Besides, the analysis of lipid and carbohydrate content in IBs indicated that the presence of Jun and Fos in J-GFP-F increased the presence of these co-contaminants while the strategy for the production of IBs using J/F-GFP showed no significant differences with GFP in lipid and carbohydrate content (Fig. 8).

Fig. 7
figure 7

Protein purity of IBs. For each construct, protein purity was determined by densiometric image-analysis, comparing the amount of the specific recombinant construct to the whole protein amount of the IB. Different letters depict statistically significant differences (p ≤0.05)

Fig. 8
figure 8

Lipid (a) and carbohydrate (b) content of for each IB construct. Different letters indicate statistically significant differences (p ≤0.01 (a) and p ≤0.05, (b))

Discussion

IBs are protein aggregates ranging at nanoscale that have been widely studied from different perspectives. Since most of the recombinant proteins form IBs when overexpressed in bacteria [3], and some of them are only produced in this insoluble format, these nanoparticles have been broadly used as a source to obtain the soluble form of a wide range of protein of interest [32]. On the other hand, by being protein nanoparticles rich in functional recombinant protein, IBs have also been explored as a new class of biomaterial with promising applications in biocatalysis, tissue engineering, and human and animal therapies [7,8,9, 11]. However, although different aggregation tags have been used to promote their formation, the impact of these aggregation domains in the quality of protein aggregates has not been addressed so far. For that, in this work, we have explored if the use of Jun and Fos LZ dimerization motifs could drive the formation of IBs in a more controlled way in terms of quality. The results obtained proved that the constructs containing Jun and Fos (Fig. 1a), although they showed some toxicity to the producer cells (Fig. 2), they were produced as IBs with aggregation ratios higher than GFP (Fig. 3), proving that the aggregation propensity can be improved by using this strategy. The enhancement of the aggregation propensity can be attributed to the formation of Jun and Fos interactions among constructs, confirming that while the Jun and Fos domains seem to explore a number of partially disordered conformations when free in solution, this does not hinder the formation of stable aggregates. Interestingly, the aggregation ratio was observed to be higher for J/F-GFP than for J-GFP-F (Fig. 3). This can be rationalized in terms of competition between intramolecular and intermolecular Jun/Fos interactions. In J-GFP-F the Jun and Fos domains are entangled due to intramolecular interactions between them (Additional file 1: Figure S1). Their association within the same construct competes with the formation of dimers, or multimers, where Jun and Fos fragments belonging to different molecules interact and bond with each other. In both cases the two interacting fragments are the same, thus their energetic is expected to be close. At room temperature we might expect the competition and coexistence of monomers, dimers, and multimers in solution. However, it is also worth noting that in J-GFP-F, each Jun unit will be close to a Fos domain; this corresponds to a high local concentration of the partner, thus favoring the intramolecular entangled J-GFP-F. On the other hand, when the Jun/Fos domains are uncoupled, independently bound to different GFPs (as in J/F-GFP), there is obviously no competition between intermolecular and intramolecular interactions.

The study of the specific physicochemical characteristics of purified GFP, J-GFP-F and J/F-GFP IBs showed that, depending on the LZ strategy used, it’s possible to modulate IBs features such as specific activity, size, protein purity and presence of contaminants such as lipids and carbohydrates. The presence of LZ has a significant impact on both the size and shape of the aggregates, with IB diameters around 250 nm and more amorphous forms, while GFP IBs have a size of 400–500 nm and a higher average surface and higher roughness (Fig. 4). However, despite the differences in size, which can be correlated with the lower production yields achieved for those constructs containing LZ (Table 1), surface rugosity and roundness, ultrastructural morphometry of J-GFP-F and J-/F-GFP IBs (Fig. 4) agree with other conventional IBs produced in E. coli [33] and other cell factories such as Lactococcus lactis [34] or Pichia pastoris [35]. Interestingly, more amorphous shape and rough surface of J-GFP-F and J/F-GFP IBs could be indicative of significant differences in protein production and nanoparticle formation. In fact, the nanoarchitectural aspect of GFP IBs shows nanostructures with higher size, compactness and smooth surface than Jun and Fos IBs, which appear as more soft particles with a rough surface, more amorphous shape, and lower mean size. These nanoscale differences can play a great role in differential functionalities of proteins forming IBs in potential therapeutical applications. Then, low size, amorphous shape, and rough surfaces increase surface/mass ratio and, therefore, potential protein release at in vitro and in vivo conditions, becoming suitable and desired morphometric characteristics for more efficiently releasing nano or micro-platforms of drug delivery systems. Moreover, although smaller, GFP proteins forming J-GFP-F nanoparticles had a specific activity significantly higher than the other constructs that we tested (Fig. 5a), which is in accordance with FTIR spectra (Fig. 6 and Additional file 1: Figure S3). In particular, after H/D exchange, the main β-sheet peak of J-GFP-F IBs was observed to be very close to that observed for the soluble GFP, indicating the maintenance of native-like conformational features. Although J/F-GFP IBs showed no differences in the quality of the protein forming the IBs when compared to GFP (Fig. 5a), both J-GFP-F and J/F-GFP proteins forming IBs presented activities closer to the soluble form than GFP produced without tags (Fig. 5b). This indicates that the presence of Jun and Fos had a positive impact on the IB protein quality. This result agrees with previous studies that have described that the recombinant protein produced or the strain used can have an impact on the conformational quality of the proteins forming IBs [36,37,38].

On the other hand, the analysis of the elements forming such aggregates showed that J-GFP-F and J/F-GFP IBs had less protein impurities (Fig. 7), which suggests that Jun and Fos sequences drive a more controlled formation of the protein nanoparticles in terms of protein composition. By contrast, the presence of lipids and carbohydrates could not be decreased by the addition of LZ (Fig. 8). The amount of both lipids and carbohydrates in J/F-GFP IBs was comparable to the levels in GFP IBs, whereas in the case of J-GFP-F was even higher (Fig. 8). This could be correlated with the LZ interactions in each construct. As previously discussed, in J-GFP-F the Jun and Fos domains of different molecules compete with intramolecular interactions between the two domains, and this could drive the formation of less compact aggregates which could contain more lipids and carbohydrates.

Overall, these results demonstrate that aggregation-seeding domains based on LZ peptide–peptide interaction can drive the formation of a specific type of IBs, improving their quality in terms of protein content (Fig. 7), and in one of the approaches increasing specific activity (Figs. 5a and 6). However, those IBs with higher specific activity (J-GFP-F) are produced through a strategy that has a negative impact in the presence of contaminants such as lipids and carbohydrates (Fig. 8). On the contrary, the strategy based on the coexpression of Jun-GFP and Fos-GFP (J/F-GFP) to form IBs, even though it did not display any increase in the protein quality (Fig. 5a), had no negative impact in the content of lipids and carbohydrates (Fig. 8), therefore showing that it could be a promising approach for the production of IBs with higher recombinant protein content and less protein impurities.

Conclusions

Altogether this study proved that the use of Jun and Fos LZ is a good strategy for the production of IBs with more appealing characteristics that might be useful for pharma or biotech applications. This is especially relevant for J/F-GFP approach, which allowed producing hybrid IBs with increased aggregation ratio and protein purity without negatively affecting their activity and lipid and carbohydrate content.

Methods

Construction of expression plasmids

The sequence encoding amino acid residues 2–238 of the enhanced GFP (EGFP) was fused downstream of sequence encoding Fos 118–210 (bFos) or Jun 257–318 (bJun) using the linker sequence encoding SGGGSGGS to construct Fos-GFP and Jun-GFP, respectively. For the Jun-GFP-Fos construct, the sequence encoding Jun 257–318 (bJun) was fused at the N-terminus, while the sequence encoding Fos 118–210 (bFos) was fused at the C-terminus, using in both cases the linker sequence encoding SGGGSGGS. The Jun-GFP-Fos construct was cloned into pET22b (AmpR) vector (pET22b-Jun-GFP-Fos), while Jun-GFP and Fos-GFP constructs were cloned in pETDuet-1 (AmpR) to co-express them (pETDuet-1-Jun-GFP/Fos-GFP). As a control, residues 2–238 of EGFP were cloned into pET22b. The DNA sequences corresponding to each gene sequence were codon optimized for its expression in Escherichia coli (GeneArt, Germany).

Preparation of E. coli competent cells

Escherichia coli BL21 (DE3) cultures were grown overnight (ON) in LB medium at 37 °C with shaking at 250 rpm. A 1/100 inoculum was done in 50 ml of LB and the culture was grown until the optical density (OD600nm) reached a value between 0.2 and 0.4. After that, cultures were centrifuged (4000×g) at 4 °C for 15 min. Pellets were resuspended in 12.5 ml of cold and sterile 50 mM CaCl2 and incubated for 45 min in an ice bath. Cells were centrifuged again as described above and resuspended in 1.25 ml of cold and sterile 50 mM CaCl2 in glycerol (15% v/v) to prepare aliquots of 200 µl, which were stored at − 80 °C. To transform the cells, 40 ng of plasmid DNA were added to 200 µl of competent cells. The mixtures were incubated on ice for 30–60 min, warmed up to 42 °C for 45 s and placed on ice for 30 s. After incubation, 800 µl of LB media were added, and transformed cells were incubated at 37 °C for 1 h. Finally, the cells were plated on LB-agar plates containing the corresponding antibiotic.

Protein production

Escherichia coli BL21 (DE3)/pET22b-GFP, E. coli BL21(DE3)/pET22b-Jun-GFP-Fos (J-GFP-F) and E. coli BL21(DE3)/pETDuet-1-Jun-GFP/Fos-GFP (J/F-GFP) ON cultures were inoculated in 50 ml of LB media with 100 µg/ml ampicillin in 200-ml flasks at an initial OD600nm = 0.05. Each culture was grown at 37 °C and 250 rpm until the OD600nm was 0.5 and 1 mM isopropyl-β-d-thiogalactoside (IPTG) was added to induce recombinant protein expression. At times 0, 1, 3, and 5 h after IPTG induction, 1 ml samples were collected for protein fractioning. These cultures were performed by triplicate.

Protein fractionation

Samples of 1 ml were harvested by centrifugation at 6000×g at 4 °C for 15 min and the pellet was resuspended in 0.5 ml phosphate buffered saline (PBS) supplemented with protease inhibitor (Complete EDTA-free, Roche, Switzerland) to prevent protein proteolysis. Then, ice-jacketed samples were disrupted by sonication (2 cycles of 1.5 min at 10% amplitude under 0.5 s cycles) (Branson Ultrasonic SA, Switzerland). These samples were centrifuged at 15,000×g and 4 °C for 15 min to separate soluble and insoluble fractions protein fractions. Samples were stored at − 80 °C.

Protein determination

Soluble and insoluble protein fractions were analyzed by denaturing SDS-PAGE (15% acrylamide). Denaturing buffer (Laemli 4×: Tris base 1.28 g, glycerol 8 ml, SDS 1.6 g, β-mercaptoethanol 4 ml, urea 9.6 g in 20 ml) was added to the insoluble and soluble fractions to a final concentration of 1× (see protein fractioning). Soluble and insoluble protein fractions were boiled for 10 and 45 min, respectively. At that time, samples were loaded onto the gel. SDS-PAGE protein bands were transferred onto PVDF membranes and identified using a commercial anti-GFP antibody (1:1000, sc-9996, Santa Cruz Biotechnology, USA), followed by an incubation with a secondary ALP-conjugated anti-mouse IgG (whole molecule) antibody (1:20,000, A4313, Sigma-Aldrich, USA). The amounts of recombinant protein were estimated by comparison with known amounts (usually ranging from 125 to 1000 ng) of T22-GFP protein [36]. Protein bands were visualized with a solution of NBT/BCIP (B6404, Sigma-Aldrich, USA), and images were obtained using a Color Image Scanner. ImageJ software was used to perform densitometric analyses of the bands.

Determination of specific fluorescence

Fluorescence intensity of the three constructs in 1 ml samples was determined in a Varian Cary Eclipse fluorescence spectrometer (Agilent Technologies, Australia) at excitation and emission wavelengths of 480 and 510 nm, respectively. Since protein amounts in 1 ml samples for each construct, fraction and time are known (see protein determination), specific fluorescence was calculated by dividing the fluorescence intensity by the amount of protein for each construct, time and fraction.

Protein purity assessment

After separation, the gels were carefully transferred to a plastic tray filled with 200 ml of distilled water and agitated at 50 rpm to remove SDS traces. The staining solution (Coomassie Brilliant BlueR-250 Staining Solution, Bio-Rad, USA) was added into the gels (and incubated for 1 h at room temperature (RT) and revealed with destaining solution (50% H2O, 40% methanol, 10% acetic acid (v/v)) until bands were clearly visible. The gel images were acquired by a Color Image Scanner and analyzed with the ImageJ software.

Protein aggregation ratio

The amount of recombinant protein of the soluble and insoluble fraction for each construct was determined as explained above. After that, the aggregation ratio was calculated by dividing the quantity of the insoluble fraction for each time, replicate and construct by the quantity of protein in the respective soluble fraction.

Purification of protein nanoparticles

Bacterial cultures were processed 3 h post-induction through a combination of mechanical and enzymatic disruption methods. Protease inhibitors (Complete EDTA-free, Roche, Switzerland), phenylmethanesulphonylfluoride (PMSF) and lysozyme were added to the culture at a final concentration of 0.4 mM (Sigma-Aldrich, USA) and 1 µg/ml (Sigma-Aldrich, USA), respectively. After 2 h of incubation at 37 °C and 250 rpm the culture was centrifuged at 6000×g and resuspended in 30 ml of PBS supplemented with protease inhibitors (Complete EDTA-free, Roche, Switzerland). Then, the mixture was ice-jacketed and sonicated for 4 cycles of 1.5 min at 10% amplitude under 0.5 s cycles (Branson Ultrasonic SA, Switzerland). After sonication, the mixture was frozen ON at − 80 °C. The mixture was thawed and Triton X-100 was added (0.4% (v/v)) and incubated for 1 h at RT. After this treatment, the mixture was frozen at − 80 °C for 2 h and then thawed for several cycles until no viable bacterial growth was detected. After that, 125 µl of Nonidet P40 (NP-40) was added and incubated for 1 h at 4 °C. Then, DNA was removed with DNAse at a final concentration of 0.6 µg/ml and 0.6 mM MgSO4 for 1 h at 37 °C. Samples were centrifuged at 15,000×g for 15 min at 4 °C. Pellets containing IBs were washed with 25 ml lysis buffer (50 mM Tris–HCl pH 8, 100 mM NaCl, 1 mM EDTA and Triton X-100 0.5% (v/v)). Finally, pellets were centrifuged at 4 °C for 15 min and 15,000×g and stored at − 80 °C until analysis. The IBs were quantified by western blot using a monoclonal anti-GFP antibody (1:1000, sc-9996; Santa Cruz Biotechnology, USA). All incubations were done under agitation.

Electron microscopy

Field Emission Scanning Electron Microscope (FESEM) was used to visualize the ultrastructural morphology (size and shape) of protein nanoparticles in a nearly native state. For that, protein samples were directly deposited over silicon wafers (Ted Pella, USA), air dried, and observed with a high resolution standard secondary electron detector through a FESEM Merlin (Zeiss, Germany) operating at 2 kV. As a quantitative morphometric measurement, the mean area of IBs for each construct was analyzed as estimator of size with the Image J software. The number of particles was 30, 9, 11, for GFP, J-GFP-F and J/F-GFP, respectively. Using mean area values, the mean diameter of each particles was calculated, and the diameter was also used as an estimator of size. The roundness of IB particles for each construct was also evaluated (n = 50/construct), as an assessment of particle shape.

Total carbohydrate and lipid analysis

The total lipid amount in IBs was determined following a sulfo-phospho-vanillin colorimetric assay. Briefly, 500–1000 µl of each sample were centrifuged at 15,000×g and 4 °C for 15 min. The supernatants were removed, and each pellet was dried by means of a vacuum lyophilizer (SpeedVac™, ThermoFisher) and dissolved in 200 µl of chloroform in a capped glass tube. The chloroform was evaporated at 63 °C in a fume hood. To each tube, 2 ml of 18 M sulfuric acid were added. The samples were then incubated for 10 min in a boiling water bath. After that, each tube was cooled down in ice for 5 min. Five milliliters of phosphoric acid-vanillin reagent (100 ml of 85% (v/v) phosphoric acid and 0.12 g of vanillin, (Sigma-Aldrich, Germany)) were added to the tubes and incubated for 15 min at 37 °C. The tubes were cooled for 15 min in ice and absorbance was measured at 530 nm. For the standard curve, a range of 10-100 μg triolein (dissolved in chloroform) was used. To determine the total content of carbohydrates present in the IBs samples, a phenol–sulfuric assay was performed. A glucose standard was prepared ranging from 0 to 150 nanomols and 150-225 µl of IBs were centrifuged at 15,000×g and 4 °C for 15 min. The supernatant was removed, and each pellet was dried with a vacuum lyophilizer. A total of 150 μl of 18 M sulfuric acid was added, followed immediately by 30 μl of 5% phenol. Samples were incubated at 90 °C for 5 min and then cooled to RT, and absorbance was measured at 490 nm. All assays were performed in triplicate.

Structure modelling

J-GFP-F and J/F-GFP were modelled by iterative threading, as implemented the I-TASSER server [26] without applying any additional restrains. Given the primary sequence of the three constructs the server allows generating possible three-dimensional models by multiple threading alignment. The quality of the models was assessed both by their C-score, which was negative in all cases, and visually. For J-GFP-F, both the first (C-score = − 1.86) and the second (C-score = − 2.36) modes were chosen as these models showed different arrangements of the Jun/Fos domains. For both J-GFP (C-score = − 1.68) and F-GFP (C-score = − 2.87) the first model was chosen.

Molecular dynamics protocol

For all the systems, the free construct was minimized, placed in a cubic box with a water layer of 0.7 nm and Na+ Cl- ions to neutralize the system, and a second minimization was performed. We used AMBER99SB-ILDN [37] force field and Simple Point Charge water. NVT and NPT equilibrations were run for 100 ps, followed by 250 ns NPT production run at 300 °K. The temperature was controlled with a modified Berendsen thermostat [38], the pressure with an isotropic Parrinello-Rahman at 1 bar. The iteration time step was set to 2 fs with the Verlet integrator and LINCS [39] constraint. Periodic boundary conditions were used. All simulations and their analysis were run as implemented in the GROMACS package [40].

FTIR analysis

IB samples were resuspended in a few hundred microliters of 20 mM of sodium phosphate buffer pH 7.4 and 2–3 µl of the suspension at 1 μg/μl were deposited in several drops on a BaF2 window and dried at RT in order to obtain a protein film. In particular, different sample concentrations were analyzed in order to obtain thin (to avoid excessive absorption) and uniform films without cracks [28]. Soluble GFP at 2 mg/ml in 20 mM sodium phosphate buffer pH 7.4, was analyzed in the form of a protein film as for IB samples. The transmission FTIR absorption spectra were then acquired by the Varian 670-IR FTIR spectrometer coupled to the Varian 610-IR infrared microscope (both from Varian, Australia) equipped with a mercury cadmium telluride nitrogen-cooled detector. The following conditions were employed: 2 cm−1 spectral resolution, 25 kHz scan speed, 512 scan co-additions, and triangular apodization. Several areas for each sample were measured to verify the reproducibility of the spectral results. Only absorption spectra with the Amide I band intensity below 0.8 were considered reliable. For H/D exchange, the protein film on the BaF2 window was rehydrated by the deposition of 8 μl of D2O around the dried film. The chamber was then tightly closed by a second window using a flat O-ring and incubated for 5 h at RT to allow H/D exchange [28, 34]. The FTIR spectra of the D2O-rehydrated samples were collected as described above.

Protein spectra were obtained after subtraction of the proper reference spectra strictly collected under the same conditions. The second derivatives [41] were calculated after spectral smoothing by the Savitsky-Golay method. Data collection and analysis were performed using the Resolutions-Pro software (Varian, Australia).

Statistical analysis

All quantitative data are presented as mean values ± standard error of the mean (x̅ ± SEM). Normality of the data was determined by a Shapiro–Wilk test. For IB roundness, the data were normalized with the following formula: sqrt(max(x + 1) - x), where x is IB roundness and max the maximal roundness value in the data. A one-way ANOVA was conducted for all quantitative results except for IB surface and specific fluorescence measures. For the latter, values were compared with an independent sample t-test. Finally, we did Post hoc comparisons using the Tukey HSD test for all data analyzed by the ANOVA method. The level of significance was set at p < 0.05. Measures were done in triplicate, except for IB mean area, diameter and roundness, as indicated above. All statistical analyses were performed using the RStudio Statistical Software (RStudio, Inc., USA).

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Abbreviations

LZ:

Leucine zipper

IB:

Inclusion body

GFP:

Green fluorescent protein

IPTG:

Isopropyl-β-d-thiogalactoside

FESEM:

Field Emission Scanning Electron Microscope

H/D:

Hydrogen/deuterium

References

  1. García-Fruitós E, Sabate R, de Groot NS, Villaverde A, Ventura S. Biological role of bacterial inclusion bodies: a model for amyloid aggregation. FEBS J. 2011;278(14):2419–27. https://doi.org/10.1111/j.1742-4658.2011.08165.x.

    Article  PubMed  CAS  Google Scholar 

  2. Choi JH, Lee SY. Secretory and extracellular production of recombinant proteins using Escherichia coli. Appl Microbiol Biotechnol. 2004;64(5):625–35. https://doi.org/10.1007/s00253-004-1559-9.

    Article  PubMed  CAS  Google Scholar 

  3. Villaverde A, Corchero JL, Seras-Fanzoso J, Garcia-Fruitós E. Functional protein aggregates: just the tip of the iceberg. Nanomedicine. 2015;10(18):2881–91.

    Article  CAS  Google Scholar 

  4. García-Fruitós E, Vázquez E, Díez-Gil C, Corchero JL, Seras-Franzoso J, Ratera I, et al. Bacterial inclusion bodies: making gold from waste. Trends Biotechnol. 2012;30(2):65–70. https://linkinghub.elsevier.com/retrieve/pii/S0167779911001685.

  5. Ventura S, Villaverde A. Protein quality in bacterial inclusion bodies. Trends Biotechnol. 2006;24(4):179–85. https://linkinghub.elsevier.com/retrieve/pii/S0167779906000527.

  6. García-Fruitós E, González-Montalbán N, Morell M, Vera A, Ferraz RM, Arís A, et al. Aggregation as bacterial inclusion bodies does not imply inactivation of enzymes and fluorescent proteins. Microb Cell Fact. 2005;4:1–6.

    Article  CAS  Google Scholar 

  7. Rinas U, Garcia-Fruitós E, Corchero JL, Vázquez E, Seras-Franzoso J, Villaverde A. Bacterial Inclusion Bodies: Discovering Their Better Half. Trends Biochem Sci. 2017;42(9):726–37. https://linkinghub.elsevier.com/retrieve/pii/S0968000417300269.

  8. de Marco A, Ferrer-Miralles N, Garcia-Fruitós E, Mitraki A, Peternel S, Rinas U, et al. Bacterial inclusion bodies are industrially exploitable amyloids. FEMS Microbiol Rev. 2019;43(1):53–72. https://academic.oup.com/femsre/article/43/1/53/5144214.

  9. Pesarrodona M, Jauset T, Díaz-Riascos ZV, Sánchez-Chardi A, Beaulieu M, Seras-Franzoso J, et al. Targeting antitumoral proteins to breast cancer by local administration of functional inclusion bodies. Adv Sci. 2019;6(18):1900849. https://doi.org/10.1002/advs.201900849.

    Article  CAS  Google Scholar 

  10. Roca-Pinilla R, López-Cano A, Saubi C, Garcia-Fruitós E, Arís A. A new generation of recombinant polypeptides combines multiple protein domains for effective antimicrobial activity. Microb Cell Fact. 2020;19(1):122. https://doi.org/10.1186/s12934-020-01380-7.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  11. Gifre-Renom L, Seras-Franzoso J, Rafael D, Andrade F, Cano-Garrido O, Martinez-Trucharte F, et al. The biological potential hidden in inclusion bodies. Pharm. 2020;12(2):157. https://doi.org/10.3390/pharmaceutics12020157.

    Article  CAS  Google Scholar 

  12. Unzueta U, Cespedes MV, Sala R, Alamo P, Sánchez-Chardi A, Pesarrodona M, et al. Release of targeted protein nanoparticles from functional bacterial amyloids: A death star-like approach. J Control Release. 2018;279:29–39. https://linkinghub.elsevier.com/retrieve/pii/S0168365918301780.

  13. Unzueta U, Seras-Franzoso J, Céspedes MV, Saccardo P, Cortés F, Rueda F, Garcia-Fruitós E, Ferrer-Miralles N, Mangues R, Vázquez E, Villaverde A. Engineering tumor cell targeting in nanoscale amyloidal materials. Nanotechnology. 2017;28(1):15102. https://doi.org/10.1088/0957-4484/28/1/015102.

    Article  CAS  Google Scholar 

  14. Seras-Franzoso J, Peternel S, Cano-Garrido O, Villaverde A, García-Fruitós E. Bacterial inclusion body purification. New York: Humana Press; 2015. p. 293–305. https://doi.org/10.1007/978-1-4939-2205-5_16.

    Book  Google Scholar 

  15. Wang X, Zhou B, Hu W, Zhao Q, Lin Z. Formation of active inclusion bodies induced by hydrophobic self-assembling peptide GFIL8. Microb Cell Fact. 2015 \;14(1):88. http://www.microbialcellfactories.com/content/14/1/88.

  16. Wu W, Xing L, Zhou B, Lin Z. Active protein aggregates induced by terminally attached self-assembling peptide ELK16 in Escherichia coli. Microb Cell Fact. 2011;10(1):9. https://doi.org/10.1186/1475-2859-10-9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  17. Slouka C, Kopp J, Spadiut O, Herwig C. Perspectives of inclusion bodies for bio-based products: curse or blessing? Appl Microbiol Biotechnol. 2019;103(3):1143–53. https://doi.org/10.1007/s00253-018-9569-1.

    Article  PubMed  CAS  Google Scholar 

  18. Rinas U, Bailey J. Protein compositional analysis of inclusion bodies produced in recombinant Escherichia coli. Appl Microbiol Biotechnol. 1992. https://doi.org/10.1007/BF00240735.

    Article  PubMed  Google Scholar 

  19. Fahnert B, Lilie H, Neubauer P. Inclusion bodies: formation and utilisation. Berlin: Springer; 2004. p. 93–142. https://doi.org/10.1007/b93995.

    Book  Google Scholar 

  20. Hakoshima T. Leucine Zippers. Encyclopedia of Life Sciences. Chichester: Wiley; 2005. https://doi.org/10.1038/npg.els.0005049.

    Chapter  Google Scholar 

  21. Porte D, Oertel-Buchheit P, Granger-Schnarr M, Schnarr M. Fos Leucine Zipper variants with increased association capacity. J Biol Chem. 1995;270(39):22721–30. https://doi.org/10.1074/jbc.270.39.22721.

    Article  PubMed  CAS  Google Scholar 

  22. Vinson C, Myakishev M, Acharya A, Mir AA, Moll JR, Bonovich M. Classification of human B-ZIP proteins based on dimerization properties. Mol Cell Biol. 2002;22(18):6321–35. https://mcb.asm.org/content/22/18/6321.

  23. Hu C-D, Kerppola TK. Simultaneous visualization of multiple protein interactions in living cells using multicolor fluorescence complementation analysis. Nat Biotechnol. 2003;21(5):539–45. http://www.nature.com/articles/nbt816.

  24. Choi S, Lee SJ, Yeom S, Kim HJ, Rhee YH, Jung H, et al. Controlled localization of functionally active proteins to inclusion bodies using Leucine Zippers. PLoS ONE. 2014;9(6):e97093. https://doi.org/10.1371/journal.pone.0097093.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. Han GH, Seong W, Fu Y, Yoon PK, Kim SK, Yeom SJ, et al. Leucine zipper-mediated targeting of multi-enzyme cascade reactions to inclusion bodies in Escherichia coli for enhanced production of 1-butanol. Metab Eng. 2017;40:41–9. https://doi.org/10.1016/j.ymben.2016.12.012.

    Article  PubMed  CAS  Google Scholar 

  26. Roy A, Kucukural A, Zhang Y. I-TASSER: a unified platform for automated protein structure and function prediction. Nat Protoc. 2010;5(4):725–38. http://www.nature.com/articles/nprot.2010.5.

  27. Ami D, Natalello A, Taylor G, Tonon G, Maria Doglia S. Structural analysis of protein inclusion bodies by Fourier transform infrared microspectroscopy. Biochim Biophys Acta. 2006;1764(4):793–9. https://linkinghub.elsevier.com/retrieve/pii/S1570963905004401.

  28. Natalello A, Doglia SM. Insoluble protein assemblies characterized by Fourier transform infrared spectroscopy. Methods in molecular biology. Clifton: Humana Press Inc; 2015. p. 347–69.

    Google Scholar 

  29. Barth A. Infrared spectroscopy of proteins. Biochim Biophys Acta - Bioenerg. 2007;1767(9):1073–101.

    Article  CAS  Google Scholar 

  30. Scheyhing CH, Meersman F, Ehrmann MA, Heremans K, Vogel RF. Temperature-pressure stability of green fluorescent protein: a Fourier transform infrared spectroscopy study. Biopolymers. 2002;65(4):244–53.

    Article  CAS  Google Scholar 

  31. Tedeschi G, Mangiagalli M, Chmielewska S, Lotti M, Natalello A, Brocca S. Aggregation properties of a disordered protein are tunable by pH and depend on its net charge per residue. Biochim Biophys Acta - Gen Subj. 2017;1861(11):2543–50.

    Article  CAS  Google Scholar 

  32. Singhvi P, Saneja A, Srichandan S, Panda AK. Bacterial inclusion bodies: a treasure trove of bioactive proteins. Trends Biotechnol. 2020;38(5):474–86. https://linkinghub.elsevier.com/retrieve/pii/S0167779919303063.

  33. Cano-Garrido O, Rodríguez-Carmona E, Díez-Gil C, Vázquez E, Elizondo E, Cubarsi R, et al. Supramolecular organization of protein-releasing functional amyloids solved in bacterial inclusion bodies. Acta Biomater. 2013;9(4):6134–42. https://linkinghub.elsevier.com/retrieve/pii/S1742706112005831.

  34. Cano-Garrido O, Sánchez-Chardi A, Parés S, Giró I, Tatkiewicz WI, Ferrer-Miralles N, et al. Functional protein-based nanomaterial produced in microorganisms recognized as safe: a new platform for biotechnology. Acta Biomater. 2016;43:230–9. https://doi.org/10.1016/j.actbio.2016.07.038.

    Article  PubMed  CAS  Google Scholar 

  35. Rueda F, Gasser B, Sánchez-Chardi A, Roldán M, Villegas S, Puxbaum V, et al. Functional inclusion bodies produced in the yeast Pichia pastoris. Microb Cell Fact. 2016;15(1):166. https://doi.org/10.1186/s12934-016-0565-9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  36. Unzueta U, Céspedes MV, Ferrer-Miralles N, Casanova I, Cedano J, Corchero JL, et al. Intracellular CXCR4 + cell targeting with T22-empowered protein-only nanoparticles. Int J Nanomed. 2012;7:4533–44.

    CAS  Google Scholar 

  37. Lindorff-Larsen K, Piana S, Palmo K, Maragakis P, Klepeis JL, Dror RO, et al. Improved side-chain torsion potentials for the Amber ff99SB protein force field. Proteins Struct Funct Bioinforma. 2010. https://doi.org/10.1002/prot.22711.

    Article  Google Scholar 

  38. Bussi G, Donadio D, Parrinello M. Canonical sampling through velocity rescaling. J Chem Phys. 2007;126(1):014101. https://doi.org/10.1063/1.2408420.

    Article  PubMed  CAS  Google Scholar 

  39. Hess B, Bekker H, Berendsen HJC, Fraaije JGEM. LINCS: a linear constraint solver for molecular simulations. J Comput Chem. 1997;18(12):1463–72. https://doi.org/10.1002/%28SICI%291096-987X%28199709%2918%3A12%3C1463%3A%3AAID-JCC4%3E3.0.CO%3B2-H.

    Article  CAS  Google Scholar 

  40. Pronk S, Páll S, Schulz R, Larsson P, Bjelkmar P, Apostolov R, et al. GROMACS 45: a high-throughput and highly parallel open source molecular simulation toolkit. Bioinformatics. 2013;29(7):845–54. https://doi.org/10.1093/bioinformatics/btt055.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  41. Susi H, Byler DM. Resolution-enhanced Fourier transform infrared spectroscopy of enzymes. Methods in enzymology. Cambridge: Academic Press; 1986. p. 290–311.

    Google Scholar 

Download references

Acknowledgements

The authors are indebted to ICTS ‘‘NANBIOSIS”, more specifically the Protein Production Platform of CIBER in Bioengineering, Biomaterials & Nanomedicine (CIBER- BBN)/IBB, at the UAB sePBioEs scientific-technical service (http://www.nanbiosis.es/unit/u1-protein-production-platform-ppp/) and to ‘‘Servei de Microscopia” (UAB). The authors also acknowledge CINECA for the availability of high performance computing resources and support.

Funding

This work was supported by an IRTA Cutting Edge Research project. We are indebted to CERCA Programme (Generalitat de Catalunya) and the European Social Fund for supporting our research. RRP received a PhD fellowship from Secretaria d’Universitats I Recerca del Departament d’Economia i Coneixement de la Generalitat de Catalunya (AGAUR), EGF a post-doctoral fellowship from INIA (DOC-INIA)), and DA a post-doc research fellow (Assegno di Ricerca) from the University of Milano-Bicocca. We acknowledge the CINECA Award N. HP10BTJPER, 2017, for the availability of high performance computing resources and support.

Author information

Authors and Affiliations

Authors

Contributions

RRP carried out most of the experiments, performed the data analysis and prepared a draft of the manuscript. SF performed the structural modeling and computational study. AN and DA carried out the FTIR analysis and ASC was in charge of electron microscopy analysis. EGF and AA directed the study and prepared the final manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Anna Arís or Elena Garcia-Fruitós.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

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.

Figure S1. Overlap between starting model (lighter shades) and final configuration (darker shades) of the J-GFP-F (two models: a and b) and J/F-GFP constructs (one model each), after 250 ns of molecular dynamics simulation. The generated models were minimized, placed in a cubic water box, minimized again, equilibrated and, for each construct, 250 ns of molecular dynamics simulation were performed. Large rearragements of the Jun/Fos domains were observed. Construct domains are color coded as follow: GFP (green), Jun (blue), Fos (red). Figure S2. (a) Representative FESEM images of the isolated IBs for each construct: GFP IBs, J-GFP-F IBs and J/F-GFP IBs. Bars size represent 200 nm. (b) Frequency distribution of IBs ultrastructural morphometry quantification for each construct: size (area (nm2) and diameter (nm)) and shape (roundness (%)). Figure S3. A) FTIR absorption spectra of the protein films. B) FTIR absorption spectra collected after re-hydration of the protein films with D2O for 5 h. GFP and J/F-GFP IBs displayed similar absorption spectra both as film and after re-hydration, while J-GFP-F IBs showed distinct spectral features. As a control, the absorption spectra of the soluble GFP are also shown. Supplementary Table 1. Statistics for the protein aggregation ratio (%) for each construct over time. (a) Aggregation ratio (%) differences between the three constructs and (b) aggregation ratio (%) differences for each construct over time. Different letters mean statistically significant difference (Post-hoc Tukey HSD (THSD) comparisons).

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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Roca-Pinilla, R., Fortuna, S., Natalello, A. et al. Exploring the use of leucine zippers for the generation of a new class of inclusion bodies for pharma and biotechnological applications. Microb Cell Fact 19, 175 (2020). https://doi.org/10.1186/s12934-020-01425-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12934-020-01425-x

Keywords