Computational Molecular Biology 2014, Vol. 4, No. 15, 1-5
http://cmb.biopublisher.ca
5
2.2 Alignment between domain region and whole
enzyme
The domain region was extracted from the whole
protein sequence and subjected to clustalw2 for
multiple sequence alignment. From this study the
conserved residues among domain region and whole
enzyme sequence were obtained. The match regions
were denoted with *, mis-match regions with . and the
gaps were denoted with -. Here the Figure 1 shows the
common residues.
2.3 Homology modeling analysis
The homologous models for both enzyme and domains
were generated. The helices are denoted red colours,
helices with yellow colours and the loops are denoted
with green colours respectively. The VDW radius, beta
factor, mass of structures, stability of the objects were
found from yasara tool. Then total number of atoms,
formal charge sum, molecular surface area and solvent
accessible surface area were found from pymol tool.
2.4 Model validation and optimization
The backbone confirmations of the structures were
generated from rampage server. 91.8% of the residues lie
in most favored region whereas 3.3% of the residues lie
in outlier region. The QMEAN score and Z-scores were
found from both anolea server and from prosa server
.
2.5 Active site analysis result
The active sites were found from castp server. The
balls with green color show the binding regions of
whole protein and the domain region. Many pockets
were found, but out of them the first volume was
selected. The pockets of domain were mainly ALA,
ARG, CYS, HIS, GLY, VAL, SER, TYU, LEU, PHE,
ILE, THR, PRO and similarly the pockets of whole
enzyme were PHE, ILE, LEU, VAL, SER, HIS, GLY,
CYS, ALA, MET, GLU, THR, ASP, TYR, ASN, GLN,
PRO, TRP, GLN respectively.
2.6 Protein-ligand interaction
The pdb files of protein, domain and the ligand were
submitted and by adding polar hydrogens to protein
and domain the interaction was studied. The binding
energy which was found between the interactions of
whole protein-ligand (FAD) was -11.23, with -14.81
intermolecular energy, 3.58 torsion energy and 13.66
internal energy. Similarly the interaction which was
performed between the domain region and ligand
(FAD) showed -10.91 binding energy, -14.49
intermolecular energy, 3.58 torsion energy and 14.04
internal energy respectively.
From the interaction it is found that the FAD is
binding to ALA6, ARG1, AR536, TYR10, ASP55,
ALA57, CYS38 at the domain part and in whole
protein the domain is binding to ARG134, ASP128,
SER131, TYR83, SER85, LEU135, ARG86, GLN136,
LEU135, LEU47.
3 Discussions
From the above study it is found that the model of whole
enzyme and the domain are suitable for future analysis as
most of the residues lie in favored region. The main
objective of this study was to find the stronger interaction
of ligand (FAD) with the domain part of enzyme and
with the whole sequence, which showed that the binding
interaction is stronger between the domain part and FAD
than in comparison to the interaction between FAD and
whole amino acid sequence
References
Ashikari et al., 2005, Cytokinin Oxidase Regulates Rice Grain Production,
science journal, vol.309, No.5735,741-745
Bilal et al., 2013, Generation of a 3D model for human cereblon using
comparative modeling, journal of bioinformatics and sequence analysis,
vol. 5, No.1, 10-15
Gao et al., 2014, CYTOKININ OXIDASE/DEHYDROGENASE4
Integrates Cytokinin and Auxin Signaling to Control Rice Crown Root
Formation, plant physiology journal, vol.165, 1035-1046
Jasim et al., 2013, Protein-Ligand interaction studies on 2, 4, 6-trisubstituted
triazine derivatives as anti-malarial DHFR agents using AutoDock,
IOSR Journal of Applied Physics, vol.4, No.5, 2278-4861
Li et al., 2013, Rice zinc finger protein DST enhances grain production
through controlling Gn1a/OsCKX2 expressio
lant biology, vol. 10,
No.8, 3167-72.
Mameaux et al., 2011, Molecular, phylogenetic and comparative genomic
analysis of the
cytokinin oxidase/dehydrogenase
gene family in the
Poaceae, plant biotechnology journal, vol.10, No.1, 67-82
Rambabu et al., 2012, Computational analysis and function prediction of a
hypothetical protein 1RW0, International Journal of Computational
Bioinformatics and In Silico Modeling, vol.1, No.5, 58-62
Sali et al., 1993, Comparative Protein Modelling by Satisfaction of Spatial
Restraints, Journal of Molecular Biology, vol.234, No.3, 779-815
Seeliger et al., 2010, Ligand docking and binding site analysis with PyMOL
and Autodock/Vina, J Comput Aided Mol Des, vol.4, 417-422
Sungh et al., 2004, Comparative modeling and analysis of 3-D structure of
Hsp 70, in Cancer irroratus, Biological Forum — 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 journal, vol. 22, 4673-468