Optimizing Granulocyte Colony-Stimulating Factor Transcript for Enhanced Expression in Escherichia coli

The human granulocyte colony-stimulating factor (G-CSF) is a hematopoietic growth factor used to prevent and treat neutropenia. G-CSF stimulates the bone marrow to produce infection-fighting granulocytes. Food and Drug Administration of the United States approved G-CSF in 1991 and its PEGylated version in 2002 as a prophylactic and therapeutic measure against neutropenia. Recombinant human G-CSF is produced in surrogate host Escherichia coli and is PEGylated at N-terminal. Besides neutropenia, G-CSF is also used in bone marrow transplantation for the mobilization and maturation of peripheral blood stem cells. Considering the requirement of producing G-CSF therapeutic in large quantities, construct designing for high expression is critical for the biopharmaceutical and industrial application. Earlier studies have employed approaches such as codon optimization, use of strong promoters, employment of protein tags, secretion signals, optimization of protein folding, etc., for increasing expression and yield of therapeutic proteins. In this study, it was observed that mRNA transcribed from the native human cDNA of G-CSF and the codon-optimized variant leads to low protein expression in E. coli. To understand the underlying reasons, the mRNA secondary structure of the 5′ end of the G-CSF transcript was analyzed. This analysis revealed the presence of stable secondary structures at the 5′ end of the G-CSF transcript, arising from the native human gene and even from the codon-optimized sequence. These secondary structures were disrupted through translationally silent mutations within the first 24 nucleotides of the transcript without affecting the protein sequence. Interestingly, through this approach, the G-CSF protein expression was increased 60 folds as compared to native G-CSF construct. We believe that these findings create a roadmap for optimization of G-CSF transcript for enhanced expression in E. coli and could be employed to increase the expression of other therapeutic proteins.

Introduction

Human granulocyte colony-stimulating factor (G-CSF) is a hematopoietic cytokine that plays a critical role in the stimulation, proliferation, mobilization, maturation and activation of granulocytes, including neutrophils (Anderlini et al., 1996; Welte et al., 1996). Neutrophils are one of the most abundant cell types amongst the leukocytes and thus play a critical role in the clearance of invading pathogens and modulate immune system homeostasis (Malech et al., 2014). A decrease in circulating neutrophils is called neutropenia, and it often leads to systemic infections and associated problems (Lyman and Kuderer, 2002). Due to its ability to induce proliferation and activation of neutrophils, G-CSF is clinically used to prevent and treat cancer chemotherapy-induced neutropenia (Crawford et al., 2004; Lyman, 2006). Over the years, G-CSF and its PEGylated forms have become the most valuable therapeutic proteins (Dale, 2002).

For clinical interventions, a recombinant form of G-CSF produced in Escherichia coli is utilized. The recombinant G-CSF differs from the native human G-CSF through the addition of N-terminal methionine for facilitating the expression in the surrogate host, and it lacks inherent glycosylation at Threonine-133 as observed in humans. However, the recombinant G-CSF possesses similar biological activity as that of the native G-CSF (Frampton et al., 1994). G-CSF is a 18.8 kDa protein consisting of one free cysteine and two disulfide bonds. Clinically used PEGylated form of G-CSF is also derived from recombinant G-CSF produced in E. coli. G-CSF’s recommended dosage regimen is 230 μg/m²/day for 2 weeks (Frampton et al., 1994), while 6 mg of PEGylated G-CSF is administered once per cycle of chemotherapy as a prophylactic and therapeutic measure against neutropenia (Rifkin et al., 2010). These dosages suggest that recombinant G-CSF needs to be produced in large quantities as a recombinant protein in E. coli. Given the requirement of large quantities of G-CSF therapeutic protein, construct designing for high expression of G-CSF is critical for managing production costs. Several approaches are employed to produce this protein in higher quantities, including codon optimization, use of strong promoters, employment of protein tags, secretion signals, etc. (Jin et al., 2011; Do et al., 2014; Zhou et al., 2016). In this study, I have employed mRNA transcript engineering for increasing the G-CSF expression levels in the E. coli expression system. This approach has resulted in a notable increase in the expression of functional G-CSF.

Materials and Methods

Materials

The details of all chemicals, molecular biology reagents, G-CSF standard, protein purification resins, columns, etc., were reported in Kumari et al. (2020). Protein concentrators were procured from Merck Millipore. RNA isolation and DNase treatment kits were procured from Qiagen. cDNA synthesis kit was procured from Bio-Rad Laboratories, Inc. Highest analytical grade reagents were used in the current study. AKTA Pure M (25M1) chromatographic system was utilized for protein purification, and analysis was performed using Unicorn 7.3 software. All cell culture reagents were from Gibco, Thermo Fisher Scientific.

G-CSF Transcript Engineering to Increase the Protein Expression in E. coli

The human G-CSF cDNA, codon-optimized gene sequence and cloning details are as provided earlier (Kumari et al., 2020). Briefly, the native human G-CSF cDNA sequence was PCR amplified using forward primer SD13 and reverse primer SD14 and cloned in pET 23a. The E. coli codon-optimized G-CSF sequence was also cloned at NdeI/HindIII sites of the expression vector pET 23a (Novagen) using forward primer SD1 and reverse primer SD2 (primer sequences are provided in Table 1). Translationally silent mutations were incorporated in the codon-optimized gene sequence to engineer the G-CSF transcript’s 5′ region. The incorporation of translationally silent mutations was confirmed by DNA sequencing. The constructs were transformed in E. coli BL21 (DE3) from Novagen.

TABLE 1

Table 1. Details of primers used in this study.

Purification of G-CSF

The G-CSF constructs used in the study were cloned in pET 23a which is a T7 RNA polymerase inducible promoter based expression vector and were transformed in the BL21 (DE3) strain of E. coli. Isopropyl β-D-1-thiogalactopyranoside (IPTG) was used to induce protein expression. The bacterial pellet was disrupted by sonication. After centrifugation, the supernatant (soluble fraction) and the pellet (insoluble fraction) were analyzed using 12.5% polyacrylamide gel electrophoresis (SDS-PAGE) and stained with Coomassie Brilliant Blue. G-CSF was purified using the method reported earlier in Kumari et al. (2020). Briefly, the G-CSF protein was expressed as inclusion bodies, which were solubilized using 2 M urea. The solubilized G-CSF protein was subjected to refolding and purified using cation exchange chromatography (CEC).

RT PCR for Different Constructs for Analyzing mRNA Expression Levels in E. coli

The RNA was extracted after IPTG induction using the RNeasy Kit from Qiagen (Cat No. 74104). The RNA extracted was also treated with RNase-Free-DNase from Qiagen (Cat. No. 79245). cDNA from different constructs were synthesized using the iScript^TM cDNA Synthesis Kit from Bio-Rad. Different primer sets covering full, 5′ and 3′ regions were used to check mRNA expression levels.

Analytical Characterization of G-CSF

Purified G-CSF was analyzed using 15% SDS-PAGE and stained with Coomassie Brilliant Blue (CBB). Western blot analysis was performed to confirm the identity of the purified G-CSF protein. After running purified recombinant G-CSF protein and the commercially available product on sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE), the gel was electro-transferred onto the PVDF membrane. The membrane was blocked with 2% milk powder or 1% BSA in 1X PBS (phosphate-buffered saline) for 1–2 h. The blocking solution was removed, and the blots were incubated for 1 h with a monoclonal antibody (mAb) against human G-CSF at 1:1000 dilution in 1X PBS-T. The intact molecular mass of G-CSF and the control was determined using liquid chromatography-mass spectrometry (LC-MS) Agilent 6550 system. Circular dichroism (CD) spectrometry was performed, as mentioned earlier on the Jasco J-815 spectropolarimeter (Kumari et al., 2020).

In vitro Biological Activity of G-CSF

The in vitro biological activity of G-CSF was determined by measurement of cell proliferation assay using the M-NFS-60 cell line (ATCC CRL- 1838). The cells were treated with various concentrations of WHO G-CSF standard and E8 G-CSF for 48 h and metabolic activity was assessed using standard XTT assay protocol. The cells were grown in RPMI-1640 medium containing 10% fetal bovine serum, 1X penicillin and streptomycin, 0.05 mM β-mercaptoethanol and 62 ng/ml of human recombinant macrophage colony-stimulating factor (M-CSF) as recommended by ATCC. M-NFS-60 cells (35,000 cells/well) were seeded into a 96-well flat-bottom plate containing 2% FBS and serum-starved for 24 h. Varying concentrations of GCSF and G-CSF standard (0.01, 0.1, 1, 10, 100, 1,000 pg/mL) were added to each well. G-CSF formulation buffer was used as a control. Each concentration was set up in triplicate. After 48 h of incubation, XTT was added to each well, and the cells were kept for incubation at 37°C for 4 h. The optical density of the 96 well plate was measured at 490 nm and at 650 nm to normalize the blank reading using an ELISA reader. Data were plotted using GraphPad Prism 6 for Windows, Version 6.05. The standard error of the mean (SEM) was calculated.

Results

Evaluation of Transcript Levels and Protein Expression From Native Human G-CSF and Codon-Optimized Variant of G-CSF

This study was initiated to engineer a novel construct that could lead to an enhanced expression of the recombinant G-CSF in E. coli. Toward this, the cDNA sequence of human G-CSF was retrieved from NCBI GenBank using accession number M13008.1. This sequence was codon-optimized and synthesized (via GenScript, United States) for expression in E. coli. Both the native and the codon-optimized sequences were PCR amplified and cloned in pET-23a. The constructs were transformed in E. coli BL21 (DE3) strain to analyze the protein expression level. Protein expression from both the constructs was analyzed using SDS-PAGE analysis. It is important to note that the expression of recombinant G-CSF using native human sequence was extremely low, and even the codon-optimized sequence results in low expression (Figure 1A). Having seen low protein expression levels, mRNA levels of the constructs with native human G-CSF gene and codon-optimized gene were analyzed using RT-PCR analysis. This analysis suggested that transcripts from both the constructs are expressed sufficiently to comparable levels (Figure 1B).

FIGURE 1

Figure 1. Low expression of G-CSF protein from constructs with native human cDNA and codon-optimized sequence despite sufficient transcription. (A) Depicts the protein expression profile from native human G-CSF and codon optimized G-CSF constructs. The constructs were overexpressed in E. coli. After sonication of IPTG induced transformed E. coli, the supernatant and pellet were analyzed on SDS-PAGE and stained with Coomassie Brilliant Blue. Hu indicates native human G-CSF, CO for codon optimized G-CSF. Control is commercially available G-CSF, which was used as positive control, the arrow head shows the G-CSF protein, S-supernatant, P-pellet and M-marker. (B) Depicts the RT-PCR profile of native human G-CSF and codon optimized G-CSF constructs. Primers (SD11 and SD2) targeting the 3′ end of the transcript were used for RT-PCR. The expected product size is 139bp.

mRNA Structural Analysis and G-CSF Transcript Engineering

Stable secondary structures at the 5′ end of the transcript are known to modulate the mRNA’s translation efficiency. Thus we aimed at studying the presence of secondary structures at the 5′ end. To analyze the structural architecture of the 5′ region of the G-CSF transcript, computational analysis of the mRNA sequence was performed. RNAfold web server was used (Lorenz et al., 2011) to analyze the secondary structures at the 5′ end of the G-CSF transcript. This analysis suggested the formation of hairpins and highly stable secondary structures at the 5′ end of the mRNA transcript (Figures 2A–C), raising the possibility that these complex structures might be hindering the translation of the transcript. To further understand the secondary structures’ stability, the thermodynamic ensemble’s free energy was analyzed for the G-CSF transcript’s first 24 nucleotides. In agreement with the low expression of recombinant G-CSF, this analysis suggested a high GC content at 5′ prime end of the transcript and negative value of free energy, both of which predict the presence of highly structured stable RNA structures that could affect the translation efficiency (Table 2). These data indicate that the structural features of mRNA transcribed from the native G-CSF cDNA and its codon-optimized version may be hindering the translation efficiency.

FIGURE 2

Figure 2. Depiction of stable secondary structure along with minimum free energy (MFE). (A) Native human G-CSF sequence; Hu (B) Codon optimized sequence; CO and (C) Engineered construct; E8. The free energy of the thermodynamic ensemble is also given below each drawing. The MFE structure is colored by base-pairing probabilities and for unpaired regions the color denotes the probability of being unpaired as depicted in the scale bar. The ViennaRNA Web Service (RNAfold web server) was used to predict the MFE structure. (D) Depicts the 5′ sequence alignment of native human G-CSF, codon optimized and engineered construct E8.

TABLE 2

Table 2. mRNA sequence of native, codon optimized and engineered G-CSF transcripts, their GC content and free energy of the thermodynamic ensemble.

Next, we utilized the mRNA’s secondary structure and envisaged translationally silent mutations (Figure 2D) to disrupt the secondary structures. The mRNA encoded by these sequences have free energy of the thermodynamic ensemble in a range of −2.46 kcal/mol to −7.88 kcal/mol, and GC content range from 45 to 75% (Table 2). Furthermore, the analysis of secondary structures of designed sequences using the RNAfold web server suggested that the mutagenesis would lead to instability in the secondary structures in the 5′ of the G-CSF transcript and decrease the base-paring probability. The free energy of the thermodynamic ensemble of the highly structured RNA was more negative, and it was reduced upon the engineering of the G-CSF transcript (Figure 2).

G-CSF Protein Expression and Purification

We selected the engineered 8 variant (referred to as E8 in this manuscript) with eight translationally silent modifications for further analyses (Figure 2D). The translationally silent mutations were introduced using PCR primers with the desired substitutions at the specific positions in the 5′ coding sequence. This engineered construct was transformed in BL21 (DE3) strain of E. coli to analyze protein expression levels. To investigate the effect of these translationally silent mutations on G-CSF transcription, RT-PCR analysis was performed. We utilized different primer set targeting the full-length transcript, or the 5′ end or the 3′ end of the transcript (Figure 3A). No significant difference in the transcript expression of the G-CSF transcript between the engineered construct and the codon-optimized construct was observed (Figure 3B).

FIGURE 3

Figure 3. Increase in protein expression of G-CSF through mRNA engineering without affecting transcription. (A) Depicts the position of different primer sets used in RT PCR. Three different primer sets were used covering full length of the G-CSF transcript; Full (SD1 & SD2 for CO and SD8 & SD2 for E8), 5′ prime end; 5′ (SD1 & SD12 for CO and SD8 & SD12 for E8) and 3′ prime end; 3′ (SD11 & SD2 for both CO and E8) of the transcript. (B) Depicts the RT-PCR profile of G-CSF codon-optimized (CO) and engineered construct (E8). Expected RT-PCR product size is 525, 134, and 139 bp respectively. (C) Depicts the protein expression profile of native G-CSF, codon-optimized G-CSF and engineered G-CSF construct, E8. The constructs were overexpressed in E. coli. After sonication of IPTG induced E. coli, the supernatant and pellet were analyzed on SDS-PAGE and stained with Coomassie Brilliant Blue. Lane identity are given at top; Hu-Native human cDNA sequence, CO-codon optimized G-CSF, E8-engineered G-CSF construct. Commercially available G-CSF was used as control. M- marker, S- supernatant and P- pellet. (D) Cation exchange chromatography profile of the G-CSF. Protein was eluted by using 1 M Tris-Cl. Parameters such as absorbance at 280 nm and conductance have been represented with blue and brown lines, respectively. (E) Bar graph depicting the yield of CEC purified G-CSF from construct using native human cDNA sequence (Hu), construct using codon-optimized gene for G-CSF expression (CO) and engineered construct (E8). Data are Mean ± SEM from three biological experiments. Students’ t-test was performed for establishing the statistical significance of the results. ****indicates a p value of <0.0001 and *** indicates p value of 0.001.

The E8 variant, along with the native G-CSF and codon-optimized constructs, were compared for the protein expression levels to determine the effect of translationally silent mutations on the translational efficiency. Interestingly, a substantial increase in the protein expression from the modified transcript was observed. The enhanced G-CSF protein expression was observed in the E8 variant (Figure 3C). To quantitate the increase in translational efficiency, G-CSF was expressed from these constructs, followed by protein purification and quantitation of purified protein. The G-CSF protein resided in the inclusion bodies and was purified using urea followed by two-step refolding and cation exchange chromatography (Figure 3D). The purified protein was quantitated using Thermo Scientific NanoDrop 2000. Interestingly, engineered construct E8 showed a 30–60 fold increase in CEC purified protein as compared to the native human G-CSF sequence (Figure 3E).

Physicochemical Characterization of E8 G-CSF

The purified protein was analyzed using SDS-PAGE (Figure 4A). Western blot analysis was performed to assess the quality and identity of G-CSF protein (Figure 4B). Commercially available G-CSF as control was also included in these gels, and data from both the analyses confirmed the purity and identity of purified G-CSF protein from E8 construct. Circular dichroism (CD) spectroscopy was performed for probing the secondary structure of G-CSF and control. G-CSF possesses a predominantly helical structure, wherein four-alpha-helices are connected through unstructured regions. We observed that the far-UV CD spectra of E8 G-CSF and control G-CSF were practically identical (Figure 4C). The spectra had the same shape, indicating alpha-helices’ predominance in both the proteins. The superimposed spectra from control G-CSF and the E8 G-CSF is suggestive of identical secondary structures. The intact molecular mass of E8 G-CSF and the control G-CSF was determined using liquid chromatography-mass spectrometry (LC-MS) and was found to be similar (Figures 4D,E).

FIGURE 4

Figure 4. Characterization of the E8 G-CSF variant. (A) Depicts the SDS-PAGE profile of eluted E8 G-CSF protein and the commercially available G-CSF as control, exhibiting the general purity of the eluted G-CSF protein. (B) Shows the western blot analysis of G-CSF proteins 1, 2, and 4 μg of protein sample were used. (C) Shows a comparison of far UV circular dichroic spectra of G-CSF and commercial available G-CSF. (D) Depicts the LC/MS prolife of control protein. (E) Depicts the LC/MS prolife of E8 G-CSF protein. Control/C stands for commercial available G-CSF protein.

In vitro Biological Activity of G-CSF

Next, we also analyzed the in vitro biological activity of G-CSF protein using a G-CSF responsive M-NFS-60 cell line. The treatment of M-NFS-60 cells to increasing concentration of G-CSF induces a dose-dependent proliferation of M-NFS-60 cells. In this study, the WHO G-CSF standard obtained from National Institute for Biological Standards and Control (NIBSC) was used as a reference. The results demonstrate that E8 possesses similar cell proliferative activity/potential as that of the G-CSF standard (Figure 5). Thus, the results highlight that the free energy value of the first 24 nucleotide thermodynamic ensemble is crucial in predicting translational efficiency. This study suggested that such an approach could also be employed to increase the expression of other therapeutic proteins, wherein even mRNA engineering of the first 24 nucleotides in transcript could result in a substantial increase in protein expression. These findings delineate a simplified and robust strategy that could increase the protein expression of other therapeutic proteins and proteins required for research applications.

FIGURE 5

Figure 5. In vitro biological activity of G-CSF using M-NFS-60 cell line. 0.035 million M-NFS-60 cells were seeded into a 96-well flat-bottom plate containing 2% FBS and serum-starved for 24 h. Varying concentrations of GCSF and G-CSF standard (0.01, 0.1, 1, 10, 100, 1,000 pg/mL) were added to each well. 48 h after providing the cytokine, cellular proliferation was analyzed using XTT assay. G-CSF WHO standard was used as a positive control. Data are Mean ± SEM.

Discussion

G-CSF is an important biotherapeutic protein. Given its market value, it needs to be produced in large quantities. In this study, it was observed that native human cDNA sequence for G-CSF and its codon-optimized variant are transcribed with reasonable efficiency in E. coli through T7 promoter. However, the protein expression for both the two transcripts was low. Analysis of the secondary structures in the 5′ end of the transcript revealed the presence of stable secondary structures with high base pairing probability. Disruption of stable secondary structures at the 5′ end of the transcript resulted in a significant increase in the protein expression. Thus, this study has outlined a strategy for optimizing transcript sequence to improve the protein expression of the clinically/industrially important proteins.

High expression of therapeutic proteins is always desired. Several studies have earlier utilized different methodologies for improving the yield of G-CSF protein. Most of the studies have focused on modification of media composition, IPTG use (Vanz et al., 2008; Boubeva et al., 2012), ethanol utilization (Mishra et al., 2020), optimization of fermentation (Kim et al., 2014), or subsequent purification strategy/methodology (Wang et al., 2008; Li et al., 2012; Vemula et al., 2015) for improving the yields of G-CSF protein. Codon optimization of G-CSF for improved expression has been extensively used for expression in yeast (Maity et al., 2016) and E. coli (Gomes et al., 2012). Literature also suggests that the optimization of AT-content of codons immediately downstream of the initiation codon for high-level G-CSF expression. The culture conditions were also optimized wherein 1% glucose was used as a supplement, and E. coli BL21 (DE3) PLysS strain was utilized. The productivity was increased 1.5 folds (Krishna Rao et al., 2008).

However, the genetic architecture of the gene encoding G-CSF or its transcript has been largely ignored. Importantly, mRNA secondary structure regulates protein expression by modulating the transcript’s functional half-life (Mauger et al., 2019). Literature is indicative of an important role of mRNA secondary structure optimization in modulating translation efficiency and thus the protein expression. One of the classical studies by Kozak delineated the influence of mRNA secondary structures on the initiation of translation by eukaryotic ribosomes. In this study, the incorporation of artificial secondary structures into the 5′ non-coding region of a chimeric mRNA encoding preproinsulin suggested that the 40s ribosomal subunits could melt hairpin structures with Gibbs energy in the range of −30kcal/mol, but secondary structures in the range of −50 kcal/mol resist melting and inhibit translation up to 85–95% (Kozak, 1984, 1986).

In the light of this literature, this study was initiated with the objective of engineering G-CSF transcript for enhanced recombinant G-CSF expression in E. coli. It was observed that the native human and E. coli codon-optimized sequences lead to low to moderate expression of G-CSF. Thus, a computational analysis of the 5′ end of G-CSF transcript from the native human sequence and the codon-optimized gene was carried out. We specifically focused on the RNA structure of the first 24 nucleotides rather than the structural prediction of entire mRNA since secondary structures in the 5′ end of the transcript play a more prominent role in regulating translation initiation. This analysis revealed the presence of stable secondary structures at the 5′ end of the transcript. These stable structures were driven by higher GC content at the 5′ end of the transcript. We utilized computational biology for rational sequence designing of several 5′ sequences to disrupt the stable secondary sequences to enable efficient translation initiation from the transcript. Toward this, we utilized analysis of GC content and free energy levels. The codon-optimized transcript was modified using translationally silent mutations. These translationally silent mutations were explicitly aimed at reducing the mRNA secondary structure complexity by replacing critical GC-rich regions with AT-rich regions. The percentage codon usage suggesting the frequency of the codon used per 100 codon information was also considered while incorporating changes. Emphasis was given on reducing the transcript structural complexity by decreasing the base pair probability in the paired region and increasing the probability of being unpaired in the unpaired region of transcript structure, irrespective of the Codon Adaptation Index (CAI) of the modified sequence (Sharp and Li, 1987). The rationale for modifying the 5′ end of the transcript was that the RNA secondary structures after initiation codon could act as a roadblock to impact the ribosome scanning during translation and impact translational efficiency negatively. The physicochemical characterization and immunoassay suggested that the G-CSF protein resulting through modification of transcript was identical to the commercially available G-CSF. Furthermore, in vitro biological activity was also studied and was comparable with WHO standard of G-CSF protein. Thus the transcript engineering involving translational silent mutations results in enhancing the protein expression and biologically active protein.

This study has suggested that the protein expression from the construct having native human cDNA for G-CSF and the codon-optimized gene for recombinant G-CSF is low to moderate due to the formation of stable secondary structures at the 5′ end of the G-CSF transcript. Translationally silent mutagenesis within the first 24 bp in the G-CSF construct could dramatically disrupt the stable secondary structures to increase translational efficiency and G-CSF yield. This study suggests that such an approach could also be employed to increase the yield of other therapeutic proteins, wherein even engineering of the first few nucleotides in the transcript could result in a substantial increase in protein yield.

Conclusion

Thus study suggests that the protein yield from the construct with native human cDNA G-CSF and codon-optimized sequence for recombinant G-CSF is low due to stable secondary structures at the 5′ end of the transcript. Translationally silent mutagenesis within the first 24 nucleotides of the G-CSF construct was employed to disrupt the stable secondary structures. This led to a dramatic increase in translational efficiency and enhanced protein expression. Interestingly, the free energy value of the first 24 bp thermodynamic ensemble is crucial in predicting translational efficiency. This study suggested that such an approach could also be employed to increase the yield of other therapeutic proteins and proteins involved for research purposes. Importantly, even engineering of the first few nucleotides in the transcript could substantially increase protein yield.

Data Availability Statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.

Author Contributions

SD conceived the idea, performed the experiments, interpreted the data, and wrote the manuscript.

Funding

SD acknowledges the funding from CSIR Fast-Track Translational Project MLP 011.

Conflict of Interest

The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The author gratefully acknowledges support from CSIR-IMTECH for providing all in-house infrastructure and facilities. The author also acknowledges technical help from Mr. Akshit V. Jain.

References

Anderlini, P., Przepiorka, D., Champlin, R., and Korbling, M. (1996). Biologic and clinical effects of granulocyte colony-stimulating factor in normal individuals. Blood 88, 2819–2825.

Google Scholar

Boubeva, R., Reichert, C., Handrick, R., Muller, C., Hannemann, J., and Borcharda, G. (2012). New expression method and characterization of recombinant human granulocyte colony stimulating factor in a stable protein formulation. Chimia 66, 281–285. doi: 10.2533/chimia.2012.281

PubMed Abstract | CrossRef Full Text | Google Scholar

Crawford, J., Dale, D. C., and Lyman, G. H. (2004). Chemotherapy-induced neutropenia: risks, consequences, and new directions for its management. Cancer 100, 228–237. doi: 10.1002/cncr.11882

PubMed Abstract | CrossRef Full Text | Google Scholar

Dale, D. C. (2002). Colony-stimulating factors for the management of neutropenia in cancer patients. Drugs 62(Suppl. 1), 1–15. doi: 10.2165/00003495-200262001-00001

PubMed Abstract | CrossRef Full Text | Google Scholar

Do, B. H., Ryu, H. B., Hoang, P., Koo, B. K., and Choe, H. (2014). Soluble prokaryotic overexpression and purification of bioactive human granulocyte colony-stimulating factor by maltose binding protein and protein disulfide isomerase. PLoS One 9:e89906. doi: 10.1371/journal.pone.0089906

PubMed Abstract | CrossRef Full Text | Google Scholar

Frampton, J. E., Lee, C. R., and Faulds, D. (1994). Filgrastim. A review of its pharmacological properties and therapeutic efficacy in neutropenia. Drugs 48, 731–760. doi: 10.2165/00003495-199448050-00007

PubMed Abstract | CrossRef Full Text | Google Scholar

Gomes, F. R., Maluenda, A. C., Tapias, J. O., Oliveira, F. L., Sa-Rocha, L. C., Carvalho, E., et al. (2012). Expression of recombinant human mutant granulocyte colony stimulating factor (Nartograstim) in Escherichia coli. World J. Microbiol. Biotechnol. 28, 2593–2600. doi: 10.1007/s11274-012-1068-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, H., Cantin, G. T., Maki, S., Chew, L. C., Resnick, S. M., Ngai, J., et al. (2011). Soluble periplasmic production of human granulocyte colony-stimulating factor (G-CSF) in Pseudomonas fluorescens. Prot. Express. Purific. 78, 69–77. doi: 10.1016/j.pep.2011.03.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, C. K., Choi, J. H., Lee, S. B., Lee, S. M., and Oh, J. W. (2014). Expression and purification of recombinant human granulocyte colony-stimulating factor in fed-batch culture of Escherichia coli. Appl. Biochem. Biotechnol. 172, 2425–2435. doi: 10.1007/s12010-013-0708-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Kozak, M. (1984). Point mutations close to the AUG initiator codon affect the efficiency of translation of rat preproinsulin in vivo. Nature 308, 241–246. doi: 10.1038/308241a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Kozak, M. (1986). Point mutations define a sequence flanking the AUG initiator codon that modulates translation by eukaryotic ribosomes. Cell 44, 283–292. doi: 10.1016/0092-8674(86)90762-2

CrossRef Full Text | Google Scholar

Krishna Rao, D. V., Rao, J. V., Narasu, M. L., and Bhujanga Rao, A. K. (2008). Optimization of the AT-content of codons immediately downstream of the initiation codon and evaluation of culture conditions for high-level expression of recombinant human G-CSF in Escherichia coli. Mol. Biotechnol. 38, 221–232. doi: 10.1007/s12033-007-9018-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumari, M., Sahni, G., and Datta, S. (2020). Development of Site-Specific PEGylated Granulocyte Colony Stimulating Factor With Prolonged Biological Activity. Front. Bioengine. Biotechnol. 8:572077. doi: 10.3389/fbioe.2020.572077

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, M., Fan, H., Liu, J., Wang, M., Wang, L., and Wang, C. (2012). High pH solubilization and chromatography-based renaturation and purification of recombinant human granulocyte colony-stimulating factor from inclusion bodies. Appl. Biochem. Biotechnol. 166, 1264–1274. doi: 10.1007/s12010-011-9512-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Lorenz, R., Bernhart, S. H., Honer, Zu Siederdissen, C., Tafer, H., Flamm, C., et al. (2011). ViennaRNA Package 2.0. Algorit. Mol. Biol. AMB 6:26. doi: 10.1186/1748-7188-6-26

PubMed Abstract | CrossRef Full Text | Google Scholar

Lyman, G. H. (2006). Risks and consequences of chemotherapy-induced neutropenia. Clin. Cornerst. 8(Suppl. 5), S12–S18. doi: 10.1016/s1098-3597(06)80054-2

CrossRef Full Text | Google Scholar

Lyman, G. H., and Kuderer, N. M. (2002). Filgrastim in patients with neutropenia: potential effects on quality of life. Drugs 62(Suppl. 1), 65–78. doi: 10.2165/00003495-200262001-00005

PubMed Abstract | CrossRef Full Text | Google Scholar

Maity, N., Thawani, A., Sharma, A., Gautam, A., Mishra, S., and Sahai, V. (2016). Expression and Control of Codon-Optimized Granulocyte Colony-Stimulating Factor in Pichia pastoris. Appl. Biochem. Biotechnol. 178, 159–172. doi: 10.1007/s12010-015-1865-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Malech, H. L., Deleo, F. R., and Quinn, M. T. (2014). The role of neutrophils in the immune system: an overview. Methods Mol. Biol. 1124, 3–10. doi: 10.1007/978-1-62703-845-4_1

CrossRef Full Text | Google Scholar

Mauger, D. M., Cabral, B. J., Presnyak, V., Su, S. V., Reid, D. W., Goodman, B., et al. (2019). mRNA structure regulates protein expression through changes in functional half-life. Proc. Natl. Acad. Sci. U S A. 116, 24075–24083. doi: 10.1073/pnas.1908052116

PubMed Abstract | CrossRef Full Text | Google Scholar

Mishra, B., Murthy, G., Sahoo, B., Uhm, S. J., and Gupta, M. K. (2020). Combinatorial ethanol treatment increases the overall productivity of recombinant hG-CSF in E. coli: a comparative study. Appl. Microbiol. Biotechnol. 104, 9135–9145. doi: 10.1007/s00253-020-10899-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Rifkin, R., Spitzer, G., Orloff, G., Mandanas, R., McGaughey, D., Zhan, F., et al. (2010). Pegfilgrastim appears equivalent to daily dosing of filgrastim to treat neutropenia after autologous peripheral blood stem cell transplantation in patients with non-Hodgkin lymphoma. Clin. Lymphoma Myeloma Leuk 10, 186–191. doi: 10.3816/CLML.2010.n.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharp, P. M., and Li, W. H. (1987). The codon Adaptation Index–a measure of directional synonymous codon usage bias, and its potential applications. Nucleic Acids Res. 15, 1281–1295. doi: 10.1093/nar/15.3.1281

PubMed Abstract | CrossRef Full Text | Google Scholar

Vanz, A. L., Renard, G., Palma, M. S., Chies, J. M., Dalmora, S. L., Basso, L. A., et al. (2008). Human granulocyte colony stimulating factor (hG-CSF): cloning, overexpression, purification and characterization. Microbial Cell Factor. 7:13. doi: 10.1186/1475-2859-7-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Vemula, S., Thunuguntla, R., Dedaniya, A., Kokkiligadda, S., Palle, C., and Ronda, S. R. (2015). Improved Production and Characterization of Recombinant Human Granulocyte Colony Stimulating Factor from E. coli under Optimized Downstream Processes. Prot. Express. Purific. 108, 62–72. doi: 10.1016/j.pep.2015.01.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, C., Wang, L., and Geng, X. (2008). High recovery refolding of rhG-CSF from Escherichia coli, using urea gradient size exclusion chromatography. Biotechnol. Prog. 24, 209–213. doi: 10.1021/bp070263y

PubMed Abstract | CrossRef Full Text | Google Scholar

Welte, K., Gabrilove, J., Bronchud, M. H., Platzer, E., and Morstyn, G. (1996). Filgrastim (r-metHuG-CSF): the first 10 years. Blood 88, 1907–1929.

Google Scholar

Zhou, Z., Dang, Y., Zhou, M., Li, L., Yu, C. H., Fu, J., et al. (2016). Codon usage is an important determinant of gene expression levels largely through its effects on transcription. Proc. Natl. Acad. Sci. U S A. 113, E6117–E6125. doi: 10.1073/pnas.1606724113

PubMed Abstract | CrossRef Full Text | Google Scholar