Fungal carbonatogenesis process mediates zinc and chromium removal via statistically optimized carbonic anhydrase enzyme

Introduction With rapid elevation in population, urbanization and industrialization, the environment is exposed to uncontrolled discharge of effluents filled with broad-spectrum toxicity, persistence and long-distance transmission anthropogenic compounds, among them heavy metals. That put our ecosystem on the verge or at a stake of drastic ecological deterioration, which eventually adversely influence on public health. Therefore, this study employed marine fungal strain Rhodotorula sp. MZ312369 for Zn2+ and Cr6+ remediation using the promising calcium carbonate (CaCO3) bioprecipitation technique, for the first time. Results Initially, Plackett–Burman design followed by central composite design were applied to optimize carbonic anhydrase enzyme (CA), which succeeded in enhancing its activity to 154 U/mL with 1.8-fold increase comparing to the basal conditions. The potentiality of our biofactory in remediating Zn2+ (50 ppm) and Cr6+ (400 ppm) was monitored through dynamic study of several parameters including microbial count, CA activity, CaCO3 weight, pH fluctuation, changing the soluble concentrations of Ca2+ along with Zn2+ and Cr6+. The results revealed that 9.23 × 107 ± 2.1 × 106 CFU/mL and 10.88 × 107 ± 2.5 × 106 CFU/mL of cells exhibited their maximum CA activity by 124.84 ± 1.24 and 140 ± 2.5 U/mL at 132 h for Zn2+ and Cr6+, respectively. Simultaneously, with pH increase to 9.5 ± 0.2, a complete removal for both metals was observed at 168 h; Ca2+ removal percentages recorded 78.99% and 85.06% for Zn2+ and Cr6+ remediating experiments, respectively. Further, the identity, elemental composition, functional structure and morphology of bioremediated precipitates were also examined via mineralogical analysis. EDX pattern showed the typical signals of C, O and Ca accompanying with Zn2+ and Cr6+ peaks. SEM micrographs depicted spindle, spherical and cubic shape bioliths with size range of 1.3 ± 0.5–23.7 ± 3.1 µm. Meanwhile, XRD difractigrams unveiled the prevalence of vaterite phase in remediated samples. Besides, FTIR profiles emphasized the presence of vaterite spectral peaks along with metals wavenumbers. Conclusion CA enzyme mediated Zn2+ and Cr6+ immobilization and encapsulation inside potent vaterite trap through microbial biomineralization process, which deemed as surrogate ecofriendly solution to mitigate heavy metals toxicity and restrict their mobility in soil and wastewater.


Background
Heavy metals are widely used in modern industrial production.However, due to the incomplete waste treatment after production process, a large amount of sewage and solid wastes containing such heavy metals are discharged into the environment, causing serious environmental pollution [1].Metal ions persist in the environs due to the bioaccumulation affinity of living organisms and their restricted ability to metabolize them into less toxic forms.Once they absorbed into the cell of living organisms, they bind to vital components, for instance nucleic acids, structural proteins and enzymes, impairing thereby numerous essential metabolic functions [2].
Several industrial processes such as textile dyeing, tanneries, metal electroplating, battery manufacturing units, mining, metallurgy, pigment/paint, galvanizing steel, paper bleaching and corrosion inhibition wastewater treatment plants are responsible for increasing Chromium (Cr) and Zinc (Zn) concentrations as particles released into the terrestrial atmosphere or as wastewater discharged into the ground or waterways.Bearing in mind, Cr (VI) is a well-known carcinogen, teratogen, mutagen, and have a toxic effect on all living systems with strong migration ability, whose toxicity and mutagenicity are exceeded that of Cr (III) by 100 times [3].Regrettably in similar way, exposure to certain levels of Zn, even within a short time, can cause severe health issues such as stomach cramps, anemia, nausea, vomiting, reducing of HDL cholesterol and pancreas damage [2,4].Hence, the World Health Organization (WHO) recommended a maximum acceptable Zn concentration in drinking water of 5.0 mg/L, and the allowable limit of Cr (VI) in natural water is 0.05 mg/L [5].Based on the above mentioned facts, the treatment of water and soil contaminated with toxic metals is a challenging issue to find technically, economic, feasible, ecofriendly and effective procedures.Several technical approaches have been proposed to control metals toxicity, among them, carbon adsorption, chemical precipitation, electrochemical treatment, ion exchange, reverse osmosis, and membrane separation.However, these conventional methods are expensive, energyintensive, and produce toxic sludge that requires special and additional handling [6].Remarkably, the biological remediation pathways are deemed being environmentally friendly to eliminate toxic metals, with unique advantages of low cost and high efficiency, especially at low concentrations [7].Interestingly, a plethora of former studies have proved the feasibility of bioremediation using microorganism to detoxify metals and transform them to less toxic forms [7][8][9].Microbial remediation mechanisms mainly include bioaccumulation, biosorption, bioleaching, biotransformation, and biomineralization, which basically rely on metalmicrobe interactions [10][11][12].Notably, the microbial mineralization is described as one of the most promising bioremediation techniques that could be executed through many microbial metabolic processes, which generates numerous biominerals.About 60 and more type of biominerals have been produced via microbial mineralization process either directly or indirectly such as phosphate, oxalate, and carbonate that could encompass numerous metal ions, like Mg + , Mn ++ , Fe +++ , Ca ++ , etc. [13,14].
Microbial induced calcium carbonate precipitation (MICP) is a unique biomineralization process by which different microbes induce bio-mineralization of CaCO 3 in the presence of dissolved calcium (Ca 2+ ) and carbonate (CO 3 2− ) ions [15].This process has been investigated as a potential method for the removal of heavy metals from contaminated water sources.It mediates the mineralization of heavy metals from the ionic soluble state into a stable precipitated solid form, thus, reduces the mobility and toxicity of those hazardous metals [16,17].The bio-precipitation of calcium carbonate could be implemented via either an autotrophic or heterotrophic pathways, while the later seems to be more popular.The autotrophic biosynthesis of CaCO 3 involves different mechanisms such as methanogenesis, aerobic photosynthesis and anaerobic photosynthesis [18].On the other hand, three main categories of microorganisms enhance bio-calcification technique heterotrophically.The first category induces the reduction of sulphate by sulphate reducing bacteria (SRB) [19].The second category comprises microorganisms that engage in nitrogen cycle by the subsequent ways: (A) ammonification of amino acids, (B) denitrification and C) urea hydrolysis [20].The third category enhances the reversible transformation of carbon dioxide (CO 2 ) to bicarbonate through carbonic anhydrase enzyme (CA) [17].
CA has been associated in MICP, CO 2 sequestration and subsequently lessens greenhouse effect.It catalyzes the reversible CO 2 hydration and dehydration of HCO 3−  through Ping-Pong mechanism (i.e., 2-stage); utilizing zinc bound hydroxide ion found in the enzymes active site.Through MICP process, HCO 3− can interact with Ca 2+ under alkaline conditions, which were precipitated as CaCO 3 crystals [15,17].Thus, this enzyme possesses the ability to accelerate CO 2 uptake, besides carbonate rock dissolution at the same time.Notably, it is anticipated that CA reaction will afford the key vital molecules for the formation of CaCO 3 when proceeding in the forward direction [21].In addition, CA has many other functions, as it has been established to be imperative to plant growth, stomatal development, responses to various stresses, as well as participating in photosynthetic light reactions, while capturing the atmospheric CO 2 [22,23].On the other hand, CA is requested for rapid processes in several organisms, commonly transport processes.For instance, it is necessary in the removal of CO 2 from lungs for the synthesis of eye secretions.Additionally, CA keeps optimum level of HCO 3− and CO 2 in the body, which both used as substrate for many enzymatic reactions.Moreover, CA has vital role in blood as it retains acid-base balance, helps in homeostasis the physiological pH, and also participates in respiration and ion transport [24].However, the majority of studies concerning with calcification technology addressed urea hydrolysis and few researches were interested in denitrification mechanism [25].Nonetheless, no study, till our knowledge, documented the recruitment of carbonic anhydrase enzyme in heavy metal removal through MICP process.
Remarkably, the employed microorganism in any application must be effective in its performance, however, other secondary properties (e.g., biosafety, proliferation requirements/conditions, formulation options, bioreactor application, recovery, etc.) are just as or even more influential.Therefore, noteworthy shed the light to the characteristic features of unicellular fungi (i.e., yeasts), which appeal scientific community for their recruiting in wide spectrum of applications, with less complications elicitation than filamentous fungi or bacteria.They are relatively abundant in various ecosystems, easily culturing, wide nutritional requirements, facultative anaerobes, short generation time relative to multicellular fungi, higher biomass yield and better adaptability with higher metabolites productivity per biomass unit.Let alone, their higher stability and complexity of genome organization, biofilm forming capacity, superior metal bioaccumulation potentials, enhanced tolerance to adverse circumstances with multiple detoxification mechanisms, which situated them as a platform for innovative applications [26].Strikingly, all previously mentioned traits of yeasts triggered them a promising biogenic factory for multifaceted applications in all aspects of biotechnology, which were symbolized by rainbow code, namely red biotechnology [27,28], yellow biotechnology [29,30], gray/ white biotechnology [31][32][33], gold biotechnology [34,35], green biotechnology [36][37][38][39].
Intriguingly, the capability of yeasts to generate calcite through MICP was recorded previously [11,12].Besides, their versatile metabolic diversity that enable them to exhibit a noticeable survival performance in extreme environments such as concrete, limestone, mortar, marble and granite were also documented [9].Despite several remediation mechanisms exerted by different yeast stains for purifying the environment from pollutants, the remediation through MICP proved its superior effectiveness.Wherein, the remediated metals entrapped inside durable CaCO 3 matrix; displaying thereby more potent trap than surface adsorption, internal accumulation or even hydroxide precipitates, which exposed to re-dissolution in water and generation of soluble anionic hydroxyl complexes [40].Additionally, Eltarahony et al. [40] found that there were no dissociation or dissolution of remediated metals sequestered in the calcareous trap under acidic pH (pH 3.8 ± 0.045), reflecting the stabilization of immobilized metals with low chance to release in ambient ecosystems even under highly acidic conditions such as acidic rain.
Nevertheless, no study to the best of our knowledge invested the capacity of yeast in remediating heavy metals through MICP process, in particular via enhanced performance of CA enzyme.While Barbero et al. [41] had successfully edited the genes encoding carbonic anhydrase and mineralization peptides in yeast, which would significantly enhance its ability to induce CaCO 3 precipitation.Based on the previous background, the present study, for the first time, aimed to eliminate Zn 2+ and Cr 6+ using the carbonate precipitation efficiency of carbonic anhydrase-producing yeast isolate.The core of this study focused on the statistical optimization to maximize the performance of carbonic anhydrase, which thereafter, enhanced MICP process and elevated Zn 2+ and Cr 6+ elimination efficacy .Finally, the removal of examined metals was confirmed through various mineralogical analyses (i.e., energy dispersive X-ray spectrometry (EDX), scanning electron microscopy (SEM), X-ray diffraction (XRD), and Fourier transform infrared spectroscopy (FTIR) of Zn 2+ and Cr 6+ enclosed in CaCO 3 deposits.

Microorganism, and cultural conditions
A water sample was collected from Marsa Allam, Red sea governerate, Egypt.Directly after sampling, isolation and screening of biomineralizing microbes were performed.Initially, the samples were inoculated in Basal media (B4) containing glucose (21) g/L, ca-acetate (18 g/L), NaNO 2 (1.5) at pH 7.0 ± 0.2, 30 ℃ [40].Out of 15 microbial isolate, one isolate was selected based on its higher carbonatogenesis capability.The strain was stored in glycerol (20%, v/v) at − 20 ℃ for the forthcoming investigations.

Morphological and cultural characterization
The morphology and dimensions of the selected isolate were determined from photomicrographs utilizing scanning electron microscopy (SEM) (JEOL JEM-1230, Japan-Faculty of Science-Alexandria University).For colony characterization, the cells were cultivated under both aerobic and anaerobic conditions on YPD agar (Yeast Extract, 10.0 g, Peptone 20.0 g, Dextrose, 20.0 g, Agar, 15.0 g) and incubated at 30 ℃. [42].About 0.5 McFarland equivalents to about 1.5 × 10 8 CFU/mL was inoculated in 250 mL Erlenmeyer flasks containing 100 mL of the YPD medium and incubated at 30 °C in an orbital shaker (STUART SI500) at 150 rpm.After the incubation period, the cells were centrifuged at 10.000 Xg for 20 min at 4 ℃ and the harvested cells were used for further tests.

Molecular identification
The selected isolate was identified using 18S rRNA gene sequencing.The genomic DNA of the selected isolate was extracted from overnight pure cultures and 18S rRNA gene was amplified using18S primers [45], then the purified PCR product was sequenced as described elsewhere [46].The phylogenetic affiliation was inquired by applying BLASTn analysis to determine the similarities with their available GenBank database sequences.Its generated sequence was submitted to the GenBank to obtain its corresponding accession number.For multiple alignment and phylogenetic tree construction, the software package MEGA-6 was employed.

Assay of carbonic anhydrase enzyme (CA)
The CA was determined in the fungal pellets that were collected by centrifugation at 10,000 Xg for 20 min.Then the pellets were dissolved in phosphate buffer (pH 7.0), mixed well and disrupted by ice cold TSE buffer (10 mM Tris HCl, 100 mM NaCl, and 1.0 mM EDTA at a pH of 7.8).The reaction was incubated for 30 min at 30 ℃ with good vortex every 10 min.The cell debris and unbroken cell were removed by centrifugation at 10,000 Xg for 3 min, while the supernatant represents the crude enzyme [47].CA activity was assessed by calculating the micromole of p-nitrophenol liberated from p-NPA (p-nitrophenyl acetate).The reaction mixture composed of (in order of addition):825 μL of Phosphate buffer (50 mM, pH 7.5), 175 μL of the substrate stock solution (p-NPA, 10 mM in isopropanol).The blank reactions contain the media without yeast culture.Then, 25 μL of an enzyme was added to initiate the reaction.The reaction mixture was incubated at 37 ℃ in a water bath for 5 min [48].The amount of p-nitrophenol released was estimated at 410 nm using a spectrophotometer [47].One unit of CA activity correlates to the quantity of enzyme that enhances the formation of 1 μmol of p-nitrophenol per minute under standard assay conditions.

Experimental design and statistical analysis
In order to optimize the concentration of all media components simultaneously and not by changing a single factor, statistical experimental designs were applied in two steps, namely, (i) Plackett-Burman Design "PBD" to determine the important factors affecting CaCO 3 precipitation and CA activity, (ii) Central Composite Design (CCD) as a type of Response Surface Methodology (RSM) to infer the exact concentrations of important factors that achieve the maximum calcium carbonate (CaCO 3 ) precipitation and CA activity.

Determination of significant parameters by Plackett-Burman design (PBD)
PBD is the design that intended for screening and identifying the controlled experimental parameters (nutritional, environmental and incubation conditions) based on their main effect on CA enzyme activity and CaCO 3 weight.PBD investigates (n) variables which are expressed at two levels, high (+) and low (−).Herein, a total of 7 (n) variants were studied at 2-level concentrations in 12 experimental matrices, as shown in Table 1.Each experiment was done in triplicate and CA activity and CaCO 3 weight were assessed as response [40].The Plackett-Burman experimental design is built on the first order model (Eq.1): where, Y is the response or dependent variable (CaCO 3 weight and CA activity); it will always be the variable we aim to predict, βo is intercept of the model, βi is the linear coefficient, and Xi is the level of the independent variable.The statistical analysis output will be used to calculate the significance of the variables depending on their nature; and their positive or negative effects on CA enzyme activity and the weight of CaCO 3 .

Central composite design (CCD)
To determine the optimal levels of the most interesting variables that were detected by PBD and to infer their interactions, RSM was applied using CCD.Five levels (− 2, − 1, 0, + 1, + 2) of four different variables (i.e., glucose, sodium nitrite, calcium acetate and inoculum size) were studied in 31 experiments as listed in Table 2. Along with each experiment both CA activity and CaCO 3 weight were determined.Each Experiment was performed in triplicate and the mean was calculated for subsequent statistical analysis [40].Considering statistical estimation, the relationship between the coded and actual values is represented by Eq. 2: where Xi is the coded value of the ith variable, Ui is the real value of the ith variable; Ui0 is the real value of the ith variable at the centre point and ΔUi is the step change of variable.The second order polynomial structured described in Eq. 3: where: Y is the predicted response; X1, X2, X3 are input variables that affect the response variable Y; β0, intercept; β1, β2 and β3 linear coefficients; β11, β22 and β33, squared or quadratic coefficients β12, β13, and β23 interaction coefficients.

Statistical analysis and verification of the model
The statistical software Minitab 14.0 (Minitab Inc., Pennsylvania, USA software) was used for establishing the experiment designs "matrices" and subsequent statistical analysis of PBD and CCD data (regression analysis and ANOVA).Further, the relationship between the response and variables was also illustrated graphically by threedimensional surface plots (3D) and two-dimensional contour plots (2D).In addition, the optimizer tool was used to predict the optimum level of experimental factors.Under predicted optimized conditions, the model was validated through comparing CA enzyme activity and the weight of CaCO 3 with that obtained from the basal conditions [40].

Minimal inhibitory concentration test
Firstly, the metal toxicity experiment was performed to decide the minimum inhibitory concentration (MIC) of Zn 2+ and Cr 6+ .Different concentrations (800, 400, 200, 100, 50, 25 and, 12.5 ppm) of Zn 2+ (ZnCl 2 ) and Cr 6+ (K 2 Cr 2 O 7 ) were examined in YPD agar media inoculated with 0.5 McFarland scale.The inoculated plates were incubated at 30 ℃ for 48 h in triplicate.Accordingly, 150 mL of optimized mineralizing media in 500 mL flasks supplemented with 50 ppm of Zn 2+ and 400 ppm of Cr 6+ (½ MIC) were inoculated by 10 8 CFU/mL.The flasks were incubated in a rotary shaker (150 rpm) at 30 ℃ for 7 days.At the same time, two controls were run in parallel; abiotic controls or negative controls which was devoid from fungal inoculum and the biotic controls contained inoculated biomineralization media without heavy metals [40].

Monitoring of chemical changes during Zn
Zn 2+ and Cr 6+ .During the biomineralization mechanism, the culture media was drawn at a constant time interval to evaluate the previous parameters.The cell number (CFU/mL) was detected by the pour plate method; the CA activity and pH were assessed as described previously.The concentrations of soluble Ca 2+ , Zn 2+ and Cr 6+ were estimated by an inductively coupled plasma optical emission spectrometer (Agilent ICP-OES 5110DVD) (Central Lab, Alexandria university).All experiments were performed in triplicate and the mean values were considered.At the end of the incubation period, all precipitates were centrifuged at 10,000 Xg for 20 min, then washed, dried and subjected to mineralogical analysis [40].The bioremediation efficiency of examined strains was calculated as a percentage representing the differences between the initial and final concentrations of Ca 2+ , Zn 2+ and Cr 6+ in the supernatant.

Mineralogical and morphological analysis
To evaluate the identity, morphology, microstructure and chemical constituents of the precipitated samples (positive control and remediated deposits with Zn 2+ and Cr 6+ ), X-ray diffraction (XRD), Energy dispersive X-ray spectroscopy (EDX), scanning electronic microscopy (SEM) and Fourier transform infrared spectroscopy (FTIR) were utilized.XRD was used to identify the precipitated minerals (Bruker MeaSrv D2-208219, Germany-Central Lab, Faculty of science, Alexandria University), which were irradiated with Cu Kα (λ = 0.15406 nm), generated at 30 kV and 30 mA at a scan rate of 2°/min for 2θ values over a wide range of Bragg angles 10° ≤ 2θ ≤ 80. Microchemical analysis of the samples was achieved using EDX analyzer combined with SEM (JEOL JSM 6360LA, Japan).The morphological features of deposits were visualized using SEM (JEOL JSM 6360LA, Japan-Advanced Technologies and New Materials Research Institute (ATNMRI) SRTA-City) at an accelerating voltage of 20 kV.To scrutinize the accompanying functional groups of CaCO 3 crystals, FTIR was performed by Shimadzu FTIR-8400 S, Japan with a resolution of 4 cm −1 .Sample preparation prior to FTIR analysis began with mixing the bioprecipitated CaCO 3 crystals with KBr followed by pulverization to powder and pressing into discs.
The spectrum was scanned with a range of 4000 and 400 cm −1 .

Screening of CaCO 3 producing microorganisms
Actually, the signal for the MICP approach was monitored through the presence of solid deposits on B4-agar medium [40].Herein, the isolate under study was selected among 15 isolates as precipitated the highest amount of CaCO 3 in B4-broth medium comparing to the others.Subsequently, the molecular characterization was performed by identifying the partial 18S rRNA gene sequence (≈500 bp), which revealed 99% sequence resemblance with all species of the genus Rhodotorula and less than that with other genera.Remarkably, the accepted taxonomic affiliation standard is considered at sequence similarity greater than 98% between reference and inquired strains as mentioned by Nunes et al. [43].Accordingly, it was identified as Rhodotorula sp.
Its nucleotide sequence was submitted to GenBank with the accession number of MZ312359.The phylogenetic position was constructed as illustrated in Fig. 1A, the members of genus Rhodotorula are affiliated to the basidiomycetic fungi and well-known for their ability to biosynthesize a diverse range of valuable biomolecules, including carotenoids, lipids, enzymes, and polysaccharides [49].However, it was not studied before as calcifying microbe, till our acquaintance.

Cultural, physiological and morphological characterization
On YPD agar, the colony appeared spherical, pigmented, flat or slightly raised with smooth borders as illustrated in Fig. 1B.It characterized by orange pigments which seemed to be carotenoid that acts as defensive mechanism against oxidative stress [50].It exhibited aerobic and anaerobic growth (i.e., facultative anaerobic) but at a faster pace in oxygen consuming conditions.As a unicellular eukaryot, it was visulaized as non-motile elongated spherical or oval cells shape under SEM, their length ranged from 1.5 to 2.8 μm and its width lied between 0.5 and 0.9 μm Fig. 1C.Besides, the examined Rhodotorula sp.MZ312369 showed various biochemical activities as illustrated in Table 3.It was capable of growing on a broad temperature range varies from 10 ℃ to 40 ℃ with optimum at 25-30 ℃, below and above this range the growth adversely affected as harmonized with Allahkarami et al. [51].Whereas, its growth in different pH ranged from 6 to 9 reflected its optimal performance at an initial pH of 7.0-8.0.Similarly, the physiological and biochemical characters of our isolate agreed with that examined by Nunes [43].

Determination of significant parameters by Plackett-Burman (PBD)
The Plackett-Burman experimental design is a useful tool for the rapid evaluation and screening of the significant nutritional and environmental parameters affecting on examined responses [52][53][54], which were CA activity and CaCO 3 weight.The influence of 7 independent factors on the overall carbonatogenic process (i.e., CA activity and the weight of CaCO 3 ) was detected in 12 experimental runs.The results exhibited a wide variation in enzyme activity from 2 to 116 U/mL and in the CaCO 3 precipitates from 0.02 to 1.13 g as shown in Table 1.Such variation reflected the leverage of nutritional and environmental factors combination on CA activity and CaCO 3 weight.It is important to find out the adequacy and significance of the model, hence, the data of CA activity and CaCO 3 weight were analyzed to declare multiple linear regression analysis and analysis of variance (ANOVA).Generally, the coefficient with lower probability P-value (≤ 0.05) indicates that the corresponding factor is significant and possesses a high effect on response (CaCO 3 weight & CA activity), while a coefficient with probability P-value exceeds 0.05 reveals the negligible effect on response.Besides, the sign of the examined variable's coefficient gives an insinuation about the influence of variable's concentration on response.Namely, the positive sign of coefficient indicates the higher response upon using the higher value/concentration of the corresponding variable and vice versa.Out of Table 4, the calculated confidence level and p-values manifested that glucose, Ca-acetate, NaNO 2 and inoculum size were the significant parameters on MICP process.The model coefficient of determination R 2 and adjusted-R 2 were 0.9788%, 0.9789% and 0.9417%, 0.9421% for CaCO 3 weight and CA activity, respectively, which implies that 97.88% and 97.89% of the variation in the data were illustrated by the model Table 4.While there were just 2.12% & 2.11% chance that could occur because of the noise.This again ensured a satisfactory adjustment and a good relation between the observed and the predicted values.As mentioned in Mojtaba and Fardin [55], the closer R 2 is to 1, the better the estimation of regression equation fits the sample data.However, the residuals were studentized and their values were in the range of ± 2, which fall in reasonable range as denoted by Anuar et al. [56].Additionally, as revealed by ANOVA, the first-order

Regression and analysis of variance (ANOVA)
Multiple regression analysis was employed to statistically analyse the data of CaCO 3 weight and CA activity as tabulated in Table 5, which also included the values of R 2 , adjusted-R 2 , the coefficient estimates, probability P-value, lack-of-fit, linear, quadratic and interactions impacts as well.As noticed, the values of R the independent variables and only 6.9% and 7.3% could not be explained in the view of models.Besides, the Adj-R 2 values were quantified as 0.87 and 0.862, which emphasized the model significance.Notably, the small difference between R 2 and adjusted-R 2 reflects the good coordination between the actual experimental values and the predicted values of both responses; thus, the models of the current study were optimal within the range of experimental factors to predict an efficient CaCO 3 weight and CA activity.In addition, the positive coefficient values pointed out that the linear effect of all variables, quadratic effect and mutual interactions effect of some factors exhibited synergistic leverage in CA activity and CaCO 3 weight (i.e., their higher values enhance CA activity and CaCO 3 weight).While the other factors, which displayed negative coefficient values signifies their higher impact on carbonatogenesis process at their negative value.For instance, the positive sign of glucose coefficient indicated a positive influence on the weight of CA activity and CaCO 3 at high concentrations, whereas with more increasing in its concentration, it inhibited the overall process, as evident from negative coefficient of the squared term.
For interactive terms, the interaction between (glucose & Ca-acetate), (glucose & inoculum size), (Ca-acetate & inoculum size) and (NaNO 2 & inoculum size) were described as insignificant as given by P-value, which exceeded 0.05.While the other relations, namely glucose & NaNO 2 and Ca-acetate & NaNO 2 were both significant.Meaning that they can act as limiting factors and any little difference in their values will alter vividly the CA activity and CaCO 3 weight as well [59].Table 5 informed that the relation influence between (glucose & ca-acetate), (glucose & inoculum size), (Ca-acetate & NaNO 2 ), (Ca-acetate & inoculum size) and (NaNO 2 & inoculum size) were positive (i.e., synergistic effect) as CA activity and CaCO 3 weight increase with increasing in both factors.While, antagonistic effect appeared to be prevailing between glucose and NaNO 2 .That means the effect of higher level of one parameter increases CaCO 3 precipitation process with lower level of another examined parameter.Generally, the regression models found in this study was highly significant as denoted by the low P-values with 0.000; the linear and quadratic effects appear to be predominant over interactive effect as observed in Table 6.In addition, as unveiled by ANOVA, Lack-of-fit assessed by 0.259 and 0.149 for CaCO 3 weight and CA activity, respectively, which indicated their insignificance and reflecting the robustness of both models and their accuracy.Moreover, the analysis of the data developed from RSM generally considers a second order polynomial equation, which explains the relation among the response variable(s) and the factors.In this equation the linear, the interactions and quadratic effects of each factor on the response variable(s) are determined.The second-order polynomial equation which defines the predicted response formulated as follows Eqs.

Graphical demonstrations of the response surface model
The three-dimensional surface (3D) and two-dimensional contour (2D) plots were generated to understand the interaction of the variables and responses (CA activity and CaCO 3 weight) and also predict the ideal level of each variable for maximal response; such plots are a graphical representation of the model equations achieved in the regression analysis [59].The response surface plot is a 3-D graph which represents the empirical functional relation among the response with the vertical axis and two factors on horizontal axes representing the coded levels of two explanatory factors, as the remaining factors being held at constant levels.The optimum values for the variables were achieved by moving along the major and minor axis of the contour.The plot was utilized to visualize how a response varied with differences in the factor [59]. Figure 2A, B, C, D depicted the CA activity and CaCO 3 weight as a function of glucose and NaNO 2 ; the surface plot was convex suggest that there are welldefined optimal variables.Additionally, as the variable ranges were suitable, the optimum lies in the design space.As depicted, the CA activity and CaCO 3 weight increased with increasing glucose, while decreasing NaNO 2 concentration and vice versa.Moreover, increasing variable concentrations resulted in decreasing the carbonatogenic parameters.Obviously, the contour-2D plot revealed significant antagonism interaction.On the other hand, surface plot and contour plot of calcium acetate and NaNO 2 highlighted the mutual interaction influence on both carbonatogenic parameters as represented in Fig. 2E, F, G, H. Wherein, the carbonatogenic parameters raised with uplifting the concentrations of both variables synchronously till arrived to the highest possible level.Then, CA activity and CaCO 3 weight started to decrease with more raising in variable values.Commonly, the contour plot of this mutual interaction display elliptical shape; expressing significant interaction [60].Meanwhile, Fig. 2I, J, K synergetic interaction symbolizing by circular shape contour plot.

The desirability function for prediction of the optimum conditions and model validation
The substantial aim of the statistical design of experiment focuses on attaining the maximum performance of carbonatogenic process, via both parameters of CA activity and CaCO 3 weight, in the terms of optimum levels of examined variables, which could be obtained by desirability function.The factor settings with maximum desirability were: glucose 21, ca-acetate 18, NaNO 2 1.5 g/L and inoculum size (25%).As predicted, CA activity and CaCO 3 weight assessed by 117.6048U/mL and 1.0785 g with desirability value recorded 0.94012 and 0.82649, respectively.For verification of such prediction, the experiments were carried out in triplicates for each trial.
The attained experimental values recorded 154 U/mL and 1.9 g with 1.8-fold increase comparing to the basal conditions; manifesting the good relation between the observed and predicted values.
It is noteworthy that The MICP process is influenced by many factors such as bacterial solution concentration, Ca 2+ concentration, carbon, and nitrogen sources [61].

Minimum inhibitory concentration
The MIC values of Rhodotorula sp.MZ312359 aganist Zn 2+ and Cr 6+ recorded 100 and 800 ppm, respectively.Our data agreed that reported by Grujic [62], who found that Rhodotorula sp. had approximate MIC values; implying its potency in bioremediation process of considerable concentrations of metals.

Monitoring of chemical changes during Zn 2+ and Cr 6+ removal
The biomineralization dynamics of the whole remediation process was achieved through detecting changes in the chemistry of bioremediation solution in parallel to the biotic control.Therefore, the biomineralization parameters representing in CA activity, CaCO 3 weight, pH fluctuation, soluble Ca 2+ concentration along with microbial count, Zn 2+ and Cr 6+ concentrations were determined as a function of time during 10 days' incubation (Fig. 3).Generally, a positive correlation between microbial growth, CA activity, CaCO 3 weight was clearly evident, which was synchronized with pH uplifting and elimination of soluble (Ca 2+ , Zn 2+ and Cr 6+ ).In the biotic control, the growth profile of Rhodotorula sp.MZ312359 displayed a typical growth phases (lag, logarithmic and stationary) with gradual increase in the cell number and carbonatogenic parameters till reach to the maximum CA activity at 144 h, which remained stable Fig. 2 Contour Plots (A, C, E, G, I, K) and Surface plots (B, D, F, H, J, L) for CaCO 3 weight and CA activity of Rhodotorula sp.MZ312359 showing the interactive effects of some variables to some extent till 180 h by recording 154.59 ± 0.193 U/ mL with cell density assessed by 11.84 × 10 8 ± 2.74 × 10 6 CFU/mL.Besides, the solution pH recorded an elevation consistently during incubation period, which assessed by 7.1 ± 0.2 at the onset of the experiment and reached to 9.58 ± 0.015 at the end.Similarly, CaCO 3 weight steadily increased within incubation period recorded the maximum value by 1.73 ± 0.1 g at 144 h and remained stable, without noticeable significant increase, till 10 days of incubation.Interestingly, all such increasing in the examined parameters was associated with Ca 2+ depletion, which reached to 95.83% removal at 144 h and was entirely exhausted by the end of incubation.
Whereas, the overall diminishing in cell density concurrently with CA activity during growth stages was observed in presence of Zn 2+ and Cr 6+ .That could be possibly attributed to the stress effect of both metals on the examined cells, which entailed acclimatization performance with them; symbolizing in slight retardation in the initial growth stages.Similar finding was reported by study performed by Mwandira et al. [63].Broadly, in remediation experiments, the maximum CA activity displayed by 9.23 × 10 7 ± 2.1 × 10 6 CFU/mL and 10.88 × 10 7 ± 2.5 × 10 6 CFU/mL were 124.84 ± 1.24 and 140 ± 2.5 U/mL in 132 h for Zn 2+ and Cr 6+ , respectively.This implied the leverage of both metals on the whole metabolic performance of the microbial cells.Simultaneously, a gradual raising in pH was noticed from 7.1 ± 0.2 to 9.45 ± 0.2 and 9.51 ± 0.2 for both remediated metals in the same order.On the other hand, there was vividly increasing in soluble Ca 2+ removal percentages as inferred by ICP-OES analysis, which assessed by 51.55 and 67.05% in Zn 2+ and Cr 6+ remediating experiments, respectively at 144 h.Notably, such removal percentage increased upon the end of incubation and reached to 78.99 and 85.06% for Zn 2+ and Cr 6+ , respectively.The exhaustion of such soluble form of Ca 2+ was harmonized with its precipitation in a solid phase as determined by CaCO 3 weight, which evaluated by 1.51 ± 0.2 and 1.62 ± 0.15 g.Interestingly, the effective metals remediation process recorded 81.9% for Zn 2+ and 93.9% for Cr 6+ were displayed within 144 h., respectively.A complete removal for both metals was observed at 168 h.It is worth mentioning that the absence of any precipitation in the non-biological (chemical) or abiotic control, reflects the influence of microorganisms on the change of the physical and chemical factors of the culture medium, thus, promoting the precipitation of CaCO 3 .Meanwhile, Mugwar and Harbottle [64] reported the removal of small concentration of Zn 2+ by Sporosarcina pasteurii through MICP within 1 week, which was in tandem with our results.

Mineralogical and morphological analysis
While studying the dynamic of Zn 2+ and Cr 6+ removal, the data unveiled their absence in soluble form and implied their precipitation via MICP process.The employment of EDX, XRD, FTIR and SEM analysis emphasized their incorporation in the bioremediated deposites of CaCO 3 crystals, confirming by such way the efficacy of Rhodotorula sp.MZ312359 in MICP-based bioremediation process by the dint of CA activity.

EDX analysis
The EDX microanalysis of the precipitated crystals in the biotic control sample expressed in Fig. 4A.The spectrum showed unique peaks at 0.277, 0.525, and 3.69 keV, which are related to the binding energy of carbon, oxygen, and calcium, sequentially [65,66].Other peaks relevant to the binding energy of Zn 2+ were distinguished by the Kα and Lα characteristics at 8630 and 10.12 keV (Fig. 4B).EDX also illustrated the presence of Cr 6+ through the distinctive emission peaks of Kα and Lα at 5.411 and 0.573, as demonstrated in Fig. 4C.These results affirmed that Zn 2+ and Cr 6+ were involved in calcareous bioremediated precipitates.Obviously, the existence of nitrogen (0.39 keV) and phosphorus (2.013 keV) peaks in significant amounts pointed out to the biological nature of the bioprecipitate.Actually, these elements are vital components of microbial cells biomolecules that compose proteins, nucleic acids, phospholipids and lipopolysaccharides [67].Generally, the calcium peaks intensities and their correlating atomic percentages, which were higher than carbon peak, may reflect higher purity in structure as implied by Caicedo-Pineda et al. [68].

Scanning electron microscope (SEM)
The morphological and textural properties of bioremediated samples as well as the control were visualized by SEM (Fig. 5).The control micrograph illustrated acicular spheres (red arrows; 13.3 ± 2.5 µm) with other spindle shaped particles (green arrows; 1.5 ± 0.4 µm).As described by Wagterveld et al. [69], such crystal shape called morning stars.On the other hand, the bioremediated Zn 2+ bioliths appeared mixed morphologies with spindle shaped-particles (green arrows) seemed being coalescent forming spherical particles ((red arrows; 5.3 ± 1.1 µm).Strikingly, larger cylindrical particles (blue) arrows also were detected with size assessed  with rough and wrinkled surface also were evident (5.1 ± 1.7 µm).Notably, square or cubic forms (violet arrows) were also observed with well-defined faces.

XRD
The heterotrophic precipitation of CaCO 3 is emphasized by XRD spectra.As shown in Fig. 6, sharp, distinctive, identifiable, and wide diffraction peaks were detected at 2θ values (hkl) of 20.9 (004), 24.8 (110), 27 (112), 32.8 (114), 43.7 (300), 50.1 (118), 55.9 (224) in the biotic control.These peaks accentuated the presence of crystalline CaCO 3 deposits in vaterite phase as corroborate with the standard JCPDS card No. 72-0506 [71].Notably, the crystallographic profile of the bioremediated samples, either Zn 2+ , or Cr 6+ , manifested the predominance of vaterite phase with minor shifting in some peaks, reduction in peaks intensities and also the disappearance of the others, comparing to the biotic control; implying the incorporation or replacement of remediated metals in vaterite matrix.Interestingly, a study conducted by Han et al. [72][73][74] empowered our finding.Despite the clear deposition of metal-CO 3 (i.e., ZnCO 3 or Cr 2 (CO 3 ) 3 was not observed in our study, their elimination process was fulfilled in stable solid phase as revealed by EDX.The obtained results might be assigned to the low crystallinity degree of both Zn 2+ or Cr 6+ related biominerals, amorphous existence of them and the replacement of Zn 2+ / Cr 6 to fewer positions of Ca 2+ in the vaterite matrix.Besides, Zn 2+ /Cr 6+ were continuously and gradual integrated/encapsulated by co-precipitation or even filled in the defect vacancy of vaterite during the uncontrolled yeast growth, which synchronize with continuous crystal nucleation/growth stages of CaCO 3 precipitates.Subsequently, Zn 2+ /Cr 6+ concentrations were too low to crystallize on vaterite matrix and their carbonate phases was entirely wrapped by outer CaCO 3 .Intriguingly, the flexibility of mycogenic vaterite structure with sufficient surface area, electronegativity and porosity were also taken in consideration for explaining our results.Therefore, it was plausible to describe such biogenic vaterite as disingenuous scavenger that not only absorbed Zn 2+ and Cr 6+ but also sequestered them in robust trap.In accordance with our results [75] documented the exact results in removing Cr(VI) via chemical precipitation without obvious change in XRD pattern.Similarly, [1] supports our finding, in which all minerals precipitated chemically, under controlled conditions of Ca 2+ -binary combining with either Cd 2+ or Zn 2+ , were only calcite peaks without any evidence of other metal-carbonate phases (i.e., otavite and smithsonite); implying the incorporation of both metals inside calcite crystal rather than adsorbed on the surface, followed by precipitating of Zn 2+ -bearing calcite and Cd 2+ -bearing calcite in lieu of end-member carbonate phase.Likewise, Hua et al. [76], highlighted that the adsorption of Cr 6+ during chemical synthesis of CaCO 3 crystals did not modify the mineral structure.In parallel, a study implemented by Tang et al. [77] assured the efficiency of decreasing the concentration of soluble Cr(VI) in co-precipitation process; suggesting its preferential incorporation into calcite lattice in the form of carbonate-bound Cr(VI) throughout the crystal growth process.Meanwhile, a ureolytic Vibrio harveyi strain selectively transformed CaCO 3 polymorph and maintained vaterite stabilization in the presence of Zn 2+ , comparing to Pb 2+ , Cd 2+ and Cr 6+ [78].In the same sense, Qiao et al. [79] found that ureolytic strains of Lysinibacillus sp., and Pseudochrobactrum sp.removed about 48% and 53% of soluble Zn 2+ (80 ppm) and their XRD results didn't match any reference XRD patterns of standard zinc carbonate, zinc hydrogen, zinc oxygen or even zinc hydroxide precipitations, although a clear and distinct EDX-signal was detected.All these scholars also harmonized with our findings.Besides, Al Disi et al. [80] ascribed the ability of oil-degrading Providencia rettgeri and Pseudomonas aeruginosa in remediating Cr 6+ , Cu 2+ , Zn 2+ and Ni 2+ through interacting of these heavy metals with MICP-products, which were calcite (82%) and brushite (18%), whereas, their enhancing/inhibiting calcite solubility based on metals type and their content.On the contrary, Penicillium chrysogenum CS1 immobilized a relatively lower amount (i.e., 34%) of Cr (VI) in solution containing 200 mg through fungal-based MICP in 12 days in the form of chromium oxide carbonate as identified by XRD analysis [81].Besides, strains of Sporosarcina luteola possessed the ability to precipitate Mn 2+ , Cd 2+ , Sr 2+ , Pb 2+ , Ba + , Zn 2+ , and Mg + in their carbonates of rhodochrosite (MnCO 3 ), otavite (CdCO 3 ), strontianite (SrCO 3 ), cerussite (PbCO 3 ), witherite (BaCO 3 ), hydrozincite (Zn 5 (CO 3 ) 2 (OH) 6 ) and hydromagnesite (Mg 5 (CO 3 )4 (OH) 2 •4H 2 O) [82].Generally, the discrepancies in results among such recent studies could be attributed to the differences in microbial remediation patterns (i.e., diverse metabolic behaviors, versatile microbial response and metal selectivity strategy) toward different metals during MICP [83].Interestingly, multiple concurrent immobilization mechanisms such as membrane surface ion exchange, adsorption, covalent binding, non-specific binding and particulate entrapment are possibly triggered along with calcite precipitation as referred by Qiao et al. [79].

Fourier transform infrared spectroscopy (FTIR)
FTIR spectroscopy is an expedient tool in providing important information about the structure, functional groups associated with any examined molecule and also differentiates between different carbonate polymorphs as well.Figure 7 manifested the existence of common bands in the rejoin of 3000-3800 cm −1 that could be ascribed to stretching vibrations of O-H groups of the adsorbed water molecule as refereed by Kadir [84].Besides, the vibration bands around 2357 cm −1 could be ascribed to carbon dioxide in the atmosphere [85].Meanwhile, the, the vibration band related to the stretching C-H functional group was found at 2919 and 2851 cm −1 as referred by Hamedi et al. [86].In addition, the peaks at the rejoin of 2396 to 2100 cm −1 might concern carbon dioxide in the atmosphere [85].However, peaks at 1633 and 1639 cm −1 would be ascribed to amide I signature of proteins as described by Saracho et al. [87].A same finding was recorded by Rodriguez-Navarro et al. [88], who referred to the long-term stability of CaCO 3 polymorphs by the virtue of tight binding between amide groups and CaCO 3 molecules.In addition, the singe for the combination of the main vibration frequencies of CO 3 2− assigned between ν1, (symmetric stretching) and ν4 (in-plane bending) was implied from the spectral peak at 1743 cm −1 [89].Interestingly, our results matched those obtained by Zain and Kadir [84].Also, the characterstic peaks of vaterite were detected at IR frequencies of 745, 870, 1087 and 1435 cm −1 , which are assigned to in-plane bending vibration (ν4), out-of plane bending modes (ν2), symmetric stretching (ν1) and asymmetric stretching vibration (ν3) of CO 3 2− , respectively [84].Notably, the shifting in the typical vaterite bands in the bioremediated samples was observed; implying the incorporation of the metals, ionic bond breaking and rearrangement process inside vaterite lattice.In addition, the appearance of extra peaks at 516, 576 and 670 cm −1 revealed the presence of metal oxide absorption bands.Generally, any peaks at the range of 400-700 cm −1 highlights the existence of metals/metal oxides as stated by Hassan et al. [90].Eventually, the tight association of fungal biomolecules (e.g., proteins, extracellular polysaccharide, glycoproteins, phospholipids, nucleic acids, etc.) with carbonate structure furnished the biominerals with higher stability and lower solubility, which in turn block the release of sequestered metals back to the remediated environment.As denoted by Li et al. [74] the biogenically synthesized CaCO 3 crystals exhibited a potent stability and less solubility than that formed under abiotic environments.
Collectively, based on all previous data, the mechanism of Zn 2+ and Cr 6+ removal mediated by CA of Rhodotorula sp.MZ312359 in MICP process could be deduced.Initially, the growth and proliferation of fungal cell was continued with oxidation of carbon sources (i.e., glucose and acetate) under oxic conditions; generating herby CO 2 (Equations.8 and 9).At this time, the interconversion of CO 2 implements by the catalysis of CA enzyme; generating dissociated ions ) (Eq. 10) [91].This continuous process would in turn elevate pH that eventually favors CaCO 3 precipitation in the presence of Ca 2+ ions through three-stage process of nucleation, growth and crystallization.Wherein, the fungal cells themselves serving as nucleation sites by the dint of their electronegativity nature.The negatively charged macromolecules such as polysaccharide, glycoproteins, lipids, proteins and lipopolysaccharide contain functional groups (e.g., carboxylic (R-COO − ), sulfonate (R-SO 3 ), and phosphate (R-PO 4 2− ), etc.) that attract the cations via attractive van der Waals forces [92].Once the nuclei of CaCO 3 precipitated on the fungal cell, in a supersaturated solution, crystal development begins via atom-by-atom addition during the second stage of crystal growth; generating larger size particles of CaCO 3 deposits (Eqs.11-13).Upon continuous fungal metabolic activity, which was synchronized with the production of CO 3 − and Ca 2+ ions and pH raising, the maturation of CaCO 3 crystals executed in the third crystallization stage.
In the concern of bioremediation process, it was accomplished within a proper time frame; producing stable solid phase, in spite of the absence of metal in the form of metal-CO 3 phase.Wherein, two suggestions governed the MICP-process driven by the examined fungus of the current study.The first one assumes both metals were remediated through the substitution of Ca 2+ with Zn 2+ and Cr 6+ , in particular with the presence of unconsumed percentage of soluble Ca 2+ that assessed by 21.01 and 14.94% for Zn 2+ and Cr 6+ , respectively, comparing to the biotic control that precipitated the soluble Ca 2+ entirely till the end of remediation experiment.In a convergent finding, Eltarahony [40] reported that Proteus mirabilis 10B eliminated Pb 2+ and Hg 2+ through Ca 2+ substitution in MICP process mediated by nitrate reductase enzyme.It is plausible to state that ionic exchange process between cations is controlled by several parameters for instance ionic radius, electronegativity, metals hydrolysis constant, ionic radius, electronegativity, atomic radius and hydrated radius [25].Remarkably, Zhu and Dittrich [25], Chada et al. [93] revealed that the divalent cations (Pb 2+ , Cd 2+ , Sr 2+ , Co 2+ and Zn 2+ ) have ionic radii that are closed to Ca 2+ , which subsequently enable their replacement and inclusion in calcite crystals.In the same extent, the removal of toxic Pb 2+ and Cr 6+ by Penicillium chrysogenum in MICP was achieved proposing their replacement with CO 3 − anion in calcite lattice [81].
The another suggestion depends on the characterstic properties of biogenic vaterite including its flexible structure with considerable surface area and sufficient porosity, which in turn boosted the incorporation of Zn 2+ and Cr 6+ inside vaterite lattice.Notable, [81] recorded the same observation.However, the electronegativity of vaterite also facilitated the gradual capturing of Zn 2+ and Cr 6+ during the development and growth stages of vaterite crystals; signifying the absorptive with sequestering capacity of vaterite.Interestingly, [75] found the same feature during the removal of Cr 6+ by chemical precipitation process.In this regard, Sdiri and Higashi [94] manifested the effective utilization of natural limestones in eliminating 10 ppm of Pb 2+ within 6 h.In the same extend, Du et al. [95], employed mollusk shells that mainly composed of calcite and aragonite as an efficient biosorbent for remediating divalent Pb, Cd and Zn.Broadly, all these scholars ascertained the biosorption capacity of CaCO 3 , unveiling the promising and potent performance of either MICP process or even its byproduct in purifying contaminated ecosystem.
Ultimately, the mineralogical characterization techniques assured the engulfment of Zn 2+ /Cr 6+ in vaterite co-precipitated remediated products as possible action of carbonic anhydrase.Moreover, they reflected not only the scavenging role of CaCO 3 in chelating soluble heavy metals but also their stabilization inside potent trap, which displayed superior efficiency than adsorption [93].Arguably, the adsorption is superficial process and based mainly on bonding specifications between adsorbents and adsorbate, which substantially depends on solution ionic strength, pH and availability of sufficient functional groups [96].Strikingly, the precipitation approach in hydroxides forms was extensively and efficiently studied in heavy metals removal, nevertheless, the tendency to form soluble anionic hydroxyl complexes through water redissolution considered being the main drawback [97,98].In general, the remediation of heavy metals through MICP deemed as effective, easy, ecofriendly and inexpensive procedure to restrict toxic metals mobility, bioavailability and their release to the surrounding milieu with facilitated separation and without additional step like coagulation, filtration, or flocculation.

Conclusion
In conclusion, the current study dedicated on maximizing carbonatogenic process, which was implemented by the catalysis of CA enzyme of Rhodotorula sp.MZ 312359 via employing statitical experimental designs (i.e., PBD and CCD) to remeidate Zn 2+ and Cr 6+ in MICP, for the first time.Such statistical approaches enhanced carbonatogenic parameters of CA activity and also CaCO 3 weight by 1.8-fold, comparing to un-optimized basal media.Under the optimized precipitating conditions, Rhodotorula sp.MZ 312359 eliminated entirely 50 and 400 ppm of Zn 2+ and Cr 6+ , respectively within 7 days of incubation.Subsequently, the mineralogical analysis including EDX, SEM, FTIR and XRD confirmed the immobilization of soluble toxic metals inside a potent trap of vaterite lattice, which deemed as a promising strategy for cost-effective, environmentally benign, sustainable and effective mean for alleviation heavy metals toxicity.

Table 1
Matrix of Plackett-Burman experimental design, variables and their levels along with experimental and predicted response values for screening of significant processes variables affecting CaCO 3 weight and CA activity

Fig. 1
Fig. 1 Neighbor-joining dendrogram (A), cultural characteristics on YPD agar plate (B) and morphological features of Rhodotorula sp.MZ312359 cells as examined by SEM analysis

Fig. 3
Fig. 3 Dynamic analysis of carbonatogenesis process associated with changes in pH, cell growth, CA activity, CaCO 3 weight and the removal of Ca 2+ , Zn 2+ and Cr 6+ .A Biotic control, B Zn 2+ -remediation experiment, C Cr 6+ -remediation experiment.The results were expressed as mean ± SEM.To adjust the scale, some parameters are tripled in factor of 10 as indicated on the figure

Table 2
Central composite design matrix of CaCO 3 weight and CA activity along with the predicted responses and variables concentrations

Table 4
Estimated effects, regression coefficients and corresponding P-values for the first order equation model of CaCO 3 weight and CA activity optimized by PBD

and Coefficients for CaCO 3 weight Estimated Effects and Coefficients for CA Enzyme
equations representing the optimum of CA activity and CaCO 3 weight as a function of the studied independent factors were expressed as follows in Eqs. 4 and 5:

Central composite design for optimizing CA activity and CaCO 3 weight
[57,58]d variables on response[57,58].The goal of the RSM trials was to get a more accurate evaluation of the optimal operating conditions for the factors screened from PBD at a five-level (− 2, − 1, 0, + 1, + 2) to achieve the maximum CA activity and CaCO 3 weight.
RSM is a mean of statistics that depends on key statistical concepts, randomization, replication and duplication.

Table 5
Estimated effects, regression coefficients and corresponding P-values of second order polynomial model for CaCO 3 weight and CA activity optimized by CCD

Table 6
ANOVA for quadratic polynomial model of calcium carbonate weight and CA activity NaNO 2 + 0.00125 Ca − acetate * inoculum size + 0.035 NaNO 2 * inoculum size.