Phylogenetic analysis of ADP-glucose pyrophosphorylase subunits reveals a role of subunit interfaces in the allosteric properties of the enzyme

A critical inter-subunit interaction for the transmission of the allosteric signal in the Agrobacterium tumefaciens ADP-glucose pyrophosphorylase

ADP-glucose pyrophosphorylase is a key regulatory enzyme involved in starch and glycogen synthesis in plants and bacteria, respectively. It has been hypothesized that inter-subunit communications are important for the allosteric effect in this enzyme. However, no specific interactions have been identified as part of the regulatory signal. The enzyme from Agrobacterium tumefaciens is a homotetramer allosterically regulated by fructose 6-phosphate and pyruvate. Three pairs of distinct subunit-subunit interfaces are present. Here we focus on an interface that features two symmetrical interactions between Arg11 and Asp141 from one subunit with residues Asp141 and Arg11 of the neighbor subunit, respectively. Previously, scanning mutagenesis showed that a mutation at the Arg11 position disrupted the activation of the enzyme. Considering the distance of these residues from the allosteric and catalytic sites, we hypothesized that the interaction between Arg11 and Asp141 is critical for allosteric signaling rather than effector binding. To prove our hypothesis, we mutated those two sites (D141A, D141E, D141N, D141R, R11D, and R11K) and performed kinetic and binding analysis. Mutations that altered the charge affected the regulation the most. To prove that the interaction per se (rather than the presence of specific residues) is critical, we partially rescued the R11D protein by introducing a second mutation (R11D/D141R). This could not restore the activator effect on kcat , but it did rescue the effect on substrate affinity. Our results indicate the critical functional role of Arg11 and Asp141 to relay the allosteric signal in this subunit interface.

Genome-Wide Identification of DOF Gene Family and the Mechanism Dissection of SbDof21 Regulating Starch Biosynthesis in Sorghum

Starch is one of the main utilization products of sorghum (Sorghum bicolor L.), the fifth largest cereal crop in the world. Up to now, the regulation mechanism of starch biosynthesis is rarely documented in sorghum. In the present study, we identified 30 genes encoding the C2-C2 zinc finger domain (DOF), with one to three exons in the sorghum genome. The DOF proteins of sorghum were divided into two types according to the results of sequence alignment and evolutionary analysis. Based on gene expressions and co-expression analysis, we identified a regulatory factor, SbDof21, that was located on chromosome 5. SbDof21 contained two exons, encoding a 36.122 kD protein composed of 340 amino acids. SbDof21 co-expressed with 15 genes involved in the sorghum starch biosynthesis pathway, and the Pearson correlation coefficients (PCCs) with 11 genes were greater than 0.9. The results of qRT-PCR assays indicated that SbDof21 is highly expressed in sorghum grains, exhibiting low relative expression levels in the tissues of roots, stems and leaves. SbDOF21 presented as a typical DOF transcription factor (TF) that was localized to the nucleus and possessed transcriptional activation activity. Amino acids at positions 182–231 of SbDOF21 formed an important structure in its activation domain. The results of EMSA showed that SbDOF21 could bind to four tandem repeats of P-Box (TGTAAAG) motifs in vitro, such as its homologous proteins of ZmDOF36, OsPBF and TaPBF. Meanwhile, we also discovered that SbDOF21 could bind and transactivate SbGBSSI, a key gene in sorghum amylose biosynthesis. Collectively, the results of the present study suggest that SbDOF21 acts as an important regulator in sorghum starch biosynthesis, exhibiting potential values for the improvement of starch contents in sorghum.

Phosphorylation of ADP-Glucose Pyrophosphorylase During Wheat Seeds Development

Starch is the dominant reserve polysaccharide accumulated in the seed of grasses (like wheat). It is the most common carbohydrate in the human diet and a material applied to the bioplastics and biofuels industry. Hence, the complete understanding of starch metabolism is critical to design rational strategies to improve its allocation in plant reserve tissues. ADP-glucose pyrophosphorylase (ADP-Glc PPase) catalyzes the key (regulated) step in the synthetic starch pathway. The enzyme comprises a small (S) and a large (L) subunit forming an S₂L₂ heterotetramer, which is allosterically regulated by orthophosphate, fructose-6P, and 3P-glycerate. ADP-Glc PPase was found in a phosphorylated state in extracts from wheat seeds. The amount of the phosphorylated protein increased along with the development of the seed and correlated with relative increases of the enzyme activity and starch content. Conversely, this post-translational modification was absent in seeds from Ricinus communis. In vitro, the recombinant ADP-Glc PPase from wheat endosperm was phosphorylated by wheat seed extracts as well as by recombinant Ca²⁺-dependent plant protein kinases. Further analysis showed that the preferential phosphorylation takes place on the L subunit. Results suggest that the ADP-Glc PPase is a phosphorylation target in seeds from grasses but not from oleaginous plants. Accompanying seed maturation and starch accumulation, a combined regulation of ADP-Glc PPase by metabolites and phosphorylation may provide an enzyme with stable levels of activity. Such concerted modulation would drive carbon skeletons to the synthesis of starch for its long-term storage, which later support seed germination.

Genetic and genomic analysis of the seed-filling process in maize based on a logistic model

Seed filling is a dynamic process that determines seed size and nutritional quality. This time-dependent trait follows a logistic (S-shaped) growth curve that can be described by a logistic function, with parameters of biological relevance. When compared between genotypes, the filling dynamics variations are explained by the differences of parameter values; as such, the parameter estimates can be considered as “traits” for genetic analysis to identify loci that are associated with the seed-filling process. We carried out genetic and genomic analysis of the seed-filling process in maize, using a recombinant inbred line (RIL) population derived from the two inbred lines with contrasting seed-filling dynamics. We recorded seed dry weight at 14 time points after pollination, spanning the early filling phases to the late maturation stages. Fitting these data to a logistic model allowed for estimating 12 characteristic parameters that can be used to meaningfully describe the seed-filling process. Quantitative trait locus (QTL) mapping of these parameters identified a total of 90 nonredundant loci. Using bulked segregant RNA-sequencing (BSR-seq) analysis, we identified eight genes that showed differential gene expression patterns at multiple time points between the extreme pools, and these genes co-localize with the mapped QTL regions. Two of the eight genes, GRMZM2G391936 and GRMZM2G008263, are implicated in starch and sucrose metabolism, and biosynthesis of secondary metabolites that are well known for playing a vital role in seed filling. This study suggests that the logistic model-based approach can efficiently identify genetic loci that regulate dynamic developing traits.

Enhancement of Heat Stability and Kinetic Parameters of the Maize Endosperm ADP-Glucose Pyrophosphorylase by Mutagenesis of Amino Acids in the Small Subunit With High B Factors

ADP-glucose pyrophosphorylase (AGPase) is an important enzyme in starch synthesis and previous studies showed that the heat lability of this enzyme is a determinant to starch synthesis in the maize endosperm and, in turn, seed yield. Here, amino acids in the AGPase endosperm small subunit with high B-factors were mutagenized and individual changes enhancing heat stability and/or kinetic parameters in an Escherichia coli expression system were chosen. Individual mutations were combined and analyzed. One triple mutant, here termed Bt2-BF, was chosen for further study. Combinations of this heat stable, 3-PGA-independent small subunit variant with large subunits also heat stable yielded complex patterns of heat stability and kinetic and allosteric properties. Interestingly, two of the three changes reside in a protein motif found only in AGPases that exhibit high sensitivity to 3-PGA. While not the 3-PGA binding site, amino acid substitutions in this region significantly alter 3-PGA activation kinetics.

Evolutionary, structural and expression analysis of core genes involved in starch synthesis

Starch is the main storage carbohydrate in plants and an important natural resource for food, feed and industrial raw materials. However, the details regarding the pathway for starch biosynthesis and the diversity of biosynthetic enzymes involved in this process are poorly understood. This study uses a comprehensive phylogenetic analysis of 74 sequenced plant genomes to revisit the evolutionary history of the genes encoding ADP-glucose pyrophosphorylase (AGPase), starch synthase (SS), starch branching enzyme (SBE) and starch de-branching enzyme (DBE). Additionally, the protein structures and expression patterns of these four core genes in starch biosynthesis were studied to determine their functional differences. The results showed that AGPase, SS, SBE and DBE have undergone complicated evolutionary processes in plants and that gene/genome duplications are responsible for the observed differences in isoform numbers. A structure analysis of these proteins suggested that the deletion/mutation of amino acids in some active sites resulted in not only structural variation but also sub-functionalization or neo-functionalization. Expression profiling indicated that AGPase-, SS-, SBE- and DBE-encoding genes exhibit spatio-temporally divergent expression patterns related to the composition of functional complexes in starch biosynthesis. This study provides a comprehensive atlas of the starch biosynthetic pathway, and these data should support future studies aimed at increasing understanding of starch biosynthesis and the functional evolutionary divergence of AGPase, SS, SBE, and DBE in plants.

Resurrecting the Regulatory Properties of the Ostreococcus tauri ADP-Glucose Pyrophosphorylase Large Subunit

ADP-glucose pyrophosphorylase (ADP-Glc PPase) catalyzes the first committed step for the synthesis of glycogen in cyanobacteria and starch in green algae and plants. The enzyme from cyanobacteria is homotetrameric (α₄), while that from green algae and plants is heterotetrameric (α₂β₂). These ADP-Glc PPases are allosterically regulated by 3-phosphoglycerate (3PGA, activator) and inorganic orthophosphate (Pi, inhibitor). Previous studies on the cyanobacterial and plant enzymes showed that 3PGA binds to two highly conserved Lys residues located in the C-terminal domain. We observed that both Lys residues are present in the small (α) subunit of the Ostreococcus tauri enzyme; however, one of these Lys residues is replaced by Arg in the large (β) subunit. In this work, we obtained the K443R and R466K mutants of the O. tauri small and large subunits, respectively, and co-expressed them together or with their corresponding wild type counterparts. Our results show that restoring the Lys residue in the large subunit enhanced 3PGA affinity, whereas introduction of an Arg residue in the small subunit reduced 3PGA affinity of the heterotetramers. Inhibition kinetics also showed that heterotetramers containing the K443R small subunit mutant were less sensitive to Pi inhibition, but only minor changes were observed for those containing the R466K large subunit mutant, suggesting a leading role of the small subunit for Pi inhibition of the heterotetramer. We conclude that, during evolution, the ADP-Glc PPase large subunit from green algae and plants acquired mutations in its regulatory site. The rationale for this could have been to accommodate sensitivity to particular metabolic needs of the cell or tissue.

A brittle-2 transgene increases maize yield by acting in maternal tissues to increase seed number

The enzyme ADP-glucose pyrophosphorylase is essential for starch biosynthesis and is highly regulated. Here, mutations that increased heat stability and interactions with allosteric effectors were incorporated into the small subunit of the isoform known to be expressed at high levels in the maize endosperm. The resulting variants were transformed into maize with expression targeted to the endosperm. Transgenes harboring the changes increased yield some 35%; however, yield enhancement occurred via an increase in seed number rather than by increased seed weight. Interestingly, seed number increase is controlled by the genotype of the plant rather than the genotype of the seed as seeds increase in number whether or not they contain the transgene as long as the maternal parent has the transgene. The transgene is however expressed in the endosperm, and the altered allosteric and stability properties initially seen in Escherichia coli expression experiments are also seen with the endosperm-expressed gene. The extent of seed number increase is positively correlated with the average daily high temperature during the first 4 days postpollination. While these results were unexpected, they echo the phenotypic changes caused by the insertion of an altered large subunit of this enzyme reported previously (Plant Cell, 24, 2012, 2352). These results call into question some of the reported fundamental differences separating starch synthesis in the endosperm vis-à-vis other plant tissues.

Lineage-Specific Evolutionary Histories and Regulation of Major Starch Metabolism Genes during Banana Ripening

Starch is the most widespread and abundant storage carbohydrate in plants. It is also a major feature of cultivated bananas as it accumulates to large amounts during banana fruit development before almost complete conversion to soluble sugars during ripening. Little is known about the structure of major gene families involved in banana starch metabolism and their evolution compared to other species. To identify genes involved in banana starch metabolism and investigate their evolutionary history, we analyzed six gene families playing a crucial role in plant starch biosynthesis and degradation: the ADP-glucose pyrophosphorylases (AGPases), starch synthases (SS), starch branching enzymes (SBE), debranching enzymes (DBE), α-amylases (AMY) and β-amylases (BAM). Using comparative genomics and phylogenetic approaches, these genes were classified into families and sub-families and orthology relationships with functional genes in Eudicots and in grasses were identified. In addition to known ancestral duplications shaping starch metabolism gene families, independent evolution in banana and grasses also occurred through lineage-specific whole genome duplications for specific sub-families of AGPase, SS, SBE, and BAM genes; and through gene-scale duplications for AMY genes. In particular, banana lineage duplications yielded a set of AGPase, SBE and BAM genes that were highly or specifically expressed in banana fruits. Gene expression analysis highlighted a complex transcriptional reprogramming of starch metabolism genes during ripening of banana fruits. A differential regulation of expression between banana gene duplicates was identified for SBE and BAM genes, suggesting that part of starch metabolism regulation in the fruit evolved in the banana lineage.

Glu-370 in the large subunit influences the substrate binding, allosteric, and heat stability properties of potato ADP-glucose pyrophosphorylase

ADP-glucose pyrophosphorylase (AGPase) is a key allosteric enzyme in plant starch biosynthesis. Plant AGPase is a heterotetrameric enzyme that consists of large (LS) and small subunits (SS), which are encoded by two different genes. In this study, we showed that the conversion of Glu to Gly at position 370 in the LS of AGPase alters the heterotetrameric stability along with the binding properties of substrate and effectors of the enzyme. Kinetic analyses revealed that the affinity of the LS^E370GSS^WT AGPase for glucose-1-phosphate is 3-fold less than for wild type (WT) AGPase. Additionally, the LS^E370GSS^WT AGPase requires 3-fold more 3-phosphogyceric acid to be activated. Finally, the LS^E370GSS^WTAGPase is less heat stable compared with the WT AGPase. Computational analysis of the mutant Gly-370 in the 3D modeled LS AGPase showed that this residue changes charge distribution of the surface and thus affect stability of the LS AGPase and overall heat stability of the heterotetrameric AGPase. In summary, our results show that LS^E370 intricately modulate the heat stability and enzymatic activity of potato the AGPase.

Mapping of quantitative trait loci for tuber starch and leaf sucrose contents in diploid potato

Most QTL for leaf sucrose content map to positions that are similar to positions of QTL for tuber starch content in diploid potato. In the present study, using a diploid potato mapping population and Diversity Array Technology (DArT) markers, we identified twelve quantitative trait loci (QTL) for tuber starch content on seven potato chromosomes: I, II, III, VIII, X, XI, and XII. The most important QTL spanned a wide region of chromosome I (42.0–104.6 cM) with peaks at 63 and 84 cM which explained 17.6 and 19.2% of the phenotypic variation, respectively. ADP-glucose pyrophosphorylase (AGPase) is the key enzyme for starch biosynthesis. The gene encoding the large subunit of this enzyme, AGPaseS-a, was localized to chromosome I at 102.3 cM and accounted for 15.2% of the variance in tuber starch content. A more than 100-fold higher expression of this gene was observed in RT-qPCR assay in plants with the marker allele AGPaseS-a1334. This study is the first to report QTL for sucrose content in potato leaves. QTL for sucrose content in leaves were located on eight potato chromosomes: I, II, III, V, VIII, IX, X and XII. In 5-week-old plants, only one QTL for leaf sucrose content was detected after 8 h of darkness; four QTL were detected after 8 h of illumination. In 11-week-old plants, 6 and 3 QTL were identified after dark and light phases, respectively. Of fourteen QTL for leaf sucrose content, eleven mapped to positions that were similar to QTL for tuber starch content. These results provide genetic information for further research examining the relationships between metabolic carbon molecule sources and sinks in potato plants.

AGPase: its role in crop productivity with emphasis on heat tolerance in cereals

AGPase, a key enzyme of starch biosynthetic pathway, has a significant role in crop productivity. Thermotolerant variants of AGPase in cereals may be used for developing cultivars, which may enhance productivity under heat stress. Improvement of crop productivity has always been the major goal of plant breeders to meet the global demand for food. However, crop productivity itself is influenced in a large measure by a number of abiotic stresses including heat, which causes major losses in crop productivity. In cereals, crop productivity in terms of grain yield mainly depends upon the seed starch content so that starch biosynthesis and the enzymes involved in this process have been a major area of investigation for plant physiologists and plant breeders alike. Considerable work has been done on AGPase and its role in crop productivity, particularly under heat stress, because this enzyme is one of the major enzymes, which catalyses the rate-limiting first committed key enzymatic step of starch biosynthesis. Keeping the above in view, this review focuses on the basic features of AGPase including its structure, regulatory mechanisms involving allosteric regulators, its sub-cellular localization and its genetics. Major emphasis, however, has been laid on the genetics of AGPases and its manipulation for developing high yielding cultivars that will have comparable productivity under heat stress. Some important thermotolerant variants of AGPase, which mainly involve specific amino acid substitutions, have been highlighted, and the prospects of using these thermotolerant variants of AGPase in developing cultivars for heat prone areas have been discussed. The review also includes a brief account on transgenics for AGPase, which have been developed for basic studies and crop improvement.

AGPase: its role in crop productivity with emphasis on heat tolerance in cereals

AGPase, a key enzyme of starch biosynthetic pathway, has a significant role in crop productivity. Thermotolerant variants of AGPase in cereals may be used for developing cultivars, which may enhance productivity under heat stress. Improvement of crop productivity has always been the major goal of plant breeders to meet the global demand for food. However, crop productivity itself is influenced in a large measure by a number of abiotic stresses including heat, which causes major losses in crop productivity. In cereals, crop productivity in terms of grain yield mainly depends upon the seed starch content so that starch biosynthesis and the enzymes involved in this process have been a major area of investigation for plant physiologists and plant breeders alike. Considerable work has been done on AGPase and its role in crop productivity, particularly under heat stress, because this enzyme is one of the major enzymes, which catalyses the rate-limiting first committed key enzymatic step of starch biosynthesis. Keeping the above in view, this review focuses on the basic features of AGPase including its structure, regulatory mechanisms involving allosteric regulators, its sub-cellular localization and its genetics. Major emphasis, however, has been laid on the genetics of AGPases and its manipulation for developing high yielding cultivars that will have comparable productivity under heat stress. Some important thermotolerant variants of AGPase, which mainly involve specific amino acid substitutions, have been highlighted, and the prospects of using these thermotolerant variants of AGPase in developing cultivars for heat prone areas have been discussed. The review also includes a brief account on transgenics for AGPase, which have been developed for basic studies and crop improvement.

AGPase: its role in crop productivity with emphasis on heat tolerance in cereals

AGPase, a key enzyme of starch biosynthetic pathway, has a significant role in crop productivity. Thermotolerant variants of AGPase in cereals may be used for developing cultivars, which may enhance productivity under heat stress. Improvement of crop productivity has always been the major goal of plant breeders to meet the global demand for food. However, crop productivity itself is influenced in a large measure by a number of abiotic stresses including heat, which causes major losses in crop productivity. In cereals, crop productivity in terms of grain yield mainly depends upon the seed starch content so that starch biosynthesis and the enzymes involved in this process have been a major area of investigation for plant physiologists and plant breeders alike. Considerable work has been done on AGPase and its role in crop productivity, particularly under heat stress, because this enzyme is one of the major enzymes, which catalyses the rate-limiting first committed key enzymatic step of starch biosynthesis. Keeping the above in view, this review focuses on the basic features of AGPase including its structure, regulatory mechanisms involving allosteric regulators, its sub-cellular localization and its genetics. Major emphasis, however, has been laid on the genetics of AGPases and its manipulation for developing high yielding cultivars that will have comparable productivity under heat stress. Some important thermotolerant variants of AGPase, which mainly involve specific amino acid substitutions, have been highlighted, and the prospects of using these thermotolerant variants of AGPase in developing cultivars for heat prone areas have been discussed. The review also includes a brief account on transgenics for AGPase, which have been developed for basic studies and crop improvement.

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.

Enhanced heterotetrameric assembly of potato ADP-glucose pyrophosphorylase using reverse genetics

ADP-glucose pyrophosphorylase (AGPase) is a key allosteric enzyme in plant starch biosynthesis. Plant AGPase is a heterotetrameric enzyme that consists of large (LS) and small subunits (SS), which are encoded by two different genes. Computational and experimental studies have revealed that the heterotetrameric assembly of AGPase is thermodynamically weak. Modeling studies followed by the mutagenesis of the LS of the potato AGPase identified a heterotetramer-deficient mutant, LS(R88A). To enhance heterotetrameric assembly, LS(R88A) cDNA was subjected to error-prone PCR, and second-site revertants were identified according to their ability to restore glycogen accumulation, as assessed with iodine staining. Selected mutations were introduced into the wild-type (WT) LS and co-expressed with the WT SS in Escherichia coli glgC(-). The biochemical characterization of revertants revealed that LS(I90V)SS(WT), LS(Y378C)SS(WT) and LS(D410G)SS(WT) mutants displayed enhanced heterotetrameric assembly with the WT SS. Among these mutants, LS(Y378C)SS(WT) AGPase displayed increased heat stability compared with the WT enzyme. Kinetic characterization of the mutants indicated that the LS(I90V)SS(WT) and LS(Y378C)SS(WT) AGPases have comparable allosteric and kinetic properties. However, the LS(D410G)SS(WT) mutant exhibited altered allosteric properties of being less responsive and more sensitive to 3-phosphoglyceric acid activation and inorganic phosphate inhibition. This study not only enhances our understanding of the interaction between the SS and the LS of AGPase but also enables protein engineering to obtain enhanced assembled heat-stable variants of AGPase, which can be used for the improvement of plant yields.

Enhanced heat stability and kinetic parameters of maize endosperm ADPglucose pyrophosphorylase by alteration of phylogenetically identified amino acids

ADP-glucose pyrophosphorylase (AGPase) controls the rate-limiting step in starch biosynthesis and is regulated at various levels. Cereal endosperm enzymes, in contrast to other plant AGPases, are particularly heat labile and transgenic studies highlight the importance of temperature for cereal yield. Previously, a phylogenetic approach identified Type II and positively selected amino acid positions in the large subunit of maize endosperm AGPase. Glycogen content, kinetic parameters and heat stability were measured in AGPases having mutations in these sites and interesting differences were observed. This study expands on our earlier evolutionary work by determining how all Type II and positively selected sites affect kinetic constants, heat stability and catalytic rates at increased temperatures. Variants with enhanced properties were identified and combined into one gene, designated Sh2-E. Enhanced properties include: heat stability, enhanced activity at 37 °C, activity at 55 °C, reduced Ka and activity in the absence of activator. The resulting enzyme exhibited all improved properties of the various individual changes. Additionally, Sh2-E was expressed with a small subunit variant with enhanced enzyme properties resulting in an enzyme that has exceptional heat stability, a high catalytic rate at increased temperatures and significantly decreased Km values for both substrates in the absence of the activator.

Insight into the 3D structure of ADP-glucose pyrophosphorylase from rice (Oryza sativa L.)

ADP-glucose pyrophosphorylase (E.C. 2.7.7.27; AGPase) is a key regulatory enzyme that catalyzes the rate-limiting step of starch biosynthesis in higher plants. AGPase consists of pair of small (SS) and large (LS) subunits thereby constituting a heterotetrameric structure. No crystal structure of the native heterotetrameric enzyme is available for any species, thus limiting the complete understanding of structure-function relationships of this enzyme. In this study, an attempt was made to deduce the heterotetrameric assembly of AGPase in rice. Homology modeling of the three-dimensional structure of the LS and SS was performed using the Swiss Model Server, and the models were evaluated and docked using GRAMM-X to obtain the stable heterodimer orientation (LS as receptor and SS as ligand) and then the heterotetrameric orientation. The initial heterotetrameric orientation was further refined using the RosettaDock Server. MD simulation of the representative heterodimer/tetramer was performed using NAMD, which indicated that the tail-to-tail interaction of LS and SS was more stable than the head-to-head orientation, and the heterotetramer energy was also minimized to -767,011 kcal mol(-1). Subunit-subunit interaction studies were then carried out using the programs NACCESS and Dimplot. A total of 57 interface residues were listed in SS and 63 in LS. The residues plotted by Dimplot were similar to those listed by NACCESS. Multiple sequence alignment of the sequences of LS and SS from potato, maize and rice validated the interactions inferred in the study. RMSD of 1.093 Å was obtained on superimposition of the deduced heterotetramer on the template homo-tetramer (1YP2), showing the similarity between the two structures.

Insight into the 3D structure of ADP-glucose pyrophosphorylase from rice (Oryza sativa L.)

ADP-glucose pyrophosphorylase (E.C. 2.7.7.27; AGPase) is a key regulatory enzyme that catalyzes the rate-limiting step of starch biosynthesis in higher plants. AGPase consists of pair of small (SS) and large (LS) subunits thereby constituting a heterotetrameric structure. No crystal structure of the native heterotetrameric enzyme is available for any species, thus limiting the complete understanding of structure-function relationships of this enzyme. In this study, an attempt was made to deduce the heterotetrameric assembly of AGPase in rice. Homology modeling of the three-dimensional structure of the LS and SS was performed using the Swiss Model Server, and the models were evaluated and docked using GRAMM-X to obtain the stable heterodimer orientation (LS as receptor and SS as ligand) and then the heterotetrameric orientation. The initial heterotetrameric orientation was further refined using the RosettaDock Server. MD simulation of the representative heterodimer/tetramer was performed using NAMD, which indicated that the tail-to-tail interaction of LS and SS was more stable than the head-to-head orientation, and the heterotetramer energy was also minimized to -767,011 kcal mol(-1). Subunit-subunit interaction studies were then carried out using the programs NACCESS and Dimplot. A total of 57 interface residues were listed in SS and 63 in LS. The residues plotted by Dimplot were similar to those listed by NACCESS. Multiple sequence alignment of the sequences of LS and SS from potato, maize and rice validated the interactions inferred in the study. RMSD of 1.093 Å was obtained on superimposition of the deduced heterotetramer on the template homo-tetramer (1YP2), showing the similarity between the two structures.

Accelerated evolution and coevolution drove the evolutionary history of AGPase sub-units during angiosperm radiation

BACKGROUND AND AIMS

ADP-glucose pyrophosphorylase (AGPase) is a key enzyme of starch biosynthesis. In the green plant lineage, it is composed of two large (LSU) and two small (SSU) sub-units encoded by paralogous genes, as a consequence of several rounds of duplication. First, our aim was to detect specific patterns of molecular evolution following duplication events and the divergence between monocotyledons and dicotyledons. Secondly, we investigated coevolution between amino acids both within and between sub-units.

METHODS

A phylogeny of each AGPase sub-unit was built using all gymnosperm and angiosperm sequences available in databases. Accelerated evolution along specific branches was tested using the ratio of the non-synonymous to the synonymous substitution rate. Coevolution between amino acids was investigated taking into account compensatory changes between co-substitutions.

KEY RESULTS

We showed that SSU paralogues evolved under high functional constraints during angiosperm radiation, with a significant level of coevolution between amino acids that participate in SSU major functions. In contrast, in the LSU paralogues, we identified residues under positive selection (1) following the first LSU duplication that gave rise to two paralogues mainly expressed in angiosperm source and sink tissues, respectively; and (2) following the emergence of grass-specific paralogues expressed in the endosperm. Finally, we found coevolution between residues that belong to the interaction domains of both sub-units.

CONCLUSIONS

Our results support the view that coevolution among amino acid residues, especially those lying in the interaction domain of each sub-unit, played an important role in AGPase evolution. First, within SSU, coevolution allowed compensating mutations in a highly constrained context. Secondly, the LSU paralogues probably acquired tissue-specific expression and regulatory properties via the coevolution between sub-unit interacting domains. Finally, the pattern we observed during LSU evolution is consistent with repeated sub-functionalization under 'Escape from Adaptive Conflict', a model rarely illustrated in the literature.

The condensed chromatin fiber: an allosteric chemo-mechanical machine for signal transduction and genome processing

Allostery is a key concept of molecular biology which refers to the control of an enzyme activity by an effector molecule binding the enzyme at another site rather than the active site (allos = other in Greek). We revisit here allostery in the context of chromatin and argue that allosteric principles underlie and explain the functional architecture required for spacetime coordination of gene expression at all scales from DNA to the whole chromosome. We further suggest that this functional architecture is provided by the chromatin fiber itself. The structural, mechanical and topological features of the chromatin fiber endow chromosomes with a tunable signal transduction from specific (or nonspecific) effectors to specific (or nonspecific) active sites. Mechanical constraints can travel along the fiber all the better since the fiber is more compact and regular, which speaks in favor of the actual existence of the (so-called 30 nm) chromatin fiber. Chromatin fiber allostery reconciles both the physical and biochemical approaches of chromatin. We illustrate this view with two supporting specific examples. Moreover, from a methodological point of view, we suggest that the notion of chromatin fiber allostery is particularly relevant for systemic approaches. Finally we discuss the evolutionary power of allostery in the context of chromatin and its relation to modularity.

Contrasted patterns of selection since maize domestication on duplicated genes encoding a starch pathway enzyme

Maize domestication from teosinte (Zea mays ssp. parviglumis) was accompanied by an increase of kernel size in landraces. Subsequent breeding has led to a diversification of kernel size and starch content among major groups of inbred lines. We aim at investigating the effect of domestication on duplicated genes encoding a key enzyme of the starch pathway, the ADP-glucose pyrophosphorylase (AGPase). Three pairs of paralogs encode the AGPase small (SSU) and large (LSU) subunits mainly expressed in the endosperm, the embryo and the leaf. We first validated the putative sequence of LSU(leaf) through a comparative expression assay of the six genes. Second, we investigated the patterns of molecular evolution on a 2 kb coding region homologous among the six genes in three panels: teosintes, landraces, and inbred lines. We corrected for demographic effects by relying on empirical distributions built from 580 previously sequenced ESTs. We found contrasted patterns of selection among duplicates: three genes exhibit patterns of directional selection during domestication (SSU(end), LSU(emb)) or breeding (LSU(leaf)), two exhibit patterns consistent with diversifying (SSU(leaf)) and balancing selection (SSU(emb)) accompanying maize breeding. While patterns of linkage disequilibrium did not reveal sign of coevolution between genes expressed in the same organ, we detected an excess of non-synonymous substitutions in the small subunit functional domains highlighting their role in AGPase evolution. Our results offer a different picture on AGPase evolution than the one depicted at the Angiosperm level and reveal how genetic redundancy can provide flexibility in the response to selection.

Statistical properties of the branch-site test of positive selection

The branch-site test is a likelihood ratio test to detect positive selection along prespecified lineages on a phylogeny that affects only a subset of codons in a protein-coding gene, with positive selection indicated by accelerated nonsynonymous substitutions (with ω = d(N)/d(S) > 1). This test may have more power than earlier methods, which average nucleotide substitution rates over sites in the protein and/or over branches on the tree. However, a few recent studies questioned the statistical basis of the test and claimed that the test generated too many false positives. In this paper, we examine the null distribution of the test and conduct a computer simulation to examine the false-positive rate and the power of the test. The results suggest that the asymptotic theory is reliable for typical data sets, and indeed in our simulations, the large-sample null distribution was reliable with as few as 20-50 codons in the alignment. We examined the impact of sequence length, the strength of positive selection, and the proportion of sites under positive selection on the power of the branch-site test. We found that the test was far more powerful in detecting episodic positive selection than branch-based tests, which average substitution rates over all codons in the gene and thus miss the signal when most codons are under strong selective constraint. Recent claims of statistical problems with the branch-site test are due to misinterpretations of simulation results. Our results, as well as previous simulation studies that have demonstrated the robustness of the test, suggest that the branch-site test may be a useful tool for detecting episodic positive selection and for generating biological hypotheses for mutation studies and functional analyses. The test is sensitive to sequence and alignment errors and caution should be exercised concerning its use when data quality is in doubt.

Characterization of the AGPase large subunit isoforms from tomato indicates that the recombinant L3 subunit is active as a monomer

The enzyme AGPase [ADP-Glc (glucose) pyrophosphorylase] catalyses a rate-limiting step in starch synthesis in tomato (Solanum lycopersicon) fruit, which undergoes a transient period of starch accumulation. It has been a generally accepted paradigm in starch metabolism that the enzyme naturally functions primarily as a heterotetramer comprised of two large subunits (L) and two small subunits (S). The tomato genome harbours a single gene encoding S and three genes for L proteins, which are expressed in both a tissue- and time-specific manner. In the present study the allosteric contributions of the different L subunits were compared by expressing each one in Escherichia coli, in conjunction with S and individually, and characterizing the resulting enzyme activity. Our results indicate different kinetic characteristics of the tomato L1/S and L3/S heterotetramers. Surprisingly, the recombinant L3 protein was also active when expressed alone and size-exclusion and immunoblotting showed that it functioned as a monomer. Subunit interaction modelling pointed to two amino acids potentially affecting subunit interactions. However, directed mutations did not have an impact on subunit tetramerization. These results indicate a hitherto unknown active role for the L subunit in the synthesis of ADP-Glc.

Probing allosteric binding sites of the maize endosperm ADP-glucose pyrophosphorylase

Maize (Zea mays) endosperm ADP-glucose pyrophosphorylase (AGPase) is a highly regulated enzyme that catalyzes the rate-limiting step in starch biosynthesis. Although the structure of the heterotetrameric maize endosperm AGPase remains unsolved, structures of a nonnative, low-activity form of the potato tuber (Solanum tuberosum) AGPase (small subunit homotetramer) reported previously by others revealed that several sulfate ions bind to each enzyme. These sites are also believed to interact with allosteric regulators such as inorganic phosphate and 3-phosphoglycerate (3-PGA). Several arginine (Arg) side chains contact the bound sulfate ions in the potato structure and likely play important roles in allosteric effector binding. Alanine-scanning mutagenesis was applied to the corresponding Arg residues in both the small and large subunits of maize endosperm AGPase to determine their roles in allosteric regulation and thermal stability. Steady-state kinetic and regulatory parameters were measured for each mutant. All of the Arg mutants examined--in both the small and large subunits--bound 3-PGA more weakly than the wild type (A(50) increased by 3.5- to 20-fold). By contrast, the binding of two other maize AGPase allosteric activators (fructose-6-phosphate and glucose-6-phosphate) did not always mimic the changes observed for 3-PGA. In fact, compared to 3-PGA, fructose-6-phosphate is a more efficient activator in two of the Arg mutants. Phosphate binding was also affected by Arg substitutions. The combined data support a model for the binding interactions associated with 3-PGA in which allosteric activators and inorganic phosphate compete directly.