Comparative structural analysis of cytokinin dehydrogenase enzymes of O.sativa, A. thaliana and Zea mays leading to predict the best enzyme  

Pranati Swain , Lambodar Behera
Central rice research institute, cuttack, odisha, India
Author    Correspondence author
Computational Molecular Biology, 2014, Vol. 4, No. 14   doi: 10.5376/cmb.2014.04.0014
Received: 03 Dec., 2014    Accepted: 26 Dec., 2014    Published: 30 Dec., 2014
© 2014 BioPublisher Publishing Platform
This is an open access article published under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Preferred citation for this article:

Swain and Behera, 2014, Comparative structural analysis of cytokinin dehydrogenase enzymes of O.sativa, A. thaliana and Zea mays leading to predict the best enzyme, Computational Molecular Biology, Vol.4, No.14, 1-9 (doi: 10.5376/cmb.2014.04.0014)


The cytokinin dehydrogenase enzymes of rice, wheat and maize were selected for this study. The enzymes of this crops show the same function i.e, leading to high grain production. The CKX1, CKX2, CKX3, CKX4, CKX5, CKX6, CKX7 and CKX8 enzymes of these three crops vary from each other on the basis of their physico-chemical analysis. A brief study has been done on all these enzymes showing the similarities, domains, secondary structures, homology models, backbone confirmations, best generated models and functions. From the study we got that Q4ADV8 of O.sativa has the best secondary structure with strong helix, turns and sheets along with highest molecular weight. A2XVN3 of O.sativa has the longest domain region, Q9FUJ1 of A.thaliana has highest numbers of total negatively charged residues (ASP+GLU) and Q9LY71 of A.thaliana has highest numbers of total positively charged residues (ARG+LYS). The Q9LTS3 of A.thaliana has highest VDW radius from geometric center. Q8LNV6 is the best model depending upon the Z-score calculated by ANOLEA. from the above study it is concluded that the enzymes of O.sativa andA.thaliana has more strong physico-chemical characteristics than in comparision to the zea mays.

Cytokinin dehydrogenase enzymes; Physico-chemical analysis; Domain region prediction; Predicting the best enzyme

Due to the growth of population the demands for more grain production are increasing day by day. The main crops harvested in world are rice, maize and wheat. These are the main source of food. Human beings absorb many proteins from these crops. These crops are blessed with such miraculous proteins which help a lot in high grain production in order to fulfill the demands. Here we have made a brief study about the cytokinin dehydrogenase enzymes of rice, maize and wheat leading to high grain production. cytokinin dehydrogenase belongs tooxidoreductasefamily and catalyzes the chemical reaction. Cytokinin dehydrogenase is otherwise known as N6-dimethylallyladenine: (acceptor) oxidoreductase, 6-N-dimethylallyladenine: acceptor oxidoreductase, OsCKX2, CKX, and cytokinin oxidase/dehydrogenase. cytokinin dehydrogenase is helpfull in the degradation of cytokininsisopentenyladenine, zeatin, and their ribosides (Chmülling et al., 2003). Gn1a gene of rice possessescytokinin oxidase/dehydrogenase (OsCKX2), an enzyme which degrades the phytohormonecytokinin (Ashikari et al., 2005). The Gn1a/OsCKX2 of Oryza sativa L, grain number 1a/Cytokinin oxidase 2 gene, encodes a cytokinin oxidase, which acts as a major quantitative trait locus contributing to grain number improvement in rice breeding practice (Lia et al., 2012). In wheat the CKX2 expression is activated by the IKU transcription factor WRKY10 directly and promotes endosperm growth (Li et al., 2013).In wheat the growth control of endosperm by CKX2 integrates genetic and epigenetic regulations. In angiosperms, cytokines are highly active in endosperm, IKU effectors coordinate environmental and physiological factors, resulting in modulation of seed size. Higher activity of the ckx3 ckx5 in wheat increase in seed yield, highlighting the relevance of sink strength as a yield factor. CKX3 and CKX5 regulate the activity of the reproductive meristems of Arabidopsis thaliana (Batrina et al., 2011). In this study we have considered the CKX1 of zea mays (Maize), CKX2 of Oryza Sativa (japonica), Arabidopsis thaliana (Mouse-ear cress), CKX3 of Oryza Sativa (japonica), Arabidopsis thaliana (Mouse-ear cress), CKX4 of Oryza Sativa (japonica), CKX5 of Arabidopsis thaliana (Mouse-ear cress), CKX6 of Arabidopsis thaliana (Mouse-ear cress) and CKX7 of Arabidopsis thaliana (Mouse-ear cress) and CKX8 of Oryza Sativa(indica). The main objective of the study is to find out their physic-chemical properties, similarities and differences among all the enzymes. However the function of all the eznzymes is same, but their properties are explained below.
Material and Methods
The cytokinin dehydrogenase enzymes of rice, wheat and maize were selected in this study due to its role in high grain production in crops.
Sequence retrieval
The cytokinin dehydrogenase enzymes of O.sativa (japonica), A.thaliana (mouse-ear cress) and zea mays (maize) were used for this study. Their amino acid sequences were retrived from uniprot (http://www., which is given in Table 1.

Table 1 The informations of all enzymes which were collected from Uniprot

Sequence alignment
The cytokinin dehydrogenase enzymes of all the plants were subjected for multiple sequence alignment using ClustalW2 ( Tools/msa/clustalw2/) analysis (Thompson et al., 1994). The alignment file was saved in clustal format and conserved residues were identified successfully. A phylogenetic tree was generated. The conserved residues of all the enzymes are shown in Figure 1.

Figure 1 the multiple sequence alignment result of all the enzyme

Physico-Chemical Characterization of all the enzymes
The amino acid sequences of all the enzymes were subjected to ProtParamtool ( protparam/) for further analysis, from where we got the physical and chemical properties of all the enzymes including molecular weight, theoretical PI, extinction coefficients, instability index, aliphatic index, grand average of hydrophobicity, total number of positively and negatively charged residues [9]. Detail information is shown in Table 2, Table 3, Table 4.
Secondary structure prediction of enzymes
The CFFSP server ( chou-fasman/) was used for secondary structure analysis where we submitted the amino acid sequences as an input and the result showed the details about helices, sheets and turns respectively (David et al., 1999). The detail information is shown in Table 5 and Figure 2-Figure 11.

Table 2 Physico-chemical characterization of enzymes

Table 3 Information about the positively and negatively charged residues of all enzymes

Table 4 Domains of proteins

Table 5 The secondary structural analysis of all enzymes

Figure 2 Q67YU0-CKX5 (A.thaliana)

Figure 3 Q4ADV8-CKX2(O.sativa)

Figure 4 Q5JLP4-CKX4(O.sativa)

Figure 5 Q8LNV6-CKX3(O.sativa)

Figure 6 Q9LTS3-CKX3 (A.thaliana)

Figure 7 Q9FUJ3-CKX2 (A.thaliana)

Figure 8 Q9FUJ1-CKX7(A.thaliana)

Figure 9 Q9LY71 (A.thaliana)

Figure 10 A2XVN3-CKX8(O.sativa)

Figure 11 Q9T0N8-CKX1 (Zea mays)

Structural alignment
As the NMR crystallographic structures of these enzymes are not available at PDB, so the tertiary structures of all the enzymes of plants were designed by using phyre2 server whichgenerate reliable protein models when other widely used methods such as PSI-BLAST cannot. Phyre2 server (http://www.sbg. generate the model using the principles and techniques of homology modeling. The target can be modelled with reasonable accuracy on a very distantly related sequence of template. The phyre2 server uses a profile-profile alignment algorithm based on each proteins position-specific scoring matrix. Phyre2 server include protein structure prediction, function prediction, domain prediction, domain boundary prediction, evolutionary classification of proteins, guiding site-directed mutagenesis and solving protein crystal structures by molecular replacement (Christie et al., 2012, Bilal et al., 2013, Singh et al., 2009). The homology models of enzymes are shown in Figure 12-Figure 21 and the detail information after visualization is shown in Table 6.

Figure 12 A2XVN3-CKX8(O.sativa)

Figure 13 Q4ADV8-CKX2(O.sativa)

Figure 14 Q5JLP4-CKX4(O.sativa)

Figure 15 Q8LNV6-CKX3(O.sativa)

Figure 16 Q9FUJ1-CKX7(A.thaliana)

Figure 17 Q9FUJ3-CKX2 (A.thaliana)

Figure 18 Q9LTS3-CKX3 (A.thaliana)

Figure 19 Q9LY71 (A.thaliana)

Figure 20 Q9T0N8-CKX1 (Zea mays)

Figure 21 Q67YU0-CKX5 (A.thaliana)

Table 6 properties of all enzymes which were obtained by using pymol visualiser

Model validation and optimization
The final tertiary structures of the enzymes were visualized usigPyMol visualization tool. Then by using PyMolvisualiser we got the number of atoms present in all the enzymes, formal charge sum, and partial charge sum, molecular and solvent accessible surface area. The backbone confirmation of all the enzymes were checked using Rampage server ( ~rapper/rampage.php) from where we got the residues lying in allowed, favoured and outlier regions [13]. Then the models were subjected to ANOLEA-swiss model from where we got the Z-score, QMEAN score. The model with least Z-score is generally selected as the best model. The detail information is shown in Table 7.

Table 7 Backbone confirmation of all enzymes

Result and Discussion
Sequence retrieval analysis
The amino acid sequences were successfully retrived from uniprot, which showed the following informations.
Multiple sequence alignment result
From the MSA analysis we got that there are many conserved residues among all the enzymes i.e, D, F, G, V, P, L, H, S, Q, W, Y, N, C, T, R, A, E, V, P, K, I. the similarities are denoted with *, dissimilarities with . and gaps with -.
Physico-chemical analysis result
From the physico-chemical analysis it is found that the A2XVN3 enzyme of O.sativa has longest domain region whereas Q9FJU3 of A.thaliana has smallest domain region. The Q67YU0 of A.thaliana has highest molecular weight, A2XVN3 enzyme of O.sativa has highest theoretical PI and aliphatic index, Q4ADV8 enzyme of O.sativa has highest total number of atoms.
Secondary structure prediction result of all enzymes
The secondary structures of all the enzymes were predicted using CFFSP server from where we could got the percentage of alpha, helix, turns. The helix is denoted with red colours, sheet is denoted with green colours, turn is denoted with blue colours. The secondary structures of all the enzymes are given from Figure 2-Figure 11.
Homology modeling results
The models which were generated by using phyre2 server were visualized by using PyMol visualization tool, which showed the following result.
Backbone confirmation result of enzymes
The final generated models of all the enzymes were submitted to RAMPAGE server from where we got the residues lying in allowed, favoured and outlier regions. The detail information is given below.
From the above study we got the functions of all the enzymes, the most common residues among them, from the physico-chemical analysis we got that A2XVN3 of O.sativa has the most longest domain region. The best model which was selected on the basis of Z-score was Q8LNV6 ofO.sativa. All the enzymes has most common domain i.e, FAD binding PCMH type. Depending upon the back-bome confirmation of enzymes the best model was Q9T0N8 of zea mays. The function of all the enzymes are same leading to high grain production, however they are having many structural and chemical differences which is represented in this study.
Ashikari et al., 2005, Cytokinin Oxidase Regulates Rice Grain Production. Science journal, vol.309, No.5735, 741-745
Batrina et al., 2011, Cytokinin regulates the activity of reproductive meristems, flower organ size, ovule formation, and thus seed yield in Arabidopsis thaliana, The Plant Cell journal, vol.23, No.1, 69-80
Bilal et al., 2013, Generation of a 3D model for human cereblon using comparative modeling.ournal of Bioinformatics and Sequence Analysis Vol. 5(1), pp. 10-15
chmülling et al., 2003, Structure and function of cytokinin oxidase/dehydrogenase genes of maize, rice, Arabidopsis and other species. Journal of Plant Research. Vol.116, No.3,241-52
Christie et al., 2012, phylogenetic analysis and Homology modelling of laccase gene of fusarium Oxysporum lcc2. International Journal of Pharma Sciences and Research (IJPSR), Vol 3, No 12
David et al., 1999, Protein Secondary Structure Prediction Based on Position-specific Scoring Matrices, journal of molecular biology, vol.292, 195-202
Lia et al., 2012, Rice zinc finger protein DST enhances grain production through controlling Gn1a/OsCKX2 expression
Li et al., 2013, Integration of epigenetic and genetic controls of seed size by cytokinin in Arabidopsis. Plant biology journal, vol.110, No.38, 15479-15484
Singh et al., 2009, Comparative modeling and analysis of 3-D structure of Hsp 70, in Cancer irroratus An International Journal, vol.1, No.2, 1-4

Thompson et al., 1994, CLUSTAL W:improving the sensitivity of progressive multiple sequence alignment through sequence weighting,position-specific gap penalties and weight matrix choice, Nucleic Acids Research, Vol. 22, No. 22, 4673-4680

Computational Molecular Biology
• Volume 4
View Options
. PDF(727KB)
. Online fPDF
Associated material
. Readers' comments
Other articles by authors
. Pranati Swain
. Lambodar Behera
Related articles
. Cytokinin dehydrogenase enzymes
. Physico-chemical analysis
. Domain region prediction
. Predicting the best enzyme
. Email to a friend
. Post a comment