Next Article in Journal
Characterization of SOD and GPX Gene Families in the Soybeans in Response to Drought and Salinity Stresses
Next Article in Special Issue
Impact of Combined Heat and Salt Stresses on Tomato Plants—Insights into Nutrient Uptake and Redox Homeostasis
Previous Article in Journal
A Functional Orange Juice Fortified with Beetroot By-Products Attenuates Hyperlipidemia and Obesity Induced by A High-Fat Diet
Previous Article in Special Issue
Interactive Effects of Melatonin and Nitrogen Improve Drought Tolerance of Maize Seedlings by Regulating Growth and Physiochemical Attributes
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Antioxidants
Volume 11
Issue 3
10.3390/antiox11030458
antioxidants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Computational Study of the Role of Secondary Metabolites for Mitigation of Acid Soil Stress in Cereals Using Dehydroascorbate and Mono-Dehydroascorbate Reductases

by
Shuvasish Choudhury
1,2,*,†,
Muhammed Khairujjaman Mazumder
2,3,†,
Debojyoti Moulick
2,†,
Parul Sharma
4,†,
Sandeep Kumar Tata
5,†,
Dibakar Ghosh
6,
Hayssam M. Ali
7,
Manzer H. Siddiqui
7,
Marian Brestic
8,
Milan Skalicky
9 and
Akbar Hossain
10,*
1
Plant Stress Biology and Metabolomics Laboratory, Department of Life Science and Bioinformatics, Assam University, Silchar 788011, India
2
Central Instrumentation Laboratory, Assam University, Silchar 788011, India
3
Department of Zoology, Dhemaji College, Dhemaji 787057, India
4
Don Bosco School, Silchar 788003, India
5
Department of Botany, Balmiki Rajniti Mahila College, Munger University, Munger 811201, India
6
Division of Agronomy, ICAR-Indian Institute of Water Management, Bhubaneswar 751023, India
7
Department of Botany and Microbiology, College of Science, King Saud University, Riyadh 11451, Saudi Arabia
8
Department of Plant Physiology, Slovak University of Agriculture, Nitra, Tr. A. Hlinku 2, 94901 Nitra, Slovakia
9
Department of Botany and Plant Physiology, Faculty of Agrobiology, Food, and Natural Resources, Czech University of Life Sciences Prague, Kamycka 129, 16500 Prague, Czech Republic
10
Department of Agronomy, Bangladesh Wheat and Maize Research Institute, Dinajpur 5200, Bangladesh
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Antioxidants 2022, 11(3), 458; https://doi.org/10.3390/antiox11030458
Submission received: 28 January 2022 / Revised: 9 February 2022 / Accepted: 18 February 2022 / Published: 25 February 2022
(This article belongs to the Special Issue Antioxidant Mechanisms in Plants)
Browse Figures
Versions Notes

Abstract

:
The present study investigates the potential ameliorative role of seven secondary metabolites, viz., ascorbate (AsA), reduced glutathione (GSH), jasmonic acid (JA), salicylic acid (SA), serotonin (5-HT), indole–3–acetic acid (IAA) and gibberellic acid (GA3), for mitigation of aluminium (Al3+) and manganese (Mn2+) stress associated with acidic soils in rice, maize and wheat. The dehydroascorbate reductase (DHAR) and mono-dehydroascorbate reductase (MDHAR) of the cereals were used as model targets, and the analysis was performed using computational tools. Molecular docking approach was employed to evaluate the interaction of these ions (Al3+ and Mn2+) and the metabolites at the active sites of the two target enzymes. The results indicate that the ions potentially interact with the active sites of these enzymes and conceivably influence the AsA–GSH cycle. The metabolites showed strong interactions at the active sites of the enzymes. When the electrostatic surfaces of the metabolites and the ions were generated, it revealed that the surfaces overlap in the case of DHAR of rice and wheat, and MDHAR of rice. Thus, it was hypothesized that the metabolites may prevent the interaction of ions with the enzymes. This is an interesting approach to decipher the mechanism of action of secondary metabolites against the metal or metalloid - induced stress responses in cereals by aiming at specific targets. The findings of the present study are reasonably significant and may be the beginning of an interesting and useful approach towards comprehending the role of secondary metabolites for stress amelioration and mitigation in cereals grown under acidic soil conditions.

1. Introduction

Soil acidity is one of the major constraints that considerably affects global agricultural productivity. More than 40% of global arable lands are considered acidic, and they usually contain exceedingly high concentrations of aluminium (Al3+) and manganese (Mn2+), along with several other associated factors [1,2,3,4,5,6]. The presence of both Al3+ and Mn2+ at remarkably high concentrations remains as a major cause for the decline in crop yield and productivity in acidic soils [7]. While in neutral and basic soil conditions, aluminium exists as non-toxic forms such as aluminium oxide or aluminosilicates; however, it mobilizes to toxic trivalent cations (Al3+) under acidic soil conditions with pH < 5 [6,8]. The primary site of Al3+ injury is the plant root system, causing root growth inhibition [1,9,10,11,12,13]. Al3+ induces the production of reactive oxygen species (ROS) such as hydrogen peroxide (H2O2) and superoxide radical (O2∙−), and imparts oxidative stress, leading to loss of plasma membrane integrity and lipid peroxidation [14,15,16,17,18,19,20,21,22]. Besides causing damage to plasma membrane, Al3+–induced ROS production also results in growth inhibition, ATP depletion, inhibition of cellular respiration, alteration of redox homeostasis, mitochondrial dysfunction and metabolic alterations [22,23]. Like other stressors (abiotic) and metal/metalloid, Al3+-induced ROS production and subsequent alteration of physiological and cellular functions have been well documented in a large variety of plant species [24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39]. Manganese (Mn2+) toxicity is recognised as a major factor affecting plant growth and metabolism in acidic and poorly drained soils [40]. Besides its toxic impacts on plants at higher concentrations, Mn2+ is also a vital element required by plants for various physiological and cellular functions [41]. Mn2+ plays a crucial role during photosynthesis and also acts as co-factor for many important enzymes [42,43,44]. Some recent studies have demonstrated that the presence of Mn2+ at higher concentrations causes a strong detrimental impact on a plant’s physiological and biochemical functions, leading to cell death [44,45,46,47]. However, in comparison to plant responses to Al3+ in acidic soils, relatively very little is known regarding the basis of Mn2+ in plants.
Monodehydroascorbate reductase (MDHAR) and dehydroascorbate reductase (DHAR) are two major enzymes of the ascorbate–glutathione (AsA–GSH) cycle [48,49,50]. The AsA–GSH cycle is responsible for the detoxification of H2O2, produced through light-dependant photosynthetic reactions. During the ascorbate peroxidase enzyme mediated reduction of H2O2 to H2O, monodehydroascorbate (MDHA) radical is formed, which is then reduced to ascorbate by MDHAR or is unequally converted to ascorbate and dehydroascorbate (DHA) [48,49]. The DHA, thus formed, is converted to ascorbate by DHAR with the expense of reduced glutathione (GSH) to form oxidized glutathione (GSSG). The roles of plant hormones and secondary metabolites in regulating abiotic stress by modulating the antioxidant pool have been widely studied [51,52,53]. These include a diverse class of chemical compounds such as gibberellins (GA3), salicylic acid (SA), jasmonic acid (JA), serotonin (5-HT), ascorbic acid (AsA), glutathione (GSH), etc. The SA (or ortho hydroxyl benzoic acid) is a phenolic compound that regulates diverse aspects of plant responses to abiotic stresses [39,54,55,56]. SA can modulate stress responses and enhances stress tolerance by regulating the antioxidative defence system to scavenge the ROS production [57]. The role of SA in alleviating heavy metal and metalloid-induced stress responses have been reported in rice, maize, legumes and Cassia tora [54,58,59,60]. JA plays important role in regulating various essential processes in plants and is also believed to have a major ameliorating role in controlling the deleterious effects of various abiotic stresses [61,62]. Auxins or indole-3-acetic acid (IAA) and GA3 are important plant growth regulators which also play a major role in plants during abiotic stresses [63]. AsA and GSH are well known to ameliorate abiotic stresses by scavenging the toxic ROS produced in cells [64,65]. They are known to function stringently in controlling the levels of ROS produced during a wide range of stresses such as drought, salinity, heavy metals, metalloids, etc. [66,67,68]. Another metabolite that emerged as an important component in regulating abiotic stress in plants is serotonin (5-hydroxytryptamine, 5-HT) [52,69,70]. However, its precise role in regulating abiotic stress including Al3+ and Mn2+ is still not clearly known.
Growth regulators and metabolites are major components of plants which influence physiological and cellular functions. Though the signalling roles of these chemical entities during abiotic stress have been largely investigated, their specific functional role during abiotic stresses is largely unknown. Computational approaches to study the interaction of ions with antioxidant enzymes and the possible role of plant secondary metabolites thereon are seldom performed [28,29,59]. In this investigation, we studied the role of some plant growth regulators and secondary metabolites in ameliorating the possible Al3+ and Mn2+-induced stress in cereal crops by targeting MDHAR and DHAR as model antioxidant enzymes using an in silico approach.

2. Materials and Methods

2.1. The Target Enzymes

The three-dimensional structure of the DHAR of rice (OsDHAR) in complex with AsA was available at the RCSB Protein Data Bank (PDB) with id 5D9W. The structure was determined by X-ray diffraction at a resolution of 1.69 Å. Likewise, the structure of MDHAR from rice (OsMDHAR) with PDB id 5JCI, determined at 1.70 Å resolution employing X-ray diffraction, was available in PDB. These structures were downloaded from the database in .pdb format. The three-dimensional structures of DHAR and MDHAR of maize and wheat were obtained from the SWISS-MODEL Repository database (https://swissmodel.expasy.org/repository/. accessed on 20 June 2020 ), which is a repository of three-dimensional structures which were modelled using the SWISS-MODEL homology modelling tool [71]. The models of DHAR from wheat (TaDHAR) and maize (ZmDHAR) with accession numbers Q84UH6 and B4FT31, respectively were obtained from this database. Models of MDHAR from maize (ZmMDHAR) and wheat (TaMDHAR), bearing accession numbers C4J4E4 and K4HRS2, respectively, were obtained accordingly. All the structures were downloaded in .pdb format.

2.2. Ligands (Metabolites and Ions)

The two-dimensional structures of Al3+ and Mn2+ were obtained from the NCBI PubChem compounds database (https://pubchem.ncbi.nlm.nih.gov/, accessed on 20 June 2020), while three-dimensional structures of AsA, JA, SA, GSH, IAA, GA3, 5-HT, MDHA and DHA were downloaded as .sdf files from this database. The three-dimensional structure of MDHA was not available at the database, and hence the two-dimensional structure was downloaded and used as input file at the https://cactus.nci.nih.gov/translate/ (accessed on 20 June 2020). The output three-dimensional structure was downloaded as an .sdf file. The physical and chemical properties of the metabolites, such as molecular weight, molecular formula, topological polar surface area (TPSA), octanol/water partition coefficient (XlogP3) and numbers of hydrogen bond donor (HBD) and acceptor (HBA) groups were obtained from the PubChem compounds database. Details of the ligands are given in Table 1. The structures of DHA and MDHA were used as reference ligands, while others were test ligands.

2.3. Molecular Docking

Molecular docking is a powerful in silico tool which determines the interactions between a protein and a compound (called ligand). The strength of their interactions, site of interaction and optimal binding orientations are also determined during docking. In the present study, molecular docking was performed using Molegro Virtual Docker 6.0 (MVD), following [72,73]. MVD is a potent docking tool which provides over 87% accuracy [74]. For OsDHAR, the docking was performed at the AsA binding site, having coordinates X: 2.47; Y: −5.84; Z: 6.72, and amino acids within a radius of 20 Å were included in the docking site. The three-dimensional structure of ZmDHAR was obtained from SWISS-MODEL database. Available cavities of the structure were determined using MVD, and docking was accordingly performed at the cavities cantered around X: 0.50; Y: 31.04 and Z: −0.75, incorporating amino acid residues within a radius of 30 Å. TaDHAR was docked at X: 11.06; Y: −0.51; Z: 9.75, incorporating amino acids within a radius of 20 Å. OsMDHAR was docked at the FAD-binding active site, with coordinates X; 28.94; Y: 7.29; Z: 17.67, including amino acids within a radius of 20 Å. ZmMDHAR was docked at X: −11.73; Y: 24.87; Z: −15.20, including amino acids within a radius of 30 Å. TaMDHAR was docked at X: 27.13; Y: 17.74; Z: 15.78, including amino acids within a radius 25 Å.
The docking parameters used in the study include the MolDock scoring function, with 10 runs each of 1500 iterations for each ligand, and five best poses were retained, in terms of MolDoc score (i.e., highest free energy of binding of the ligand and the receptor), and they were selected for further analysis. For Al3+ and Mn2+, the five best poses were used for the analysis. Furthermore, the Rerank score and hydrogen bond score for each ligand with each target were recorded.

2.4. Visualization of the Interactions

The docking poses obtained following the docking using MVD were imported into the workspace to analyse the binding orientations and interactions. The electrostatic surfaces of the docked ligands or poses at the active site were generated for each ligand. The electrostatic surfaces of the best pose of all the docked secondary metabolites and the ions at the active site of the enzymes were individually developed. Following this, the surfaces were observed for possible overlap for each ion and metabolite to see if the surfaces of ions overlap with those of the secondary metabolites. Further, the amino acid residues of the active site of the receptors involved in different types of interactions (stearic, electrostatic and hydrogen bonding) with different ligands were visualised. In addition, energy maps were generated for each receptor/target to see the sites which favour these interactions and the ligands which interact at these sites.

2.5. Statistical Analysis

To identify the properties of the ligands which affect the docking scores, a correlation analysis was performed by determining Pearson’s correlation coefficient. The physical and chemical properties of the ligands considered for the analysis include molecular weight, TPSA, XlogP3, HBD and HBA groups present in the ligands, following Mazumder et al. [72,73]. These properties were analysed against MolDoc, Re-rank and Hydrogen bond scores. However, because docking scores are negative values, and smaller values represent better inhibition, all the docking scores were converted to positive values before the statistical analysis, following Choudhury et al. [59].

3. Results

3.1. Docking Scores against DHAR

All the ligands showed negative docking scores with the DHAR targets in rice, wheat and maize (Table 2). Against OsDHAR, JA shows the best MolDock score, while GSH showed the highest scores against ZmDHAR and TaDHAR. Compared to Al3+, these docking scores are 3.61-fold, 2.51-fold and 4.38-fold higher for these metabolites with the targets. Similarly, the scores of these metabolites with the three targets are higher by 3.46-fold, 2.51-fold and 4.12-fold, respectively, compared to Mn2+. However, AsA showed the highest hydrogen bonding scores against all the targets.

3.2. Docking Scores against MDHAR

Against MDHAR of rice, wheat and maize, GSH showed the best docking (MolDock) scores, followed by GA3 (Table 3). The docking scores of GSH were found to be higher than the scores of Al3+ and Mn2+ by 3.07-fold, 3.17-fold and 4.36-fold respectively, against the MDHAR of rice, maize and wheat. Against all the three targets, GSH also shows the highest hydrogen bond scores, indicating that hydrogen bonding is critical in the interactions and docking scores.

3.3. Docking Orientation and Overlapping Surfaces

Against OsDHAR, the ions bind at three sites of the enzyme, while all the metabolites bind to a single active site. On the other hand, in the case of ZmDHAR, all the ions bind to the same site while the metabolites bind at two pockets, and there is no common binding pocket. Interestingly, against TaDHAR, the metabolites bind to two pockets, Al3+ binds to a single site and Mn2+ binds to two sites. One of the sites of the binding of Mn2+ falls within one of the binding sites of the metabolites, while the other site and the site of binding of Al3+ are different from those of metabolites (Figure 1).
On developing electrostatic surfaces of the individual metabolites and the five best poses of the ions on OsDHAR, it was found that only one of the three poses of the ions partially overlap with each of the surfaces of the metabolites (Figure 2). In the case of ZmDHAR, there is no such overlapping of the surfaces (Figure S1). In the case of TaDHAR, there is partial overlap of the electrostatic surfaces of MDHA, DHA, AsA, SA, 5-HT and IAA with that of Al3+, and as also revealed from the docking poses, there is overlap of the surfaces of one of the poses of Mn2+ with MDHA, DHA, AA, SA, 5-HT and IAA. Here, the surfaces of GSH and GA3 do not overlap with those of the ions (Figure 3).
In OsMDHAR, the ions (Al3+ and Mn2+) bind to the same pocket of the enzyme as that of the metabolites with three different poses of each ion. However, out of the three poses of the two ions, the surfaces of only one of the poses overlaps with the surfaces of the metabolites (Figure 4). In the case of ZmMDHAR, the metabolites bind to three different pockets and the ions bind to an isolated single pocket. Thus, there is no overlap in the electrostatic surfaces (Figure S1). In the case of TaMDHAR, all the metabolites bind to a single elongated pocket of the enzyme, and the ions bind to another isolated single pocket, and thus no overlap in the electrostatic surfaces is seen (Figure S1).

3.4. Energy Map

With OsDHAR and ZmDHAR, the metabolites interact mainly in the electrostatic, stearic and HBD favourable regions, and not in the HBA region. However, the ions fall in the electrostatic and stearic favourable regions. In the case of TaDHAR, the metabolites dock similarly at electrostatic, stearic, HBA and HBD interaction favourable regions, while the ions interact at electrostatic interaction favourable regions. The metabolites interact with all the favourable regions in the case of OsMDHAR, while the ions dock at the electrostatic and stearic interaction favourable regions. In the case of ZmMDHAR and TaMDHAR, the metabolites interact more at the stearic favourable regions compared to other regions (Figure 5).

3.5. Ligand Properties Affecting the Interactions

Against OsDHAR, MolDock and Rerank scores are positively correlated to molecular weight, while hydrogen bond scores were found to be positively correlated to TPSA, HBD and HBA, and negatively correlated to XlogP3, albeit weakly. In the case of ZmDHAR, while the MolDock score has no correlation with any of the ligand properties studies, the hydrogen bond score was found to be negatively correlated to XlogP3, and the Rerank score was positively correlated to molecular weight, HBD and TPSA. Again, with TaDHAR, both MolDock and Rerank scores were not correlated to any of the properties of the ligands, while the hydrogen bond score was found to be positively correlated to HBA and TPSA, and negatively correlated to XlogP3 (Table 4).
Thus, the results indicate that the hydrogen bond scores were negatively correlated to XlogP3 for all the targets, while TPSA and HBA are positively correlated to the hydrogen bond scores for OsDHAR and TaDHAR. In case of all the MDHAR targets, MolDock scores were found to be mildly positively correlated to molecular weights, and the Rerank score was positively correlated only to the molecular weight in the case of TaMDHAR. Interestingly, with the hydrogen bonding scores of all the targets, while XlogP3 was negatively correlated, HBD, HBA and TPSA were found to be positively correlated. The correlations, both positive and negative, in the case of the hydrogen bond scores were stronger (Table 5).

3.6. Interactions of Different Metabolites

Interactions between the ligands and the active site of OsDHAR. (A) AsA; (B) JA; (C) SA; (D) 5-HT; (E) GA3; (F) GSH; (G) IAA; (H) DHA and (I) MDHA are shown in Figure 6. The dotted lines represent hydrogen bonding and stearic interactions between the ligands and the amino acid residues shown.
Figure 7 represents the interactions of different metabolites such as (A) MDHA; (B) DHA; (C) AsA; (D) GSH; (E) SA; (F) JA; (G) 5-HT; (H) GA3 and (I) IAA with TaDHAR. The amino acids shown are the ones with which the metabolites form hydrogen bonding and stearic interactions.
Figure 8 highlights the interactions of different metabolites (A) Co-crystallized FAD; (B) MDHA; (C) GSH; (D) AsA; (E) JA; (F) 5-HT; (G) SA; (H) GA3 and (I) IAA with the active site amino acid resides of OsMDHAR. The amino acids shown are the ones with which the metabolites form hydrogen bonds and stearic interactions. The interactions of ZmDHAR, ZmMDHAR and TaMDHAR are provided in Figures S2–S4, respectively.

4. Discussion

In the present study, molecular docking was performed between different secondary metabolites and two metal ions (Al3+, Mn2+) with the DHAR and MDHAR of rice, maize and wheat. Here, the enzymes were taken as models, and the potential of the metabolites in conferring protection against the ions were analysed. When a metabolite binds to an enzyme, it may prevent interaction and/or binding of another. This depends on the affinities, which are measured as docking scores. Docking scores represent the amount of energy liberated when a compound docks/binds with an enzyme. A ligand having a higher docking (a more negative) score may potentially prevent another which has a lower (a less negative) docking score. The validity and accuracy of the present modelling study is revealed from the fact that co-crystallized ligands were found to bind to the same active site with the same orientation when docked using MVD. Additionally, the metal/metalloid ions showed no hydrogen bond scores (Table 2), which further validates the modelling study. Furthermore, compared to other similar modelling algorithms, MVD has a better accuracy [74]. It is known that several plant secondary metabolites, including AsA, GSH, GA3, JA, SA, 5-HT and IAA, confer protection against metal/metalloid-induced abiotic stresses [75,76,77,78,79,80,81,82,83,84,85,86,87,88]. The present study is an attempt to understand whether these metabolites may prevent the interactions of Al3+ and Mn2+ with the two critical anti-oxidant enzymes (DHAR and MDHAR) of rice, maize and wheat.
In the present study, it was found that the docking scores of the metal ions were lower compared to the metabolites in the case of all the targets (Table 2). This is largely because the ions provide lesser surfaces for interactions and binding, compared to the metabolites. However, among the five best docking poses of the ligands with the OsDHAR, both the ions bind at three different locations, and out of them only one of the poses overlaps partially with the electrostatic surfaces of the metabolites (Figure 2). In the case of TaDHAR, while Al3+ binds at one, Mn2+ binds at two locations. The electrostatic surfaces of Al3+ and one of the poses of Mn2+ overlap with the surfaces of the metabolites (Figure 3). This indicates that the metabolites may potentially interfere with the interaction of the ions. However, no such overlapping occurs in the case of ZmDHAR (Figure S1), which indicates that the metabolites may not confer protection against these ions. In the case of OsMDHAR, both the metal ions bind at the active site of the enzyme at three locations, and out of them the surface of only one ion overlaps with all the metabolites, indicating a possible mode of protection (Figure 4). However, the other two MDHARs show no such interactions and overlaps (Figures S1, S3 and S4). Since the docking scores of the metal ions are lower than the secondary metabolites, and because they have overlapping electrostatic surfaces in the case of OsDHAR, TaDHAR and OsMDHAR, it is argued that the metabolites may prevent the ions from interacting with the enzymes. Thus, the results demonstrate that the toxicity of the ions may be reduced or prevented in the presence of higher cellular quantities of the metabolites. Similar findings have recently been reported by Choudhury et al. [54] in rice using DHAR as a model. Thus, the present study hypothesizes a novel mechanism of protection by plant secondary metabolites against metal-induced abiotic stress. The study is novel as it hypothesizes a novel mechanism of amelioration of metal-induced abiotic stress by plant secondary metabolites in cereals.
However, the major limitation of the present study is that the structures of DHAR and MDHAR of maize and wheat were not available, and their three-dimensional models were obtained from the SWISS-MODEL database. This may be one of the factors which influenced the overlapping of the electrostatic surfaces. It is assumed that, with actual structures, a better modelling result may be obtained. Nevertheless, the findings are crucial in understanding the mechanism of action of plant metabolites in protecting metal-induced abiotic stress, oxidative stress in particular, in cereals.

5. Conclusions

Abiotic stresses, such as Al3+ and Mn2+ stress, affect the antioxidant defence mechanisms in plants and alters the cellular redox homeostasis. This not only affects growth and productivity, but maintaining the redox homeostasis remains a major challenge. The present study demonstrated that the stressors such as Al3+ and Mn2+ associated with acidic soils may potentially interact with the two critical antioxidant enzymes of the AsA-GSH cycle, DHAR and MDHAR, and possibly affect their activities. However, secondary metabolites, including AsA, GSH, JA, SA, IAA, GA3 and 5-HT may prevent interaction of the metal ions with OsDHAR, TaDHAR and OsMDHAR. This is hypothesized to be affected by stronger interactions of the metabolites with the active sites of the enzymes. Because the findings hypothesize a novel mechanism in understanding the role of plant metabolites against metal-induced abiotic stress in cereals, it is highly significant. Thus, increasing the cellular levels of these metabolites through metabolic and genetic engineering, and plant breeding approaches, is hypothesized to ameliorate abiotic stress caused by the two metal ions.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/antiox11030458/s1, Figure S1: Electrostatic surfaces of the metabolites and ions at ZmDHAR (A,B), ZmMDHAR (C,D) and TaMDHAR (E,F). In (A–D), surfaces of all the metabolites are shown in white, and Al3+ in yellow and Mn2+ in red. The figure shows that there is no overlap between the metabolites and ions; Figure S2: Interactions of different metabolites with the active site amino acid residues of ZmDHAR with (A) DHA; (B) GSH; (C) AsA; (D) JA; (E) 5-HT; (F) SA; (G) GA3; (H) IAA. The dotted lines represent hydrogen bonding and stearic interactions between the ligands and the amino acid residues shown; Figure S3: Interactions of different metabolites (A) MDHA; (B) DHA; (C) AsA; (D) GSH; (E) SA; (F) JA; (G) 5-HT; (H) GA3; (I) IAA with different amino acids of the active site of ZmMDHAR. The amino acids shown are the ones with which the metabolites form hydrogen bonds and stearic interactions; Figure S4: Interactions of the metabolites (A) MDHA; (B) DHA; (C) AsA; (D) GSH; (E) SA; (F) JA; (G) 5-HT; (H) GA3; (I) IAA with the active site residues of TaMDHAR. The amino acids shown are the ones with which the metabolites form hydrogen bonds and stearic interactions.

Author Contributions

Conceptualization: S.C., M.K.M., D.M., P.S., S.K.T. and D.G.; methodology: S.C., M.K.M., D.M., P.S., S.K.T. and D.G.; software: S.C., M.H.S., A.H. and D.M., validation: S.C., M.K.M., D.M., P.S., S.K.T. and D.G.; formal analysis: S.C. and D.G.; investigation: S.C., M.K.M., D.M., P.S., S.K.T. and D.G.; resources: M.K.M., S.C., H.M.A. and M.H.S.; data curation: S.C., D.M., A.H., H.M.A. and D.G.; writing—original draft: S.C., M.K.M., D.M., P.S., S.K.T. and D.G.; writing—review and editing: M.B., M.S., H.M.A., M.H.S. and A.H.; visualization: S.C., M.K.M., D.M., P.S., S.K.T. and D.G.; supervision: D.G., M.K.M., S.C. and S.K.T.; project administration: M.K.M., H.M.A. and M.H.S.; funding acquisition: H.M.A., M.H.S. and A.H. All authors have read and agreed to the published version of the manuscript.

Funding

The study was funded by Plant Stress Biology and Metabolomics Laboratory Central Instrumentation Laboratory (CIL), Assam University, Silchar-788 011, India and Department of Biotechnology, Government of India (DBT - Research Associateship). The study was also supported by the Researchers Supporting Project number (RSP-2021/347), King Saud University, Riyadh, Saudi Arabia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the relevant data is available in the article and its Supplementary Materials.

Acknowledgments

The authors would like to extend their sincere appreciation to the Researchers Supporting Project number (RSP-2021/347), King Saud University, Riyadh, Saudi Arabia.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Foy, C.D. Physiological effects of hydrogen, Al and manganese toxicities in acid soil. In Soil Acidity and Liming; Pearson, R.W., Adams, F., Eds.; American Society of Agronomy: Madison, WI, USA, 1984; pp. 57–97. [Google Scholar]
  2. Vaughan, D.; Ord, B.G. Influence of phenolics acids on the sodium, calcium and chloride contents of Pisum sativum under axenic conditions. Soil Biol. Biochem. 1991, 23, 1191–1193. [Google Scholar] [CrossRef]
  3. Hafner, H.; Ndunguru, B.J.; Bationo, A.; Marschner, H. Effect of nitrogen, phosphorus and molybdenum application on growth and symbiotic N2-fixation of groundnut in an acid sandy soil in Niger. Fertz. Res. 1992, 31, 69–77. [Google Scholar] [CrossRef]
  4. Baziramakenga, R.; Leroux, G.D.; Simard, R.R. Effects of benzoic and cinnamic acids on membrane permeability of soybean roots. J. Chem. Ecol. 1995, 21, 1271–1285. [Google Scholar] [CrossRef]
  5. Barceló, J.; Poschenrieder, C.; Vázquez, M.D.; Gunsé, B. Aluminum phytotoxicity: A challenge for plant scientists. Fertz. Res. 1996, 43, 217–223. [Google Scholar] [CrossRef]
  6. Kidd, P.S.; Proctor, J. Why plant grows poorly on very acid soil: Are ecologists missing the obvious? J. Exp. Bot. 2001, 52, 791–799. [Google Scholar] [CrossRef] [Green Version]
  7. Xia, J.; Yamaji, N.; Ma, J.F. Plasma membrane-localized transporter for aluminium in rice. Proc. Natl. Acad. Sci. USA 2010, 107, 18381–18385. [Google Scholar] [CrossRef] [Green Version]
  8. Kochian, L.V.; Piñeros, M.A.; Hoekenga, O.A. The physiology, genetics and molecular biology of plant aluminium resistance and toxicity. Plant Soil 2005, 274, 175–195. [Google Scholar] [CrossRef]
  9. Delhaize, E.; Craig, S.; Beatson, C.D.; Bennet, R.J.; Jagadish, V.C.; Randall, P.J. Aluminium tolerance in wheat (Triticum aestivum L.) (I. Uptake and distribution of aluminium in root apices). Plant Physiol. 1993, 103, 685–693. [Google Scholar] [CrossRef] [Green Version]
  10. Ryan, P.R.; Kochian, L.V. Interaction between aluminium toxicity and calcium uptake at root apex in near—Isogenic lines of wheat (Triticum aestivum L.) differing in aluminium tolerance. Plant Physiol. 1993, 102, 975–982. [Google Scholar] [CrossRef]
  11. Sivaguru, M.; Horst, W.J. The distal part of the transition zone is the most aluminium sensitive apical root zone of maize. Plant Physiol. 1998, 116, 155–163. [Google Scholar] [CrossRef] [Green Version]
  12. Kochian, L.V. Cellular mechanism of aluminium toxicity and resistance in plants. Ann. Rev. Plant Physiol. Plant Mol. Biol. 1995, 46, 237–260. [Google Scholar] [CrossRef]
  13. Ma, J.F. Syndrome of aluminium toxicity and diversity of aluminum resistance in higher plants. Int. Rev. Cytol. 2007, 264, 225–252. [Google Scholar] [PubMed]
  14. Ono, K.; Yamamoto, Y.; Hachiya, A.; Matsumoto, H. Synergistic inhibition of growth by aluminium and iron of tobacco (Nicotiana tabacum L.) cells in suspension culture. Plant Cell Physiol. 1995, 36, 115–125. [Google Scholar]
  15. Yamamoto, Y.; Hachiya, A.; Matsumoto, H. Oxidative damage to membranes by a combination of aluminium and iron in suspension-cultured tobacco cells. Plant Cell Physiol. 1997, 38, 1333–1339. [Google Scholar] [CrossRef] [Green Version]
  16. Yamamoto, Y.; Kobayashi, Y.; Matsumoto, H. Lipid peroxidation is an early symptom triggered by aluminium, but not the primary cause of elongation inhibition in pea roots. Plant Physiol. 2001, 12, 199–208. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Yamaguchi, Y.; Yamamoto, Y.; Matsumoto, H. Cell death process initiated by combination of aluminium and iron in suspension cultured tobacco cells. Soil Sci. Plant Nutr. 1999, 45, 647–657. [Google Scholar] [CrossRef]
  18. Kawano, T.; Kadono, T.; Furuichi, T.; Muto, S.; Lapeyrie, F. Aluminium induced distortion in calcium signaling involving oxidative burst and channel regulation in tobacco BT-2 cells. Biochem. Biophys. Res. Commun. 2003, 308, 35–42. [Google Scholar] [CrossRef]
  19. Exley, C. The pro-oxidant activity of aluminium. Free Radic. Biol. Med. 2004, 236, 380–387. [Google Scholar] [CrossRef]
  20. Pereira, J.F.; Zhou, G.; Delhaize, E.; Richardson, T.; Zhou, M.; Ryan, P.R. Engineering greater aluminium resistance in wheat by over-expressing TaALMT1. Ann. Bot. 2010, 106, 205–214. [Google Scholar] [CrossRef] [Green Version]
  21. Matsumoto, H.; Motoda, H. Aluminium toxicity recovery process in root apices. Possible association with oxidative stress. Plant Sci. 2012, 185–186, 1–8. [Google Scholar] [CrossRef]
  22. Choudhury, S.; Sharma, P. Aluminum stress inhibits root growth and alters physiological and metabolic responses in chickpea (Cicer arietinum L.). Plant Physiol. Biochem. 2014, 85, 63–70. [Google Scholar] [CrossRef]
  23. Yamamoto, Y.; Kobayashi, Y.; Devi, S.R.; Rikiishi, S.; Matsumoto, H. Aluminium toxicity is associated with mitochondrial dysfunction and the production of reactive oxygen species in plant cells. Plant Physiol. 2002, 128, 63–72. [Google Scholar] [CrossRef] [PubMed]
  24. Moulick, D.; Samanta, S.; Sarkar, S.; Mukherjee, A.; Pattnaik, B.K.; Saha, S.; Awasthi, J.P.; Bhowmick, S.; Ghosh, D.; Samal, A.C.; et al. Arsenic contamination, impact and mitigation strategies in rice agro-environment: An inclusive insight. Sci. Total Environ. 2021, 800, 149477. [Google Scholar] [CrossRef] [PubMed]
  25. Moulick, D.; Chowardhara, B.; Panda, S.K. Agroecotoxicological aspect of arsenic (As) and cadmium (Cd) on Field crops and its mitigation: Current status and future prospect. In Plant-Metal Interactions; Springer: Cham, Switzerland, 2019; pp. 217–246. [Google Scholar]
  26. Hossain, A.; Pramanick, B.; Bhutia, K.L.; Ahmad, Z.; Moulick, D.; Maitra, S.; Ahmad, A.; Aftab, T. Emerging roles of osmoprotectant glycine betaine against salt-induced oxidative stress in plants: A major outlook of maize (Zea mays L.). In Frontiers in Plant-Soil Interaction; Academic Press: Cambridge, MA, USA, 2021; pp. 567–587. [Google Scholar]
  27. Hossain, A.; Ahmad, Z.; Moulick, D.; Maitra, S.; Bhadra, P.; Ahmad, A.; Garai, S.; Mondal, M.; Roy, A.; Sabagh, A.E.; et al. Jasmonates and Salicylates: Mechanisms, Transport and Signalling During Abiotic Stress in Plants. In Jasmonates and Salicylates Signaling in Plants; Springer: Cham, Switzerland, 2021; pp. 1–29. [Google Scholar]
  28. Mazumder, M.K.; Moulick, D.; Choudhury, S. Iron (Fe 3+)-mediated redox responses and amelioration of oxidative stress in cadmium (Cd 2+) stressed mung bean seedlings: A biochemical and computational analysis. J. Plant Biochem. Biotechnol. 2021, 31, 49–60. [Google Scholar] [CrossRef]
  29. Mazumder, M.K.; Sharma, P.; Moulick, D.; Tata, S.K.; Choudhury, S. Salicylic acid ameliorates zinc and chromium-induced stress responses in wheat seedlings: A biochemical and computational analysis. Cereal Res. Commun. 2021, 1–12. [Google Scholar] [CrossRef]
  30. Moulick, D.; Ghosh, D.; Santra, S.C. Evaluation of effectiveness of seed priming with selenium in rice during germination under arsenic stress. Plant Physiol. Biochem. 2016, 109, 571–578. [Google Scholar] [CrossRef]
  31. Moulick, D.; Santra, S.C.; Ghosh, D. Seed priming with Se alleviate As induced phytotoxicity during germination and seedling growth by restricting As translocation in rice (Oryza sativa L cv IET-4094). Ecotoxicol. Environ. Saf. 2017, 145, 449–456. [Google Scholar] [CrossRef]
  32. Moulick, D.; Santra, S.C.; Ghosh, D. Seed priming with Se mitigates As-induced phytotoxicity in rice seedlings by enhancing essential micronutrient uptake and translocation and reducing As translocation. Environ. Sci. Pollut. Res. 2018, 25, 26978–26991. [Google Scholar] [CrossRef]
  33. Chowardhara, B.; Borgohain, P.; Saha, B.; Awasthi, J.P.; Moulick, D.; Panda, S.K. Phytotoxicity of Cd and Zn on three popular Indian mustard varieties during germination and early seedling growth. Biocatal. Agric. Biotechnol. 2019, 21, 101349. [Google Scholar] [CrossRef]
  34. Saha, B.; Chowardhara, B.; Kar, S.; Devi, S.S.; Awasthi, J.P.; Moulick, D.; Tanti, B.; Panda, S.K. Advances in heavy metal-induced stress alleviation with respect to exogenous amendments in crop plants. In Priming and Pretreatment of Seeds and Seedlings; Springer: Singapore, 2019; pp. 313–332. [Google Scholar]
  35. Moulick, D.; Samanta, S.; Saha, B.; Mazumder, M.K.; Dogra, S.; Panigrahi, K.C.; Banerjee, S.; Ghosh, D.; Santra, S.C. Salinity Stress Responses in Three Popular Field Crops Belonging to Fabaceae Family: Current Status and Future Prospect. In The Plant Family Fabaceae; Springer: Singapore, 2020; pp. 519–541. [Google Scholar]
  36. Meriga, B.; Reddy, B.K.; Rao, K.R.; Raddy, L.A.; Kavi Kishor, P.B. Aluminium-induced production of oxygen radicals, lipid peroxidation and DNA damage in seedlings of rice (Oryza sativa). J. Plant Physiol. 2004, 161, 63–68. [Google Scholar] [CrossRef]
  37. Dipierro, N.; Mondelli, D.; Paciolla, C.; Brunetti, G.; Dipierro, S. Changes in ascorbate system in the response of pumpkin (Cucarbita pepo L.) roots to aluminium stress. J. Plant Physiol. 2005, 162, 529–536. [Google Scholar] [CrossRef] [PubMed]
  38. Schuch, M.W.; Cellini, A.; Masia, A.; Marino, G. Aluminium induced effects on growth, morphogenesis and oxidative stress reactions in in vitro cultures of quince. Sci. Hortic. 2010, 125, 151–158. [Google Scholar] [CrossRef]
  39. Wang, L.J.; Li, S.H. Salicylic acid induced heat or cold tolerance in relation to Ca2+ homeostasis and antioxidant system in young grape plants. Plant Sci. 2006, 170, 685–694. [Google Scholar] [CrossRef]
  40. Horst, W.J. The physiology manganese toxicity. In Manganese in Soil and Plants; Graham, R.D., Hannam, R.J., Uren, N.C., Eds.; Kluwer Academic Publishers: Amsterdam, The Netherlands, 1998; pp. 175–188. [Google Scholar]
  41. Doncheva, S.; Poschenrieder, C.; Stoyanova, Z.; Georgieva, K.; Velichkova, M.; Bareceló, J. Silicon amelioration of manganese toxicity in Mn sensitive and Mn tolerant maize varieties. Environ. Exp. Bot. 2009, 65, 189–197. [Google Scholar] [CrossRef]
  42. Sharma, A. Manganese in cell metabolism of higher plants. Bot. Rev. 1991, 57, 117–149. [Google Scholar]
  43. Gonzalez, A.; Steffen, K.L.; Lynch, J.P. Light and excess manganese. Plant Physiol. 1998, 118, 493–504. [Google Scholar] [CrossRef] [Green Version]
  44. Huang, H.; Zhao, Y.; Xu, Z.; Zhang, W.; Jiang, K. Physiological responses of Broussonetia papyrifera to manganese stress, a candidate plant for Phytoremediation. Ecotoxicol. Environ. Saf. 2019, 181, 18–25. [Google Scholar] [CrossRef]
  45. Braun, H.P. Effect of manganese toxicity on the proteome of the leaf apoplast in cowpea. Plant Physiol. 2003, 133, 1935–1946. [Google Scholar]
  46. Delhaize, E.; Gruber, B.D.; Pittman, J.K.; White, R.G.; Leung, H.; Miao, Y.; Jiang, L.; Ryan, P.R.; Richardson, A.E. A role for AtMTP11 gene of Arabidopsis in manganese transport and tolerance. Plant J. 2007, 51, 198–210. [Google Scholar] [CrossRef]
  47. Huang, H.; Zhao, Y.; Xu, Z.; Ding, Y.; Zhang, W.; Wu, L. Biosorption characteristics of a highly Mn(II)-resistant Ralstonia pickettii strain isolated from Mn ore. PLoS ONE 2018, 13, e0203285. [Google Scholar] [CrossRef]
  48. Noctor, G.; Foyer, C.H. Ascorbate and glutathione: Keeping active oxygen under control. Ann. Rev. Plant Physiol. Plant Mol. Biol. 1998, 49, 249–279. [Google Scholar] [CrossRef] [PubMed]
  49. Jimenez, A.; Hernandez, J.A.; del Rio, L.A.; Sevilla, F. Evidence for the presence of the ascorbate-glutathione cycle in mitochondria and peroxisomes of pea leaves. Plant Physiol. 1997, 114, 275–284. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Meyer, A.J.; Brach, T.; Marty, L.; Kreye, S.; Rouhier, N.; Jacquot, J.P.; Hell, R. Redox sensitive GFP in Arabidopsis thaliana is a quantitative biosensor for the redox potential of the cellular glutathione redox buffer. Plant J. 2007, 52, 973–986. [Google Scholar] [CrossRef] [PubMed]
  51. Bari, R.; Jones, J.D.G. Role of plant hormones in plant defense responses. Plant Mol Biol. 2009, 69, 473–488. [Google Scholar] [CrossRef] [PubMed]
  52. Zhang, J.; Li, D.; Wei, J.; Ma, W.; Kong, X.; Rengel, Z.; Xiao, L.; Xia, G. Melatonin alleviates aluminium induced root growth inhibition by interfering with nitric oxide production in Arabidopsis. Environ. Exp. Bot. 2018, 161, 157–165. [Google Scholar] [CrossRef]
  53. Zhu, C.Q.; Zhang, J.H.; Sun, L.M.; Zhu, L.F.; Abliz, B.; Hu, W.J.; Zhong, C.; Bai, Z.G.; Sajid, H.; Cao, X.C.; et al. Hydrogen sulfide alleviates aluminium toxicity via decreasing apoplast and symplast Al contents in rice. Front. Plant Sci. 2018, 9, 294. [Google Scholar] [CrossRef] [Green Version]
  54. Choudhury, S.; Panda, S.K. Role of salicylic acid in regulating cadmium (Cd) induced oxidative stress in Oryza sativa L. roots. Bulg. J. Plant Physiol. 2004, 30, 95–110. [Google Scholar]
  55. Cao, Y.; Zhang, Z.W.; Xue, L.W.; Du, J.B.; Shang, J.; Xu, F.; Yuan, S.; Lin, H.-H. Lack of salicylic acid in Arabidopsis protects plants against moderate salt stress. Z. Nat. C J. Biosci. 2009, 64, 231–238. [Google Scholar] [CrossRef]
  56. Nazar, R.; Umar, S.; Khan, N.A. Exogenous salicylic acid improves photosynthesis and growth through increase in ascorbate-glutathione metabolism and S assimilation in mustard under salt stress. Plant Signal. Behav. 2015, 10, e1003751. [Google Scholar] [CrossRef] [Green Version]
  57. Ma, X.; Zheng, J.; Zhang, X.; Hu, Q.; Qian, R. Salicylic acid alleviates the adverse effects of salt stress on Dianthus superbus (Carophyllaceae) by activating photosynthesis, protecting morphological structure and enhancing the antioxidant system. Front. Plant Sci. 2017, 8, 600. [Google Scholar] [CrossRef]
  58. Chen, J.; Zhu, C.; Li, L.P.; Sun, Z.Y. Effects of exogenous salicylic acid on growth and H2O2—Metabolizing enzymes in rice seedlings under lead stress. J. Environ. Sci. 2007, 19, 44–49. [Google Scholar] [CrossRef]
  59. Choudhury, S.; Moulick, D.; Mazumder, M.K. Secondary metabolites protect against metal and metalloid stress in rice: An in silico investigation using dehydroascorbate reductase. Acta Physiol. Plant 2020, 43, 3. [Google Scholar] [CrossRef]
  60. Zengin, F. Exogenous treatment with salicylic acid alleviating copper toxicity in bean seedlings. Proc. Natl. Acad. Sci. India Sect. B Biol. Sci. 2014, 84, 749–755. [Google Scholar] [CrossRef]
  61. Yoon, J.Y.; Hamayun, M.; Lee, S.K.; Lee, I.J. Methyl jasmonate alleviated salinity stress in soybean. J. Crop Sci. Biotechnol. 2009, 12, 63–68. [Google Scholar] [CrossRef]
  62. Clarke, S.M.; Cristescu, S.M.; Miersch, O.; Harren, F.J.; Wastenack, C.; Mur, L.A. Jasmonates act with salicylic acid to confer basal thermotolerance in Arabidopsis thaliana. New Phytol. 2009, 182, 175–187. [Google Scholar] [CrossRef] [PubMed]
  63. Sharma, A.; Zheng, B. Molecular responses during plant grafting and its regulation by auxins, cytokinins, and gibberellins. Biomolecules 2019, 9, 397. [Google Scholar] [CrossRef] [Green Version]
  64. Hare, P.D.; Cress, W.A.; Van, S.J. Dissecting the role of osmolyte accumulation during stress. Plant Cell Environ. 1998, 21, 535–553. [Google Scholar] [CrossRef]
  65. Asada, K. The water—Water cycle in the chloroplast: Scavenging active oxygen species and dissipation of excess photons. Ann. Rev. Plant Physiol. Mol. Biol. 1999, 50, 601–639. [Google Scholar] [CrossRef]
  66. Hernandez, J.A.; Jimenez, A.; Mullineaux, P.M.; Sevilla, F. Tolerance of pea (Pisum sativum L.) to long term salt stress is associated with induction of antioxidant defenses. Plant Cell Environ. 2000, 23, 853–862. [Google Scholar] [CrossRef]
  67. Shao, H.B.; Jiang, S.Y.; Li, F.M.; Chu, L.Y.; Zhao, C.X.; Shao, M.A.; Zhao, X.; Li, F. Some advances in plant stress physiology and their implications in the systems biology era. Biointerfaces 2007, 54, 33–36. [Google Scholar] [CrossRef]
  68. Sreenivasulu, N.; Grimma, B.; Wobusa, U.; Wescheka, W. Differential responses of antioxidant compounds to salinity stress in salt-tolerant and salt-sensitive seedlings of foxtail millet (Setaria itilica). Physiol. Plant. 2000, 109, 435–442. [Google Scholar] [CrossRef]
  69. Mukherjee, S. Novel perspective of molecular crosstalk mechanisms of serotonin and melatonin in plants. Plant Physiol. Biochem. 2018, 132, 33–45. [Google Scholar] [CrossRef]
  70. Ke, Q.; Ye, J.; Wang, B.; Ren, J.; Yin, L.; Deng, X.; Wang, S. Melatonin mitigates salt stress in wheat seedlings by modulating polyamine metabolism. Front. Plant Sci. 2018, 9, 914. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Bienert, S.; Waterhouse, A.; de Beer, T.A.; Tauriello, G.; Studer, G.; Bordoli, L.; Schwede, T. The SWISS-MODEL Repository—New features and functionality. Nucleic Acids Res. 2017, 45, D313–D319. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Mazumder, M.K.; Borah, A.; Choudhury, S. Inhibitory potential of plant secondary metabolites on anti-Parkinsonian drug targets: Relevance to pathophysiology, and motor and non-motor behavioural abnormalities. Med. Hypotheses 2020, 137, 109544. [Google Scholar] [CrossRef] [PubMed]
  73. Mazumder, M.K.; Choudhury, S.; Borah, A. An in silico investigation on the inhibitory potential of the constituents of Pomegranate juice on antioxidant defense mechanism: Relevance to neurodegenerative diseases. IBRO Rep. 2019, 6, 153–159. [Google Scholar] [CrossRef] [PubMed]
  74. Thomsen, R.; Christensen, M.H. MolDock: A new technique for high-accuracy molecular docking. J. Med. Chem. 2006, 49, 3315–3321. [Google Scholar] [CrossRef]
  75. Ahmad, P.; Abass Ahanger, M.; Nasser Alyemeni, M.; Wijaya, L.; Alam, P.; Ashraf, M. Mitigation of sodium chloride toxicity in Solanum lycopersicum L. by supplementation of jasmonic acid and nitric oxide. J. Plant Interact. 2018, 13, 64–72. [Google Scholar] [CrossRef] [Green Version]
  76. Ahmad, P.; Alyemeni, M.N.; Ahanger, M.A.; Egamberdieva, D.; Wijaya, L.; Alam, P. Salicylic acid (SA) induced alterations in growth, biochemical attributes and antioxidant enzyme activity in faba bean (Vicia faba L.) seedlings under NaCl toxicity. Russ. J. Plant Physiol. 2018, 65, 104–114. [Google Scholar] [CrossRef]
  77. Ahmad, P.; Alyemeni, M.N.; Vijaya, L.; Alam, P.; Ahanger, M.A.; Alamri, S.A. Jasmonic acid alleviates negative impacts of cadmium stress by modifying osmolytes and antioxidants in faba bean (Vicia faba L.). Arch. Agron. Soil Sci. 2017, 63, 1889–1899. [Google Scholar] [CrossRef]
  78. Ahmad, P.; Nabi, G.; Ashraf, M. Cadmium-induced oxidative damage in mustard [Brassica juncea (L.) Czern. & Coss.] plants can be alleviated by salicylic acid. S. Afr. J. Bot. 2011, 77, 36–44. [Google Scholar]
  79. Ahmad, P.; Sarwat, M.; Sharma, S. Reactive oxygen species, antioxidants and signaling in plants. J. Plant Biol. 2008, 51, 167–173. [Google Scholar] [CrossRef]
  80. Mukherjee, S.; David, A.; Yadav, S.; Baluška, F.; Bhatla, S.C. Salt stress-induced seedling growth inhibition coincides with differential distribution of serotonin and melatonin in sunflower seedling roots and cotyledons. Physiol. Plant 2014, 152, 714–728. [Google Scholar] [CrossRef] [PubMed]
  81. Viehweger, K. How plants cope with heavy metals. Bot. Stud. 2014, 55, 35. [Google Scholar] [CrossRef] [Green Version]
  82. Byeon, Y.; Choi, G.H.; Lee, H.Y.; Back, K. Melatonin biosynthesis requires N-acetylserotonin methyltransferase activity of caffeic acid O-methyltransferase in rice. J. Exp. Bot. 2015, 66, 6917–6925. [Google Scholar] [CrossRef] [Green Version]
  83. Bali, S.; Kaur, P.; Sharma, A.; Ohri, P.; Bhardwaj, R.; Alyemeni, M.N.; Wijaya, L.; Ahmad, P. Jasmonic acid-induced tolerance to root-knot nematodes in tomato plants through altered photosynthetic and antioxidative defense mechanisms. Protoplasma 2017, 255, 471–484. [Google Scholar] [CrossRef]
  84. Bali, S.; Kaur, P.; Kohli, S.K.; Ohri, P.; Thukral, A.K.; Bhardwaj, R.; Wijaya, L.; Alyemeni, M.N.; Ahmad, P. Jasmonic acid induced changes in physio-biochemical attributes and ascorbate-glutathione pathway in Lycopersicon esculentum under lead stress at different growth stages. Sci. Total Environ. 2018, 645, 1344–1360. [Google Scholar] [CrossRef]
  85. Kohli, S.K.; Handa, N.; Sharma, A.; Gautam, V.; Arora, S.; Bhardwaj, R.; Alyemeni, M.N.; Wijaya, L.; Ahmad, P. Combined effect of 24-epibrassinolide and salicylic acid mitigates lead (Pb) toxicity by modulating various metabolites in Brassica juncea L. seedlings. Protoplasma 2018, 255, 11–24. [Google Scholar] [CrossRef]
  86. Kohli, S.K.; Handa, N.; Sharma, A.; Gautam, V.; Arora, S.; Bhardwaj, R.; Alyemeni, M.N.; Ahmad, P. Interaction of 24-epibrassinolide and salicylic acid regulates pigment contents, antioxidative defense responses, and gene expression in Brassica juncea L. seedlings under Pb stress. Environ. Sci. Pollut. Res. 2018, 25, 15159–15173. [Google Scholar] [CrossRef]
  87. Mir, M.A.; John, R.; Alyemeni, M.N.; Alam, P.; Ahmad, P. Jasmonic acid ameliorates alkaline stress by improving growth performance, ascorbate glutathione cycle and glyoxylase system in maize seedlings. Sci. Rep. 2018, 8, 2831. [Google Scholar] [CrossRef]
  88. Mir, M.A.; Sirhindi, G.; Alyemeni, M.N.; Alam, P.; Ahmad, P. Jasmonic acid improves growth performance of soybean under nickel toxicity by regulating nickel uptake, redox balance, and oxidative stress metabolism. J. Plant Growth Regul. 2018, 37, 1195–1209. [Google Scholar] [CrossRef]
Antioxidants 11 00458 g001 550
Figure 1. Docking poses of the metabolites and ions at the active sites of the six targets. The protein chain is shown in the secondary structures.
Figure 1. Docking poses of the metabolites and ions at the active sites of the six targets. The protein chain is shown in the secondary structures.
Antioxidants 11 00458 g001
Antioxidants 11 00458 g002 550
Figure 2. Overlapping electrostatic surfaces of metabolites and metals at the active site of OsDHAR. (A) AsA; (B) JA; (C) SA; (D) 5-HT; (E) GA3; (F) GSH; (G) IAA and (H) DHA with Mn2+; and (I) AsA; (J) JA; (K) SA; (L) 5-HT; (M) GA3; (N) GSH; (O) IAA and (P) DHA with Al3+. The spherical purple and red spheres represent the surfaces of Al3+ and Mn2+ ions, respectively.
Figure 2. Overlapping electrostatic surfaces of metabolites and metals at the active site of OsDHAR. (A) AsA; (B) JA; (C) SA; (D) 5-HT; (E) GA3; (F) GSH; (G) IAA and (H) DHA with Mn2+; and (I) AsA; (J) JA; (K) SA; (L) 5-HT; (M) GA3; (N) GSH; (O) IAA and (P) DHA with Al3+. The spherical purple and red spheres represent the surfaces of Al3+ and Mn2+ ions, respectively.
Antioxidants 11 00458 g002
Antioxidants 11 00458 g003 550
Figure 3. Electrostatic surfaces of the ions and metabolites and their overlaps in the case of TaDHAR: (A) MDHA; (B) DHA; (C) AsA; (D) GSH; (E) SA; (F) 5-HT; (G) GA3; (H) IAA with Al3+. (I) MDHA; (J) DHA; (K) AsA; (L) GSH; (M) SA; (N) 5-HT; (O) GA3; (P) IAA with Mn2+. The dark blue and yellow spheres represent the electrostatic surfaces of Al3+ and Mn2+.
Figure 3. Electrostatic surfaces of the ions and metabolites and their overlaps in the case of TaDHAR: (A) MDHA; (B) DHA; (C) AsA; (D) GSH; (E) SA; (F) 5-HT; (G) GA3; (H) IAA with Al3+. (I) MDHA; (J) DHA; (K) AsA; (L) GSH; (M) SA; (N) 5-HT; (O) GA3; (P) IAA with Mn2+. The dark blue and yellow spheres represent the electrostatic surfaces of Al3+ and Mn2+.
Antioxidants 11 00458 g003
Antioxidants 11 00458 g004 550
Figure 4. Electrostatic surfaces of the ligands at the active site of OsMDHAR. Overlapping electrostatic surfaces of docked ligands and metal ions. (A) MDHA; (B) GSH; (C) AsA; (D) JA; (E) 5-HT; (F) SA; (G) GA3 and (H) IAA with Al3+; and (I) MDHA; (J) GSH; (K) AsA; (L) JA; (M) 5-HT; (N) SA; (O) GA3 and (P) IAA with Mn2+. The dark blue spheres represent the surfaces of Al3+ and the green spheres represent Mn2+.
Figure 4. Electrostatic surfaces of the ligands at the active site of OsMDHAR. Overlapping electrostatic surfaces of docked ligands and metal ions. (A) MDHA; (B) GSH; (C) AsA; (D) JA; (E) 5-HT; (F) SA; (G) GA3 and (H) IAA with Al3+; and (I) MDHA; (J) GSH; (K) AsA; (L) JA; (M) 5-HT; (N) SA; (O) GA3 and (P) IAA with Mn2+. The dark blue spheres represent the surfaces of Al3+ and the green spheres represent Mn2+.
Antioxidants 11 00458 g004
Antioxidants 11 00458 g005 550
Figure 5. Energy maps of the receptors showing favourable surfaces for different types of interactions. The docked poses of the ligands including metabolites and ions are also shown.
Figure 5. Energy maps of the receptors showing favourable surfaces for different types of interactions. The docked poses of the ligands including metabolites and ions are also shown.
Antioxidants 11 00458 g005
Antioxidants 11 00458 g006 550
Figure 6. Interactions between the ligands and the active site of OsDHAR. (A) AsA; (B) JA; (C) SA; (D) 5-HT; (E) GA3; (F) GSH; (G) IAA; (H) DHA and (I) MDHA. The dotted lines represent hydrogen bonding and stearic interactions between the ligands and the amino acid residues shown.
Figure 6. Interactions between the ligands and the active site of OsDHAR. (A) AsA; (B) JA; (C) SA; (D) 5-HT; (E) GA3; (F) GSH; (G) IAA; (H) DHA and (I) MDHA. The dotted lines represent hydrogen bonding and stearic interactions between the ligands and the amino acid residues shown.
Antioxidants 11 00458 g006
Antioxidants 11 00458 g007 550
Figure 7. Interactions of different metabolites (A) MDHA; (B) DHA; (C) AsA; (D) GSH; (E) SA; (F) JA; (G) 5-HT; (H) GA3 and (I) IAA with TaDHAR. The amino acids shown are the ones with which the metabolites form hydrogen bonding and stearic interactions.
Figure 7. Interactions of different metabolites (A) MDHA; (B) DHA; (C) AsA; (D) GSH; (E) SA; (F) JA; (G) 5-HT; (H) GA3 and (I) IAA with TaDHAR. The amino acids shown are the ones with which the metabolites form hydrogen bonding and stearic interactions.
Antioxidants 11 00458 g007
Antioxidants 11 00458 g008 550
Figure 8. Interactions of different metabolites (A) Co-crystallized FAD; (B) MDHA; (C) GSH; (D) AsA; (E) JA; (F) 5-HT; (G) SA; (H) GA3 and (I) IAA with the active site amino acid resides of OsMDHAR. The amino acids shown are the ones with which the metabolites form hydrogen bonds and stearic interactions.
Figure 8. Interactions of different metabolites (A) Co-crystallized FAD; (B) MDHA; (C) GSH; (D) AsA; (E) JA; (F) 5-HT; (G) SA; (H) GA3 and (I) IAA with the active site amino acid resides of OsMDHAR. The amino acids shown are the ones with which the metabolites form hydrogen bonds and stearic interactions.
Antioxidants 11 00458 g008
Table
Table 1. Details of the compounds/ligands used in the study. XlogP3: Octanol-water partition coefficient; HBD: Number of hydrogen bond donor acceptor; HBA: Number of hydrogen bond acceptor; TPSA: Topological polar surface area.
Table 1. Details of the compounds/ligands used in the study. XlogP3: Octanol-water partition coefficient; HBD: Number of hydrogen bond donor acceptor; HBA: Number of hydrogen bond acceptor; TPSA: Topological polar surface area.
LigandPubChem IDMol. FormulaMol. Wt (g/mol)XlogP3HBDHBATPSA (Å2)
Salicylic acid338C7H6O3138.122.32357.5
IAA802C10H9NO2175.181.42253.1
Serotonin5202C10H12N2O176.210.23262
GA36466C19H22O6346.40.236104
GSH124886C10H17N3O6S307.33−4.568160
DHA440667C6H6O6174.11−126101
Jasmonic acid5281166C12H18O3210.271.61354.4
MDHA53262277C6H7O6175.12−12590.8
Ascorbic acid54670067C6H8O6176.12−1.646107
Aluminium ion104727Al3+26.981538-000
Manganese ion27854Mn2+54.93804-000
Table
Table 2. Docking scores (MolDock, Rerank and Hydrogen bond) of DHAR of rice, maize and wheat with different ligands, including ions.
Table 2. Docking scores (MolDock, Rerank and Hydrogen bond) of DHAR of rice, maize and wheat with different ligands, including ions.
LigandsOsDHARZmDHARTaDHAR
MolDockRerankHBondMolDockRerankHBondMolDockRerankHBond
SA−53.61−42.70−4.22−65.25−57.66−8.19−58.23−52.20−6.82
IAA−87.50−69.73−5.15−84.28−60.29−4.08−80.75−63.68−4.90
5-HT−72.01−58.53−4.99−92.70−50.79−8.91−79.74−64.88−2.50
GA3−84.02−70.38−3.46−99.20−78.81−5.72−92.83−58.92−5.62
GSH−91.60−75.41−5.74−103.79−85.35−8.92−96.60−82.38−9.31
JA−93.23−72.72−4.60−96.34−18.03−5.00−84.62−65.49−4.62
AsA−62.38−56.41−8.61−72.86−67.44−9.82−70.58−63.03−13.44
Al3+−25.81−30.590.00−41.36−34.430.00−22.05−21.930.00
Mn2+−26.98−30.410.00−41.37−34.420.00−23.44−21.660.00
Table
Table 3. Docking scores (MolDock, Rerank and Hydrogen bond) of MDHAR of rice, maize and wheat with different ligands, including ions.
Table 3. Docking scores (MolDock, Rerank and Hydrogen bond) of MDHAR of rice, maize and wheat with different ligands, including ions.
LigandsOsMDHARZmMDHARTaMDHAR
MolDockRerankHBondMolDockRerankHBondMolDockRerankHBond
SA−69.24−56.19−5.00−71.05−59.15−4.82−70.20−59.84−5.78
IAA−102.96−83.96−4.81−100.38−79.58−5.73−98.33−81.03−3.17
5-HT−102.67−84.13−6.06−99.78−78.83−4.12−103.73−84.33−5.87
GA3−116.56−78.70−6.36−117.03−60.76−7.58−119.83−95.04−5.81
GSH−126.49−78.36−12.95−125.39−91.09−15.93−141.58−122.03−18.32
JA−111.14−91.62−4.44−109.56−89.51−7.06−114.21−92.53−6.59
AsA−81.00−66.94−8.07−88.97−80.97−13.70−89.68−78.91−12.48
Al3+−41.24−41.050.00−39.50−31.690.00−32.45−32.750.00
Mn2+−41.24−41.050.00−39.51−31.700.00−32.45−32.750.00
Table
Table 4. Correlation analysis of docking scores of the DHAR of rice, maize and wheat with different ligand properties. NC: Not correlated.
Table 4. Correlation analysis of docking scores of the DHAR of rice, maize and wheat with different ligand properties. NC: Not correlated.
Ligand PropertyOsDHARZmDHARTaDHAR
MolDockRerankHbondMolDockRerankHbondMolDockRerankHbond
Mol. Wt0.610.64NCNC0.51NCNCNCNC
XlogP3NCNC−0.51NCNC−0.58NCNC−0.58
HBDNCNC0.60NC0.72NCNCNCNC
HBANCNC0.73NCNCNCNCNC0.68
TPSANCNC0.74NC0.53NCNCNC0.58
Table
Table 5. Correlation analysis of docking scores (MolDock, Rerank and Hydrogen bond) of the MDHAR of rice, maize and wheat with different ligand properties. NC: Not correlated.
Table 5. Correlation analysis of docking scores (MolDock, Rerank and Hydrogen bond) of the MDHAR of rice, maize and wheat with different ligand properties. NC: Not correlated.
Ligand PropertyOsMDHARZmMDHARTaMDHAR
MolDockRerankHbondMolDockRerankHbondMolDockRerankHbond
Mol. Wt0.57NCNC0.54NCNC0.590.54NC
XlogP3NCNC−0.95NCNC−0.87NCNC−0.93
HBDNCNC0.73NCNC0.79NCNC0.79
HBANCNC0.83NCNC0.85NCNC0.83
TPSANCNC0.90NCNC0.88NCNC0.89
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Choudhury, S.; Mazumder, M.K.; Moulick, D.; Sharma, P.; Tata, S.K.; Ghosh, D.; Ali, H.M.; Siddiqui, M.H.; Brestic, M.; Skalicky, M.; et al. A Computational Study of the Role of Secondary Metabolites for Mitigation of Acid Soil Stress in Cereals Using Dehydroascorbate and Mono-Dehydroascorbate Reductases. Antioxidants 2022, 11, 458. https://doi.org/10.3390/antiox11030458

AMA Style

Choudhury S, Mazumder MK, Moulick D, Sharma P, Tata SK, Ghosh D, Ali HM, Siddiqui MH, Brestic M, Skalicky M, et al. A Computational Study of the Role of Secondary Metabolites for Mitigation of Acid Soil Stress in Cereals Using Dehydroascorbate and Mono-Dehydroascorbate Reductases. Antioxidants. 2022; 11(3):458. https://doi.org/10.3390/antiox11030458

Chicago/Turabian Style

Choudhury, Shuvasish, Muhammed Khairujjaman Mazumder, Debojyoti Moulick, Parul Sharma, Sandeep Kumar Tata, Dibakar Ghosh, Hayssam M. Ali, Manzer H. Siddiqui, Marian Brestic, Milan Skalicky, and et al. 2022. "A Computational Study of the Role of Secondary Metabolites for Mitigation of Acid Soil Stress in Cereals Using Dehydroascorbate and Mono-Dehydroascorbate Reductases" Antioxidants 11, no. 3: 458. https://doi.org/10.3390/antiox11030458

APA Style

Choudhury, S., Mazumder, M. K., Moulick, D., Sharma, P., Tata, S. K., Ghosh, D., Ali, H. M., Siddiqui, M. H., Brestic, M., Skalicky, M., & Hossain, A. (2022). A Computational Study of the Role of Secondary Metabolites for Mitigation of Acid Soil Stress in Cereals Using Dehydroascorbate and Mono-Dehydroascorbate Reductases. Antioxidants, 11(3), 458. https://doi.org/10.3390/antiox11030458

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop