Structural Comparison, Substrate Specificity, and Inhibitor Binding of AGPase Small Subunit from Monocot and Dicot: Present Insight and Future Potential
ADP-glucose pyrophosphorylase (AGPase) is the first rate limiting enzyme of starch biosynthesis pathway and has been exploited as the target for greater starch yield in several plants. The structure-function analysis and substrate binding specificity of AGPase have provided enormous potential for understanding the role of specific amino acid or motifs responsible for allosteric regulation and catalytic mechanisms, which facilitate the engineering of AGPases. We report the three-dimensional structure, substrate, and inhibitor binding specificity of AGPase small subunit from different monocot and dicot crop plants. Both monocot and dicot subunits were found to exploit similar interactions with the substrate and inhibitor molecule as in the case of their closest homologue potato tuber AGPase small subunit. Comparative sequence and structural analysis followed by molecular docking and electrostatic surface potential analysis reveal that rearrangements of secondary structure elements, substrate, and inhibitor binding residues are strongly conserved and follow common folding pattern and orientation within monocot and dicot displaying a similar mode of allosteric regulation and catalytic mechanism. The results from this study along with site-directed mutagenesis complemented by molecular dynamics simulation will shed more light on increasing the starch content of crop plants to ensure the food security worldwide.
Starch is a basic constituent of the human and animal diet. It is an important carbohydrate considered as one of the primary energy sources for plants and a very important raw material for industrial processes. In many different plant species it has been demonstrated that ADP-glucose pyrophosphorylase (AGPase) (EC 126.96.36.199) is one of the major enzymes for starch biosynthesis. The overall crop yield potential is greatly influenced by the enzyme which modulates the photosynthetic efficiency in source tissues and determines the level of starch storage in sink tissues . Combined participation of AGPase, starch synthase, and branching enzyme is solely responsible for biosynthesis of starch in plant [2, 3]. In starch biosynthesis, AGPase is the first regulatory allosteric enzyme which converts ATP and glucose-1-phosphate (Glc1P) to adenosine-5′-diphosphoglucose (ADPGlc) and inorganic pyrophosphate (PPi) [4–8] (see Figure 1).
Mutant analysis and transgenic plant provide strong evidences of the allosteric properties of AGPase in controlling the rate of starch biosynthesis in higher plants [9–13]. In most cases the regulation of AGPase depends on the ratio of 3-phosphoglyceric acid and inorganic phosphate (3PGA/Pi) showing a direct correlation between the concentration of 3-PGA and starch accumulation and an inverse correlation between Pi concentration and the starch content . Although the overall kinetic mechanism of AGPase appears to be similar in bacteria and higher plants, their quaternary structures differ from each other . Bacterial AGPases are composed of four identical subunits (α) to form α4 homotetramer whereas plant AGPases are heterotetramer of two different yet evolutionarily related subunits containing a pair of identical small (SS or α) and identical large subunits (LS or β) to form α2β2 heterotetramer [14–17]. The two subunits vary in their molecular weight and genetic origin and are encoded by two different genes [16, 18]. Primary sequence analysis of LS and SS of AGPase has shown considerable sequence homology, suggesting a common evolutionary origin . Several researchers have reported that SS of AGPase has both catalytic and regulatory functions whereas LS has only regulatory function [19–24]. The hypothesis is well supported by former reports, reflecting SS is capable of forming a homotetramer with catalytic properties whereas LS is incompetent of forming an oligomeric structure with catalytic activities [1, 14, 20, 25]. In contrast, Kavakli et al. 2001  and Hwang et al. 2006 and 2008 [27, 28] suggested that the LS may bind to substrate ATP as well as glucose-1 phosphate and may allow the LS to interact in tandem with the catalytic SS influencing the net catalysis. In addition, specific regions of both the LS and the SS were found to be important for enzyme stability and subunit association . Study on chimeric maize/potato small subunits reflects a polymorphic motif of 55-amino acid region between the residues 322–376 plays a critical role during the interaction with LS and contributes to the overall stability of the enzyme . All these reports suggest that both the subunits are of equal importance for the catalysis and allosteric regulation of the enzyme.
Due to the difficulty of obtaining AGPase in stable form neither the LS nor the heterotetrameric AGPase (α2β2) atomic resolution structure from plant species has been solved yet. In 2005 Jin et al.  reported the first atomic resolution structure of AGPase SS from Solanum tuberosum. The crystal structure of SS was found in a homotetrameric form. Since then not a single crystallographic structure of AGPase has been reported. Although the AGPase gene(s) offers an attractive tool for engineering crop plants to enhance the yield potential of starch content, the understanding of structure-function relationships and the unique substrate specificity of AGPase has remained elusive.
In the absence of experimental three-dimensional structures, comparative modeling of protein is considered as one of the most accurate methods of model building and is often considered fundamental for understanding their function . This approach provides reasonable result based on the assumption that the tertiary structure of two proteins will be similar if they share high percentage of sequence similarity . It is widely being used when there is a clear relationship of homology between the target protein sequences and at least with an experimental (XRD or NMR) protein structure. Comparative structural analysis coupled with docking study has been immensely used for understanding the structure function relationship, mode of enzyme substrate interaction, and key residues involved in interaction without requiring further biochemical or immunological data [33–48]. Comparative modeling of AGPase has previously been studied for understanding the structure function relationship and to investigate the subunit interaction for enzyme activity [28, 49–58].
In the present study, we report a comparative structure analysis of AGPase SS from different monocot and dicot crop plants based on the available atomic resolution structure of S. tuberosum AGPase SS (PDB ID: 1YP2). A detailed structural comparison of both monocot and dicot AGPase SS along with their specificity towards substrate (ATP) and inhibitor (sulphate) binding has been elucidated. The mode of interactions of the SS of AGPases with sulphate inhibitor is studied with the aid of molecular docking. Detailed structural comparison of AGPase SS and the key amino acid residues involved in substrate and inhibitor binding from the selected crop species will highlight the important structural aspects of AGPase SS and may provide insights into the enzyme’s catalytic mechanism and understanding of the inhibitor binding specificity.
2. Materials and Methods
2.1. Computational Resources
All steps in this research were carried computationally on a Xeon, 2.13 GHz server equipped with the windows server 2003 environment. Preparation of three-dimensional structures, structure refinement, superimpositions, and docking were performed in Discovery Studio (DS3.5) (Accelrys, San Diego, CA, USA).
2.2. Sequence Analysis
Fasta formatted amino acid sequence of AGPase SS from three monocot crop plants, that is, Oryza sativa ssp. japonica, Hordeum vulgare, and Triticum aestivum, and six dicot crop plants, that is, Arabidopsis thaliana, Solanum lycopersicum, Beta vulgaris, Vicia faba, Cicer arietinum, and Brassica napus, was retrieved from the UniProtKB (http://www.uniprot.org/help/uniprotkb) database of ExPaSy. Primary structural study of the protein was done by computing various Physicochemical properties such as molecular weight, isoelectric point, instability index, aliphatic index, and grand average hydropathy (GRAVY) using ProtParam tool (http://web.expasy.org/protparam/) . The secondary structure of AGPase SS was predicted from its primary amino acid sequence using CONCORD (http://helios.princeton.edu/CONCORD)  secondary structure prediction server. This is an accurate secondary structure prediction method that incorporates seven popular secondary structure prediction methods, namely, PSIPRED, DSC, GOR, Predator, Prof, PROFphd, and SSpro, for predicting the consensus out of them. The disordered regions of AGPase SS were predicted by protein disordered metaprediction server (metaPrDOS) (http://prdos.hgc.jp/meta/) .
2.3. Domain Analysis and Linker Prediction
Protein domain boundaries and architecture knowledge is essential for understanding and characterising of protein function. Detection of protein domain and architecture in the absence of three-dimensional structure benefits many areas of protein science, such as protein engineering and protein structure prediction . Putative conserved domain, family, and superfamily possessed by AGPase SS were predicted based on sequence similarity search with its closest orthologous family members. Different tools and databases, namely, InterProScan (http://www.ebi.ac.uk/Tools/pfa/iprscan/) , Proteins Families Database (Pfam) (http://pfam.sanger.ac.uk/) , NCBI Conserved Domains Database (NCBI-CDD) (http://www.ncbi.nlm.nih.gov/Structure/cdd/cdd.shtml) , and SMART (http://smart.embl-heidelberg.de/)  server, were used for performing the task.
Overall functionality and efficiency of multiple domain proteins are affected by linker sequences. Cooperation and interaction between domains are affected by linker sequences which are flexible in 3D space, nonglobular, unstructured, or low complexity segment . The linker sequence joining the discrete domains of AGPase SS was inferred manually.
2.4. Multiple Sequence Alignment
ClustalW  was used to construct a multiple sequence alignment of monocot and dicot AGPase SS along with their closest structural homologue S. tuberosum AGPase SS to have a better knowledge on conservation and variation of different amino acids at sequence level. ESPript (Easy Sequencing in PostScript) (http://www.ipbs.fr/ESPript)  was used for rendering the result which facilitates the rapid visualisation of sequence alignment via PostScript output. It produces a synthesis of sequence and structural information by reading secondary structure files such as that created by the program DSSP .
2.5. Comparative Modeling
2.5.1. Template Identification, Model Building, and Refinement
Full length amino acid sequences of both monocot and dicot AGPase SS were subjected to BLASTP analysis against PDB in order to find suitable templates for comparative modeling. Blosum-62 matrix was used with a default threshold -value of 10 and inclusion threshold value of 0.005. Template was selected based on high level of sequence identity, query coverage, and alignment quality which promises a more reliable and good quality model. GeneSilico MetaServer (https://genesilico.pl/meta2)  which uses a consensus approach to predict the template for model building was also employed to have better confidence and to ensure the sensitivity and accuracy of template selection for AGPase SS.
Theoretical three-dimensional modeling of AGPase SS was built using create homology model module of DS3.5. Initially twenty different models for each protein were built and were ranked according to their normalized discrete optimized protein energy (DOPE) scores and the model with the lowest DOPE score was selected for further validation. Modrefiner server (http://zhanglab.ccmb.med.umich.edu/ModRefiner/)  which is a high resolution protein structure refinement algorithm at atomic level was used to refine the modeled 3D structure of AGPase SS of both monocot and dicot closer to their native state. To increase the compatibility score of each residue, the target models were further refined by loop modeling and side chain refinement using DS3.5 ab-initio loop prediction algorithm Looper  and ChiRotor  for refining protein side-chain conformations. Looper generates a set of low energy conformations for the specified loop region and ChiRotor systematically search side-chain conformation and scores based on their CHARMM  energy. Out of five different models generated by Looper and ChiRotor, the best model was selected based on the lowest DOPE score. Refined models were subjected to energy minimisation by DS3.5 with the minimisation protocol. The minimisation protocol employs the steepest descent and conjugate gradient methods of minimisation algorithms used with a generalized born implicit solvent model. Parameters of a distance-dependent dielectric constant = 1, nonbonded radius of 14 Å, CHARMM force field, spherical electrostatic cutoff, and the steepest descent algorithm were used to remove close van der Waals contacts for maximum steps of 2000 with 0.1 minimising RMS gradient.
2.6. Structural Assessment and Refinement
Generated models were tested for quality by both geometric and energetic means. PROCHECK , ERRAT , and VERIFY3D  tools, which are embedded in structure analysis and validation server (SAVES) (http://nihserver.mbi.ucla.edu/SAVES/), were used for validation of the modeled proteins. The PROCHECK provides an idea of the stereo chemical quality of the protein. It analyses the Ramachandran plot quality, peptide bond planarity, nonbonded interactions, main chain hydrogen bond energy, Cα chiralities, and overall factor. ERRAT checks the overall quality factor of the protein and was used to check the statistics of nonbonded interactions between different atom types. VERIFY3D was used to access the compatibility of the atomic models with its own amino acid sequence. A high VERIFY3D profile score indicates the high quality of protein model. MetaMQAPII (https://genesilico.pl/toolkit/unimod?method=MetaMQAPII)  which uses the result of VERIFY3D, PROSA, BALA, ANOLEA, PROVE, TUNE, REFINER, and PROQRES was used to identify global and/or local accuracy in models. Furthermore, secondary structure, solvent accessibility, and depth within the structure are also analysed by MetaMQAPII and are being used to assess the deviations of C-α atoms for a given model together with linear regression. Then, the global accuracy of the model is calculated based on these predictions and expressed as Global Distance Test_Total Score (GDT_TS). This indicator uses a set of distance thresholds of 1 Å, 2 Å, 4 Å, and 8 Å to find the average of the percentages of matched residue pairs between the C-α atoms in the model and the experimentally determined structure.
After each loop and side chain refinement steps, the above model quality assessment programs were employed to check the error at each residue in the protein. This process was repeated iteratively until the most geometrically and energetically stable structural conformation was attained.
To investigate how well the modeled structure matches the X-ray data of template protein, RMSD (root mean square deviation) between equivalent Cα and backbone atom pairs (target and template) was calculated by structural superimposition using superimpose module of DS3.5. STRIDE (http://webclu.bio.wzw.tum.de/stride/)  server was used to distinguish secondary structural elements of the predicted three-dimensional models from their atomic coordinates.
2.6.1. Calculation of Noncovalent Interactions
Noncovalent interactions are weak electromagnetic interactions between atoms or molecules and help in understanding many chemical and biological phenomena. These interactions are critical in maintaining the three-dimensional structure of protein and hold the key to understand the molecular basis of stability and functions of protein. Various noncovalent interactions such as hydrogen bonds, ionic interactions, hydrophobic interactions, disulphide bonds, aromatic-aromatic interactions, aromatic-sulphur interactions, and Cation-pi interactions of both monocot and dicot AGPase SS were calculated and compared with the template protein using protein interactions calculator (PIC) server . PIC server computes various interactions as mentioned above within a protein upon submission of the three- dimensional coordinate set of the protein.
Interaction between enzyme and its substrate provides an accurate picture of the interacting amino acid residues between the substrate and the active site. Different binding site prediction methods were employed for finding the binding site amino acid residues of both monocot and dicot AGPase SS. Binding sites of these models were selected based on the ligand-binding pocket of the template. DS3.5 binding site prediction module was also employed to predict the binding sites amino acid residues and functional residues which identifies the binding sites based on eraser and flood-filling algorithm . MetaPocket2.0 server (http://projects.biotec.tu-dresden.de/metapocket/)  was also used for predicting binding sites of these models. It employs a consensus method by combining the results of four different methods (LIGSITEcs, PASS, Q-SiteFinder, and SURFNET) to improve the prediction success rate. The potential ligand binding sites were generated using a probe radius of 5.0 Å and the binding site having highest -score was considered for further investigation.
Previous study by Jin et al., 2005,  reported that sulphate molecule acts as an inhibitor of S. tuberosum AGPase and later Boehlein et al., 2010,  explored the sulphate ion binding sites within the SS homotetramer of S. tuberosum to probe the allosteric binding sites of the maize endosperm AGPase. In addition, sulphate shows its structural similarity to all known allosteric regulators of higher plant AGPases which contain one or more phosphate moieties. For these reasons, we are using sulphate as an inhibitor to probe the allosteric binding sites of both monocot and dicot AGPase SS.
For ligand-protein interaction, both protein and ligand molecules were optimized using the “prepare protein and ligands tool” of DS3.5 which adds hydrogen ions to the protein and adds charges and hydrogen and applies force field to the ligand based on the CHARMM force field.
In this study, CDOCKER  module of DS3.5 was used to carry out the docking analysis. It is a grid based docking which uses CHARMM molecular simulation program to dock ligands within the active site of receptors. Prepared ligand molecule was docked into the active site of both monocot and dicot AGPase SS to elucidate its binding affinity towards the inhibitor molecule which in turn provides an insight into the allosteric regulation of the protein. The binding affinity of the ligand molecule into the active site of the protein was calculated based on the consensus scoring scheme of CDOCKER ENERGY, CDOCKER_Interaction Energy, Ligscore1_Dreiding, LigScore2_Dreiding, PLP1, PLP2, Jain, PMF, and PMF4 implemented in the protein-ligand interaction module of DS3.5.
3. Results and Discussions
3.1. Sequence Analysis
Single gene encoding the AGPase SS from monocot (Oryza sativa ssp. japonica, Hordeum vulgare, and Triticum aestivum) and dicot (Arabidopsis thaliana, Solanum lycopersicum, Beta vulgaris, and Brassica napus) crop plants was selected for the present study. In addition, two different SS encoding genes of dicot crop plant Vicia faba (AGPP, AGPC) and Cicer arietinum (CagpS1, CagpS2) were also considered for the present investigation. AGPP and AGPC genes of Vicia faba will be termed as Vicia faba1 and Vicia faba2 and CagpS1 and CagpS2 genes of Cicer arietinum will be termed as Cicer arietinum1 and Cicer arietinum2, respectively, for simplicity. The UniProtKB accession number, organism name, gene name, number of amino acids, subunit structure, and subcellular location of AGPase SS belonging to monocot and dicot crop plants are shown in Table 1.
Primary sequence analysis of AGPase SS signifies that molecular weight of the selected sequences falls between 52 and 57 KDa. Theoretical isoelectric point, low GRAVY indices, and high aliphatic index suggest that both monocot and dicot AGPase SS under consideration are slightly acidic in nature, possess high affinity towards water, and are thermostable at wide range of temperatures. Protein with instability index less than 40 is predicted as stable or else unstable . Computed instability index suggests that except Oryza sativa (P15280), Solanum lycopersicum (Q42882), Vicia faba2 (P52417), and Cicer arietinum1 (Q9AT06), all other AGPase SS sequences under consideration are stable proteins. Computed physicochemical properties of AGPase SS from monocot and dicot crop plants are reported in Table 2.
Comparison of the predicted secondary structure statistics of both monocot and dicot AGPase SS along with their structural homologue using CONCORD reveals that random coils dominated among secondary structure elements followed by strands and helices which are shown in Table 3.
Proteins often in their native states have regions with very flexible and unstable structures treated as disordered regions which are involved in many biological processes such as regulation, signaling, and cell cycle control [87, 88]. During the interaction with ligands, it is often observed that disordered regions transit to order where the flexibility of the region provides high specificity and low affinity towards multiple partners . Results from metaPrDOS reveal that the N-terminal amyloplast target sequence of approximately 80 residues and approximately five residues at the end of C-terminal region of AGPase SS is falling in the disordered region. It is also observed that the disordered regions of AGPase SS are prone to charged residues, low sequence complexity regions, and residues involved in phosphorylation.
3.2. Domain Analysis and Linker Prediction
Results from different domain prediction tools confirm a consensus prediction. Both monocot and dicot SS of AGPase are composed of an N-terminal catalytic ADP-glucose pyrophosphorylase domain and a C-terminal left-handed parallel beta helix domain. The N-terminal catalytic domain is approximately 200 residues and is structurally similar to Rossmann fold, typically present in nucleotide-binding domains. The C-terminal domain is composed of approximately 105 residues and is involved in cooperative allosteric regulation and oligomerisation. Domain positions of both the domains in the sequence are shown in Table 3. Smart analysis reflects that none of the SS of AGPase under investigation possess any transmembrane domain and they appeared to be soluble proteins.
Discrete domains are often associated with multiple functions of protein where domains are connected by inter-domain linkers. They keep the domains apart and provide great extent of flexibility to move individually. This phenomenon is a part of their catalytic function. The linker regions of AGPase SS were manually delineated and the amino acids propensities and order in linkers were examined. The multiple sequence alignment of the linker region (data not shown) shows a high percentage of sequence identity between monocot and dicot AGPase SS. Polar charged and uncharged hydrophilic residues and nonpolar hydrophobic residues are in equal propensities in the linker region. The number of proline residues is not adequate for rigidity of the linker which keeps apart the discrete domains present in AGPase SS.
3.3. Multiple Sequence Alignment
The alignment shows a high percentage (more than 90%) of sequence identity between monocot and dicot crop plants under investigation. Secondary structure information of S. tuberosum AGPase SS was used for analysing the conservation of secondary structure elements which in turn provides the functional characteristics. Comparison of secondary structure elements (α helices and β strands) with S. tuberosum AGPase SS shows a high percentage of sequence conservation throughout the alignment which suggests a similar functionality of this protein in both monocot and dicot. Analysis of both the domain architecture and boundaries in the multiple alignments shows a high percentage of sequence conservation among all the AGPase SS along with their structural homologue. This reflects the conservation of domains throughout the evolutionary period and suggests their conserved role in substrate and inhibitor binding. The consensus sequence along with the template and the secondary structure elements is represented in Figure 2.
3.4. Comparative Modeling
BLAST search results reveal that high resolution (2.11 Ǻ) crystal structure of S. tuberosum AGPase SS (PDB ID: 1YP2_A) is the suitable template for model building. It has a high sequence identity, query coverage, less -value, and a high level of alignment quality with query sequences. Table 4 shows the query coverage, -value, and sequence identity of the query sequences against their template protein S. tuberosum AGPase SS.
The crystal structure of S. tuberosum AGPase SS consists of 442 amino acids, lacking the amyloplast target sequence. Prior to modeling, the amyloplast target sequence of both monocot and dicot AGPase SS at the N-terminal end was removed to exclude the random coil fragment at this region and to achieve a better global superimposition with the template. Based on single-template approach, 20 different models were generated and the model having the lowest DOPE score is selected for model refinement and validation.
3.5. Model Quality Assessment
Different model quality assessment programs were employed for validating the refined and optimised models. The PROCHECK analysis shows that a high percentage of residue’s and angles are in the favoured region of Ramachandran plot contributing to the correctness of the models which is evident from Table 5 (see Supplementary Figure S1 available online at http://dx.doi.org/10.1155/2014/583606). Overall factor is also in acceptable range as shown in Table 5 indicating that the designed models are of good quality and acceptable. An acceptable value of -factor in Procheck is between 0 and −0.5 with the best model displaying values close to zero . VERIFY3D program assessed the packing quality of each residue of the model where the compatibility of the model residues with their environment is assessed by a score function. Residues with a score over 0.2 (cut-off score > 0) should be considered reliable. It is evident from Table 5 that the score of both monocot and dicot models maximally lies above 0.2 which corresponds to acceptable side chain environment. A very high ERRAT score (Table 5) contributes to the acceptance of the models. Global distance total test score (GDT_TS)/RMSD score predicted by MetaMQAPII server is well above the cut-off score suggesting the acceptability of the models (Table 5). An ideal model has GDT_TS score over 59 and a RMSD around 2.0 Å. Model quality assessment score predicted by different programs is shown in Table 5.
The above model quality assessment programs check the stereo chemical quality, nonbonded interactions of the residues, the compatibility of the side chain environment, packing quality, and the energy profile of the predicted AGPase SS models and their result signify the high quality, reliability, and acceptability of the proposed models.
Computed RMSD of the Cα and backbone atom pairs for all the models are very low (Table 5). Low RMSD values indicate that the generated models are reasonably good and share high structural similarity and a common folding pattern with the template. Atomic coordinates of the models are deposited in Protein Model Database (PMDB) and can be accessed at http://mi.caspur.it/PMDB using PMDB ID: PM0079235-244 and PM0079249.
3.5.1. Detailed Structural Comparison of AGPase Small Subunit
Theoretical three-dimensional models of AGPases SS belonging to different monocot and dicot species were analysed extensively to have a wide spectrum on the three-dimensional structure and the role of key residues responsible for catalytic and inhibitory function. All the predicted structures of AGPase SS (Supplementary Figure S2) are composed of 15 helices (~23%), 10 310 helices (2.3%), and 29-30 strands (~27.1%). Both monocot and dicot AGPase SS possess a loop and two gamma turns. The loop has a strand of 11 amino acid residues positioning Tyr217-Ser227 and a strand of 5 residues positioning Val260-Leu264 interconnected by a 32-residue long loop. The gamma loop forming residues of both monocot and dicot subunits are Gly40-Ala41-Asn42 and Val200-Asp201-Thr202. Structural comparison of modeled proteins shows strong conservation of secondary structural elements among them. Hydrogen bonding is most prominent among the beta turn forming residues thereby stabilising the connecting link of helices and strands. The comparative structural analysis of all the AGPase SS in the present investigation shares the same architecture. Structural superimposition of all the AGPase SS with its template shows that the secondary structure elements are superposed well and the key residues of allosteric regulation, that is, Arg32 equivalent to Arg41 in the template, Arg44/Arg53, Lys60/Lys69, His75/His84, His125/His134, Gln305/Gln314, Arg307/Arg316, Lys395/Lys404, and Lys432/Lys441, are fully conserved and are allosterically significant (Figure 3).
Tyrosine 135 which is a key residue within the active site of S. tuberosum AGPase SS for allosteric regulation was substituted to Asn126 in both monocot and dicot AGPase SS. Previous study on this enzyme reports that GXGXRL loop, PAVP motif, and residue equivalent to Arg33 in S. tuberosum SS play a key role for ATP binding which has been demonstrated by mutagenesis study [21, 91]. Structural superimposition shows that the strong conservation of GXGXRL loop positioning 20–25, PAVP motif positioning 35–38, and conservation of Arg24 (equivalent to Arg33 of S. tuberosum AGPase SS) in both monocot and dicot SS of AGPases firmly reflects the similar mode of action in this family of enzymes. Previous study by Jin et al. (2005)  concluded that metal-mediated catalytic mechanism is also used by AGPase. Residues equivalent to Asp145 and Asp280 (in S. tuberosum AGPase SS) chelate the metal ion and play a crucial role in metal-mediated catalytic mechanism. In several organisms the absolute need of a metal ion for AGPase has been biochemically demonstrated [6, 17, 92]. Taken together we have tried to check the binding specificity of these residues in our models to have a better confidence about the importance and involvement of these residues in metal-mediated catalytic mechanism across the family. Structural superimposition of our modeled proteins with the template protein reveals that Asp136 and Asp271 of both monocot and dicot equivalent to Asp145 and Asp280 of S. tuberosum AGPase SS are strongly conserved and follow common folding pattern with a very low RMSD difference which may prompt for metal-mediated catalytic mechanism. Comparison of secondary structure elements of AGPase SS from its three-dimensional coordinate using STRIDE (uses the three-dimensional structure for prediction of their secondary structure elements) (Table 6) reveals similar statistics to that of secondary structures from the primary amino acid sequence predicted by CONCORD server. This signifies the accuracy and reliability of the secondary structure elements of AGPase SS, assigned through homology modeling.
3.5.2. Analysis on Noncovalent Interactions of Both Monocot and Dicot AGPase Small Subunit along with Their Structural Homologue
Collective effort of various noncovalent interactions determines the overall structure and behavior of proteins. An extensive hydrogen bond network (Main Chain-Main Chain/Main Chain-Side Chain/Side Chain-Side Chain) is observed in both monocot and dicot AGPase SS, which is parallel to the hydrogen bond network formed by their structural homologue. As protein folding and function is significantly contributed by hydrogen bonds, it is conceivable that both monocot and dicot AGPase SS share a common structural and functional similarity with that of the template protein. When positively charged amino acids and aromatic amino acids are in close proximity, they form the Cation-pi interactions which are significant to the protein structure. Arginine (R) is dominant over lysine (K) and the order of likeliness to participate in cation-pi interactions of aromatic amino acids is Tryptophan (W) followed by Tyrosine (Y) followed by Phenylalanine (F). In both monocot and dicot AGPase SS equal numbers of energetically significant cation-pi interactions were observed which is parallel to the cation-pi interactions of the template protein.
3.5.3. Docking Analysis
The enzyme-substrate complex provides greater insight into the interactions and structural complementarities between the substrate and the active site of a protein. Assuming that the inhibitor binding modes are similar in the target and the template protein structure, active sites of the molecules were delineated based on the ligand-binding pocket of the template protein. The binding site module of DS3.5 and MetaPocket2.0 server also predicted the consensus binding site cavities and carries the similar conserved residues. Superposition (overlapping of Cα atoms) of modeled AGPase SS with their structural homologue in DS3.5 reveals that the spatial position and orientation of substrate and inhibitor binding site are highly conserved (Figures 4(a) and 4(b) and Figures 5(a) and 5(b)) and shows a very low RMSD within the substrate and inhibitor binding site which is evident from Table 7.
Three sulphate molecules were successively docked into the active sites of both monocot and dicot AGPase SS to elucidate their structural and functional relevance in terms of inhibitor binding. Docking of inhibitor into the active site of modeled subunit reveals that Arg32, Arg44, Lys395, and Lys432 are directly involved in the interaction with the first sulphate molecule of both monocot and dicot AGPase SS by strong hydrogen and hydrophobic bonding, except for Beta vulgaris where Lys395 is missing. Tyr135 which is an aromatic, partially hydrophobic, amino acid was substituted to polar Asn126 in both monocot and dicot and plays a major role in the interaction with the second sulphate molecule along with Lys60 and His125 via hydrogen bonding. The third sulphate molecule binds with Arg44, Gln305, and Arg307 in both monocot and dicot AGPase SS, whereas H75 also shows its binding efficiency for sulphate in Cicer arietinum1, Beta vulgaris, Solanum lycopersicum, and Vicia faba2. Comparative docking study of sulphate inhibitor into the binding cavity of both monocot and dicot AGPase SS reflects a similar mode of binding specificity. Strong conservation and involvement of similar residues in the interaction convincingly suggest a similar inhibitor binding mechanism in both monocot and dicot AGPase SS in relation to the template protein. Stick representation of the docked complex of both monocot and dicot AGPase SS are shown in Figure 6.
3.5.4. Comparison of Electrostatic Surface Potential of Binding Cavity from Both Monocot and Dicot AGPase Small Subunit
The electrostatic surface potential energy of the modeled protein was calculated using the Adaptive Poisson-Boltzmann solver (APBS) package via PDB2PQR web portal (http://kryptonite.nbcr.net/pdb2pqr/) . This state-of-the-art suite offers the computation of Poisson-Boltzmann electrostatic calculations on biomolecules. Energy refinement force field (AMBER) is used to assign atomic charges and ionic radii to the three- dimensional coordinate protein file. At pH 7.0, protonation state to the proteins is assigned by PROPKA. Solvent dielectric constant of 78.54 and a low dielectric constant of 2.0 were assigned to the protein. Most positively charged and most negatively charged surfaces were coloured using −5 and +5 settings where is Boltzman Constant, is temperature, and is the charge of electron. Charge distribution and patches within the binding cavity of both monocot and dicot crop plants were studied extensively. Total electrostatic energies of both monocot and dicot AGPase SS range between kJ/mol and kJ/mol which is parallel to the total electrostatic energy of the template protein ( kJ/mol). The isocontour representation data fascinatingly displayed a similar topology with a similar electrostatic potential distribution pattern within the binding site of both monocot and dicot AGPase SS along with their template protein (Figure 7). This similar distribution of the electrostatic potential energy of both monocot and dicot AGPase SS substantiates similar structural organisation and perfectly correlates with their similar mode of binding activity.
The AGPase enzyme plays a pivotal role in starch biosynthesis for both photosynthetic and nonphotosynthetic plant tissues. Catalytic activity and allosteric regulation of this enzyme have significantly contributed to the overall yield potential of many crop plants. In the present study the theoretical three-dimensional models of AGPase SS from three different monocot and six different dicot crop plants were constructed for understanding the structure function relationship, substrate, and inhibitor binding specificity. The models were validated and further used for docking analysis with sulphate inhibitor. Superimposition of models with the crystal structure of the SS of S. tuberosum AGPase (PDB ID: 1YP2) shows a relatively low RMSD difference indicating a high structural similarity among the subunits. Structural superimposition of both monocot and dicot SS along with the template protein reflects strong conservation of secondary structure elements, common folding patterns, and similarity in the domains. GXGXRL loop and PAVP motif positioning 20–25 and 35–38, respectively, and Arg24 (equivalent to Arg33 of S. tuberosum AGPase SS) which play a crucial role for ATP binding are strongly conserved. It also follows a common folding pattern with a very low RMSD value in both monocot and dicot SS of AGPases reflecting similar mode of ATP binding in this family of enzymes. Structural superimposition of our modeled proteins with the template protein reveals that the aspartic acid positioning 136 and 271 (equivalent to Asp145 and Asp280 of S. tuberosum) of both monocot and dicot is strongly conserved and may act for metal mediated catalytic mechanism. Previous Alanine-scanning mutagenesis study performed by Boehlein et al., 2010,  suggests that mutating any of the Arg residues that are equivalent to S. tuberosum sulfate binding Arg residues drastically alters the maize AGPase overall allosteric properties. It clearly suggests the involvement of these Arg residues in allosteric effectors binding. Our docking study reveals the participation of similar Args (Arg32, Arg44, Arg307, and Arg361) residues in sulfate binding. Additionally Lys60, His75, His125, Asp126, Gln305, Lys395, and Lys432 are the key residues responsible for inhibitor binding and are allosterically significant in both monocot and dicot AGPase SS. The charge distribution patterns displayed a similar topology with a similar electrostatic potential distribution pattern within the binding site of both monocot and dicot AGPase SS and signify their similar mode of binding activity.
It is one of the major concerns about the analysis of the monocot AGPases because it is known that monocot plants like barley AGPases are insensitive to 3PGA and Pİ regulation. In addition to the barley endosperm AGPase, pea embryos and wheat grains and AGPases from nonphotosynthetic tissues have also been reported to be insensitive or weakly affected by 3PGA and Pİ regulation [94–96]. However, our findings are parallel with the result of Doan et al., 1999,  who suggested that barley endosperm AGPase SS is sensitive to 3-phosphoglycerate activation and is inhibited by inorganic phosphate although the active heterotetramer is allosterically unregulated. Therefore it can by hypothesised that SS of AGPases from both monocot and dicot under present investigation has similar substrate and inhibitor binding mechanism in this family of enzyme and the LS may be playing a crucial role for showing different variant of allosterically regulated and unregulated AGPases.
The active AGPase enzyme is a heterotetramer of two SS and two LS. Several studies have indicated that maximum stability, catalytic, and regulatory properties of the enzyme from higher plants result from the synergy of both SS and LS. Thus for understanding the actual mechanism of the active heterotetrameric AGPase and the dynamics involved in SS and LS assembly, their interaction and stability should be thoroughly studied. As the mode of action and dynamics of this protein depend on the delicate equilibrium between the protein and its native environment, further structural characterisation is bound to shed more light on the detailed mode of action of this family of protein in different crop species.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
The authors thankfully acknowledge the financial support for Agri-Bioinformatics Promotion Program by Bioinformatics Initiative Division, Department of Information Technology, Ministry of Communications & Information Technology, Government of India, New Delhi, and Biotechnology Information System Network, Department of Biotechnology, Government of India. The authors are grateful to Assam Agricultural University, Jorhat, Assam, for providing the necessary facilities and constant support to carryout research endeavor. Also, authors are thankful to Assam University, Silchar, for support and guidance.
Supplementary Figure S1. Ramachandran plot of both monocot and dicot AGPase small subunit. (i) Oryza sativa ssp. japonica AGPase small subunit, (ii) Hordeum vulgare AGPase small subunit, (iii) Triticum aestivum AGPase small subunit, (iv) Arabidopsis thaliana AGPase small subunit, (v) Solanum lycopersicum AGPase small subunit, (vi) Beta vulgaris AGPase small subunit, (vii) Vicia faba1 AGPase small subunit, (viii) Vicia faba2 AGPase small subunit, (ix) Cicer arietinum1 AGPase small subunit, (x) Cicer arietinum2 AGPase small subunit, (xi) Brassica napus AGPase small subunit .
Supplementary Figure S2. Theoritical three dimensional models of both monocot and dicot AGPase small subunit. (i) Oryza sativa ssp. japonica AGPase small subunit, (ii) Hordeum vulgare AGPase small subunit, (iii) Triticum aestivum AGPase small subunit, (iv) Arabidopsis thaliana AGPase small subunit, (v) Solanum lycopersicum AGPase small subunit, (vi) Beta vulgaris AGPase small subunit, (vii) Vicia faba1 AGPase small subunit, (viii) Vicia faba2 AGPase small subunit, (ix) Cicer arietinum1 AGPase small subunit, (x) Cicer arietinum2 AGPase small subunit, (xi) Brassica napus AGPase small subunit
P. R. Salamone, T. W. Greene, I. H. Kavakli, and T. W. Okita, “Isolation and characterization of a higher plant ADP-glucose pyrophosphorylase small subunit homotetramer,” FEBS Letters, vol. 482, no. 1-2, pp. 113–118, 2000.View at: Publisher Site | Google Scholar
C. Martin and A. M. Smith, “Starch biosynthesis,” Plant Cell, vol. 7, no. 7, pp. 971–985, 1995.View at: Publisher Site | Google Scholar
J. Preiss, “Biology and molecular biology of starch synthesis and its regulation,” Plant Molecular and Cellular Biology, vol. 7, pp. 59–114, 1991.View at: Google Scholar
J. Preiss, “Bacterial glycogen synthesis and its regulation,” Annual Review of Microbiology, vol. 38, pp. 419–458, 1984.View at: Publisher Site | Google Scholar
C. J. Slattery, I. H. Kavakli, and T. W. Okita, “Engineering starch for increased quantity and quality,” Trends in Plant Science, vol. 5, no. 7, pp. 291–298, 2000.View at: Publisher Site | Google Scholar
M. A. Ballicora, A. A. Iglesias, and J. Preiss, “ADP-glucose pyrophosphorylase: a regulatory enzyme for plant starch synthesis,” Photosynthesis Research, vol. 79, no. 1, pp. 1–24, 2004.View at: Publisher Site | Google Scholar
L. A. Kleczkowski, “Is leaf ADP-glucose pyrophosphorylase an allosteric enzyme?” Biochimica et Biophysica Acta, vol. 1476, no. 1, pp. 103–108, 2000.View at: Publisher Site | Google Scholar
A. Tiessen, J. H. M. Hendriks, M. Stitt et al., “Starch synthesis in potato tubers is regulated by post-translational redox modification of ADP-glucose pyrophosphorylase: a novel regulatory mechanism linking starch synthesis to the sucrose supply,” The Plant Cell, vol. 14, no. 9, pp. 2191–2213, 2002.View at: Publisher Site | Google Scholar
C. Y. Tsai and O. E. Nelson, “Starch-deficient maize mutant lacking adenosine diphosphate glucose pyrophosphorylase activity,” Science, vol. 151, no. 3708, pp. 341–343, 1966.View at: Publisher Site | Google Scholar
L. C. Hannah and O. E. Nelson Jr., “Characterization of ADP glucose pyrophosphorylase from shrunken 2 and brittle 2 mutants of maize,” Biochemical Genetics, vol. 14, no. 7-8, pp. 547–560, 1976.View at: Publisher Site | Google Scholar
T. P. Lin, T. Casper, C. R. Sommerville, and J. Preiss, “A starch-deficient mutant of Arabidopsis thaliana with low ADP-glucose pyrophosphorylase activity lacks one of the two subunits of the enzyme,” Plant Physiology, vol. 88, pp. 1175–1181, 1988.View at: Google Scholar
B. Muller-Rober, U. Sonnewald, and L. Willmitzer, “Inhibition of the ADP-glucose pyrophosphorylase in transgenic potatoes leads to sugar-storing tubers and influences tuber formation and expression of tuber storage protein genes,” The EMBO Journal, vol. 11, no. 4, pp. 1229–1238, 1992.View at: Google Scholar
D. M. Stark, K. P. Timmerman, G. I. Barry, J. Preiss, and G. M. Kishore, “Regulation of the amount of starch in plant tissues by ADP glucose pyrophosphorylase,” Science, vol. 258, no. 5080, pp. 287–292, 1992.View at: Publisher Site | Google Scholar
A. A. Iglesias, G. F. Barry, C. Meyer et al., “Expression of the potato tuber ADP-glucose pyrophosphorylase in Escherichia coli,” The Journal of Biological Chemistry, vol. 268, no. 2, pp. 1081–1086, 1993.View at: Google Scholar
T. H. Haugen, A. Ishaque, and J. Preiss, “Biosynthesis of bacterial glycogen. Characterization of the subunit structure of Escherichia coli B glucose 1 phosphate adenyltransferase (EC 188.8.131.52),” The Journal of Biological Chemistry, vol. 251, no. 24, pp. 7880–7885, 1976.View at: Google Scholar
T. W. Okita, P. A. Nakata, J. M. Anderson, J. Sowokinos, M. Morell, and J. Preiss, “The subunit structure of potato tuber ADPglucose pyrophosphorylase,” Plant Physiology, vol. 93, no. 2, pp. 785–790, 1990.View at: Publisher Site | Google Scholar
M. A. Ballicora, A. A. Iglesias, and J. Preiss, “ADP-glucose pyrophosphorylase, a regulatory enzyme for bacterial glycogen synthesis,” Microbiology and Molecular Biology Reviews, vol. 67, no. 2, pp. 213–225, 2003.View at: Publisher Site | Google Scholar
B. J. Smith-White and J. Preiss, “Comparison of proteins of ADP-glucose pyrophosphorylase from diverse sources,” Journal of Molecular Evolution, vol. 34, no. 5, pp. 449–464, 1992.View at: Publisher Site | Google Scholar
K. Ball and J. Preiss, “Allosteric sites of the large subunit of the spinach leaf ADPglucose pyrophosphorylase,” Journal of Biological Chemistry, vol. 269, no. 40, pp. 24706–24711, 1994.View at: Google Scholar
M. A. Ballicora, M. J. Laughlin, Y. Fu, T. W. Okita, G. F. Barry, and J. Preiss, “Adenosine 5′-diphosphate-glucose pyrophosphorylase from potato tuber: significance of the N terminus of the small subunit for catalytic properties and heat stability,” Plant Physiology, vol. 109, no. 1, pp. 245–251, 1995.View at: Publisher Site | Google Scholar
T. W. Greene, S. E. Chantler, M. L. Kahn, G. F. Barry, J. Preiss, and T. W. Okita, “Mutagenesis of the potato ADPglucose pyrophosphorylase and characterization of an allosteric mutant defective in 3-phosphoglycerate activation,” Proceedings of the National Academy of Sciences of the United States of America, vol. 93, no. 4, pp. 1509–1513, 1996.View at: Publisher Site | Google Scholar
M. A. Ballicora, Y. Fu, N. M. Nesbitt, and J. Preiss, “ADP-glucose pyrophosphorylase from potato tubers. Site-directed mutagenesis studies of the regulatory sites,” Plant Physiology, vol. 118, no. 1, pp. 265–274, 1998.View at: Publisher Site | Google Scholar
M. J. Laughlin, J. W. Payne, and T. W. Okita, “Substrate binding mutants of the higher plant ADP-glucose pyrophosphorylase,” Phytochemistry, vol. 47, no. 4, pp. 621–629, 1998.View at: Publisher Site | Google Scholar
I. H. Kavakli, J. S. Park, C. J. Slattery, P. R. Salamone, J. Frohlick, and T. W. Okita, “Analysis of allosteric effector binding sites of potato ADP-glucose pyrophosphorylase through reverse genetics,” The Journal of Biological Chemistry, vol. 276, no. 44, pp. 40834–40840, 2001.View at: Publisher Site | Google Scholar
P. R. Salamone, I. H. Kavakli, C. J. Slattery, and T. W. Okita, “Directed molecular evolution of ADP-glucose pyrophosphorylase,” Proceedings of the National Academy of Sciences of the United States of America, vol. 99, no. 2, pp. 1070–1075, 2002.View at: Publisher Site | Google Scholar
I. H. Kavakli, T. W. Greene, P. R. Salamone, S. Choi, and T. W. Okita, “Investigation of subunit function in ADP-glucose pyrophosphorylase,” Biochemical and Biophysical Research Communications, vol. 281, no. 3, pp. 783–787, 2001.View at: Publisher Site | Google Scholar
S. K. Hwang, S. Hamada, and T. W. Okita, “ATP binding site in the plant ADP-glucose pyrophosphorylase large subunit,” FEBS Letters, vol. 580, no. 28-29, pp. 6741–6748, 2006.View at: Publisher Site | Google Scholar
S. K. Hwang, Y. Nagai, D. Kim, and T. W. Okita, “Direct appraisal of the potato tuber ADP-glucose pyrophosphorylase large subunit in enzyme function by study of a novel mutant form,” The Journal of Biological Chemistry, vol. 283, no. 11, pp. 6640–6647, 2008.View at: Publisher Site | Google Scholar
J. M. Cross, M. Clancy, J. R. Shaw et al., “A polymorphic motif in the small subunit of ADP-glucose pyrophosphorylase modulates interactions between the small and large subunits,” The Plant Journal, vol. 41, no. 4, pp. 501–511, 2005.View at: Publisher Site | Google Scholar
X. S. Jin, M. A. Ballicora, J. Preies, and J. H. Geiger, “Crystal structure of potato tuber ADP-glucose pyrophosphorylase,” The EMBO Journal, vol. 24, no. 4, pp. 694–704, 2005.View at: Publisher Site | Google Scholar
F. Melo and E. Feytmans, “Assessing protein structures with a non-local atomic interaction energy,” Journal of Molecular Biology, vol. 277, no. 5, pp. 1141–1152, 1998.View at: Publisher Site | Google Scholar
R. G. Bodade, S. D. Beedkar, A. V. Manwar, and C. N. Khobragade, “Homology modeling and docking study of xanthine oxidase of Arthrobacter sp. XL26,” International Journal of Biological Macromolecules, vol. 47, no. 2, pp. 298–303, 2010.View at: Publisher Site | Google Scholar
S. Witz, P. Panwar, M. Schober et al., “Structure-function relationship of a plant NCS1 member—homology modeling and mutagenesis identified residues critical for substrate specificity of PLUTO, a nucleobase transporter from arabidopsis,” PLoS ONE, vol. 9, Article ID e91343, 2014.View at: Google Scholar
B. Kamaraj and R. Purohit, “In-silico analysis of Betaine Aldehyde Dehydrogenase2 of Oryza sativa and significant mutations responsible for fragrance,” Journal of Plant Interactions, vol. 8, pp. 321–333, 2013.View at: Publisher Site | Google Scholar
R. Xu, S. Zhang, J. Huang, and C. Zheng, “Genome-wide comparative in silico analysis of the RNA helicase gene family in Zea mays and Glycine max: a comparison with arabidopsis and Oryza sativa,” PLoS ONE, vol. 8, Article ID e78982, 2013.View at: Google Scholar
M. Lakshmanan, B. Mohanty, and D. Y. Lee, “Identifying essential genes/reactions of the rice photorespiration by in silico model-based analysis,” Rice, vol. 6, pp. 1–5, 2013.View at: Google Scholar
R. Kolodny and M. Kosloff, “From protein structure to function via computational tools and approaches,” Israel Journal of Chemistry, vol. 53, no. 3-4, pp. 147–156, 2013.View at: Publisher Site | Google Scholar
F. H. Wallrapp, J. J. Pan, G. Ramamoorthy et al., “Prediction of function for the polyprenyl transferase subgroup in the isoprenoid synthase superfamily,” Proceedings of the National Academy of Sciences of the United States of America, vol. 110, no. 13, pp. E1196–E1202, 2013.View at: Publisher Site | Google Scholar
G. L. Gac, C. Ka, R. Joubrel et al., “Structure-function analysis of the human ferroportin iron exporter (SLC40A1): effect of hemochromatosis type 4 disease mutations and identification of critical residues,” Human Mutation, vol. 34, pp. 1371–1380, 2013.View at: Google Scholar
N. Chaturvedi, V. K. Singh, and P. N. Pandey, “Computational identification and analysis of arsenate reductase protein in Cronobacter sakazakii ATCC BAA-894 suggests potential microorganism for reducing arsenate,” Journal of Structural and Functional Genomics, vol. 14, no. 2, pp. 37–45, 2013.View at: Publisher Site | Google Scholar
S. Kashyap, “In Silico modeling and functional interpretations of Cry1Ab15 toxin from Bacillus thuringiensis BtB-Hm-16,” BioMed Research International, vol. 2013, Article ID 471636, 10 pages, 2013.View at: Publisher Site | Google Scholar
K. Sarma, B. Dehury, J. Sahu et al., “A comparative proteomic approach to analyse structure, function and evolution of rice chitinases: a step towards increasing plant fungal resistance,” Journal of Molecular Modeling, vol. 18, no. 11, pp. 4761–4780, 2012.View at: Publisher Site | Google Scholar
D. K. Rai and E. Rieder, “Homology modeling and analysis of structure predictions of the bovine rhinitis B virus RNA dependent RNA polymerase (RdRp),” International Journal of Molecular Sciences, vol. 13, no. 7, pp. 8998–9013, 2012.View at: Publisher Site | Google Scholar
B. Ganguly and S. Prasad, “Homology modeling and functional annotation of bubaline pregnancy associated glycoprotein 2,” Journal of Animal Science and Biotechnology, vol. 3, article 13, 2012.View at: Publisher Site | Google Scholar
A. G. Thorsell, W. H. Lee, C. Persson et al., “Comparative structural analysis of lipid binding START domains,” PLoS ONE, vol. 6, no. 6, Article ID e19521, 2011.View at: Publisher Site | Google Scholar
J. R. Manning, M. A. Bailey, D. C. Soares, D. R. Dunbar, and J. J. Mullins, “In silico structure-function analysis of pathological variation in the HSD11B2 gene sequence,” Physiological Genomics, vol. 42, no. 3, pp. 319–330, 2010.View at: Publisher Site | Google Scholar
S. W. Chen and J. L. Pellequer, “Identification of functionally important residues in proteins using comparative models,” Current Medicinal Chemistry, vol. 11, no. 5, pp. 595–605, 2004.View at: Publisher Site | Google Scholar
J. A. Peterson and S. E. Graham, “A close family resemblance: the importance of structure in understanding cytochromes P450,” Structure, vol. 6, no. 9, pp. 1079–1085, 1998.View at: Publisher Site | Google Scholar
C. Dawar, S. Jain, and S. Kumar, “Insight into the 3D structure of ADP-glucose pyrophosphorylase from rice (Oryza sativa L.),” Journal of Molecular Modeling, vol. 19, no. 8, pp. 3351–3367, 2013.View at: Publisher Site | Google Scholar
W. Wang and J. Messing, “Analysis of ADP-glucose pyrophosphorylase expression during turion formation induced by abscisic acid in Spirodela polyrhiza (greater duckweed),” BMC Plant Biology, vol. 12, article 5, 2012.View at: Publisher Site | Google Scholar
M. Petreikov, M. Eisenstein, Y. Yeselson, J. Preiss, and A. A. Schaffer, “Characterization of the AGPase large subunit isoforms from tomato indicates that the recombinant L3 subunit is active as a monomer,” Biochemical Journal, vol. 428, no. 2, pp. 201–212, 2010.View at: Publisher Site | Google Scholar
I. Baris, A. Tuncel, N. Ozber, O. Keskin, and I. H. Kavakli, “Investigation of the interaction between the large and small subunits of potato ADP-glucose pyrophosphorylase,” PLoS Computational Biology, vol. 5, no. 10, Article ID e1000546, 2009.View at: Publisher Site | Google Scholar
N. Georgelis, J. R. Shaw, and L. C. Hannah, “Phylogenetic analysis of ADP-glucose pyrophosphorylase subunits reveals a role of subunit interfaces in the allosteric properties of the enzyme,” Plant Physiology, vol. 151, no. 1, pp. 67–77, 2009.View at: Publisher Site | Google Scholar
A. Tuncel, I. H. Kavakli, and O. Keskin, “Insights into subunit interactions in the heterotetrameric structure of potato ADP-glucose pyrophosphorylase,” Biophysical Journal, vol. 95, no. 8, pp. 3628–3639, 2008.View at: Publisher Site | Google Scholar
N. Georgelis and L. C. Hannah, “Isolation of a heat-stable maize endosperm ADP-glucose pyrophosphorylase variant,” Plant Science, vol. 175, no. 3, pp. 247–254, 2008.View at: Publisher Site | Google Scholar
D. Kim, S. K. Hwang, and T. W. Okita, “Subunit interactions specify the allosteric regulatory properties of the potato tuber ADP-glucose pyrophosphorylase,” Biochemical and Biophysical Research Communications, vol. 362, no. 2, pp. 301–306, 2007.View at: Publisher Site | Google Scholar
M. A. Ballicora, E. D. Erben, T. Yazaki et al., “Identification of regions critically affecting kinetics and allosteric regulation of the Escherichia coli ADP-glucose pyrophosphorylase by modeling and pentapeptide-scanning mutagenesis,” Journal of Bacteriology, vol. 189, no. 14, pp. 5325–5333, 2007.View at: Publisher Site | Google Scholar
C. M. Bejar, X. Jin, M. A. Ballicora, and J. Preiss, “Molecular architecture of the glucose 1-phosphate site in ADP-glucose pyrophosphorylases,” The Journal of Biological Chemistry, vol. 281, no. 52, pp. 40473–40484, 2006.View at: Publisher Site | Google Scholar
E. Gasteiger, C. Hoogl, A. Gattiker et al., “Protein identification and analysis tools on the ExPASy server,” in The Proteomics Protocols Handbook, J. M. Walker, Ed., pp. 571–607, Humana Press, New Jersey, NJ, USA, 2005.View at: Google Scholar
Y. Wei, J. Thompson, and C. A. Floudas, “CONCORD: a consensus method for protein secondary structure prediction via mixed integer linear optimization,” Proceedings of the Royal Society A, vol. 468, no. 2139, pp. 831–850, 2012.View at: Publisher Site | Google Scholar | MathSciNet
T. Ishida and K. Kinoshita, “Prediction of disordered regions in proteins based on the meta approach,” Bioinformatics, vol. 24, no. 11, pp. 1344–1348, 2008.View at: Publisher Site | Google Scholar
L. Kong and S. Ranganathan, “Delineation of modular proteins: domain boundary prediction from sequence information,” Briefings in bioinformatics, vol. 5, no. 2, pp. 179–192, 2004.View at: Publisher Site | Google Scholar
E. Quevillon, V. Silventoinen, S. Pillai et al., “InterProScan: protein domains identifier,” Nucleic Acids Research, vol. 33, no. 2, pp. W116–W120, 2005.View at: Publisher Site | Google Scholar
M. Punta, P. C. Coggill, R. Y. Eberhardt et al., “The Pfam protein families database,” Nucleic Acids Research, vol. 40, no. 1, pp. D290–D301, 2012.View at: Publisher Site | Google Scholar
A. Marchler-Bauer, C. Zheng, F. Chitsaz et al., “CDD: conserved domains and protein three-dimensional structure,” Nucleic Acids Research, vol. 41, no. 1, pp. D348–D352, 2013.View at: Publisher Site | Google Scholar
I. Letunic, T. Doerks, and P. Bork, “SMART 7: recent updates to the protein domain annotation resource,” Nucleic Acids Research, vol. 40, no. 1, pp. D302–D305, 2012.View at: Publisher Site | Google Scholar
J. C. Wootton, “Non-globular domains in protein sequences: automated segmentation using complexity measures,” Computers and Chemistry, vol. 18, no. 3, pp. 269–285, 1994.View at: Publisher Site | Google Scholar
J. D. Thompson, D. G. Higgins, and T. J. Gibson, “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, pp. 4673–4680, 1994.View at: Publisher Site | Google Scholar
P. Gouet, E. Courcelle, D. I. Stuart, and F. Métoz, “ESPript: analysis of multiple sequence alignments in PostScript,” Bioinformatics, vol. 15, no. 4, pp. 305–308, 1999.View at: Publisher Site | Google Scholar
W. Kabsch and C. Sander, “Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features.,” Biopolymers, vol. 22, no. 12, pp. 2577–2637, 1983.View at: Publisher Site | Google Scholar
M. A. Kurowski and J. M. Bujnicki, “GeneSilico protein structure prediction meta-server,” Nucleic Acids Research, vol. 31, no. 13, pp. 3305–3307, 2003.View at: Publisher Site | Google Scholar
D. Xu and Y. Zhang, “Improving the physical realism and structural accuracy of protein models by a two-step atomic-level energy minimization,” Biophysical Journal, vol. 101, no. 10, pp. 2525–2534, 2011.View at: Publisher Site | Google Scholar
V. Z. Spassov, P. K. Flook, and L. Yan, “LOOPER: a molecular mechanics-based algorithm for protein loop prediction,” Protein Engineering, Design and Selection, vol. 21, no. 2, pp. 91–100, 2008.View at: Publisher Site | Google Scholar
V. Z. Spassov, L. Yan, and P. K. Flook, “The dominant role of side-chain backbone interactions in structural realization of amino acid code. ChiRotor: a side-chain prediction algorithm based on side-chain backbone interactions,” Protein Science, vol. 16, no. 3, pp. 494–506, 2007.View at: Publisher Site | Google Scholar
B. Brooks, R. E. Bruccoleri, B. D. Olafson, D. J. States, S. Swaminathan, and M. Karplus, “CHARMM: a program for macromolecular energy, minimization, and dynamics calculations,” Journal of Computational Chemistry, vol. 4, pp. 187–217, 1983.View at: Google Scholar
R. A. Laskowski, M. W. MacArthur, D. S. Moss, and J. M. Thornton, “PROCHECK-a program to check the stereochemical quality of protein structures,” Journal of Applied Crystallography, vol. 26, pp. 283–291, 1993.View at: Google Scholar
C. Colovos and T. O. Yeates, “Verification of protein structures: patterns of nonbonded atomic interactions,” Protein Science, vol. 2, no. 9, pp. 1511–1519, 1993.View at: Publisher Site | Google Scholar
R. Luthy, J. U. Bowie, and D. Eisenberg, “Assesment of protein models with three-dimensional profiles,” Nature, vol. 356, no. 6364, pp. 83–85, 1992.View at: Publisher Site | Google Scholar
M. Pawlowski, M. J. Gajda, R. Matlak, and J. M. Bujnicki, “MetaMQAP: a meta-server for the quality assessment of protein models,” BMC Bioinformatics, vol. 9, article 403, 2008.View at: Publisher Site | Google Scholar
D. Frishman and P. Argos, “Knowledge-based protein secondary structure assignment,” Proteins: Structure, Function and Genetics, vol. 23, no. 4, pp. 566–579, 1995.View at: Publisher Site | Google Scholar
K. G. Tina, R. Bhadra, and N. Srinivasan, “PIC: protein interactions calculator,” Nucleic Acids Research, vol. 35, no. 2, pp. W473–W476, 2007.View at: Publisher Site | Google Scholar
C. M. Venkatachalam, X. Jiang, T. Oldfield, and M. Waldman, “LigandFit: a novel method for the shape-directed rapid docking of ligands to protein active sites,” Journal of Molecular Graphics and Modelling, vol. 21, no. 4, pp. 289–307, 2003.View at: Publisher Site | Google Scholar
Z. Zhang, Y. Li, B. Lin, M. Schroeder, and B. Huang, “Identification of cavities on protein surface using multiple computational approaches for drug binding site prediction,” Bioinformatics, vol. 27, no. 15, Article ID btr331, pp. 2083–2088, 2011.View at: Publisher Site | Google Scholar
S. K. Boehlein, J. R. Shaw, L. C. Hannah, and J. D. Stewart, “Probing allosteric binding sites of the maize endosperm ADP-glucose pyrophosphorylase,” Plant Physiology, vol. 152, no. 1, pp. 85–95, 2010.View at: Publisher Site | Google Scholar
G. Wu, D. H. Robertson, C. L. Brooks, and M. Vieth, “Detailed analysis of grid-based molecular docking: a case study of CDOCKER—a CHARMm-based MD docking algorithm,” Journal of Computational Chemistry, vol. 24, no. 13, pp. 1549–1562, 2003.View at: Publisher Site | Google Scholar
K. Guruprasad, B. V. B. Reddy, and M. W. Pandit, “Correlation between stability of a protein and its dipeptide composition: a novel approach for predicting in vivo stability of a protein from its primary sequence,” Protein Engineering, vol. 4, no. 2, pp. 155–161, 1990.View at: Publisher Site | Google Scholar
H. J. Dyson and P. E. Wright, “Intrinsically unstructured proteins and their functions,” Nature Reviews Molecular Cell Biology, vol. 6, no. 3, pp. 197–208, 2005.View at: Publisher Site | Google Scholar
V. N. Uversky, C. J. Oldfield, and A. K. Dunker, “Showing your ID: intrinsic disorder as an ID for recognition, regulation and cell signaling,” Journal of Molecular Recognition, vol. 18, no. 5, pp. 343–384, 2005.View at: Publisher Site | Google Scholar
A. K. Dunker, J. D. Lawson, C. J. Brown et al., “Intrinsically disordered protein,” Journal of Molecular Graphics and Modelling, vol. 19, no. 1, pp. 26–59, 2001.View at: Publisher Site | Google Scholar
J. Kyte and R. F. Doolittle, “A simple method for displaying the hydropathic character of a protein,” Journal of Molecular Biology, vol. 157, no. 1, pp. 105–132, 1982.View at: Publisher Site | Google Scholar
J. Sivaraman, V. Sauvé, A. Matte, and M. Cygler, “Crystal structure of Escherichia coli glucose-1-phosphate thymidylyltransferase (RffH) complexed with dTTP and Mg2+,” The Journal of Biological Chemistry, vol. 277, no. 46, pp. 44214–44219, 2002.View at: Publisher Site | Google Scholar
A. A. Iglesias and J. Preiss, “Bacterial glycogen and plant starch biosynthesis,” Biochemical Education, vol. 20, no. 4, pp. 196–203, 1992.View at: Publisher Site | Google Scholar
N. A. Baker, D. Sept, S. Joseph, M. J. Holst, and J. A. McCammon, “Electrostatics of nanosystems: application to microtubules and the ribosome,” Proceedings of the National Academy of Sciences of the United States of America, vol. 98, no. 18, pp. 10037–10041, 2001.View at: Publisher Site | Google Scholar
C. Hylton and A. M. Smith, “The rb mutation of peas causes structural and regulatory changes in ADP glucose pyrophosphorylase from developing embryos,” Plant Physiology, vol. 99, no. 4, pp. 1626–1634, 1992.View at: Publisher Site | Google Scholar
M. Olive, R. J. Ellis, and W. W. Schuch, “Isolation and nucleotide sequences of cDNA clones encoding ADP-glucose pyrophosphorylase polypeptides from wheat leaf and endosperm,” Plant Molecular Biology, vol. 12, no. 5, pp. 525–538, 1989.View at: Publisher Site | Google Scholar
H. Weber, U. Heim, L. Borisjuk, and U. Wobus, “Cell-type specific, coordinate expression of two ADP-glucose pyrophosphorylase genes in relation to starch biosynthesis during seed development of Vicia faba L.,” Planta, vol. 195, no. 3, pp. 352–361, 1995.View at: Publisher Site | Google Scholar
D. N. P. Doan, H. Rudi, and O. A. Olsen, “The allosterically unregulated isoform of ADP-glucose pyrophosphorylase from barley endosperm is the most likely source of ADP-glucose incorporated into endosperm starch,” Plant Physiology, vol. 121, no. 3, pp. 965–975, 1999.View at: Publisher Site | Google Scholar