Heterologous, Expression, and Characterization of Thermostable Glucoamylase Derived from Aspergillus flavus NSH9 in Pichia pastoris

A novel thermostable glucoamylase cDNA without starch binding domain (SBD) of Aspergillus flavus NSH9 was successfully identified, isolated, and overexpressed in Pichia pastoris GS115. The complete open reading frame of glucoamylase from Aspergillus flavus NSH9 was identified by employing PCR that encodes 493 amino acids lacking in the SBD. The first 17 amino acids were presumed to be a signal peptide. The cDNA was cloned into Pichia pastoris and the highest expression of recombinant glucoamylase (rGA) was observed after 8 days of incubation period with 1% methanol. The molecular weight of the purified rGA was about 78 kDa and exhibited optimum catalytic activity at pH 5.0 and temperature of 70°C. The enzyme was stable at higher temperature with 50% of residual activity observed after 20 min at 90°C and 100°C. Low concentration of metal (Mg++, Fe++, Zn++, Cu++, and Pb++) had positive effect on rGA activity. This rGA has the potential for use and application in the saccharification steps, due to its thermostability, in the starch processing industries.


Introduction
Glucoamylase (1,4--D-glucan glucohydrolase, EC 3.2.1.3, GA) is an exoacting enzyme that yields -D-glucose from the nonreducing ends of starch and related oligo-and polysaccharide chains by hydrolyzing -1,4 and -1,6 linkages [1]. This enzyme is able to completely hydrolyze starch if incubated for extended periods of time and hence called the saccharifying enzyme. Glucoamylase is an important group of enzymes in starch processing in the food industries, as it is used for the production of glucose and fructose syrup from liquefied starch [2]. It is also used in baking, juice, beverage, pharmaceuticals, confectionery, and many fermented foodstuffs industries for commercial production [3]. It is the most important type of industrial enzyme because of its widespread uses together with -amylases and debranching enzymes, in the saccharification of starch to yield soluble sugars, which are widely used by many food industries and in the production of bioethanol [4].
From a structural point of view, glucoamylase enzymes constitute the bulk of family GH15 of glycoside hydrolases [5]. Glucoamylase (GA) can be derived mostly in microorganisms and also in animals and plants. A large number of microbes, including bacteria, yeast, and fungi are capable of producing glucoamylase. Filamentous fungi constitute the major source of all microorganisms. However, the exclusive 2 BioMed Research International productions of glucoamylase in the enzyme industry are achieved by Aspergillus niger [6], Aspergillus oryzae, and Rhizopus oryzae [7].
Fungi often produce more than one form of enzyme that may have different molecular weights, amino acid composition, glycosylation, ability to digest soluble and raw starches, and stability [8]. The majority of fungal glucoamylase is multidomain consisting of a catalytic domain at the N-terminus and a starch binding domain (SBD) at the C-terminus [9], which is separated by an O-glycosylated linker region rich in serine and threonine residues. But with exception, glucoamylases with a single domain structure are also available in A. niger [10], Rhizopus oryzae [11], and A. oryzae [12]. Two different forms of glucoamylase, with and without starch binding domain, are produced from one gene or two different genes. For example, A. niger and A. awamori var. kawachi produce two forms of glucoamylase (GA) from the same gene [13], while A. oryzae produces two GA from two different genes [14].
Due to industrial importance, glucoamylase encoding gene has been cloned, characterized, and heterologously expressed which derived the Aspergillus species, A. niger [15], A. awamori var. kawachi [16], and A. oryzae [12]. El-Abyad and his coworker [17] screened 21 species for amylolytic activity and found that A. flavus was among the most active amylase producer. It is also reported that A. flavus NSH9 was the best amylase producer among fourteen isolates of Aspergillus sp. [18]. Although glucoamylase from A. flavus has been characterized [18,19], little is known about its molecular information. The genetic sequence, molecular information, substrate specificity, and biochemical characteristic of glucoamylase from A. flavus are not yet well known. On the other hand, Pichia pastoris is one of the most used eukaryotic systems for the production of recombinant proteins. It is easy to culture and can reach high cell densities and has strong and regulated promoters and ability to secrete proteins as well as introducing posttranslational modifications [20,21]. Here, we described the isolation, cloning, and sequence analysis of glucoamylase gene (GA) derived from Aspergillus flavus and its expression in Pichia pastoris.

Maintenance and Culturing of Aspergillus Strain.
A newly fungal strain of Aspergillus flavus NSH9 with starch degrading properties was previously isolated from the sago humus and kept as collection at the Department of Molecular Biology, Faculty of Resource Science and Technology, Universiti Malaysia Sarawak, Malaysia. For enzyme induction, the actively growing fungal mycelium was transferred from potato dextrose agar (PDA) plate to a culture medium containing (in g L −1 ) 20 g raw sago starch, 3 g KH 2 PO 4 , 1 g (NH 4 ) 2 SO 4 , 0.5 g MgSO 4 ⋅7H 2 O, and 4 g of yeast extract. The mycelia were collected after 4 days of growth in the culture media at room temperature of 28 ∘ C with shaking at 150 rpm and then subsequently used to extract the total RNA.

Extraction of Total RNA and DNA and First-Strand cDNA
Synthesis. Total RNA was extracted by using TRIzol Reagent according to the protocol as described by Schumann et al. [22] and the genomic DNA was extracted according to the method of Cubero et al. [23]. Total RNA was purified by being treated with DNAse and reverse transcriptions were carried out in 20 L reaction mixture containing the following: 0.1-5 g DNA free RNA as a template, 0.5 g oligo (dt) 18 , 4 L of 5x reaction buffer, 2 L 10 mM dNTP Mix, 20 U Riblock RNase Inhibitor, and 40 U Moloney Murine Leukemia Virus (M-MuLV) reverse transcriptase (Thermo Scientific). The reaction was performed at 37 ∘ C for 60 min and at 70 ∘ C for 5 min. The product was either used immediately in the subsequent PCR analysis or stored at −80 ∘ C.

Isolation and
Subcloning of Glucoamylase cDNA. Glucoamylase gene from Aspergillus flavus NSH9 was isolated by using GAAsp-F and GAAsp-R primers that were designed based on sequences of A. flavus NRRL3357 (glucoamylase precursor putative mRNA, XM 002374964.1) obtained from the National Center for Biotechnology Information (NCBI). The PCR was carried out in 50 L reaction mixture containing the following: 5 L 10x Pol Buffer A (1x final concentration), 5 mM MgCl 2 , 0.2 mM dNTP, 0.1-0.5 M each of the forward and reverse gene specific primers, 0.5 g cDNA template, and 1.25 U Taq DNA polymerase (EURx, Gdansk Poland). The full length of A. flavus NSH9 glucoamylase cDNA and glucoamylase gene were amplified by PCR, using the following oligonucleotide primers: GAAsp F (5 -CGCATGCGGAACAACCTTCTTT-3 ) and GAAsp R (5 -CTACCACGACCCAACAGTTGG-3 ). The PCR reaction was carried out using the following conditions: denaturation at 94 ∘ C for 2 min for one cycle, followed by denaturation at 94 ∘ C for 45 sec, annealing at 55 ∘ C for 45 sec, and elongation at 72 ∘ C for 1.45 min for 35 cycles, and a final elongation was at 72 ∘ C for 5 min. The purified PCR product was then cloned into pGEMT-Easy Vector (Promega) as described by the manufacturer protocol and transformed into E. coli XL1-Blue. Purified plasmid was run on agarose gel and sequenced using T7 and SP6 primers.

Cloning of Glucoamylase cDNA into Yeast Expression
Vector. The pPICZ C (Invitrogen) was selected as expression vector and used for the cloning of glucoamylase cDNA without signal peptide sequences. A putative signal sequence was predicted by the Signal-3L program [24], and glucoamylase cDNA without the signal peptide was amplified using the following list of primers: forward primer (5 -AGCGAATTCAGATTACAAG GACGACGACGAT-AAGCCGTCCTTCCCTATCCAT-3 ) and reverse primer (5 -ACGACTCTAGAGCCCACGACCCAACAGTTGG-3 ). The designed primers contained an EcoR1 site and Xbal site in the forward and reverse primers (bold and underline), respectively. In addition, the forward primer was designed to include FLAG tag containing 24 nucleotides (DYKDDDDK) and the reverse primer lacking in the native stop codon. The primers were designed to be in frame with the C-terminal 6x His tag located in the pPICZ C. PCR was performed as follows: 1 cycle at 94 ∘ C for 2 min followed by 35 cycles of the sequence at 94 ∘ C for 45 sec, 62 ∘ C for 45 sec and 72 ∘ C for 1.45 min, and final extension 72 ∘ C for 5 min. After purification, the PCR product was digested with EcoR1 and Xba1 and ligated into pPICZ C, producing the P. pastoris expression plasmid pPICZ C GA. The recombinant expression plasmid containing glucoamylase cDNA (pPICZ C GA) was then transformed into Escherichia coli XL1-Blue. The insert in pPICZ C GA was confirmed by PCR, restriction enzyme digestion, and sequencing analysis for correct orientation.

Transformation and Selection of Glucoamylase
Transformant. The pPICZ C GA carrying the insert was linearized by SacI and transformed into P. pastoris GS115 (Invitrogen) competent cells. The competent cells were prepared by using the Pichia EasyComp Transformation kit according to the instruction guidelines (Invitrogen). Linearized vector pPICZ C, without insert, was also transformed into P. pastoris GS115 and used as a control. After transformation, Pichia cells were plated onto yeast extract peptone dextrose sorbitol (YPDS) agar plate containing 150 g mL −1 zeocin and incubated for 3 to 8 days at 30 ∘ C. Transformed colonies were confirmed by colony PCR using 5 AOX1 and 3 AOX1 primers and glucoamylase gene specific primers. Putative Mut + colony was selected based on PCR using AXO1 primers. Several Mut + colonies from YPDS agar plate were selected and independently blotted onto Buffered Methanol-complex Medium (BMMY) (containing 1% (w/v) yeast extract, 2% (w/v) peptone, 1.34% (w/v) YNB, 4 × 10 −5 % (w/v) biotin, 0.5% (v/v) methanol, and 100 mM potassium phosphate buffer pH 6.0) agar plate containing 1% (w/v) soluble starch. Glucoamylase producing transformants were identified by the formation of decolorization zone or clear halo around the Pichia colonies after the addition of iodine solution.

Expression of Recombinant Glucoamylase (rGA). The
Mut + GS115 single colony from YPDS agar plates was inoculated into 25 mL Buffered Glycerol-complex Medium ((BMGY) containing 1% (w/v) yeast extract, 2% (w/v) peptone, 1.34% (w/v) YNB, 4 × 10 −5 % (w/v) biotin, 1% (v/v) glycerol, and 100 mM potassium phosphate buffer pH 6.0) in 100 mL conical flask and incubated at 29 ∘ C with shaking of 230 rpm until the culture reached an OD 600 of 2-4 (approximately after 24 hours). After that, cells were harvested by centrifugation at 3000 ×g for 5 min at room temperature, resuspended to an OD 600 of 1.0 in Buffered Methanol-complex Medium ((BMMY) containing 1% (w/v) yeast extract, 2% (w/v) peptone, 1.34% (w/v) YNB, 4 × 10 −5 % (w/v) biotin, 0.5% (v/v) methanol, and 100 mM potassium phosphate buffer pH 6.0) (40 mL) in 250 mL conical flask, and incubated for 10 days at 29 ∘ C with shaking of 230 rpm. Approximately, 1 mL of the culture was withdrawn at every 24-hour intervals till 10-day period of incubation and added with absolute methanol into the culture to make the final concentration of 0.5% (v/v) of methanol. Based on the best production of recombinant glucoamylase, a Pichia cell was scaled up by increasing the methanol concentration to up to 5% (v/v) of final concentration. Pichia pastoris GS115 containing vector without the insert was also expressed as the above procedure and used as a control (Invitrogen 2010).

Electrophoresis of the Recombinant Glucoamylase.
Nondenaturing PAGE and complete SDS-PAGE were both used to check on the glucoamylase enzyme activity and to estimate the molecular weight of the rGA expressed. Electrophoresis was carried out by using a 12% (w/v) acrylamide resolving gel and 6% (w/v) acrylamide stacking gel, as described by Laemmli [25] using the Mini-Protean III electrophoresis system (Bio-Rad, Richmond, CA, USA). The proteins were stained with Coomassie Brilliant Blue R-250. As for amylase activity staining of the recombinant glucoamylase (samples without heating and beta-mercaptoethanol on gel) was performed on a PAGE gel after washing the gel for 30 min with 0.1 M sodium acetate buffer (pH 5.0) containing 0.2% (v/v) Triton X-100 to remove the SDS in gel. After washing the gel was washed again with 0.1 M sodium acetate buffer (pH 5.0) containing 1% (w/v) soluble starch for 30 min. The recombinant glucoamylase activity appeared as clear band on dark blue background after the staining of iodine solution (0.15% (v/v) of iodine and 0.5% (w/v) of potassium iodide). Desired band in the complete SDS-PAGE analysis after the destaining from Coomassie Brilliant Blue R-250 was cut out and subjected to peptide sequencing. The analysis of the protein was done by MALDI-TOF/TOF mass spectrometry using a 5800 Proteomics Analyzer. The spectra generated were evaluated and analyzed to identify the protein of interest using Mascot sequence matching software (Matrix Science) with MSPnr100 Database. This analysis was done by Proteomics International Pty Ltd., Broadway, Nedlands, Western Australia.

Recombinant Glucoamylase Enzyme Activity and Protein
Assay. Glucoamylase activity was determined by transferring a volume of 0.5 mL of the dilute glucoamylase enzyme followed by the addition of 0.5 mL of 1% (w/v) soluble starch solution in 0.1 M sodium acetate buffer (pH 5.0) and incubated at 55 ∘ C for 30 min. The released glucose was measured using 3,5-dinitrosalicylic acid (DNS) reagent with some modification [26] and measured at 540 nm. The enzyme activities were calculated using a calibration curve prepared using D-glucose as standard. One unit of glucoamylase activity was defined as the amount of enzyme that released 1 mol of glucose equivalent per minute from soluble starch under the assay condition. The protein content was determined by the method of Bradford [27] with bovine serum albumin as standard. The specific activity of recombinant glucoamylase was taken as unit of mg −1 protein.

Purification and Characterization of Recombinant
Glucoamylase (rGA) 2.9.1. Purification. Recombinant glucoamylase was purified by using an anti-FLAG M2 affinity gel according to the manufacturer's instruction (Sigma Aldrich).

Effect of pH and
Temperature. The optimum pH for activity was determined by measuring the glucoamylase activity at 55 ∘ C for 30 min using various types of buffers.
2.9.3. pH Stability and Thermostability of rGA. For pH stability, the recombinant glucoamylase was dispersed (1 : 1) into 0.1 M buffer solution pH 3.0 to 9.0 and incubated at 25 ∘ C for 24 h [28] and after that an aliquot was used to determine the remaining activity at optimum pH and temperature. Thermal stability of the recombinant glucoamylase was determined by incubating in 0.1 M sodium acetate buffer pH 5.0, at 60-100 ∘ C for 90 min, and also enzyme was incubated in water at 70 ∘ C to observe the stability at neutral pH. Time course aliquots were withdrawn, cooled down in ice bath, and assayed under standard conditions at optimum pH (pH 5.0) and temperature (70 ∘ C).

Effect of Substrate Concentrations.
Enzyme activity was assayed in the reaction mixture containing different amount of soluble starch (0.25-2.0% (w/v)) at optimum pH 5.0 and temperature at 70 ∘ C with under standard assay conditions. The data were plotted (according to Lineweaver-Burk) to determine max and m values.

Effect of Metal Ions.
In order to investigate the effect of metal ions (Na + , K + , Cu 2+ , Fe 2+ , Pb 2+ , Ca 2+ , Mg 2+ , Cd 2+ , and Zn 2+ ) on the activity of recombinant glucoamylase, an enzyme aliquot was incubated with a metal ion solution (the final concentration of 1 mM and 5 mM and the final volume after addition substrate) for 30 minutes at 37 ∘ C. After that the remaining activity was measured under standard assay condition with optimum pH (pH 5.0) and temperature (70 ∘ C). The relative glucoamylase enzyme activity which was assayed in the absence of metal ions was taken as 1.00.
2.9.6. Effect of Different Substrate. Different substrates such as soluble starch, amylopectin, glycogen, and gelatinized raw sago starch were used to investigate the hydrolyzing capacity of the recombinant glucoamylase produced.
2.9.7. Bioinformatics. Nucleotide sequences of cDNA of glucoamylase were used as the query sequences in the BLAST searches using the BLAST X program. Searches of protein sequences were done in the NCBI (Protein BLAST). The phylogenic tree of glucoamylase from Aspergillus flavus NSH9 was constructed by using MEGA6 software [29]. For this the alignment of amino acid sequences was performed using the CLUSTALW program, and the evolutionary tree was calculated from the manually adjusted alignment with the neighbor-joining clustering method and the bootstrapping procedure. The number of bootstrap trails used was 1000 and implemented in the MEGA6 program. The phylogenic tree of glucoamylase from Aspergillus flavus NSH9 was constructed by using MEGA6 software [29] ( Figure 1). The results of BLAST between A. flavus NSH9 glucoamylase and other glucoamylases showed about 98-99% identical of A. flavus NRRL3357 (XP 002375005.1, putative precursor), A. oryzae 3.042 (EIT73393.1), and A. parasiticus SU1 (KJK64719.1). It was also found to be 70% identical with glucoamylase precursor of A. terreus NH2624 (XP 0012151-58.1). About 57% to 60% identical were also observed from A. awamori (AEI58995.1), A. niger (AIY23067.1), A. ficuum (AAT58037.1) and Neurospora crassa OR74A (EAA27730.2), and A. oryzae (BAA00841.1) and putative glucan 1,4 -glucosidase was observed from A. flavus NRRL3357 (XP 002384-946.1). In the phylogenic tree ( Figure 1) and NCBI protein BLAST, it was found that glucoamylase from A. flavus NSH9 matched with the two types of glucoamylase from A. oryzae and A. flavus NRRL3357. It is assumed that A. flavus NSH9 also produced two glucoamylases, one that is lacking starch binding domain and the other with the starch binding domain.

Identification and Analysis of Glucoamylase
The amino acid sequences of glucoamylase deduced from the nucleotide sequences were compared with those of other glucoamylases (gene accession numbers p69327 and p69328) retrieved from the Universal Protein Resources  Knowledgebase (UniPort) and GeneBank [30] database. Based on the analysis of the data set, it is observed that the two amino acids, aspartic acid at the 203 position and glutamic acids at 206 position, were responsible for catalytic activity of glucoamylase of A. flavus NSH9 and tryptophan at 147 position was responsible for the binding.

Transformation, Expression of Glucoamylase cDNA,
and Purification of rGA. Linearized recombinant plasmid, pPICZ C GA, was transformed into P. pastoris GS115 by Pichia EasyComp Transformation kit (Invitrogen). This method produces low transformation efficiency with only 8 colonies observed after 4 days of incubation on YPDS agar plate. Five recombinant colonies were selected for colony PCR using AXO1 primers. The amplification results showed that positive recombinants produced two close bands for Mut + in GS115 (2.1 kb and 2.2 kb fragments) using 5 AXO1 and 3 AXO1 primers. Meanwhile the control yeast transformed with pPICZ C produced a 600 bp band. It was observed that recombinant GS115 (Mut + ) produced recombinant glucoamylase on BMMY plates after 4 days of growth and also produced decolourization zone/clear halo around the Pichia colonies with iodine solution compared to control ( Figure 2).
Heterologous expression of glucoamylase was achieved under the transcriptional control of the AOX1 promoter in P. pastoris, and the expression was induced by the supplement of methanol. After eight days of induction, GS115 strain had the highest glucoamylase activity of 8.24 U mL −1 at 1% (v/v) methanol (Figure 3), and the recombinant glucoamylase expression level was 0.145 mg mL −1 in the culture supernatant. At higher concentration of methanol, glucoamylase expression level was detected to decrease, and no expression was observed when induced with 5% (v/v) of methanol.

Molecular Weight Determination by SDS-PAGE.
The molecular weight of the rGA was estimated to be approximately 78 kDa when analyzed on SDS-PAGE (Figure 4). This value was higher than the 55.1 kDa calculated from the deduced amino acid sequences of rGA. The mature rGA had 514 amino acids which was calculated from the Kex2 signal cleavage site to 6x his tag or stop codon of the recombinant plasmid pPICZ C GA. Based on the results  of molecular weight in the SDS-PAGE and the calculated weight of rGA, the expressed recombinant glucoamylase in the culture was monomeric and might be glycosylated. The rGA of Aspergillus flavus NSH9 expressed in Pichia pastoris was confirmed by four methods via its ability to decolorize or form halo around the Pichia colonies after flooding with iodine solution on BMMY plate (Figure 2), zymogram analysis on SDS-PAGE (Figure 4), enzyme assay, and MALDI-TOF/TOF mass spectrometry analysis of peptide after digestion of protein by trypsin (data not shown).

Effects of pH, Temperature, and Kinetics.
The pH optimum of rGA for glucoamylase activity was pH 5.0 (Figure 5(a)), and rGA was highly stable over a broad pH range; about 80% of the residual activity was retained after being incubated at pH 3.0 to 9.0 and at 25 ∘ C for 24 hours  ( Figure 5(c)). The optimum temperature of rGA was 70 ∘ C, and the value was significantly higher from all values (Figure 5(b)). The rGA was thermally more stable at 70 ∘ C, when incubated at neutral pH (water), compared to pH 5.0 with 0.1 M sodium acetate buffer ( Figure 5(d)). Approximately, 90% rGA residual activity was observed after 45 min at 70 ∘ C at neutral pH and about 40% at pH 5.0 at the same condition. About 90% rGA residual activity was observed after 15 min of incubation at 60 ∘ C to 90 ∘ C, but at 100 ∘ C the rGA residual activity was about 59% after 15 min (Table 1). About 50% rGA residual activity was observed after incubation of 20 min at 90 ∘ C and 100 ∘ C ( Table 1). The values of Michaelis-Menten constant, m and maximum velocity, and max of the recombinant glucoamylase were determined using Lineweaver-Burk plot and recorded as 5.84 mg mL −1 and 153.85 U mg −1 protein, respectively, for soluble starch.

Effects of Metal Ions and Different
Substrates. Metal ions such as Na + , K + , Ca 2+ and Cd 2+ , ions had no significant effect on the recombinant glucoamylase activity at concentration of 1 mM, but higher concentration of 5 mM Na + had negative effect whilst Ca 2+ had positive effect on enzyme activity ( Table 2). The recombinant glucoamylase enzyme was significantly stimulated by Cu 2+ , Fe 2+ , Pb 2+ , Zn 2+ , and Mg 2+ , at 1 mM concentration of ions, but significantly inhibited by Cu 2+ , Pb 2+ , and Zn 2+ at 5 mM concentration ( Table 2). The results revealed that K + and Cd 2+ had no effect on enzyme activity at different concentrations, and rGA activity was significantly stimulated by only Mg 2+ ion at different concentration. In this work, it was observed that low concentration of metal ions had a positive effect on enzyme activity in all ions tested. The relative activity of the recombinant glucoamylase was observed on 1% (w/v) of different starches at optimum temperature (70 ∘ C) and pH (5.0). The recombinant glucoamylases produced were more active on the breakdown of gelatinized sago starch as compared to soluble starch, amylopectin, and glycogen (Table 3).

Discussion
The glucoamylase gene from Aspergillus flavus NSH9 consisted of two introns and the size of these introns (<100 bp) is typical for most fungal introns [31]. The splice sites of these introns correspond to the fungal consensus sequences, GT -AG boundaries, and internal consensus for lariat formation [32]. From protein BLAST (NCBI), it was determined that this glucoamylase cDNA belonged to the glycoside hydrolases 15 family excluding the starch binding domain (SBD). The single domain structure of glucoamylase (glucoamylase without starch binding domain) was also observed in many fungi such as in A. niger [10], Rhizopus oryzae [11], and A. oryzae [12], which was the similar to the glucoamylase from A. flavus NSH9 in this study. Two different forms of glucoamylase, with and without starch binding domain are produced from one gene or two different genes, depending on mRNA splicing [15], site of transcription initiation, posttranslational processing, and limited proteolysis [4]. For example, A. niger  and A. awamori var. kawachi produced two GA with and without SBD from the same gene due to proteolytic cleavage of SBD from large GA1 [13], but for A. oryzae from two different genes [14]. According to Hata et al. [14], they purified two glucoamylases of A. oryzae from solid-state culture (S-GA) and from submerged culture (L-GA). They concluded that the two glucoamylases have different molecular weight; and all of the enzymatic characteristics of the two glucoamylases were similar except the thermal stability and the raw starch digestibility (L-GA could digest the raw starch). The S-GA has a higher amount of carbohydrate than L-GA, and they assumed that the higher glycosylation of S-GA increased its thermal stability and pH stability, similar finding also observed in this study. Basically in the industry, they used the glucoamylase having the starch binding domain like L-GA due to digestibility of raw starch. In this study the glucoamylase that discussed was similar to S-GA from A. oryzae having low raw starch digestibility [14]. In the expression study, increasing the concentration of methanol above 1% had negative effect on expression, may be at higher concentrations of methanol causes stress to the cells, decreases cell viability, and facilitates cell lysis and proteolytic degradation [33]. Heterologous expression of glucoamylase in this study was comparatively higher compared to glucoamylases expressed in Clostridium thermosaccharolyticum [34], Sulfolobus solfataricus [35], and Thermomyces lanuginosus [36]. Our finding was in parallel with previous reported studies with molecular weights of GA in fungi to be within a range of 25-112 kDa [9]. Many researchers reported that the molecular weight of recombinant glucoamylase was higher than the estimated molecular weight and it depended on the glycosylation site within the sequence, which supported this study [21,[37][38][39]. The study did not measure the degree of glycosylation of the recombinant glucoamylase and deglycosylation where this will give a clearer picture of the recombinant enzyme properties and become one of the limitations of this study.
It has been reported that fungal glucoamylases are active at acidic pH [9], and the maximum catalytic activity at pH 5.0 determined in our reinforces many other studies in A. niger (pH 4.8) [40], A. tubingensis (4.5) [37], and F. solani (pH 4.5) [41]. The optimum catalytic activity of the recombinant glucoamylase at pH 4.0 to 5.0 expressed in P. pastoris was observed by many researchers [21,38]. A higher optimum temperature of recombinant glucoamylase activity reported in this study was similar to the recombinant glucoamylase from A. tubingensis expressed in S. cerevisiae [37] and from Rhizomucor pusillus expressed in Pichia pastoris [38]. The rGA from A. flavus NSH9 had a higher optimum temperature compared to recombinant glucoamylase from Chaetomium thermophilum [21] and A. awamori [28]. Meanwhile, pH stability of rGA reported in this study was almost similar to a number of researchers [9,18,19,39,42]. Glucoamylase that shows high stability over a wide pH range for relatively long period time indicates its suitability for industrial applications [9,19].
As for the production of glucose in the starch processing industry, it involves two stage processing steps: the liquefaction step at a high temperature (95-105 ∘ C) by thermostable alpha amylase and saccharification step by glucoamylase with additional debranching enzymes [43] are normally performed in the industry. Currently, glucoamylases derived from Aspergillus sp. are commonly used in starch industry for the saccharification process [39,43]. But the glucoamylase from fungal sources are active at acidic condition and also less thermostable, and the industrial saccharification takes a long time to produce the desired yield [4], for this it also required pH and temperature adjustments from liquefaction to saccharification steps [39]. The reported recombinant glucoamylase (rGA) from A. flavus NSH9 can be used directly after liquefaction process for a complete hydrolysis of the dextrin into -D-glucose, due to its thermostability. Many researchers reported that rapid denaturation of glucoamylase occurred at temperatures above 60 ∘ C from F. solani [41], and glucoamylase derived from A. citrinus [44] was very unstable at temperatures above 55 ∘ C. The recombinant glucoamylase from A. flavus NSH9 has lower m value compared to Acremonium sp. [45] and R. oryzae [46], respectively, but higher than F. solani [41] and A. niger NCIM [47] for soluble starch. Lower m value for starch may be due to a higher number of interactions between the active site of the enzyme and the substrate molecule, resulting in an increased affinity of the enzyme towards the substrate [47]. Our results also observed that low concentration of metal ions exerts positive effect on recombinant enzyme activity. Our results supports the findings reported by Bhatti et al. [41] for some metal such as Mg 2+ and Ca 2+ and inhibitory effect of Cu 2+ and Zn 2+ at higher concentration as reported by Koç and Metin [19]. There are many methods, which will improve the stability of rGA from Aspergillus flavus NSH9, such as immobilization, modification, and protein engineering [48]. Immobilization of enzyme enhances their properties for efficient utilization in industrial process [48]. Magnetically (magnetic nanoparticles, MNPs) immobilized enzymes have advantages for commercial application due to their ease of separation and reusability [48].

Conclusion
In this work, a novel thermostable glucoamylase cDNA without the substrate binding domain (SBD) of Aspergillus flavus NSH9 has been successfully isolated, characterized, and overexpressed in Pichia pastoris GS115 producing biologically active recombinant glucoamylase (rGA). These finding opens up for the possibility of upscaling and mass production of the recombinant thermostable glucoamylase in an easier manner. Subsequently, this rGA has the potential for use and application in the saccharification steps, due to its thermostability in the starch processing industries.