Role of Anthocleista vogelii in serum antioxidant defence system in cadmium-induced oxidative stress in Wistar rats

Abstract

Background

Cadmium (Cd) toxicity, which runs across the food chain, is chiefly regulated by in vivo antioxidant defence system or through antioxidant supplementation of biological systems predisposed to this environmental stressor. The present study was designed to examine the role of Anthocleista vogelii leaves in Cd-induced oxidative stress in the serum of Wistar rats through the application of response surface methodology (RSM) and biomonitoring of selective responses: malondialdehyde (MDA), catalase (CAT), superoxide dismutase (SOD), glutathione-s-transferase (GST) and peroxidase (POD) activities, respectively. The cold macerated plant leaves were subjected to fractionation process using methanol-hexane-chloroform (3:2:1 v/v) solvent system such that n-hexane fraction with ample antioxidant levels in terms of total phenolic content (TPC) and total flavonoid content (TFC) among others at p <  0.05 was selected for the study. The study employed central composite design (CCD) with twenty experimental “runs” of male Wistar rats for twenty-eight days, following a week of acclimatization, where n-hexane fraction of A. vogelii (NFAV), cadmium chloride (CdCl₂) and body weights of rats were considered input factors in the study.

Results

The study generated five quadratic models, which differed significantly at p <  0.05 for MDA levels as well as CAT, SOD, GST and POD activities in the sera of Wistar rats. The study revealed that exposure to Cd toxicity caused a marked increase (p <  0.05) in serum MDA levels, but a significant inhibition (p <  0.05) of serum SOD, CAT, GST and POD activities. However, Cd interaction with NFAV showed marked amelioration of Cd-induced oxidative stress, which was confirmed by significant decrease in serum MDA levels, but significant increase in serum SOD, CAT, GST and POD activities at p <  0.05 via the response surface plots. The study also confirmed the reliability and adequacy of the models for accurate prediction of the responses since R-squared (R²) values obtained were greater than 90%.

Conclusion

It was inferred from the present study that the adequacy of the models validated the potency of A. vogelii leaves graphically in the amelioration of Cd-induced oxidative stress in the serum of Wistar rats. Hence, the plant was considered a rich source of bioactive compounds with significant antioxidant properties.

1 Background

The post-discovery of Cadmium (Cd) in 1817 by a German chemist “Friedrich Stromeyer” raised environmental concerns about the toxic nature of the metal on biological systems, which was increased by industrialization and other anthropogenic activities. Cd is a confirmed environmental toxicant known to compromise the defence line of antioxidants in favour of pro-oxidation of biological systems by reactive oxygen species (ROS) due to overexposure. These species initiate lipid peroxidation of biomolecules, especially polyunsaturated fatty acids (PUFAs), and with subsequent loss of biological functions of membrane [1,2,3]. The toxicity initiated by Cd was reported to be associated with its protracted biological half-life, the potential capacity to interact and replace some metallic cofactors in their enzyme forms and indirect induction of free radicals via Fenton reaction [2,3,4]. The description of Cd toxicity is historically a function of its relative abundance in soil, which is brought about by transportation and infiltration along the food chain. Food, being a primary source of Cd exposure to consumers, shows that upon ingestion, it is absorbed through the gastrointestinal tract (GIT), transported and distributed through the blood to other tissues, where it bioaccumulates and elicits deleterious effects. Reports have shown that Cd impairs cellular structure, with consequential onsets of diverse pathologies, which include renal dysfunction, liver problems, erectile dysfunction and many more [2, 5,6,7,8]. The onset of these pathologies in tissues may be recognized by high expression of metallothionein (MT) and synthesis of other sulphur-rich proteins. In light of the aforesaid, reports have shown that liver and kidney tissues are fundamental centres of Cd toxicity [8,9,10].

The fundamental cellular mechanism of Cd toxicity in biological systems can be traced to several factors. Cuypers et al. [9] reported that Cd has a high affinity for thiol-rich molecules, especially enzymes or structural proteins, with implicative metabolic alterations. Again, zinc (Zn), Iron (Fe) and other cationic micronutrients, which are bivalent in nature are considered essential cofactors of metalloenzymes thereby competing with a bivalent Cd metal under physiological setting by redox displacement reaction or interference with the uptake of these nutrients so as to elicit its toxic action on tissues [2]. Another pivotal mechanism of Cd-induced damage in biological systems is oxidative stress, a phenomenon characterized by excessive synthesis of reactive oxygen species (ROS) with a consequential alteration on the prooxidative-antioxidative pool thereby favouring the prooxidants, and resulting in cellular damage and signal transduction [2]. The excessive generation of ROS triggers the formation of malondialdehyde (MDA), a biomarker of oxidative stress, through lipid peroxidation process with deleterious onsets of diverse pathologies [11].

Cd is considered as a non-fenton metal, owing to the fact that it lacks the direct capacity to initiate synthesis of reactive oxygen species (ROS). Cd-enhanced synthesis of ROS during exposure was reiterated in literature to be associated with Fenton mechanism. Thus, over-exposure to this toxicant may result in the overall compromise of the antioxidant defence system through excessive production of ROS with aftermath effect of oxidative stress on biological systems [2, 3, 12]. The compromise may be a consequence of alteration of enzymatic antioxidant activities during a prolonged exposure, which includes the distortion of catalase (CAT), glutathione-s-transferase (GST), superoxide dismutase (SOD), glutathione reductase (GR) and glutathione peroxidase (GPX) activities respectively, as reported in literature [13, 14]. The responsibility of empirically based improvement of Cd-induced oxidative stress may be associated with some secondary metabolites with antioxidant properties inherent in plants. The antioxidative function of these intrinsic metabolites of plants can be actualized, if their concentrations are adequate to prevent and downregulate undue synthesis of free radicals, for the enhancement of in vivo antioxidant defence system [15]. However, much is yet to be known about plant-rich metabolites with antioxidant properties, their synthesized amounts and mechanisms of actions respectively. Reports have empirically supported the numeric estimation and qualitative presence of phenolics and flavonoids with antioxidant properties in Anthocleista vogelii [16,17,18].

A. vogelii, a known genus among fifteen acknowledged forms of anthocleista genera, is reported to have traditional usefulness in the management of various pathologies such as sexually transmitted disease (STDs), metabolic disorders, hypertension, infertility and many more [19, 20]. It belongs to the Gentaniaceae family and is popularly called “Cabbage tree”. The ethnobotanical description of A. vogelii involves 6–20 m height, 0.15–0.55 m stem width, 0.4–1.5 m leaf length with 0.24–0.45 m width [21]. The study carried out by Anyanwu et al. [21] showed the significant role played by the plant in obese Wistar rats. The results succinctly showed a marked reduction in body weights, aspartate aminotransferase (AST), alanine aminotransferase (ALT) and alkaline phosphatase (ALP) activities as well as a marked significant increase in antioxidant status of obese animal models relative to control. Anyanwu and his colleagues hereby inferred that A. vogelii exhibits high antioxidant efficacy during supplementation, which was attributed to its rich, intrinsic phenolic and flavonoid contents respectively. To this end, the rich antioxidant activity of A. vogelii regardless of the aerial parts was reported to be characterized, and this may be attributed to intrinsic bioactive compounds with free hydroxyl groups notably; 1-hydroxy-3,7-dimethoxyxanthone, 1-hydroxy-3,7,8-trimethoxyxanthone, 1,8-dihydroxy-7,8-dimethoxyxanthone, 7α-hydroxysitosterol, stigmasterol, secologanin, and sitosterol-3-O-β-glucopyranoside [19, 22, 23].

Therefore, the novelty of the present study is aimed at identifying the potentiating antioxidant role of the plant in ameliorating Cd-induced oxidative stress in the serum of male Wistar rats through the application of response surface methodology (RSM). The application of the research tool (RSM) would help in identifying the antioxidant role of the plant as well as the generation of quadratic models for useful and accurate prediction of responses employed in the biomonitoring of Cd-induced oxidative stress in the serum of Wistar rats.

2 Methods

2.1 Chemical reagents

All chemical reagents used in the present study were analytical grades, and they were purchased from Sigma Aldrich, which is presently known as “Merck”. The chemicals include: anhydrous aluminium trichloride (AlCl₃), ascorbic acid, iron (iii) chloride (FeCl₃), cadmium chloride (CdCl₂), chloroform, folin-ciocalteu reagent, gallic acid, methanol, petroleum ether (40–60 °C), potassium ethanoate (CH₃COOK), rutin, sodium trioxocarbonate (iv) (Na₂CO₃), trichloroacetic acid, thiobarbiturate, hydrogen peroxide (H₂O₂), adrenaline, chloro-2,4-dinitrobenzene (CDNB), glutathione (GSH), pyrogallol, dipotassium hydrogen tetraoxophosphate (v) (K₂HPO₄) and potassium dihydrogen tetraoxophosphate (v) (KH₂PO₄) respectively.

2.2 Experimental animal models

Twenty male Wistar rats with average weight of 150.00 ± 8.11 g were procured from the animal house unit of Biochemistry Department of the University of Benin, Benin City, Nigeria, and housed in metabolic cages for the investigational period. They were allowed to acclimatize for a week, fed with pelleted livestock feeds, free access to water ad libitum and 12–12 h light/darkness cycle respectively. Thereafter, the animals were categorized into 20 experimental runs by means of central composite design (CCD) with two replications in a total of 60 experimental rats, and subjected to treatment for 28 days. The use of rat models in the present study was subject to approval by the ethical committee of the Delta State University, Abraka, Nigeria, responsible for the laboratory principles of handling and care of experimental animals in accordance with the European Community Act of 1986.

2.3 The plant of experimental interest

After harvest from the environment of Western Delta University, Oghara, Nigeria, the plant was identified and authenticated as “Anthocleista vogelii” at the Department of Botany, University of Benin, Benin City, Nigeria, by Dr. H.A. Akinnibosun. A sample of the plant was deposited in his herbarium with a voucher designation, “UBHa0258”.

2.4 Preparation of crude sample and its fractionation

Exactly 3.14 kg were dried at room temperature 27.0 ± 2.0 °C, milled into fine powdery form of uniform sizes by means of a warring mechanical blender, filtered and extracted by cold maceration technique for 4 days using 70% methanol. The extraction process was consecutively followed in two repetitions for absolute yield of 34.09% (1070.42 g). This was obtained by means of a rotary evaporator and stored at 4 °C in air-tight container for further use. Furthermore, 230.0 g of crude sample was fractionated by thin-layer chromatography (TLC) aided by a precoated silica gel plate and vacuum liquid chromatography (VLC). The fatty component of the extract was removed using petroleum ether (40–60 °C) type and the impenetrable portion was fractionated using methanol-hexane-chloroform, (3:2:1 v/v), solvent system. Upon two repetitions, the respective fractions were further vacuum-dried to a thick brownish-green sticky paste using a rotary evaporator and stored at 4 °C for further use. Apiamu et al. [16] revealed qualitatively the presence of phenolics and flavonoids in the leaves of A. vogelii. In light of this, the bioactive compounds with antioxidant potentials were comparatively quantified in the respective fractions to identifying the fraction required for the in vivo study. The scheme shown in Fig. 1 explains the extraction process of A. vogelii leaves for the study.

Fig. 1

The order of extraction and fractionation of A. vogelii leaves

2.5 Evaluation of in vitro total phenolic content (TPC)

This was estimated following the procedure reported by Ibrahim et al. [24] with minor changes. A 6.5-ml mixture containing 1.0 ml of test sample, 0.5 ml of Folin-ciocalteu reagent and 5 ml of double distilled water (ddH₂O) was incubated for 5 min at room temperature. To the mixture, 1.0 ml of 0.005 M Na₂CO₃ was added and subjected to 60 min incubation in the dark. Similarly, the assay procedure was done in each case, using 50% ethanol and gallic acid (GA) as negative and positive controls, and the respective absorbances were read in triplicate determinations with UV-VIS spectrophotometer double beam (Labtech-2802). The intrinsic TPC of test sample (A. vogelii) was assessed in relation to GA standard calibration curve (y = 0.0018x + 0.4562; R² = 0.9168, 0–1000 μg/ml) and expressed as mg GAE/g of dried sample.

2.6 Evaluation of in vitro total flavonoid content (TFC)

This was estimated following the procedure reported by Chang et al. [25] with minor changes. A 5.0-ml mixture containing 0.5 ml of test sample, 0.1 ml of AlCl₃ (10% w/v), 0.1 ml of 1.0 M CH₃COOk solution and 4.3 ml of ddH₂O was subjected to 30 min incubation. The assay procedure was repeated with 50% ethanol and rutin (RU) as negative and positive controls, where triplicate absorbances were read in each case using UV-VIS spectrophotometer double beam (Labtech-2802). The intrinsic TFC of test sample (A. vogelii) was assessed in relation to RU standard calibration curve (y = 0.1255x + 0.0067; R² = 0.966, 0–1000 μg/ml) and expressed as mg RUE/g of dried sample.

2.7 Experimental design using RSM

Prior to the design, the study recognized the report of Ayakeme et al. [26] and Virk et al. [27], stating that 15.0 mg CdCl₂ with reference to body weight of rats, when administered orally, resulted in the induction of oxidative stress. Again, the design measured three input factors or process variables (body weight of rats, plant extract and Cd dose rates respectively) and response variables such as serum MDA levels and serum enzymatic antioxidant: catalase (CAT), superoxide dismutase (SOD), glutathione-s-transferase (GST) and peroxidase (POD) activities respectively (Table 1). The RSM-CCD was employed in the design of the experiment with 20 sets of experimental runs (N), defined as; N = 2^k + 2^k + n_o, where k reflects the number of process variables and repetitions with eight factorial, six axial and six centre points. Thus, the general regression model for a second-degree polynomial is defined hereunder.

$$ Y={\alpha}_0+{\alpha}_1\mathrm{A}+{\alpha}_2\mathrm{B}+{\alpha}_3\mathrm{C}+{\alpha}_{11}{\mathrm{A}}^2+{\alpha}_{12}{\mathrm{B}}^2+{\alpha}_{13}{\mathrm{C}}^2+{\alpha}_{21}\mathrm{AB}+{\upalpha}_{22}\mathrm{AC}+{\upalpha}_{23}\mathrm{BC} $$

(1)

Table 1 Quantification and TLC analysis of bioactive compounds with antioxidant properties in A. vogelii

Such that; Y designates a response variable, α₀, α_1–3, α_11–13 and α_21–23 depict fixed intercept term, fixed linear coefficients, fixed-squared coefficients and fixed coefficients of interactions of the regression model. In the study, A (plant extract dose rate), B (Cd dose rate) and C (body weight of rats) represent the process variables, and a mixture of this sort; AB, AC and BC is an indication of significant interaction of the variables.

2.8 Optimization of enzymatic antioxidant defence system

The optimal activities of antioxidant enzymes in serum were actualized through the application of numerical optimization approach (NOA) on the experimental variables. Thus, plant extract and Cd dose rates were set at 200 mg/kg and 15 mg/kg, while the body weights of rats were left in range (100–200 g). The optimization process functions primarily to identifying the significant role played by A. vogelii in the improvement of Cd-induced oxidative stress in Wistar rats via minimization of MDA levels.

2.9 Experimental sampling

The rats were anaesthetized using chloroform in a desiccator, and whole blood samples were collected by heart puncture into EDTA bottles using hypodermal syringes and centrifuged at 2500 g for 15 min. The serum obtained from the centrifugation process was properly refrigerated at 4 °C for further use.

2.10 Biochemical analysis

Oxidative stress is monitored through the measurement of malondialdehyde (MDA) levels formed during lipid peroxidation process, following the assay procedure described by Iqbal et al. [28]. The assay procedure described by Aebi [29] and Luck [30] was employed in the assessment of serum CAT activity on the basis of H₂O₂ breakdown into H₂O and O₂, respectively. Serum SOD activity was evaluated following the assay method described by Misra and Fridovich [31] via the formation of adrenochrome from auto-oxidation of adrenaline. Serum GST activity was analysed following the method described by Habig et al. [32] through complexation of GSH with Chloro - 2, 4-dinitrobenzene (CDNB) with slight alteration. The assay method described by Reddy et al., [33] was employed in the analysis of serum POD activity and serum total protein was evaluated by the method of Tietz [34].

2.11 Empirical data analysis

The Design-Expert Software 10.0.7 (Stat Ease Inc., Minneapolis, USA) was employed in the design of the present study, mathematical modelling, graph plottings, optimization process and statistical analysis. In each case, the quality of developed models was evaluated through analysis of variance (ANOVA) embedded in the software, regression coefficient (R² value), predicted R² and adjusted R² values, respectively. The validation of optimum conditions was actualized using the software with consideration of mean values of triplicate result for respective responses, where models at p <  0.05 were considered significant.

3 Results

3.1 The choice of fraction for RSM-CCD study

The reduction of in vivo antioxidants’ activities during over-exposure to Cd toxicity may result in the compromise of their defensive role; hence, in vitro supplementation strategy (ISS) may be suggestive during non-affirmative exposure of biological systems to this environmental stressor. Here, the crude leaf extract, methanol and hexane fractions assessed differ significantly (p <  0.05) in relation to their TPC and TFC. However, the chloroform fraction was discarded since it qualitatively showed the absence of phenolic and flavonoid compounds. On the basis of their content, the n-hexane fraction showed an ample amount of TPC and TFC with marked significant difference in relation to the other fractions (Table 1). Therefore, the n-hexane fraction was selected for the in vivo study since researchers reported that polyphenols are invaluable bioactive compounds exhibiting antioxidant activity and capacity to chelate metal ions with toxicity tendencies [23].

3.2 Computation of mean results of responses for modelling

On the basis of the RSM-CCD, the mean values of responses, which include serum biomarker of oxidative stress (MDA levels) and enzymatic antioxidant (CAT, SOD, GST and POD) activities were computed and subjected to RSM analysis for modelling and optimization processes (Table 2). The analysis revealed that there were no outliers in the mean values of the respective responses.

Table 2 Experimental RSM-CCD and results of responses from the present study

3.3 Model fitting

The application of the RSM analysis generated response quadratic models for biomarker of oxidative stress (MDA levels) and enzymatic antioxidant (CAT, SOD, GST and POD) activities in the serum of Wistar rats under combination of twenty experimental runs with variable experimental conditions, as indicated in Eqs. (2)–(6), to evaluate the relationship between process variables (A-plant extract at 0–200 mg/kg b.wt of rats, B-CdCl₂ at 0–15 mg/kg b.wt of rats, and C-body weight changes of rats) and each response variable with the elimination of non-significant terms;

$$ {Y}_{\mathrm{MDA}\ \left(\mathrm{mmoles}/\mathrm{mg}\ \mathrm{protein}\right)}=+37.63-19.18\mathrm{A}+12.49\mathrm{B}-9.16\mathrm{AB}+3.91{\mathrm{C}}^2 $$

(2)

$$ {Y}_{\mathrm{CAT}\ \left(\mathrm{units}/\mathrm{mg}\ \mathrm{protein}\right)}=+36.34+8.47\mathrm{A}-29.72\mathrm{B}+17.21{\mathrm{B}}^2 $$

(3)

$$ {Y}_{\mathrm{SOD}\ \left(\mathrm{units}/\mathrm{mg}\ \mathrm{protein}\right)}=+53.97+23.18\mathrm{A}+6.73\mathrm{B}-7.36\mathrm{AB}-5.19\mathrm{BC}-3.46{\mathrm{C}}^2 $$

(4)

$$ {Y}_{\mathrm{GST}\ \left(\mathrm{units}/\mathrm{mg}\ \mathrm{protein}\right)}=+36.04+7.31\mathrm{A}+5.67\mathrm{AB}-11.99{\mathrm{B}}^2 $$

(5)

$$ {Y}_{\mathrm{POD}\ \left(\mathrm{units}/\mathrm{mg}\ \mathrm{protein}\right)}=+62.54+16.37\mathrm{A}-34.03\mathrm{B}+29.18{\mathrm{A}}^2 $$

(6)

In Table 3, the analysis of variance (ANOVA) for serum MDA levels showed that the quadratic model differs significantly (F = 105.48, p <  0.0001) and the lack of fit test was considered not significant (p = 0.9121), indicating the suitability of the model for adequate prediction. Also, the ANOVA for serum CAT, SOD, GST and POD activity indicated quadratic models with marked significant differences at high F and low p values, and the insignificant lack of fit tests for these responses further confirmed the reliability of the models for useful predictions during Cd-induced oxidative stress, where A. vogelii played a fundamental role of an antioxidant. Again, the regression coefficient (R²) of 0.9916, 0.9802, 0.9712, 0.9290 and 0.9921 for serum MDA levels, CAT, SOD, GST and POD activities suggested that 99.16, 98.02, 97.12, 92.90 and 99.21% of the variations observed in experimental data could be accounted for by the quadratic models.

Table 3 Analysis of variance (ANOVA) for serum oxidative-antioxidative indicators in relation to input variables

3.4 Graphical analysis

The adequacy and precision of the quadratic models for the responses were also checked by diagnostic (predicted vs. actual) plots. The proximity of the actual and predicted values to the regression line indicates the reliability of the models (Fig. 2). Consequently, the high correlation could be explained by the high R² values as well as the closeness of the predicted and adjusted R² values associated with the generated models in Table 3.

Fig. 2

Diagnostics plots of actual vs. predicted values for a serum MDA levels, b CAT activity, c SOD activity, d GST activity, and e POD activity

The 3-D response surface plots of different shapes generated by the models were employed to describe the relationship between the process variables and the responses under these experimental conditions (Figs. 3, 4, 5, 6 and 7). In each case, two process variables with varying degrees of interactive effects on the responses were considered at fixed state of the third process variable. Thus, increasing dose rate of CdCl₂ showed a significant induction (p < 0.05) of oxidative stress, characterized by increasing levels of MDA in serum of Wistar rats (Fig. 3). However, the interaction between CdCl₂ and A. vogelii (n-hexane fraction) extract indicated a significant improvement on Cd-induced oxidative stress in the serum of Wistar rats at p < 0.05. Comparatively, the interaction between body weights of rats and plant extract or CdCl₂ showed similar effects on serum MDA levels with marked significant difference (p < 0.05) characterized by rising contour bridges. Also, the study revealed a significant inhibitory effect on the activities of serum enzymatic antioxidants (CAT, SOD, GST and POD) during exposure to Cd toxicity, but a co-treatment with plant extract revealed a significant increase (p < 0.05) in the activities of antioxidant enzymes, which were evidently established by the rising contour bridges in the response surface graphs, describing the extent of interaction between the process variables (Figs. 3, 4, 5, 6 and 7).

Fig. 3

3-D surface plots for serum MDA levels under the influences of process variables

Fig. 4

3-D surface plots for serum CAT activity under the influence of process variables

Fig. 5

3-D surface plots for serum SOD activity under the influence of process variables

Fig. 6

3-D surface plots for serum GST activity under the influence of process variables

Fig. 7

3-D surface plots for serum POD activity under the influence of process variables

4 Discussion

4.1 Appraisal of RSM-CCD in the present study using n-hexane fraction for the treatment of Cd-induced oxidative stress in Wistar rats

The protracted environmental hazards initiated by cadmium (Cd) compounds on biological systems, in view of several gradations of exposure, are inevitably encompassing in the areas of toxicological investigations. It was then reiterated in literature that Cd readily accumulates along the food chain to toxic levels with detrimental effects on biological systems as a result of its prolonged biological half-life, pervasive and soluble nature [3, 35]. To this end, Cd was considered a known toxicant, whose underlying mechanism of absorption and cellular toxicity was principally reported as a concern to researchers in view of its deleterious effects to microorganisms, plants, and animals [2, 11, 36]. Although, in vivo enzymatic and non- enzymatic antioxidants may offer protection against Cd-induced oxidative stress through a mop-up of excess ROS, but over-exposure to this toxicant was reported to compromise the antioxidant defence system, which seemingly supports lipid peroxidation process in biological systems [10, 37]. Therefore, the use of medicinal plants rich in bioactive compounds with antioxidant activity were highly suggested as an alternative route of amelioration of undue generation of ROS in living organisms by exposure to Cd toxicity [15]. In the present study, this underscored the need for the use of A. vogelii with antioxidant potentials in the amelioration of Cd-induced oxidative stress.

The principal focus of toxicologists in the field of life sciences encompasses the amelioration and/or prevention of excessive generation of ROS brought about by a toxicant and subsequent optimization of the antioxidant defence system of vulnerable biological systems. To actualize this objective in the present study with adequate precision and accuracy, a set of synchronized experimental designs and analyses termed “Response Surface Methodology, (RSM)” encrypted with central composite design (CCD) was employed. The critical statistical analysis, which was done in terms of analysis of variance (ANOVA), indicated quadratic models for serum responses (Table 3). The validation of these regression models for adequate prediction of responses was reported to be a function of its regression coefficient (R²), lack of fit testing and other instruments of measurement [38]. Thus, serum MDA, CAT, SOD, GST and POD quadratic models presented R² values of 0.9916, 0.9896, 0.9712, 0.9290 and 0.9921 with insignificant lack of fit at 95% confidence interval (i.e. p > 0.05), which means that the model, in each case, consequently accounts for 99.16, 98.96, 97.12, 92.90 and 99.21% of data generated for the responses in Table 3. The analysis indicated significant quadratic models for the respective responses (< 0.0001), and adequate to denoting the mathematical relationship between variables and to measuring the degree of variability in response variables accounted for by the experimental process variables and their specific interactions [39]. The quadratic models were acknowledged good and adequate for the prediction of the responses, owing to their high R² values, low p values and high F values, respectively [38, 40]. It was reported that the magnitude of adjusted and predicted R² values can be applied as a measure of proximity of experimental and predicted values [40, 41]. In the present study, the adjusted and predicted R² values for the respective responses were in close conformity, with indication of high correlation between the experimental and predicted data (Table 3). Thus, it was evidently shown in the diagnostic plots of serum MDA levels, CAT, SOD, GST and POD activities (Fig. 2). Kanmani et al. [40] and Ramanan et al. [42] reported the numerical importance of adequate precision to determining which experimental variable generates larger signal in relation to noise during the experiment. It was then inferred that adequate precision of experimental responses greater than four was considered beneficial. Therefore, the adequate precisions; 36.692, 33.136, 24.512, 10.507 and 39.045 for MDA, CAT, SOD, GST and POD quadratic models were suggestive of sufficient signal to noise ratio for adequate prediction (Table 3).

The statistical appraisal of mathematical models is a function of significance and involvement of each process variable using p values [43]. Also, it was stated that the bigger and smaller the magnitudes of F and p values, the more significant is the corresponding coefficients of process variable [44, 45]. Thus, this necessitated the statistical considerations of the linear, square and interaction coefficients of the respective process variable terms in the modelling course, where the insignificant terms (p > 0.05) were eliminated (Eqs. 2–6). The regression models of these sorts were typified graphically in terms of 3-D response surface plots to mean interaction between two process variables at a given time with fixed level of another process variable, assessment of optimal levels and their effects on the response variables [40, 46, 47]. Again, experimental studies on contour lines revealed that it might be saddle, rising bridges or elliptical mounds in nature [40, 48]. Thus, the plot of Cd dose rate versus plant extract (Fig. 3) is a rising bridge, indicating moderate significant interaction between the two process variables. Optimum MDA levels, as biomarker of oxidative stress, was observed from the plot at 15 mg/kg Cd dose rate and zero threshold of the plant extract, which validated the reports report of Ayakeme et al. [26] and Virk et al. [27] that oral administration of 15 mgCdCl₂/kg b.wt of rats induced oxidative stress. The reduction of MDA levels became significant (p < 0.05) at maximum dose rate of plant extract (200 mg/kg). Therefore, the graphical description of response surface plots shown in Fig. 3 identified the effects of the process variables on the biomarker of oxidative stress. Rao et al. [47] reported that the curvature of response surfaces is a function of well-defined optimal conditions for the variables, but not much variation between single parameters and optimized values was obtained when the response surfaces are flat and symmetrical in proximity to optimal levels. Figures 3, 4, 5, 6 and 7 showed the rising bridges of interactions between process variables, especially Cd and plant extract, on enzymatic antioxidant activities. The inhibitory effect on serum CAT activity was significant at 12.0 mg/kg of Cd, but was alleviated through significant increase (p < 0.05) in the dose rate of leaf extract from 100 to 200 mg/kg (Fig. 4). The rising bridges or saddle points obtained for serum SOD activity indicated that the enzyme activity was inhibited using 6–15 mg/kg CdCl₂ for 28 days, which was absolutely reversed in the presence of the plant extract at peak dose rate of 200 mg/kg (Fig. 5). There was significant interaction between Cd and plant extract, owing to the nature of the contour plots and resulted in maximum serum GST activity. Maximum GST activity (36.03 mU/mg protein) was obtained for 9.0 mg Cd/kg b.wt and 200 mg/kg b.wt of plant extract (Fig. 6). Serum POD activity was significantly inhibited when Cd dose rate was increased from 6 to 15 mg/kg with minimum inhibitory effect at 12.0 mg/kg of Cd dose rate. However, the inhibitory effect was improved via a significant increase (p < 0.05) in the dose rate of leaf extract (100–200 mg/kg) (Fig. 7). The graphical outcomes of the present study conformed with previous reports, showing that exposure to Cd results in elevated levels of MDA, a derivative of lipid peroxidation and biomarker of oxidative damage [10, 49, 50]. Again, the results proposed that the leaf extract of A. vogelii ameliorated Cd-initiated oxidative stress through enhancement of enzymatic antioxidant defence system in male Wistar rats (Figs. 3, 4, 5, 6 and 7). This was attributed to the antagonistic effects brought about by intrinsic phenolic and flavonoid contents of the plant with antioxidant potentials, as characterized by Alaribe, et al. [22], Anyanwu et al. [19] and Valentão et al. [23], on Cd-induced oxidative damage in the serum of Wistar rats.

4.2 Confirmation of models

The essence of the verification study is to show that the developed fitted quadratic models can adequately and correctly be employed in the prediction of the responses’ functions, as it pertains to identifying the role of A. vogelii in Cd-induced oxidative stress in Wistar rats. Table 4 clearly elucidates that the application of design-expert-enhanced models predicted serum MDA levels to be 37.63 mmol/mg protein and serum CAT, SOD, GST and POD activities to be 36.34, 53.97, 36.04 and 62.54 units/mg protein, respectively. Therefore, the application of the set conditions in Table 4 in the laboratory for 28 days confirmed the reliability of the models for the prediction of the responses.

Table 4 Validation of response surface quadratic models

5 Conclusion

The culmination of the present study using RSM-CCD clearly shows successful modelling and optimization of experimental variables for accurate prediction of the responses, and the quadratic models enhanced by 3-D response surface plots established the substantial role played by A. vogelii leaves in the improvement of Cd-induced oxidative stress in the sera of Wistar rats, which validated the ample antioxidant property of the plant.