Modeling and analysis of soybean (Glycine max. L) Cu/Zn, Mn and Fe superoxide dismutases

Modeling and analysis of soybean (Glycine max. L) Cu/Zn, Mn and Fe superoxide dismutases

V. Ramana Gopavajhula^I; K. Viswanatha Chaitanya^I; P. Akbar Ali Khan^II; Jilani P. Shaik^II; P. Narasimha Reddy^II; Mohammad Alanazi^II

^IDepartment of Biotechnology, GITAM University, Visakhapatnam, India

^IIDepartment of Biochemistry, College of Science, King Saud University, Riyadh, Kingdom of Saudi Arabia

^{Correspondence}

ABSTRACT

Superoxide dismutase (SOD, EC 1.15.1.1) is an important metal-containing antioxidant enzyme that provides the first line of defense against toxic superoxide radicals by catalyzing their dismutation to oxygen and hydrogen peroxide. SOD is classified into four metalloprotein isoforms, namely, Cu/Zn SOD, Mn SOD, Ni SOD and Fe SOD. The structural models of soybean SOD isoforms have not yet been solved. In this study, we describe structural models for soybean Cu/Zn SOD, Mn SOD and Fe SOD and provide insights into the molecular function of this metal-binding enzyme in improving tolerance to oxidative stress in plants.

Keywords: amino acid analysis, model evaluation, molecular modeling, phylogenetic analysis, superoxide dismutase (SOD).

Introduction

Crop plants are frequently exposed to a variety of abiotic, biotic and xenobiotic stresses that cause injury, limit their growth and adversely affect their productivity. The most common result of such stress is the production of toxic reactive oxygen species (ROS). Increased levels of ROS such as superoxide (O2•^-), hydrogen peroxide (H₂O₂) and the hydroxyl radical (•OH) cause irreparable damage to cellular components such as DNA, proteins and lipids that require additional defense mechanisms (Blokhina et al., 2003). Plant responses to ROS toxicity involve the coordinated action of enzymatic and non-enzymatic antioxidant defense systems (Pallavi and Rama Shankar, 2005; Pallavi et al., 2012). Among enzymatic defenses, superoxide dismutase (SOD, EC 1.15.1.1) is the most important enzyme because of its distinct ability to neutralize superoxide anions by dismutating them into O₂ and H₂O₂. SOD is synthesized by all aerobic organisms and also by some air-tolerant and obligate anaerobic organisms (Fink and Scandalios, 2002).

SODs are metalloproteins that occur in four isoforms: copper zinc SOD (Cu/Zn SOD), manganese SOD (Mn SOD), nickel SOD (Ni SOD) and iron SOD (Fe SOD), all of which are highly stable because of the β-barrel structure and low content of α-helix strands (Renu and Sabarinath, 2004). Almost all eukaryotic organisms synthesize Mn and Cu/Zn SOD whereas Fe SOD is exclusive to plants (Kwang-Hyun et al., 2006). Ni SOD was recently reported in Streptomyces griseus and S. coelicolor (Garcia-Hernández et al., 2002). In the Arabidopsis thaliana genome, three Cu/Zn SODs (cytosolic, chloroplast and peroxisomal), one Mn SOD and three Fe SODs (Fe SOD1, 2 and 3) genes have been reported, with Fe SOD2 and 3 being chloroplastic. Chloroplast Cu/Zn SOD is located in the thylakoid membranes and Fe SOD2 and 3 are located in the stroma (Myouga et al., 2008).

Soybean (Glycine max. L) is an important legume plant that is widely cultivated for its protein and oil and is considered a miracle crop. Although the role of SOD as an antioxidant enzyme under stress conditions has been studied in different varieties of soybean (Chaitanya et al., 2009; Hassan et al., 2011), the genes that control the expression of these isoforms have not been identified. In this study, we undertook a molecular, structural and phylogenetic analysis of soybean SOD isoforms (Cu/Zn, Mn and Fe SOD) based on homology modeling using A. thaliana SODs.

Materials and Methods

Identification of EST sequences

The sequences of genes corresponding to SOD isoenzymes in Arabidopsis, namely, Mn SOD (AT3G10920), Fe SOD (AT4G25100) and Cu/Zn SOD distributed in chloroplasts (ATG12520), cytosol (AT1G08830) and peroxisomes (AT5G18100), were used as the query sequences in BLAST searches to identify the corresponding genes in soybean. These nucleotide sequences were assessed for homology at the protein level and used for phylogenetic analysis.

Phylogenetic analysis of sequences

All sequence alignments were done using ClustalW (Aiyar, 2000). Phylogenetic trees were plotted using MEGA 4.0 software with the UPGMA method (Tamura et al., 2007). Soybean genes that separated along with Arabidopsis genes were pooled and their amino acid pattern was analyzed by constructing pretty boxes using Boxshade.

Secondary structure prediction of SOD proteins

The secondary structures of the soybean SOD isoenzymes were predicted using the PSIPRED online server based on the retrieved sequences. PSIPRED incorporates two feed-forward neural networks that analyze the data generated as an output from PSI-BLAST (position-specific iterated BLAST) (Altschul et al., 1997). Validation of the procedure and performance using PSIPRED yielded an average Q3 score of 76.5%.

3D structure

The 3D structure models for soybean and Arabidopsis SODs were developed using 3D LigandSite. This software was used to the exact binding site of metal ions in the amino acid sequences. Dompred software was used to predict the domains and their boundaries for a given protein sequence.

Quaternary structure prediction

The quaternary structures of SOD proteins in soybean and Arabidopsis were predicted using the protein interfaces and surfaces tool PISA. Assemblies that could form crystals were determined by identifying the sets that represented the solutions indicated in the headings of the appropriate table (see Results). The highest values in the assemblies were considered to be the most appropriate. The MM size indicated the number of macromolecular monomeric units in that particular assembly and corresponded to an oligomeric or multimeric state. A formula was obtained to indicate the chemical composition of the assembly and denoted the number of different monomeric units. The stability of an assembly, i.e., its tendency to dissociate in solution, was also determined. The solvation free energy (∆G^int) was calculated as the difference in solvation energies of the isolated and assembly structures and indicated the free energy gain (kcal/M) during the formation of an assembly. The free energy of dissociation (∆G^diss) represented the free energy difference between the associated and dissociated states. Assemblies with ∆G^diss > 0 were thermodynamically more stable because positive values were included in external energy use during the dissociation of an assembly.

Model evaluation

The dihedral angles ϕ vs. ψ of amino acid residues in the protein structures were visualized and analyzed with Ramachandran plots (Ramachandran et al., 1963). The evaluation of models predicted in silico is essential in order to avoid errors resulting from trivial and non-trivial mistakes. To avoid ambiguities and to improve accuracy, the predicted SOD models were evaluated using the ProSA and VADAR web servers. For a specific PDB structure, ProSA calculates the overall quality score and validates a low resolution structure for approximate models using C-alpha atoms of the input structure. The output provides a z-score for the model that indicates the overall model quality; this value was determined from the plot during prediction.

Results and Discussion

ROS produced by plants are eliminated by antioxidant defense systems that enhance the tolerance of plants to environmental stress (Min-Lang et al., 2012). In view of the increasing interest in the molecular modeling of the various isoforms of SOD, in this study we investigated the structure of soybean SOD isoforms and examined their phylogenetic relationships.

Phylogenetic analysis

The availability of information from various genome-sequencing projects, cDNA libraries and EST libraries offers the possibility of complementing investigations of gene function in vivo with parallel phylogenetic analyses of multigene families to address their evolution within and across species (Vincentz et al., 2003). Within families, the protein structure and catalytic residues that determine the substrate specificity are generally conserved. Bioinformatics tools are thus useful for the functional analysis of related proteins (Henrissat et al., 2001). However, many sequence-based families are polyspecific, i.e., they include genes that encode proteins with different functions. This reflects gene duplication and evolutionary divergence, with the acquisition of new protein functions (Emanuele et al., 2004). In the present study, the phylogenetic relationships of soybean SOD genes were evaluated with respect to Arabidopsis SOD genes by using the Maximum Composite Likelihood (MCL) approach implemented in MEGA (Tamura et al., 2007).

Phylogenetic analysis of the soybean and Arabidopsis open reading frames (ORFs) provided information on the evolutionary ancestry of all the SOD groups. This analysis showed that SODs segregated into two major clusters, with cytosolic, chloroplast, peroxisomal Cu/Zn in one cluster and Mn SOD and Fe SOD in another. In this tree, soybean SOD TC332577 segregated with chloroplast Cu/Zn SOD, TC287018 with peroxisomal Cu/Zn SOD, TC282951 with the genes of Arabidopsis cytosolic Cu/Zn SOD, TC278165 with Mn SOD and TC278336 with Fe SOD (Figure 1). Arabidopsis and soybean SODs were grouped with the same branch lengths, i.e., 0.0870 for Mn SOD, 0.1257 for Fe SOD, 0.0748 for chloroplast Cu/Zn SOD, and 0.0912 and 0.1675 for cytosolic and peroxisomal clusters, respectively. This homology in grouping reflected the strong similarities in the gene patterns of these two plants.

However, there was a subtle difference in the branch lengths of the major groups. Cu/Zn SOD segregated into a major group whereas the peroxisomal enzyme of both plants grouped together with a branch length of 0.1675. Chloroplast and cytosolic enzymes grouped together with branch lengths of 0.1465 and 0.1301, respectively; they were joined to the peroxisomal enzymes via branch lengths of 0.0714 and 0.0177, respectively. The difference between these branch lengths was < 0.075. Similarly, Mn SOD and Fe SOD grouped with branch lengths of 0.4680 and 0.4294, respectively; these two major groups were linked by a branch length difference of 0.32. The UPGMA tree showed that, the difference between two branch lengths at each cluster was < 0.5 and in some cases almost zero. This result shows that these SODs are closely related to each other and that the sequences retrieved were accurate. The UPGMA tree showed that the SOD genes identified here are important and deserve further investigation.

Boxshade analysis

The comparison of homologous protein sequences is the most effective means of identifying common active sites or binding domains. Comparative studies of protein sequences allow the functional relationships among proteins to be determined and are particularly important for homology searches and threading methods in structure prediction. The alignment of multiple protein sequences is a powerful tool for grouping proteins into families and allows subsequent analysis of evolutionary issues (Balasubramanian et al., 2012). In the present study, the pattern of conserved amino acids in the soybean and Arabidopsis SOD protein sequences was studied using the Boxshade server, which split the sequences into two clusters with Cu/Zn SOD forming one cluster and Mn SOD and Fe SOD forming the second cluster (Figure 2). Fe and Mn SODs showed high similarities in sequence and structure. Rice Fe and Mn SODs also share high homology in their amino acid sequences. Mn SOD is the only form of SOD that is essential for the survival of aerobic life and plants. Mn SODs share 65% sequence similarity with each other (Youxiong et al., 2012). The degree of homology was also high among Cu/Zn SOD genes compared to that of Fe and Mn SODs. Our results indicated that the protein sequences from Cu/Zn SOD of chloroplasts, cytosol and peroxisomes had a greater number of conserved amino acid sequences than Mn and Fe SODs. The subcellular and phylogenetic distribution of SODs showed that all three SOD isoforms co-exist only in plants (Bowler et al., 1994). Comparative sequence analysis of the three SOD isoforms suggests that Fe SODs and Mn SODs are more efficient than Cu/Zn SODs, and that Fe and Mn SODs most probably arose from common ancestral enzymes, whereas Cu/Zn SODs evolved separately in eukaryotes (Smith and Dolittle, 1992).

Secondary structure analysis

The secondary structure predictions for soybean and Arabidopsis Cu/Zn, Fe and Mn SOD proteins showed that Mn SODs had a long chain length consisting of α-helices and β-strands (Figure 3). Helices were absent in chloroplast, cytosolic and peroxisomal Cu/Zn SODs of soybean and Arabidopsis and their secondary structures were identical (Table 1). Soybean and Arabidopsis SOD proteins had a similar number of domains but their locations differed. The binding sites of the SOD proteins also differed in both plants (Table 2). The heterogen counts in the SOD genes of both plants were similar with respect to the type of heterogens present in SOD proteins (Table 3).

3D structure analysis

Proteins are complex chemical entities with a large number of variable atoms and a convoluted topology that make their description complicated (Ingale and Chikhale, 2010). The 'indescribable nature' of proteins also makes the quality of an experimentally determined protein structure very difficult to assess. The rapid increase in the number of genomes being sequenced and in the number of genes being deposited in databases means there is a need to identify the protein functions involved in protein interactions that form the basis of defining protein groups. In this study, three-dimensional models of soybean and Arabidopsis Cu/Zn, Mn and Fe SOD proteins were predicted using the software 3D LigandSite. The resulting models displayed excellent global and local stereochemical properties (Figure 4). Blue colored residues were predicted to be part of the binding. Residue conservation was calculated using the Jensen-Shannon divergence score (Capra and Singh, 2008). The ligands that formed the cluster were used to predict the metal ions shown in the space-filling format (Wass et al., 2010). There was a marked distinction between the metal binding sites and normal sites without space filling that enabled us to locate the coding regions exactly. The structural symmetry of the Cu/Zn SOD groups was identical.

Quaternary structural analysis

Quaternary structure plays an important role in defining protein function by facilitating allosterism and cooperativity in the regulation of ligand binding (Matthew et al., 2001). In this study, the quaternary structures of soybean and Arabidopsis SOD proteins were predicted using PISA by considering the assembly that provided the maximum structure size with good stability as being the best (Table 4). Our results indicated that soybean and Arabidopsis chloroplast Cu/Zn SODs were tetramers, whereas peroxisomal and cytosolic Cu/Zn SODs were dimers; Mn SOD was a tetramer and Fe SOD was a monomer. The Biomol stability value for chloroplast Cu/Zn SOD was predicted to be 9 while that for Mn SOD was 1. All of the predicted structures were stable because of their positive ∆G^Diss values. There were more surface area values than buried surface values for all of the structures, indicating that these proteins had fewer folds. There was considerable similarity in the gene pattern of the soybean and Arabidopsis SOD enzymes, particularly with respect to protein structure. This similarity indicated a uniform evolutionary gene bank that was largely undisturbed (few mutations), as seen in our preliminary protein structural analysis.

Despite similarities in the secondary structures of cytosolic and peroxisomal Cu/Zn SOD there were many differences between the corresponding genes in both plants. The amino acid patterns were identical in Cu/Zn SOD compared to Mn and Fe SOD. The Fe SOD structure contained more helices, as indicated by the quality index, with more omega aberrations. Cu/Zn SODs showed more homology compared to other models with no helices in their structures. Model evaluation revealed the accuracy of the predicted models and suggested possible errors that were trivial when compared to the overall quality of the proposed structure. Quaternary structural analysis showed that all of the structures are thermodynamically stable, with dissociation energies > 0.

Model evaluation

The Z-score value, a measure of model quality that predicts the total energy of the structure (Wiederstein and Sippl, 2007), was predicted for soybean and Arabidopsis SODs using the PROSA server (Figure 5). The Z-score values for chloroplast, cytosolic and peroxisomal Cu/Zn SODs were -6.33, -7.05 and -6.93, respectively; the corresponding values for Mn SOD and Fe SOD were -7.59 and -7.92, respectively.

The stereochemical quality and exactness of the predicted soybean and Arabidopsis SOD proteins were analyzed using residue-by-residue geometry and the overall geometry of the protein structures was analyzed with Ramachandran plots (Ramachandran et al., 1963). These plots help visualize the dihedral angles ϕ vs. ψ of amino acid residues in proteins. In a polypeptide, the main chain N-C α and C α-C bonds are relatively free to rotate. These rotations are represented by the torsion angles φ and ψ, respectively. Figure 6 shows the quantitative evaluation of soybean and Arabidopsis protein structures done using Ramachandran plots provided by the VADAR web server. Soybean and Arabidopsis Mn SODs shared > 90% structural similarity with an equal number of residues in favored, allowed and outlier regions. There was little variation in the general and proline residue numbers in the αR, αL and β regions of both structures. The prominence of residues in the αR and β regions suggested that the structure was rigid with more right helices. Soybean chloroplast Cu/Zn SOD showed higher quality compared to the Arabidopsis protein, with fewer residues in the outlier region and more residues in the allowed region. The absence of residues in the αL region in both structures was an interesting feature of this protein and indicated that the structure had no left helices. Chloroplast Cu/Zn SOD of both plants had 86% of residues in the favored region (this value was lower than in previously reported SOD structures). More residues were observed in the allowed region of soybean (9%) and Arabidopsis (8.5%) SOD. There were more outliers in the Arabidopsis structure (5%) than in soybean (4%). Cytosolic Cu/Zn SODs from both plants had no residues in their outlier regions. These soybean and Arabidopsis SODs also had an equal number of residues in the favored and allowed regions (95% and 4.5%, respectively). In both structures, three proline residues were distributed in the αR and β regions. There was a subtle difference in the sparsely populated αL region, where Arabidopsis had one general and one proline while in soybean both of the residues were generally in the core region. All of the Cu/Zn structures had fewer residues in the αR region compared to the β region. The Fe SOD and Mn SOD structures had good clustering of residues and a greater number of helices. The structural quality of chloroplast Cu/Zn SOD of both plants was lower than in the remaining structures. In all of the structures, the αL region was less much less populated, i.e., few or no residues.

The soybean Fe SOD structure had 95% of its residues in expected regions and a negligible proportion (1.3%) in the outlier region; the corresponding values for Arabidopsis were 93% and 2.4%, respectively. However, Arabidopsis had 5% of its residues in allowed regions whereas soybean had 3.6%. Soybean peroxisomal Cu/Zn SOD had < 89% of its residues in the favored region. The number of residues in allowed and outlier regions was also high, indicating structural aberrations whereas in Arabidopsis no residues are observed in the outlier region and ~94% of residues were in the favorable region, indicating the quality of the structure. Similar modeling and Ramachandran plot analyses to those described here have been used in the structural and functional analysis of spinach antioxidant proteins and the models were evaluated by computational tools (Sahay and Shakya, 2010).

Conclusion

Proteins are ubiquitous molecules that are involved in numerous crucial functions in organisms. Proteins accomplish their functions by positioning specific amino acids at target sites. Knowledge of the structural arrangement of amino acids is very important for understanding the molecular mechanisms by which proteins perform their functions. SOD is an important antioxidant enzyme that provides the first line of defense against ROS toxicity. The accurate and reliable molecular structural analysis of SOD isoenzymes is important for understanding their function in response to oxidative stress. In this study, structural models of soybean Cu/Zn, Mn and Fe SOD were analyzed and compared with those for Arabidopsis. These analyses provided insights into the molecular function of SOD isoenzymes with respect to their interactions with different cellular organelles. Further studies are in progress to understand the possible SOD gene interactions that may improve our understanding of the role of SODs in minimizing ROS toxicity.

Acknowledgments

The project was supported by the Research Center, College of Science, King Saud University, Saudi Arabia.

Internet Resources

Received: May 21, 2012

Accepted: February 10, 2013

Associate Editor: Adriana S. Hemerly

ERRATUM

V. Ramana Gopavajhula, K. Viswanatha Chaitanya, P. Akbar Ali Khan P, Jilani P. Shaik, P. Narasimha Reddy and Mohammad Alanazi. Modeling and analysis of soybean (Glycine max. L) Cu/Zn, Mn and Fe superoxide dismutases. Genetics and Molecular Biology 36(2), 225-236, 2013

The correct affiliation for the last four authors is:

Department of Biochemistry, College of Science, King Saud University, Riyadh, Kingdom of Saudi Arabia.

• Aiyar A (2000) The use of CLUSTAL W and CLUSTAL X for multiple sequence alignment. Meth Mol Biol 132:221-241.
• Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W and Lipman DJ (1997) Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res 25:3389-3402.
• Balasubramanian J, Shahul Hammed MK, Tamilselvan R and Vijayakumar N (2012) Artificial neural network: A forecast in pharmaceutical science. Nerve 1:7-12.
• Blokhina O, Virolainen E and Fagerstedt KV (2003) Antioxidants, oxidative damage and oxygen deprivation stress: A review. Ann Bot 91:179-194.
• Bowler C, Van Camp W, Van Montagu M and Inze D (1994) Superoxide dismutases in plants. Crit Rev Plant Sci 13:199-213.
• Capra J and Singh M (2008) Characterization and prediction of residues determining protein functional. Bioinformatics 24:1473-1480.
• Chaitanya KV, Spandana MS, Srinivas D and Kumar AL (2009) Antioxidative responses of soybean (Glycine max L) seedlings to salinity stress. J Plant Biol 36:83-87.
• Emanuele L, Yi L and McQueen-Mason SJ (2004) Phylogenetic analysis of the plant endo-β-1,4-glucanase gene family. J Mol Evol 58:506-515.
• Fink RC and Scandalios JG (2002) Molecular evolution and structure-function relationships of the superoxide dismutase gene families in angiosperms and their relationship to other eukaryotic and prokaryotic superoxide dismutases. Arch Biochem Biophys 399:19-36.
• Garcia-Hernández M, Berardini TZ, Chen G, Crist D, Doyle A, Huala E, Knee E, Lambrecht M, Miller N, Mueller LA, et al. (2002) TAIR: A resource for integrated Arabidopsis data. Funct Integr Genomics 2:239-253.
• Hassan M, Farrokh D, Jahanfar D, Ghorban N and Davood H (2011) Effects of water deficit stress on seed yield and antioxidants content in soybean (Glycine max L.) cultivars. Afr J Agric Res 6:1209-1218.
• Henrissat B, Coutinho PM and Davies GJ (2001) A census of carbohydrate-active enzymes in the genome of Arabidopsis thaliana Plant Mol Biol 47:55-72.
• Ingale AG and Chikhale NJ (2010) Prediction of 3D structure of paralytic insecticidal toxin (ITX-1) of Tegenaria agrestis (hobo spider). J Data Mining Genom Proteom 1:102-104.
• Kwang-Hyun B, Daniel ZS, Peng L and Xianming C (2006) Molecular structure and organization of the wheat genomic manganese superoxide dismutase gene. Genome 49:209-218.
• Matthew DG and Mark SH (2001) Quaternary structure of rice non-symbiotic hemoglobin. J Biol Chem 276:6834-6839.
• Min-Lang C, Li-Jen L, Jia-Hong L, Zin-Huang L, Yuan-Ting H and Tse-Min L (2012) Modulation of antioxidant defense system and NADPH oxidase in Plucheaindica leaves by water deficit stress. Bot Stud 53:45-54.
• Myouga F, Hosoda C, Umezawa T, Iizumi H, Kuromori T, Motohashi R, Shono Y, Nagata N, Ikeuchi M and Shinozaki K (2008) A heterocomplex of iron superoxide dismutases defends chloroplast nucleoids against oxidative stress and is essential for chloroplast development in Arabidopsis Plant Cell 20:3148-3162.
• Pallavi S and Rama Shankar D (2005) Lead toxicity in plants. Braz J Plant Physiol 17:35-52.
• Pallavi S, AmbujBhushan J, Rama Shanker D and Pessarakli M (2012) Reactive oxygen species, oxidative damage, and anti-oxidative defense mechanism in plants under stressful conditions. J Botany 217037.
• Ramachandran GN, Ramakrishnan C and Sasisek-haran V (1963) Stereochemistry of polypeptide chain configurations. J Mol Biol 7:95-99.
• Renu KC and Sabarinath S (2004) Heat-stable chloroplastic Cu/Zn superoxide dismutase in Chenopodiummurale Biochem Biophys Res Commun 320:1187-1192.
• Sahay A and Shakya M (2010) In silico analysis and homology modeling of antioxidant proteins of spinach. J Proteomics Bioinfo 3:148-154.
• Smith MW and Doolittle RF (1992) A comparison of evolutionary rates of the 2 major kinds of superoxide-dismutase. J Mol Evol 34:175-184.
• Tamura K, Dudley J, Nei M and Kumar S (2007) MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software ver. 4.0. Mol Biol Evol 24:1596-1599.
• Vincentz M, Bandeira-Kobarg C, Gauer L, Schlogl P and Leite A (2003) Evolutionary pattern of angiosperm bZIP factors homologous to the maize Opaque2 regulatory protein. J Mol Evol 56:105-116.
• Wass MN, Kelley LA and Sternberg MJ (2010) 3DLigandSite: Predicting ligand-binding sites using similar structures. Nucleic Acids Res 38:469-473.
• Wiederstein and Sippl MJ (2007) ProSA-web: Interactive web service for the recognition of errors in three-dimensional structures of proteins. Nucleic Acids Res 35:407-410.
• Youxiong Q, Jinxian L, Liping X, Jinrong G and Rukai C (2012) Molecular cloning and expression analysis of an Mn superoxide dismutase gene in sugarcane. Afr J Biotechnol 11:552-560.
• Basic Local Alignment Search Tool (BLAST), http://compbio.dfci.harvard.edu/cgi-bin/tgi/Blast/index.cgi (January 4, 2012).
  » link
• Boxshade (version 3.21, written by K. Hofmann and M. Baron), http://www.ch.embnet.org/software/BOX_form.html (January 4, 2012).
• Boxshade: For multiple sequence analysis http://sourceforge.net/projects/boxshade/ (January 5, 2012).
• Dompred, Server designed to predict putative protein domains http://bioinf.cs.ucl.ac.uk/dompred (January 5, 2012).
• 3D LigandSite, Prediction of ligand binding sites: http://www.sbg.bio.ic.ac.uk/3dligandsite/ (January 5, 2012).
• Expasy, Translation of a nucleotide (DNA/RNA) sequence: http://web.expasy.org/translate/ (January 7, 2012).
• Open Reading Frame Finder (ORF Finder), http://www.ncbi.nlm.nih.gov/projects/gorf/ (January 5, 2012).
• Protein Interfaces, Surfaces and Assemblies (PISA), http://www.ebi.ac.uk/msd-srv/prot_int/cgi-bin/piserver (January 5, 2012).
• Protein Structure Prediction Server (PSIPRED), http://bioinf.cs.ucl.ac.uk/psipred/ (January 14, 2012).
• Protein Structure Analysis (ProSA), https://prosa.services.came.sbg.ac.at/prosa.php (January 18, 2012).
• Volume, Area, Dihedral Angle Reporter (VADAR), http://redpoll.pharmacy.ualberta.ca/vadar (January 18, 2012).

Publication Dates

• Publication in this collection
  17 May 2013
• Date of issue
  2013

History

• Received
  21 May 2012
• Accepted
  10 Feb 2013