key: cord-0057262-66sa1xqs
authors: Ekpenyong, Maurice; Asitok, Atim; Antigha, Richard; Ogarekpe, Nkpa; Ekong, Ubong; Asuquo, Marcus; Essien, Joseph; Antai, Sylvester
title: Bioprocess Optimization of Nutritional Parameters for Enhanced Anti-leukemic L-Asparaginase Production by Aspergillus candidus UCCM 00117: A Sequential Statistical Approach
date: 2021-03-09
journal: Int J Pept Res Ther
DOI: 10.1007/s10989-021-10188-x
sha: ba452ddedfaefe899767ea881aedaad241fb2601
doc_id: 57262
cord_uid: 66sa1xqs

Sequential optimization of bioprocess nutritional conditions for production of glutaminase-near-free L-asparaginase by Aspergillus candidus UCCM 00117 was conducted under shake flask laboratory conditions. Catalytic and anti-cancer activities of the poly-peptide were evaluated using standard in vitro biochemical methods. Medium nutrients were selected by one-factor-at-a-time (OFAT) approach while Plackett–Burman design (PBD) screened potential factors for optimization. Path of steepest ascent (PSA) and response surface methodology (RSM) of a Min-Run-Res V fractional factorial of a central composite rotatable design (CCRD) were employed to optimize factor levels towards improved enzyme activity. A multi-objective approach using desirability function generated through predictor importance and weighted coefficient methodology was adopted for optimization. The approach set optimum bioprocess conditions as 49.55 g/L molasses, 64.98% corn steep liquor, 44.23 g/L asparagine, 1.73 g/L potassium, 0.055 g/L manganese and 0.043 g/L chromium (III) ions, at a composite desirability of 0.943 and an L-asparaginase activity of 5216.95U. The Sephadex-200 partially-purified polypeptide had a specific activity of 476.84 U/mg; 0.087U glutaminase activity, 36.46% yield and 20-fold protein purification. Anti-cancer activity potentials of the catalytic poly-peptide were dose-dependent with IC(50) (µg/mL): 4.063 (HL-60), 13.75 (HCT-116), 15.83 (HeLa), 11.68 (MCF-7), 7.61 (HepG-2). The therapeutic enzyme exhibited 15-fold more cytotoxicity to myeloid leukemia cell line than to normal (HEK 238 T) cell. Optimum temperature and pH for activity were within physiological range. However, significant interactions between exposure time and levels of each of temperature and pH made interpretations of residual enzyme activities difficult. The manganese-dependent L-asparaginase from Aspergillu s candidus UCCM 00117 is recommended for further anticancer drug investigations.

Cancers are among the topmost causes of global death and arise from uncontrolled proliferation of certain cells, tissues and organs including breast, lung, blood, liver, pancreas, stomach, colon and cervix, after exposure to certain environmental, nutritional or habitual agents (Sultana et al. 2019; Ervina et al. 2020) . Certain types of cancers can also develop genetically and are heritable within genetically-related populations (Sultana et al. 2019) . Cancerous cells characteristically have a requirement for asparagine having lost their ability to express the gene that encodes asparagine synthetase (Sarquis et al. 2004 ). In human carcinomas therefore, the primary target of anti-cancer compounds is the central sensor for nutrient and energy availability which regulates cell cycle, growth and proliferation (Luengo et al. 2017; Zhang et al. 2018 ). Cancer treatments have largely proceeded through surgery, radiotherapy, immunotherapy and chemotherapy employing diverse agents including inorganic (Jaros et al. 2019; Chen et al. 2020) and organic compounds or a combination of both (Bokas et al. 2020) . Peptide molecules from different biological sources including plants (de Oliveira et al. 2015; Eswaraiah et al. 2020 ), animals (Elrayess et al. 2020 ) and microorganisms (Sereena and Sebastian 2020; Harish et al. 2020 ) have demonstrated precision in target identification and activity, and have therefore received tremendous attention as key components of anti-cancer chemotherapy in recent years.

L-asparaginases (EC 3.5.1.1) are amidohydrolases which primarily catalyze the hydrolytic cleavage of L-asparagine into L-aspartic acid and ammonia oftentimes with secondary catalytic cleavage of L-glutamine to L-glutamic acid and ammonia using glutaminase activity (Avramis 2012; Cachumba et al. 2016 ). The enzyme is frequently included in modern treatment regime for lymphoblastic leukemia and other lymphomas, and in acrylamide reduction protocol in the food industry (Kukurova et al. 2009; Pieters et al. 2011; Shakambari et al. 2019) . L-asparaginase is frequently sourced from plants (Oza et al. 2010; Mohamed et al. 2015) and microorganisms (Sahu et al. 2007; Hassan et al. 2018) . Foremost among microorganisms that produce L-asparaginases include species of Escherichia and Erwinia however; other bacteria including Bacillus velezensis , Enterobacter cloacae (Husain et al. 2016) , Streptomyces rochei (El-Naggar and El-Shweihy 2020) and Geobacillus thermodenitrificans (Nadeem et al. 2019 ) have also been reported. Species of the genera Aspergillus, Penicillium, Fusarium (Sarquis et al. 2004) and Trichosporon (Ashok et al. 2019 ) dominate literature as fungal L-asparaginase producers. Aspergillus species most reported include Aspergillus niger (Vala et al. 2018) , Aspergillus terreus and Aspergillus tamari (Baskar and Sahadevan 2012; Sarquis et al. 2004) .

Applications of bacterial and fungal L-asparaginases, especially from species of Bacillus and Aspergillus (Prabavathy and Phoebe 2020) , in anti-cancer research have been documented. However, the difficulty in obtaining sufficient plant L-asparaginases to meet its high demand and the allergic reactions that attend microbial L-asparaginases during clinical applications owing to secondary glutaminase co-activity are significant drawbacks in anti-cancer chemotherapy (Cecconello et al. 2020) . Remedial attempts at toxicity reduction during chemotherapy have included the use of pegylated forms of the product (Appel et al. 2008 ) as well as sourcing for producing strains with low glutaminase activity (Mahajan et al. 2012) . Currently, there is an upsurge in reports on glutaminase free L-asparaginase (El-Naggar et al. 2018 ) which presupposes that anti-cancer activity of L-asparaginases may not require the associated glutaminase activity (Chan et al. 2014 ). However, complete elimination of glutaminase activity from L-asparaginase may render the medication therapeutically ineffective (Mahajan et al. 2012) . A workable remedy may well be the production of L-asparaginases with low glutaminase activity and optimizing fermentation conditions that enhance their production.

Microbial fermentations have been at the forefront of production of L-asparaginases by solid-state and submerged methods (Shakambari et al. 2019) . However, yields of the product have been low and with poor production economics (Farag et al. 2015) . Sourcing for high-yielding strains, improving them through genetic and metabolic manipulations; optimization of their bioprocess conditions of media, operations and downstream processing; remain reliable approaches towards yield improvement (Baskar and Sahadevan 2012; Yingling and Zhengfang 2013) . Nutritional conditions for microbial fermentations have been reported to account for 30-40% of production economics of bioprocesses; sources of carbon and nitrogen being the most critical. Selection of significant nutrients for inclusion in microbial fermentation media and optimizing their levels are subjects that are vigorously pursued by industries (Shakambari et al. 2019) . A sequence of reported techniques in that regard include one-factor-at-a-time (OFAT) approach that does the primary selection (Long et al. 2018 ) followed by a two-level factorial design like Plackett-Burman design (PBD) that establishes significance of selected factors (Mechmeche et al. 2019; Gururaj et al. 2020) . A third procedure in the optimization process is path of steepest ascent (PSA) which moves significant coefficients from the first-order model of PBD close to their optimum with subsequent resetting of center points for response surface methodology (RSM) (Huang 2019) . The RSM is a mathematical approximation method that employs special experimental design methods in sequential fashion to vary and combine levels for significant interaction effects of variables to bring about an improved condition of target response (Kowalski et al. 2005) .

The present study reports the exploration of the sequence of optimization experiments conducted to enhance the fermentative production of a glutaminase near-free poly-peptide with L-asparaginase activity by a new addition to Aspergillus species in the quest for low-toxicity L-asparaginases. We are not aware of a prior report on L-asparaginase production by a strain of Aspergillus candidus. A description of the spectrum of anti-cancer activity of the therapeutic catalytic poly-peptide and its prospects for safety and stability are also reported.

The producing mold, Aspergillus candidus UCCM 00117, earlier isolated by Asitok and Ekpenyong (2019) and sharing 100% sequence homology with Aspergillus candidus JN-YG-3-2 (MG554228.1), was retrieved from the University of Calabar Collection of Microorganisms, UCCM (www.wfcc.info/ccinf o/colle ction /by_id/652.). The mold was reactivated in glucose-peptone broth for 72 h at room temperature (28 ± 2 °C). Next, the culture was diluted 1000-fold and then plated on Czapek-Dox agar for 72 h at room temperature. Viability and L-asparaginase-producing potential were checked by transferring one agar plug (6 mm) of culture from Czapek-Dox agar onto freshly prepared L-asparagine agar medium composed (g/L) of Na 2 HPO 4 6; KH 2 PO 4 3; MgSO 4 ·7H 2 O 0.2; FeSO 4 ·7H 2 O 0.05, NaCl 0.5; CaCl 2 0.5; KCl 0.5, 1% asparagine and supplemented with 0.25% phenol red. Development of pink coloration between 36 and 60 h was indicative of retained viability and activity for L-asparagine hydrolysis. The reactivated culture was transferred onto potato dextrose agar (PDA) on Petri dish and maintained in a refrigerator until required for further studies.

The preferred extraneous carbon source to drive maximum L-asparaginase activity by the study mold was selected by the traditional one-factor-at-a-time (OFAT) approach. The method involved varying the source of carbon and holding every other parameter constant (Ekpenyong et al. 2017a ). Glucose (GLU), mannose (MAN) and fructose (FRU) as monosaccharides, sucrose (SUC), maltose (MAL) and lactose (LAC) as disaccharides and agro-industrial waste products including molasses (MOL), cassava processing effluent (CPE) and rice processing effluent (RPE) were the carbon substrates screened. Sugars were purchased from Sigma-Aldrich (Sigma-Aldrich Inc., St. Louis, MO, USA) while effluents were obtained from local industries in Lagos and Calabar, Nigeria and stored at -20 °C until required.

The minimal medium for screening contained (g/L) Na 2 HPO 4 6; KH 2 PO 4 3; MgSO 4 ·7H 2 O 0.2; FeSO 4 ·7H 2 O 0.05, NaCl 0.5; CaCl 2 0.5; KCl 0.5 (Asitok and Ekpenyong 2019) and was supplemented with 1 mL trace mineral solution composed ( (Ekpenyong et al. 2017b ). The fully-constituted medium was supplemented with 1% (w/v) L-asparagine (Merck) and pH adjusted to 5.8 using 1 N HCl. Medium was dispensed into 250 mL Erlenmeyer flasks at 20% (v/v) and carbon sources (except the effluents) incorporated at 1% (w/v). The effluents were retrieved from the freezer and allowed to thaw under ambient laboratory conditions and subsequently boiled to facilitate removal of un-dissolved solids. Upon cooling, effluent was centrifuged at 10,000 rpm for 20 min and supernatant dispensed as 20% (v/v) into 250 mL Erlenmeyer flasks without adjusting pH (Nitschke et al. 2004) . A control flask without extraneous carbon where asparagine served as carbon source was also set up. All arrangements were made in triplicates and flasks sterilized by autoclaving at 121 °C for 15 min.

Upon cooling, one agar plug (6 mm) of reactivated and equilibrated mold culture was aseptically added to each flask. Flasks were plugged with sterile cotton wool and incubated at 30 °C on a rotary shaker agitating at 150 rpm for 72 h. Fermentation broth was first filtered with Whatman No.1 filter paper to remove mycelia and then centrifuged at 8,000 rpm for 15 min. Supernatant was twice filtered with 0.45 µm and 0.22 µm membrane filter (Millipore, USA) and L-asparaginase activity quantified for each flask by the Nesslerization protocol of Imada et al. (1973) using Tris-HCl buffer at pH 7.2. One unit of L-asparaginase activity (U) was defined as the amount of enzyme required to liberate one micromole of ammonia by L-asparagine hydrolysis in one minute under the study conditions. A one-way analysis of variance (oneway ANOVA) in GraphPad Prism 8 software (GraphPad, San Diego, CA, USA) was used to compare performances of the carbon sources and significant mean differences were separated by Tukey HSD test using 95% confidence interval.

Extraneous nitrogen sources were also screened by the OFAT approach. Nitrogen sources included ammonium chloride (NH 4 Cl), potassium nitrate (KNO 3 ), ammonium nitrate (NH 4 NO 3 ), bean processing effluent (BPE), corn steep liquor (CSL), peptone (PEP) and tryptone (TYP). Organic and inorganic nitrogen sources were added at 1% (w/v) to minimal medium supplemented with 1% (w/v) L-asparagine and selected carbon substrate. Effluent substrates were prepared as described in the section for carbon screening and used as also described. A control flask without extraneous nitrogen was also prepared. Experimental set-up, sterilization, inoculation, incubation, harvest, L-asparaginase activity assay and statistical analyses were as described under carbon screening.

To select most suitable inoculum size (spore density) for the fermentation, inoculum sizes were screened by the OFAT approach by inoculating 3% (v/v) spore suspension of different inoculum sizes ranging from 10 3 to 10 9 spore-forming units per milliliter (sfu/mL) into minimal medium containing selected extraneous carbon and nitrogen substrates. Working spore concentrations were prepared by the spectrophotometric method as described in Ekpenyong et al. (2020a) and experimental set up, incubations, harvest, determinations of enzyme activity and statistical analyses were as described under carbon screening section.

Metal ions including Na + , K + , Mg 2+ , Ca 2+ , Zn 2+ , Mn 2+ , Ni 2+ , Cu 2+ , Co 2+ , Fe 2+ , Fe 3+ , Pb 2+ , Cr 3+ and Mo 5+ were screened for their effects on L-asparaginase production. All metal salts were purchased as chlorides from Sigma-Aldrich except ferrous and magnesium salts which were purchased as hydrated sulfates from Fisher Scientific (Fisher Scientific UK Ltd, Loughborough, UK). A complete minimal medium with all salts in place as detailed under carbon screening section plus selected carbon, nitrogen and inoculum size served as control. Effect of metal ion on L-asparaginase activity was scored on the basis of L-asparaginase activity quantified in the absence of the respective metal. Higher L-asparaginase activity indicated inhibition in the presence of the metal while lower values indicated enhancement of activity. All metal ions were tested at 1 mM concentration and experimental set-up, incubations, harvest and determinations were as described for carbon screening. Effects of metals on L-asparaginase activity were compared by one-way ANOVA and significant means separated by Dunnett test comparing metal performance to the control at 95% confidence interval.

Significant reduction in L-asparaginase activity selected respective metal for further studies.

The Plackett-Burman design (PBD) was employed to screen a combination of 12 major and minor nutrients obtained from OFAT experiments to develop an appropriate fermentation medium for production of L-asparaginase with maximum enzyme activity, Y. This 2 2 factorial design completely ignores possible two-way factor interactions but focuses on their main effects on the response(s). The variables included in this study design were molasses (MOL), corn steep liquor (CSL), inoculum volume (INV), asparagine (ASP), Na + , K + , Mg 2+ , Zn 2+ , Mn 2+ , Fe 2+ , Co 2+ and Cr 3+ . All 12 factors were tested at 2 levels; high (+ 1) and low (− 1) and the design made with MINITAB 17 (MiniTab Inc, Pennsylvania, USA). The design comprised 25 experimental runs having been augmented with 5 center points to enable estimation of experimental error and test the fit of the first-order model. The factor names and their levels as employed in PBD are presented in Table DT5 of Mendeley Data (Ekpenyong et al. 2020b ). Only predictors with significant (p < 0.05) main effects selected by backward elimination were employed to build the first-order regression model for further studies. The general fitted first-order model developed by least squares method is given as Eq. 1:

where y is the predicted response, β 0 the intercept, β i the coefficients of predictors (significant) x i , i = 1, 2,…k.

The path of steepest ascent (PSA) is an intermediary gradient protocol in RSM optimization experiment employed to move identified significant predictors of a process or system close to the optimum through sequential experimentation (Kowalski et al. 2005) . The goal of the PSA is to identify a new operating region for significant process improvement. It finds a suitable approximation in the first-order model from the PBD which provides the coefficients of significant factors to direct operating region towards PSA subject to a given distance ρ, from the center point.

Since a one unit (using coded units) steepest ascent point away from the center point is desired to maximize the response, y, the coordinates of the significant factor levels were determined as follows:

(1)

The slope of the model was calculated using the largest absolute regression coefficient, β j in the model and denoted Δx j as step size 1 ii. Other step sizes were then calculated using Eq. 2 as follows:

iii. Then, the coded step sizes were converted to natural levels using the relationship between coded and natural levels as given in Eq. 3

where dX i is the change in natural level of predictor variables, Δxi, change in coded level of predictor variables, S, the slope with respect to each predictor and ρ the natural level change (proportion) in the predictor, j, with largest absolute coefficient, β.

Therefore, using the fitted first-order model, the PSA approach moved levels of significant predictors away from the original design center along the path of steepest ascent towards maximum L-asparaginase activity. The levels of significant variables in the PSA experiment that yielded the highest L-asparaginase activity were employed as center points in RSM experiment.

The matrix for the response surface method (RSM) of optimization was obtained through central composite rotatable design (CCRD) using Design Expert 12 software (Stat Ease Inc., Minnesota, USA), with each numerical factor evaluated at 5 levels. Based on the number of significant factors selected by PBD, a special kind of fractional factorial design called Min-Run-Res V type factorial design was adopted for RSM. The design comprised 6 center points for lack-of-fit (LoF) test and 34 non-center points made up of 12 axial (-1.56508 and + 1.56508) and 22 factorial (− 1 and + 1) points; making a total of 40 experimental runs. The predictors were coded as x 1 , x 2 , x 3 , x 4 , x 5 and x 6 corresponding to molasses, corn steep liquor, asparagine, potassium, manganese and chromium. Molasses was composed (%) of dry solids 82.1; Sucrose 50.7; Invert sugar 0.4, Raffinose 1.1; Ash 12.0; Total nitrogen 1.7; Purity quotient of 63.1 and pH of 7.5. Corn steep liquor was composed, in percent weight (wt%), of crude protein 47.2; ash (oxide) 15.3; lactic acid 24.8, phytic acid 7.8, fat 0.33; reducing sugar 2.8, water 32, pH 4.1. The actual levels of the factors were varied based on the steepest ascent steps earlier developed in the PSA

section as follows: X 1 = (44.61)/4.935, X 2 = (67.97)/2.995, X 3 = (40.05)/4.175, X 4 = (1.88)/0.146, X 5 = (50)/5.000, X 6 = (47.42)/4.57. Details are as presented in Table DT9 of the Mendeley Dataset (Ekpenyong et al. 2020b) .

The requirement to report specific activity of the enzyme and the yield coefficient of the bioprocess made the inclusion of biomass concentration (g/L) and total protein (mg) as additional response variables necessary, at this stage of the study. Therefore, three response variables were determined in the RSM experiments namely biomass concentration (Y 1 ), total protein (Y 2 ) and L-asparaginase activity (Y 3 ). Biomass concentration was determined by the conventional oven dry weight technique (Banerjee et al. 1993) while total protein and L-asparaginase activity were determined by the methods of Bradford (1976) and Imada et al. (1973) respectively. All 40 completely randomized experiments were conducted in triplicates and results reported as mean responses. Data was subjected to multiple regression analysis using least squares to build/fit models. The scalar notation of the general quadratic function for RSM is given below;

where y is the response, β 0 , β i , β ii and β ij are the regression coefficients for the constant, linear (x i ), quadratic (x i 2 ) and interaction terms (x i x j ) of the k factors respectively and ɛ the error term of response variable computation. Only significant (p < 0.05) predictors of each response was included in the final model. Adequacy of models was evaluated by examining model diagnostic plots, adjusted and predicted r 2 values, LoF F and p-values and percent coefficient of variation (%CV). Significant model predictors were evaluated by t-test of their coefficients and model curvatures presented as contour (2-D), surface (3-D) plots.

Since the three responses of biomass concentration (Y 1 ), total protein (Y 2 ) and L-asparaginase (Y 3 ) were measured during the same batch fermentation set-up under the same conditions; mono-objective optimization of responses became impracticable. This research therefore undertook a multi-objective optimization approach of the responses using the desirability function provided by response optimizer in Design Expert. The major thrusts of the research were to maximize specific activity of L-asparaginase (U/ mg); calculated as the ratio of L-asparaginase activity to the amount of total protein in the culture supernatant, as well as its yield coefficient (Y p/x ) calculated as ratio of L-asparaginase concentration to biomass concentration and expressed

as g/g. To achieve this therefore, a medium formulated to maximize L-asparaginase activity while minimizing biomass concentration and total protein would be ideal. The desirability function is a mathematical protocol that recommends a set of operational conditions that meet all set goals with the best desirability (Sperb et al. 2018 ). The composite desirability of a multi-response function was defined as the geometric mean of the desirability of individual responses and is given by the expression;

where D is the composite desirability, n is the number of responses, d 1 , d 2 , d 3 and d n are desirability for individual responses. The equation could be weighted to reflect the importance of the responses in the desirability measure where the weight, w i satisfies the function 0 < 1 < w i , with w 1 + w 2 + w 3 + … + w n = 1, as follows:

Triplicate experiments were set up in 1 L Erlenmeyer flasks containing media (20% v/v) formulated according to significant factor settings recommended by the multi-objective numerical optimization. The three responses were measured as earlier described and specific L-asparaginase activity and yield coefficient calculated. Results from confirmation experiments were compared with those obtained in the numerical optimization and a less than 5% difference in L-asparaginase activities was accepted as confirmatory for optimal settings.

The sterile crude L-asparaginase broth was gradually brought to 40% saturation with ammonium sulfate (enzyme grade) and the precipitate removed by centrifugation at 12,000 rpm for 10 min. The resulting supernatant was then brought to 80% saturation with the salt and further purification conducted as described in Asitok and Ekpenyong (2019) . The resulting fraction was dialyzed against 0.05 M Tris-HCl buffer (pH 7.2) at 4 °C for 24 h with 12 changes. Total protein (Bradford 1976 ), L-asparaginase and glutaminase activities (Imada et al. 1973) were determined from the concentrated dialyzed supernatant.

The dialyzed enzyme fraction obtained from the fractionation experiment was loaded into a di-ethyl-amino-ethyl (DEAE)-cellulose 10 cm column (Sigma-Aldrich) and further purification conducted as described by Sahu et al. (2007) . Total protein and activities of L-asparaginase and glutaminase were determined from the pooled concentrated dialyzed fraction as described earlier.

The concentrated fraction of the crude enzyme from the cellulose column purification step was applied onto Sephadex G-200 gel and further purification performed as described in Husain et al. (2016) . The eluted fractions were pooled into a composite sample, dialyzed as earlier described and concentrated with 0.05 M sucrose and protein concentration, L-asparaginase and L-glutaminase activities of the partiallypurified enzyme determined as described earlier.

The in vitro anti-cancer activity of the partially-purified L-asparaginase in this study was investigated using the MTT-based cytotoxicity assay of Skehan et al. (1990) . Human myeloid leukemia (HL-60), colorectal carcinoma (HCT-116), cervical epithelial carcinoma (HeLa), breast cancer (MCF-7) and hepatocellular carcinoma (HepG-2) cell lines were employed to evaluate the therapeutic potentials of the enzyme. Toxicity to normal cells was tested against non-tumor human embryonic cell line (HEK 283 T) to evaluate the selective toxicity of the enzyme. Enzyme dosages were two-fold dilutions from enzyme stock solution ranging from 1.148 to 1174.898 µg/mL. Mycoplasma sterility tests, viability and subsequent cytotoxicity studies using the sulforhodamine B (SRB) assay were as described in Asitok and Ekpenyong (2019) and Ekpenyong et al. (2020a) . Experiments were conducted in triplicates and results, expressed as % cell viability, were presented as means of triplicate determinations. Cell viability was calculated using the relationship below:

where A570T is the absorbance at 570 nm of cells treated with L-asparaginase and A570C is the absorbance at 570 nm

of untreated (control) cells. Data was analyzed by non-linear regression using the sigmoidal four-parameter logistic (4PL) model presented as Eq. 8.

where x is the independent variable and y, the dependent variable; the 4 estimated parameters include a; the minimum value that can be obtained (i.e. what happens at 0 dose), d; the maximum value that can be obtained (i.e. what happens at infinite dose), c; the point of inflection (i.e. the point on the sigmoidal curve half-way between a and d) called IC 50 or EC 50 and b; the Hill's slope of the curve.

Selectivity index (SI) of the enzyme was calculated as a ratio of IC 50 of non-tumor cell line to those of cancer cell lines and data was analyzed by two-way ANOVA.

The temperature and pH effects on L-asparaginase activity earlier reported in Asitok and Ekpenyong (2019) were repeated to confirm the optimum temperature (5-70 °C) and pH (3-12) for maximum activity of the enzyme. Reported values were means of triplicate determinations ± standard deviation. One way analysis of variance (one-way ANOVA) was used to compare data and significant mean differences were separated by Tukey HSD at p < 0.05. Furthermore, the thermal stability of the enzyme was evaluated by first exposing the enzyme in 50 mM Tris-HCl buffer (pH 7.5) to different temperatures including 40, 50, 60, 70 and 80 °C for 0, 15, 30, 45, 60, 75, 90 min. The enzyme was allowed to refold on ice for 15 min and then residual activity determined. Similarly, stability of the enzyme to increasing pH levels was evaluated by incubating enzyme in different buffers of varying pH values (Iboyo et al. 2017 ) ranging from 4.5 to 10.5 at 4 °C for 0, 6, 12, 18, 24 and 30 h. Residual activity of the enzyme was again determined after incubation as described for thermal stability. Data for enzyme stability studies were analyzed by two-way ANOVA in GraphPad Prism 8 to compare significant means and identify possible interactions at p < 0.05.

Inhibition or enhancement of enzyme activity by metal ions was investigated by incubating the enzyme reaction mixture with 1 mM concentration of each metal ion (Na + , K + , Mg 2+ , Ca 2+ , Ba 2+ , Zn 2+ , Ni 2+ , Fe 2+ , Fe 3+ , Cr 3+ , Mn 2+ , Co 2+ , Cu 2+ and Mo 5+ ). Relative activities were calculated using assay without metal as control. One-way analysis of variance was used to compare data for significant mean

differences at p < 0.05. Range finding experiments for significant activity enhancing metal ions were conducted using 0.5, 1.0, 1.5, 2.0, 2.5 and 3.0 mM concentration. All experiments were conducted in triplicates and relative activities calculated and reported as mean value ± standard deviation.

Frequently in microbial fermentations, the carbon substrate that drives cellular metabolism differs from the one that enhances target metabolite synthesis (Abalos et al. 2002) . Among the 9 carbon substrates screened by OFAT, only sucrose, cassava processing effluent and molasses significantly (p < 0.05) enhanced L-asparaginase activity above control (Fig. 1a) (2014) had earlier reported the preference of black strap molasses as carbon source for L-asparaginase production by Acinetobacter radioresistence PS14. Production of L-asparaginase on agro-industrial wastes have been suggested as a sustainable approach towards improved production economics of biotechnological products. The OFAT experiment to select appropriate extraneous nitrogen source for improved L-asparaginase activity using molasses as extraneous carbon substrate revealed that tryptone (1387 ± 62.18 U) and corn steep liquor (1602 ± 93.34 U) were the only nitrogen sources which significantly improved the activity of the enzyme (Fig. 1b) . Corn steep liquor was selected with a slight edge over tryptone (p = 0.0188 < 0.05) as best enhancing nitrogen source while urea significantly (p < 0.0001) inhibited L-asparaginase activity, allowing only an activity of 873.8 ± 90.87 U. Inhibition of L-asparaginase activity by urea had earlier been reported by Li et al. (2018) to proceed through sustained denaturation of thermo-stable L-asparaginase from Pyrococcus yayanosii CH1 expressed in Bacillus subtilis 168. Production of L-asparaginase by Aspergillus niger on corn steep liquor had earlier been reported by Zia et al. (2013) . One of the remarkable things about using agro-industrial wastes as substrates for fermentation is the multiple roles they can play. In the present study it served as nitrogen source in combination with asparagine and carbon source in combination with beet molasses and asparagine.

The results presented as Fig. 1c show that L-asparaginase activity increased with increasing spore density from 740.4 ± 82.35 U at 10 3 spore forming units per milliliter (sfu/mL) to 1742 ± 87.16 U when spore density was 10 7 sfu/mL. Tukey multiple comparisons test of a one-way ANOVA revealed that L-asparaginase activity reported at 10 7 sfu/mL was significantly higher than those at 10 6 sfu/mL (1460 ± 72.01 U) below it and 10 8 sfu/mL (1329 ± 99.98 U) above it. A similar range of inoculum size was reported by Zia et al. (2013) for L-asparaginase production by Aspergillus niger.

Finally, out of 14 metal ions screened by OFAT for their effects on L-asparaginase activity, results, presented as Fig. 1d , show that L-asparaginase activity reduced from 1732 ± 39.76U in the control to 1135 ± 44.45U in the absence of Mn 2+ and 1237 ± 41.00 U in the absence of Cr 3+ suggesting significant (p < 0.0001) requirement for these ions in the fermentation medium. However, L-asparaginase activity significantly increased to 2049 ± 51.83U in the absence of Ni 2+ and a little less in the absence of Fig. 1 One-factor-at-a-time (OFAT) plots to select most enhancing a carbon substrate b nitrogen source c spore density d metal ions for enhanced L-asparaginase production Pb 2+ , Cu 2+ and Fe 3+ suggesting significant inhibition of L-asparaginase production or activity in their presence. This is the first report on L-asparaginase activity enhancement by trivalent chromium. However, among all metals studied by El-Naggar et al. (2018) on their effects on Streptomyces brollosae NEAE-115 L-asparaginase, maximal enzyme activity was observed in the presence of manganese. Dunnett's multiple comparisons test showed that incorporation or non-incorporation of Na + , K + and Mo 5+ ions did not significantly influence L-asparaginase activity. The dataset for these experiments is available at Mendeley Data (Ekpenyong et al. 2020b) . The report by El-Naggar et al. (2018) also supports the inhibitory effects of nickel and trivalent iron on L-asparaginase production.

In this study, PBD was applied to screen 12 factors from the OFAT experiments and the design matrix with the experimental and predicted responses of all combinations of factors is presented in Table 1 . Highest experimental and predicted L-asparaginase activities of 1948.86 and 1974.0U respectively, were obtained with the factor combination in run 22 highlighted in yellow. The full regression model is presented as Eq. 9 and incorporates all 12 predictors; significant and non-significant. Their corresponding coefficient table with t-test result has been published as dataset in Mendeley Data (Ekpenyong et al. 2020b ). where, Ŷ = predicted L-asparaginase activity, MOL = molasses, CSL = corn steep liquor, INV = inoculum volume, ASP = asparagine, Na + = sodium ions, K + = potassium ions, Mg 2+ = magnesium ions, Zn 2+ = zinc ions, Mn 2+ = manganese (II) ions, Co 2+ = cobalt (II) ions, Fe 2+ = Iron (II) ions and Cr 3+ = chromium (III) ions. The model was significant with an r 2 of 0.9747, adjusted r 2 of 0.9493, predicted r 2 of 0.8785, non-significant lack-of-fit test of F (8, 12) = 1.64, p = 0.333 > 0.05. However, to achieve the primary objective of PBD viz screening, the backward elimination method was adopted to reduce the model to include only significant (p < 0.05) predictors as in Eq. 10 below:

The reduced model was also significant (p < 0.0001) at F = 59.73, r 2 of 0.9522, adjusted r 2 of 0.9362, predicted r 2 of 0.9037, a non-significant lack-of-fit test of F (14, 18) = 2.02, p = 0.260 > 0.05. The ANOVA table for the modified model is presented as Table 2 and shows that Mn 2+ and molasses made the most contributions towards enhanced L-asparaginase activity as indicated by their large F-values of 85.22 and 82.95 respectively. A careful examination of the predicted r 2 values reveal that the reduced model gives better prediction of L-asparaginase activity (pred. r 2 = 90.37%) than the full model (pred. r 2 = 87.85%). The diagnostic plots for the second model are given as Fig. 2 . The Pareto chart (Fig. 2a) Eq. 9 as important variables for L-asparaginase production while Fig. 2b reveals their levels of importance through main effects plot. The residual plots (Fig. 2c) confirm the adequacy of the reduced model to explain the events leading to the activity of L-asparaginase observed during fermentation.

Results of the path of steepest ascent (PSA) experiments are reported in Table 3 and Fig. 3 . The table summarizes the steps taken to arrive at the steepest ascent step sizes used to develop the steepest ascent experiments. The primary objective of the PSA was to identify plausible center points for the predictors suggested by the PBD first-order model for use in RSM. Table 3 assigned codes to the six predictors as X 1 , X 2 , X 3 , X 4 , X 5 and X 6 to correspond to molasses, corn steep liquor, asparagine, K + , Mn 2+ and Cr 3+ respectively. In consideration of the predictor coefficients in Eq. 9, the predictors were varied by factors of 4.935 g/L, 2.995%, 4.175 g/L, 0.146 g/L, 5 mg/L and 4.57 mg/L respectively. With these step changes, experiment 7 yielded the highest L-asparaginase activity of 3272 ± 42.64 U beyond which the activity took a descent (Fig. 3) . The path of steepest ascent clearly identified the values of predictors which should form the center points in RSM as (X 1 , X 2 , X 3 , X 4 , X 5 , X c Diagnostic plots to test adequacy of the first-order regression model for significant predictors towards L-asparaginase production Table 3 Path of steepest ascent (PSA) optimization of significant predictors for enhanced L-asparaginase activity during fermentation by Aspergillus candidus UCCM 00117

Exp. Experiment, L-ASPNase act L-asparaginase activity, x 1x 6 model codes for predictors, X 1 -X 6 Actual predictor codes N/B: a Molasses; b Corn steep liquor; c Asparagine; d Potassium ions; e Manganese (II) ions; f Chromium (III) ions; g Predictor levels at center points in PBD; h Original step units of predictors in the PBD matrix; i Coefficients of significant predictors in the first-order model from PBD; j Steepest ascent coded unit step distance; k coded ascent natural unit step distance; l Standard error of means of triplicate determinations of L-asparaginase activity; m steepest ascent step change 

In addition to L-asparaginase activity (Y 3 ), two more responses were measured during the RSM experimentation namely biomass concentration (Y 1 ) and total protein (Y 2 ) to enable the determination of specific activity of the enzyme and yield coefficient of the process. The design matrix, experimental and predicted responses are presented in Table 4 . The table shows that maximum biomass Table 4 Design codes, experimental and predicted values of responses of a surface methodology

x 1 , x 2 , x 3 , x 4 , x 5 and x 6 = model codes for molasses, corn steep liquor, asparagine, potassium ions, manganese (II) ions and chromium (III) ions respectively; eY 1 = Experimental biomass concentration (g/L); pY 1 = Predicted biomass concentration (g/L); eY 2 = experimental total protein (mg); pY 2 = Predicted total protein (mg); eY 3 = experimental L-asparaginase activity (U); pY 3 = Predicted L-asparaginase activity (U); Values of responses are means of triplicate determinations.

x 1 x 2 x 3

x 4

x 5 x 6 eY 1 pY 1 eY 2 pY 2 eY 3 pY 3 concentration, Y 1 of 87.91 g/L (predicted as 87.34 g/L) was obtained in run 5 where molasses, corn steep liquor, asparagine and Cr 3+ were at their high levels (+ 1), with K + and Mn 2+ at their low levels (− 1). This run also accounted for highest total protein, Y 2 of 65.39 mg (predicted as 65.11 mg) indicating strong significant positive correlation (r = 0.8747, p < 0.0001) between the two responses. Highest L-asparaginase activity (Y 3 ) of 4924.78 U (predicted as 4927.74 U) was recorded for run 1 where corn steep liquor, asparagine and K + were at their low levels while molasses, Mn 2+ and Cr 3+ were at their high levels. Conversely, lowest L-asparaginase activity of 2735.98 U (predicted as 2735.3 U) was observed in run 34 where all predictors were held at their center points (0) except asparagine which was supplied at its lowest or axial level (− 1.56508). Analysis of variance (ANOVA) for biomass concentration (Y 1 ) model is presented as a reduced model containing only significant predictors of the response, Y 1 obtained by backward elimination method (Table 5) . From the table, it is observed that the modified model was significant (F = 82.17, adjusted r 2 = 0.9753, predicted r 2 = 0.9269, p < 0.0001) and adequate to explain variations about the data within the experimental region as evaluated by its non-significant lack-of-fit (LoF) test value of F (15, 20) = 0.805, p = 0.6617 > 0.05. The diagnostic plot of actual versus predicted Y 1 plot (Fig. 4a) further confirms the adequacy of the model. Coefficients of the predictors were employed to build the reduced second-order regression function for biomass concentration as follows:

Analysis of variance (ANOVA) for total protein (Y 2 ) model is also presented as a reduced quadratic function that discriminated against non-significant predictor terms through backward elimination ( Table 6 ). Out of the 27 terms required to build the full model, 5 terms with non-significant (p > 0.05) F-value were eliminated. The reduced model was significant at F = 269.35, adjusted r 2 = 0.9934, predicted r 2 = 0.9803, p < 0.0001. The non-significant LoF test value of F (12, 17) = 1.97, p = 0.2352 > 0.05 indicated adequacy of the model. The goodness-of-fit, r 2 of experimental versus predicted Y 2 plot (Fig. 4b) was also in reasonable agreement (11) Y 1 =63.88 + 7.64X 1 + 4.25X 2 + 1.73X 4 − 2.03X 5 + 7.82X 6 + 1.11X 1 X 3 + 1.83X 1 X 4 . 4 Experimental versus predicted response surface methodology plots to evaluate the adequacy of a biomass concentration b total protein c L-asparaginase activity second-order regression models with the model r 2 value of 0.9971, thus validating the adequacy of the model. Model coefficient terms were used to build the second-order polynomial equation for Y 2 as follows:

The ANOVA for L-asparaginase activity regression model (Y 3 ) is presented as a full quadratic function in Table 7 . All 27 predictors of the second-order function of 6 factors were significant for Y 3 and the model's F-value of 649.96 was also significant (p < 0.0001). The difference between adjusted and predicted goodness-of-fit, r 2 of 0.9978 and 0.9909 respectively, was not more than 0.2 suggesting that the model was adequate for predictions (Fig. 4c) . The LoF (12) Y 2 =55.27 + 7.27X 1 + 3.45X 2 − 1.57X 3 + 2.07X 4 + 8.17X 6 + 2.62X 1 X 2 + 1.52X 1 X 4 − 1.13X 1 X 5

Fischer-ratio statistic of F (7, 12) = 0.3965 was not significant (p = 0.8701 > 0.05) suggesting that whatever component of the proposed quadratic model that may have lacked fit was not important enough to alter the predicted L-asparaginase activity value. The model is presented as Eq. 13 below.

The L-asparaginase activity model predictor coefficients inform about the magnitude of contributions of the predictors to the model while the signs indicate the direction of the contribution or effect of the predictor. Equation 13 reveals that asparagine (X 3 ) and Mn 2+ (X 5 ) made the largest linear contribution to predicted L-asparaginase (13) activity. Their contributions were positive suggesting that they enhanced the activity of the enzyme. Corn steep liquor and K + , on the other hand, had the largest negative coefficients in the model suggesting they largely inhibited L-asparaginase activity. The largest positive interaction coefficient is that for corn steep liquor and asparagine (X 2 X 3 ) which enhanced L-asparaginase activity. Contrariwise, the interaction between molasses and corn steep liquor (X 1 X 2 ) and that between K + and Mn 2+ (X 4 X 5 ) had the largest negative coefficients and so largely reduced the activity of the enzyme. All the quadratic terms of the model contributed towards improved L-asparaginase activity with K + and Mn 2+ making the highest contributions to the improvement. This position is the complete opposite for Y 1 and Y 2 models where quadratic terms of their models only reduced the respective responses. Interaction model predictors for all three responses are presented as surface (3-D) plots in Fig. 5 . Figure 5a reveals that highest predicted concentration of fungal biomass would be 378.69 g/L if levels of corn steep liquor and Cr 3+ were adjusted to 0.995 (70.96%) and 0.992 (51.95 mg/L) in natural units respectively, with the remaining 4 predictors held at their center points. Highest interaction term contribution to total protein was also by corn steep liquor and Cr 3+ when their levels were adjusted to 0.984 and 0.996 respectively (Fig. 5b) corresponding to 70.92% and 51.97 mg/L respectively, in natural units. This lends credence to the strong significant positive relationship between biomass 5 Surface plots of most significant two-way factor interactions towards maximum a biomass concentration b total protein and c L-asparaginase activity concentration and total protein (r = 0.8747). Figure 5c is the surface plot for most significant two-way interaction of predictors towards enhanced L-asparaginase activity. The figure shows that L-asparaginase activity will be as high as 4124.6 U if levels of both asparagine and Mn 2+ would be adjusted to 0.988 which would correspond to 44.22 g/L and 54.99 mg/L respectively, in natural units while all other factors are held at their center points.

Biomass concentration (Y 1 ) and total protein (Y 2 ) are by themselves significant predictors of L-asparaginase concentration, and by extension its activity (Y 3 ). The total activity of an enzyme is not nearly as important as its specific activity; calculated as a ratio of total enzyme activity to the amount of total protein. Similarly, the yield coefficient (yield of product on biomass), Yp/x, of a microbial fermentation Contour and surface plots of multi-objective optimal solutions for fermentation responses a A 4-in-1 contour plot of asparagine/ Mn 2+ interaction towards minimum biomass concentration, minimum total protein and maximum L-asparaginase. b-i Surface plot for mini-mum biomass concentration. b-ii Surface plot for minimum total protein. b-iii Surface plot for maximum L-asparaginase activity. c Surface plot for the composite desirability for maximum L-asparaginase activity product is calculated in relation to the concentration of biomass that synthesized it. Since the two goals of specific activity and yield coefficient would only be enhanced at high L-asparaginase activity and concentration respectively, it became imperative to search for predictor levels that would maximize L-asparaginase concentration and its activity, while keeping the denominators at minimum levels.

The goals of the multi-objective optimization were therefore set to minimize both Y 1 and Y 2 while maximizing Y 3 . Secondly, all the factors and responses were held at equal weight of 1 but with variable scaling of the importance of the three responses. A scale of 1indicated least importance while 5 indicated most importance. Results showed that a scaling of importance of 1, 3 and 5 for Y 1 , Y 2 and Y 3 respectively was most appropriate since it produced the highest composite desirability for the study. The desirability of each response was calculated as 0.840 for Y 1 , 0.887 for Y 2 and 1.000 for Y 3 (Fig. 6a) . A composite desirability, D, solved from Eq. 6 was calculated as 0.943 by adjusting levels of predictors to minimize Y 1 and Y 2 while maximizing Y 3 . The optimal solution given by the multi-objective optimization was therefore 0.99998, -0.99999, 0.999253, -0.999508, 0.991215, -0.969278 for X 1 , X 2 , X 3 , X 4 , X 5 , X 6 respectively (Fig. 6b) . These levels corresponded to 49.55 g/L molasses, 64.98% corn steep liquor, 44.23 g/L asparagine, 1.73 g/L K + , 0.055 g/L Mn 2+ and 0.043 g/L Cr 3+ in natural units. These conditions resulted in a biomass concentration of 343.655 g/L, total protein of 219.78 mg and L-asparaginase activity of 5216.945 U. The contour, surface and desirability plots in Fig. 7 are graphical representations of the optimized solutions for the responses described.

The triplicate confirmation experiments set up to verify the optimal settings by the multi-objective function towards maximum L-asparaginase activity produced 344.76 ± 23.74 g/L fungal biomass, 218.53 ± 19.59 mg total protein and an L-asparaginase activity of 5187.65 ± 184.38U.

These results gave a yield coefficient, Yp/x, of 15.05 U/(g/L) and a specific activity of 23.74 U/mg. When these results were compared to the suggestions by the multi-objective function, it was found that none of the corresponding differences was greater than 5% indicating that the optimal conditions suggested by the optimizer were in reasonable agreement with real life situations.

Results of purification, activity and fold determinations of the L-asparaginase obtained by fermentation of Aspergillus candidus UCCM 00117 in the optimized medium is summarized in Table 8 . The table shows that enzyme purity, expressed in terms of fold, increased with each purification step until a final enzyme with ~ 20-fold purity was obtained.

On the contrary, L-asparaginase yield (%) or recovery decreased with every step of the purification process until a final yield of 36.46% of the enzyme, corresponding to 1.45 mg of a final total protein of 3.97 mg was obtained. L-asparaginase yield by Streptomyces rochei subsp. chromatogens NEAE-K reported by El-Naggar and El-Shweihy (2020) from a similar optimization study was 30%, although with lower fold of purification. The L-glutaminase activity of the enzyme also decreased from an initial activity of 82.131U in the sterile broth to a final activity of 0.087U, when L-asparaginase activity was 1891.64U with a specific activity of 476.48 U/mg. Since potent L-asparaginases in anti-cancer drug formulations possess a range of untoward reactions due to glutaminase activity, Aspergillus candidus UCCM 00117L-asparaginase holds promise for very minimal toxicity during therapy. Literature is replete with reports on microbial production and anti-cancer activities of glutaminase-free L-asparaginase (Hussain et al. 2016; El-Naggar et al. 2018; Ashok et al. 2019) . However, research is not very clear on the role of glutaminase activity during L-asparaginase treatment of human cancers. Glutaminase activity of L-asparaginase may not be required for anti-cancer activity in asparaginase 

Results of the anti-cancer activity of the partially-purified L-asparaginase presented in Fig. 8 reveal a dosedependent therapeutic effect where % cell viability decreased with increasing enzyme concentration. The human myeloid leukemia (HL-60) cell line was most . 9 Activity and stability plots of temperature, pH and metal ion effects on L-asparaginase activity. a-i Effect of temperature on L-asparaginase activity. a-ii Effect of temperature on stability of L-asparaginase activity. b-i Effect of pH on L-asparaginase activity.

b-ii Effect of pH on stability of L-asparaginase activity. c-i Effect of metal ions on L-asparaginase activity. c-ii Concentration range finding plot of L-asparaginase activity enhancing metal ions sensitive to the enzyme with 100% cytotoxicity occurring at maximum cytotoxic concentration (MCC) between 18.37 and 36.73 µg/mL after 48 h exposure (Table 9) . Half-maximal inhibitory concentration (IC 50 ) of the metabolite against HL-60 cell line was 4.063 µg/mL while those for HCT-116, HeLa, MCF-7 and HepG2 cell lines were 13.75, 15.83, 11.68 and 7.610 µg/mL respectively. Alrumman et al. (2019) reported IC 50 of 17.02, 14.55 and 11.66 µg/mL for their purified L-asparaginase against HCT-116, MCF-7 and HepG-2 cell lines respectively. Their results are in reasonable agreement with the findings of this study considering the profile likelihood of 9.214-17.99 µg/mL and 7.062-16.54 µg/mL for HCT-116 and MCF-7 respectively (Ekpenyong et al. 2020b ). The four parameter logistic models for all test cell lines were significant (Table 9) including that for non-tumor human embryonic kidney cell line (HEK-238 T). The selectivity index (SI); a measure of discrimination of anti-cancer agent against normal cells; determined for HL-60 cell was ~ 15 suggesting that the therapeutic enzyme was 15-fold more cytotoxic against the myeloid leukemia cell than to normal cell ( Table 9 ). The index was greater than 2 for all tested cell lines indicating low toxicity of the compound (Suffness and Pezzuto 1990) . El-Naggar and El-Shweihy (2020) reported similar success (SI > 2) with Streptomyces rochei L-asparaginase against HepG-2 and MCF-7 cell lines. Anti-cancer substances with IC 50 less than 30 µg/mL have been reported to hold promise for cancer drug development (de Oliveira et al. 2016) . Aspergillus candidus UCCM 00117 L-asparaginase appears to be a potent anti-cancer molecule for the treatment of a number of human carcinomas especially blood and liver cancers.

Results of investigation into optimum temperature and pH for L-asparaginase activity are presented in Figs. 9a and b respectively. Figure 9a -i reveals that optimum temperature of activity for the study enzyme was 40 °C (5402 ± 13.29 U). However, Tukey HSD multiple comparisons analysis of the one-way ANOVA informed that there was no statistically significant difference between activity of the enzyme at 35 °C and at 40 °C. Since the human physiological temperature is 37 °C which lies between 35 and 40 °C, the L-asparaginase in this study is considered suitable for treatment of different human carcinomas. This finding agrees with those of El-Naggar et al. (2018) who reported a temperature optimum of 37 °C for Streptomyces brollosae NEAE-115 L-asparaginase. Results of a two-way ANOVA of the thermal stability of the enzyme indicated significant interaction between duration of exposure and temperature thus making interpretation of residual activity results difficult as shown in Fig. 9a -ii. This implies that the residual activity is a function of both the temperature and duration of exposure. The pH optimum for the enzyme ranged from 7 to 9, accommodating the typical tight range of the human physiological pH of 7.35-7.45, and suggesting once more the suitability of the enzyme for in-vivo therapeutic applications (Fig. 9b-i) .

Once again, pH stability result was difficult to interpret in isolation from exposure time since the interaction effect was significant (Fig. 9b-ii) . Strangely, this is the only report that considers the statistical implication of physicochemical parameters on stability of enzyme activity. Relative activity of the enzyme with regard to metal inhibition or enhancement was also studied and the result is presented in Fig. 9c . Figure 9c -i shows that Ba 2+ , Ni 2+ , Fe 3+ and Cu 2+ were potent inhibitors of L-asparaginase activity while K + , Mg 2+ , Mn 2+ and Co 2+ enhanced it. No significant difference existed between relative activities of the remaining cations and that of the control. Highest enhancement (168.94%) was mediated by Mn 2+ at 1 mM which is in agreement with the findings of El-Naggar et al. (2018) who reported a relative activity of 145.15% and Li et al. (2018) who reported mediation of 190% relative L-asparaginase activity of Pyrococcus yayanosii CH1 L-asparaginase expressed in Bacillus subtilis 168 by Mn 2+ . Maximum enhancement of our L-asparaginase activity by the metal occurred at 2 mM with relative activity of 184.17% (Fig. 9c-ii) suggesting that the assay medium for this L-asparaginase may require some re-constitution by manganese boosting.

Sequential statistical optimization involving OFAT, PBD, PSA and CCRD-RSM coupled with multi-objective optimization selected and optimized levels of molasses, corn steep liquor, asparagine, K + , Mn 2+ and Cr 3+ as significant medium nutrients for enhanced L-asparaginase production by Aspergillus candidus UCCM 00117. Utilization of two renewable substrates (molasses and corn steep liquor) for fermentation may contribute significantly towards a sustainable bioprocess with attractive production economics. A combination of ammonium sulfate fractionation, column chromatography and dialysis recovered substantial amount of the therapeutic enzyme (36.46%) with 20-fold purity. The high specific activity (476.48 U/mg) of the Mn 2+ -dependent L-asparaginase, coupled with low glutaminase activity (0.087 U), low half-maximal inhibitory concentration (< 20 µg/mL) and high selectivity index (> 3) offers potential for further anti-cancer drug investigation and development.

Enugu who facilitated the completion of this study without regard to the constraints of COVID-19 pandemic

Author contributions Conceptualization

Funding The authors did not receive support from any organization for the submitted work

Data Availability The datasets generated during and/or analyzed during the current study are available in the Mendeley Data repository

Utilization of response surface methodology to optimize the culture media for the production of rhamnolipids by Pseudomonas aeruginosa AT10

Production and anticancer activity of an L-Asparaginase from Bacillus licheniformis iolated from the Red Sea

Pharmacokinetic, pharmacodynamic and intracellular effects of PEGasparaginase in newly diagnosed childhood acute lymphoblastic leukemia: results from a single agent window study

Microbes producing L-asparaginase free of glutaminase and urease isolated from extreme locations of Antarctic soil and moss

Production of the anti-leukemic therapeutic enzyme, L-asparaginase, by a brackish sediment strain of Aspergillus candidus

Asparaginases: biochemical pharmacology and modes of drug resistance

Spectrophotometric determination of mycelial biomass

Optimization of culture conditions and bench-scale production of L-asparaginase by submerged fermentation of Aspergillus terreus MTCC 1782

Mechanisms of the antitumor activity of low molecular weight heparins in pancreatic adenocarcinomas

A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding

Current applications and different approaches for microbial L-asparaginase production

Asparaginase: an old drug with new questions

The glutaminase activity of L-asparaginase is not required for anticancer activity against ASNS-negative cells

Glutaminase activity of L-asparaginase contributes to durable preclinical activity against acute lymphoblastic leukemia

Systemic treatment of breast cancer with leptomeningeal metastases using bevacizumab, etoposide and cisplatin (BEEP Regimen) significantly improves overall survival

Study of the cytotoxic activity of Styrax camporum extract and its chemical markers, egonol and homoegonol

Production of L-asparaginase by some local bacterial strains using agro-industrial residues and evaluation of its anti-tumor activity

Response surface modeling and optimization of major medium variables for glycolipopeptide production

Plackett-Burman design and response surface optimization of medium trace nutrients for glycolipopeptide biosurfactant production

Statistical and artificial neural network approaches to modeling and optimization of fermentation conditions for production of a surface/bioactive glyco-lipo-peptide

Data for: bioprocess optimization of nutritional parameters for enhanced anti-leukemic L-asparaginase production from Aspergillus candidus UCCM 00117: a sequential statistical approach

Purification, characterization and immunogenicity assessment of glutaminase-free L-asparagaianase from Streptomyces brollosae NEAE-115

Bioprocess development for L-asparaginase production by Strepromyces rochei, purification and in-vitro efficacy against various human carcinoma cell lines

Bio-selective hormonal breast cancer cytotoxic and antioxidant potencies of Melia azedarach L. wild type leaves

Identification of bioactive compounds in leaf extract of Avicennia alba by GC-MS analysis and evaluation of its invitro anticancer potential against MCF-7 and HeLa cell lines

Optimization of production of anti-tumor L-asparaginase by free and immobilized marine Aspergillus terreus

2020) Production, purification and application of cutinase in enzymatic scouring of cotton fabric isolated from Acinetobacter baumannii AU10

Production enhancement of an anti-coagulant trypsin inhibitor from Oceanimonas sp. BPMS22 and its anti-cancer activity

Purification, characterization and anticancer activity of L-asparaginase produced by marine Aspergillus terreus

Application of Agaricus bisporus industrial wastewater to produce the biomass of Pichia burtonii

Purification and characterization of glutaminase free L-asparaginase from Enterobacter cloacae: in-vitro evaluation of cytotoxic potential against human myeloid leukemia HL-60 cells

Selection of Enterobacter cloacae strain POPE6 for fermentative production of extracellular lipase on palm kernel oil processing effluent

Asparaginase and glutaminase activities of microorganisms

Light-stable polypyridine silver (I) complexes of 1,3,5-triaza-7-phosphaadamantane (PTA) and 1,3,5-triaza-7-phosphaadamantane-7-sulfide (PTA=S): significant antiproliferative activity of representative examples in aqueous media

A modified path of steepest ascent for split-plot experiments

Effect of L-asparaginase on acrylamide mitigation in a fried dough pastry model

Simultaneous cell disruption and semi-quantitative activity assays for high throughput screening of thermostable L-asparaginases

Optimization of fermentation conditions for an Escherichia coli strain engineered using the response surface method to produce a novel therapeutic DNA vaccine for rheumatoid arthritis

Targeting metabolism for cancer therapy

Efficient production of L-asparaginase from Bacillus licheniformis with low glutaminase activity: optimization, scale-up and acrylamide degradation studies

Production of bioactive peptides in tomato seed protein isolate fermented by water kefir culture: optimization of the fermentation conditions

Purification and characterization of asparaginase from Phaseolus vulgaris seeds. Evid Based Complement Altern Med

Enhanced production of glutaminase free L-asparaginase by marine Bacillus velezensis and cytotoxic activity against breast cancer cell lines

Studies on the recombinant production in E. coli and characterization of pharmaceutically important thermostable L-asparaginase from Geobacillus thermodenitrificans

A novel L-asparaginase with low glutaminase coactivity is highly efficacious against both T-and B-cell acute lymphoblastic leukemias in-vivo

Selection of microorganisms for biosurfactant production using agroindustrial wastes

Anticancer properties of highly purified L-asparaginase from Withania somnifera L. against acute lymphoblastic leukemia

L-Asparaginase treatment in acute lymphoblastic leukemia: a focus on Erwinia Asparaginase

Cytotoxic activity of L-asparaginase isolated from endophytic Aspergillus nomius of Justicia adhatoda on A549 cell lines

Partial purification and anti-leukemic activity of L-asparaginase enzyme of the actinomycete strain LA-29 isolated from the estuarine fish, Mugi cephalus (Linn.)

Production of L-asparaginase by filamentous fungi

Evaluation of anticancer and antihemolytic activity of Azurin, a novel bacterial protein from Pseudomonas aeruginosa SSj

L-Asparaginases: a promising biocatalyst for industrial and clinical applications

New colorimetric cytotoxicity assay for anticancer drug screening

Simultaneous production of biosurfactants and lipases from Aspergillus niger and optimization by response surface methodology and desirability functions

Assays related to cancer drug discovery

Complementary and alternative approach to heal cancer: a review

Characterization of L-asparaginase from marine-derived Aspergillus niger AKV-MKBU, its antiproliferative activity and bench scale production using industrial waste

Application of an integrated statistical design for optimization of culture condition for ammonium removal by Nitrosomonas europaea

Anti-cancer drug sensitivity assay with quantitative heterogeneity testing using single-cell Raman spectroscopy

Production of L-asparaginase from Aspergillus niger using agro-industrial by-products in submerged fermentation process

Acknowledgements Authors appreciate the technical assistance by Staff of Research Laboratory, Calabar and Hematology and Oncology

Conflict of interest The authors have no conflicts of interest to declare that are relevant to the content of this article.