Investigation of the Interaction between the Large and Small Subunits of Potato ADP-Glucose Pyrophosphorylase

ADP-glucose pyrophosphorylase (AGPase), a key allosteric enzyme involved in higher plant starch biosynthesis, is composed of pairs of large (LS) and small subunits (SS). Current evidence indicates that the two subunit types play distinct roles in enzyme function. Recently the heterotetrameric structure of potato AGPase has been modeled. In the current study, we have applied the molecular mechanics generalized born surface area (MM-GBSA) method and identified critical amino acids of the potato AGPase LS and SS subunits that interact with each other during the native heterotetrameric structure formation. We have further shown the role of the LS amino acids in subunit-subunit interaction by yeast two-hybrid, bacterial complementation assay and native gel. Comparison of the computational results with the experiments has indicated that the backbone energy contribution (rather than the side chain energies) of the interface residues is more important in identifying critical residues. We have found that lateral interaction of the LS-SS is much stronger than the longitudinal one, and it is mainly mediated by hydrophobic interactions. This study will not only enhance our understanding of the interaction between the SS and the LS of AGPase, but will also enable us to engineer proteins to obtain better assembled variants of AGPase which can be used for the improvement of plant yield.     

Investigation of subunit function in ADP-glucose pyrophosphorylase

ADP-glucose pyrophosphorylase (AGPase), a key regulatory enzyme in higher plant starch biosynthesis, is composed of a pair of large and small subunits (alpha(2)beta(2)). Current evidence suggests that the large subunit has primarily a regulatory function, while the small subunit has both regulatory and catalytic roles. To define the structure-function relationship of the large subunit (LS), the LS of potato AGPase was subjected to chemical mutagenesis and coexpressed with the wild-type (WT) small subunit (SS) cDNA in an AGPase defective Escherichia coli strain. An LS mutant (M143) was isolated, which accumulated very low levels of glycogen compared to the WT recombinant AGPase, but maintained normal catalytic activity when assayed under saturating conditions. Sequence analysis revealed that M143 has a single amino acid change, V463I, which lies adjacent to the C-terminus. This single mutation had no effect on the Km for ATP and Mg(2+), which were similar to the WT enzyme. The K(m) for glucose 1-P, however, was sixfold higher than the WT enzyme. These results suggest that the LS plays a role in binding glucose 1-P through its interaction with the SS.     

Generation,characterization, and heterologous expression of wild-type and up-regulated forms of Arabidopsis thaliana leaf ADP-glucose pyrophosphorylase

ADP-glucose pyrophosphorylase (AGPase), a key enzyme in starch biosynthesis of higher plants, consists of a pair of regulatory large (LS) and catalytically small (SS) subunits. In plants, these subunits are coded by multiple genes resulting in the formation of tissue-specific enzyme forms, which are differentially regulated during plant growth and development. Some AGPase isoforms differ in catalytic and regulatory properties as well as intracellular location. In an effort to gain a better understanding of the role of the leaf AGPase in carbon partitioning and its effect on plant productivity, the Arabidopsis leaf AGPase containing the mature forms of the SS and LS was expressed in a heterologous expression system and characterized enzymatically. The Arabidopsis recombinant AGPase had kinetic values for 3-phosphoglyceric acid, glucose-1-phosphate and Mg(2+) similar to those of the native enzyme. As the N-terminus of the LS has been suggested to be involved in enzyme function, the length of the N-terminal region was extended or shortened. Of the five modified LSs analyzed, only the T5 form lacking six residues of the mature N-terminus was able to form detectable levels of enzyme activity, indicating that the N-terminal region is critical for enzyme function. Two up-regulatory LS mutations that allosterically activate the potato enzyme, a stem isoform, were introduced into the corresponding Arabidopsis LS sequences and co-expressed with wild-type SS. Both modified enzymes showed up-regulatory properties, indicating that these specific residue changes were also operational in the leaf isoform.     

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.     

Allosteric regulation of the higher plant ADP-glucose pyrophosphorylase is a product of synergy between the two subunits

The higher plant ADP-glucose pyrophosphorylase (AGPase) is a heterotetramer consisting of two regulatory large subunits (LSs) and two catalytic small subunits (SSs). To further characterize the roles of these subunits in determining enzyme function, different combinations of wildtype LS (LWT) and variant forms (LUpReg1, LM345) were co-expressed with wildtype SS (SWT) and variant forms (STG-15 and Sdevo330) and their enzyme properties compared to those measured for the heterotetrameric wildtype enzyme and SS homotetrameric enzymes. Analysis of the allosteric regulatory properties of the various enzymes indicates that although the LS is required for optimal activation by 3-phosphoglyceric acid and resistance to Pi, the overall allosteric regulatory and kinetic properties are specified by both subunits. Our results show that the regulatory and kinetic properties of AGPase are not simply due to the LS modulating the properties of the SS but, instead, are a product of synergistic interaction between the two subunits.     

Isolation and characterization of cDNA clones encoding ADP-glucose pyrophosphorylase (AGPase) large and small subunits from chickpea (Cicer arietinum L.)

Four cDNA clones encoding two large subunits and two small subunits of the starch regulatory enzyme ADP-glucose pyrophosphorylase (AGPase) were isolated from a chickpea (Cicer arietinum L.) stem cDNA library. DNA sequence and Southern blot analyses of these clones, designated CagpL1, CagpL2 (large subunits) and CagpS1 and CagpS2 (small subunits), revealed that these isoforms represented different AGPase large and small subunits. RNA expression analysis indicated that CagpL1 was expressed strongly in leaves with reduced expression in the stem. No detectable expression was observed in seeds and roots. CagpL2 was expressed moderately in seeds followed by weak expression in leaves, stems and roots. Similar analysis showed that CagpS1 and CagpS2 displayed a spatial expression pattern similar to that observed for CagpL2 with the exception that CagpS1 showed a much higher expression in seeds than CagpS2. The spatial expression patterns of these different AGPase subunit sequences indicate that different AGPase isoforms are used to control starch biosynthesis in different organs during chickpea development.     

Duplications and functional divergence of ADP-glucose pyrophosphorylase genes in plants

Background  

ADP-glucose pyrophosphorylase (AGPase), which catalyses a rate limiting step in starch synthesis, is a heterotetramer comprised of two identical large and two identical small subunits in plants. Although the large and small subunits are equally sensitive to activity-altering amino acid changes when expressed in a bacterial system, the overall rate of non-synonymous evolution is ~2.7-fold greater for the large subunit than for the small subunit. Herein, we examine the basis for their different rates of evolution, the number of duplications in both large and small subunit genes and document changes in the patterns of AGPase evolution over time.     

Catalytic implications of the higher plant ADP-glucose pyrophosphorylase large subunit

ADP-glucose pyrophosphorylase, a key regulatory enzyme of starch biosynthesis, is composed of a pair of catalytic small subunits (SSs) and a pair of catalytically disabled large subunits (LSs). The N-terminal region of the LS has been known to be essential for the allosteric regulatory properties of the heterotetrameric enzyme. To gain further insight on the role of this region and the LS itself in enzyme function, the six proline residues found in the N-terminal region of the potato tuber AGPase were subjected to scanning mutagenesis. The wildtype and various mutant heterotetramers were expressed using our newly developed host-vector system, purified, and their kinetic parameters assessed. While P(17)L, P(26)L, and P(55)L mutations only moderately affected the kinetic properties, P(52)L and P(66)L gave rise to significant and contrasting changes in allosteric properties: P(66)L enzyme displayed up-regulatory properties toward 3-PGA while the P(52)L enzyme had down-regulatory properties. Unlike the other mutants, however, various mutations at P(44) led to only moderate changes in regulatory properties, but had severely impaired catalytic rates, apparent substrate affinities, and responsiveness to metabolic effectors, indicating Pro-44 or the LS is essential for optimal catalysis and activation of the AGPase heterotetramer. The catalytic importance of the LS is further supported by photoaffinity labeling studies, which revealed that the LS binds ATP at the same efficiency as the SS. These results indicate that the LS, although considered having no catalytic activity, may mimic many of the catalytic events undertaken by the SS and, thereby, influences net catalysis of the heterotetrameric enzyme.     

A low-starch barley mutant,risø 16, lacking the cytosolic small subunit of ADP-glucose pyrophosphorylase,reveals the importance of the cytosolic isoform and the identity of the plastidial small subunit

To provide information on the roles of the different forms of ADP-glucose pyrophosphorylase (AGPase) in barley (Hordeum vulgare) endosperm and the nature of the genes encoding their subunits, a mutant of barley, Ris? 16, lacking cytosolic AGPase activity in the endosperm was identified. The mutation specifically abolishes the small subunit of the cytosolic AGPase and is attributable to a large deletion within the coding region of a previously characterized small subunit gene that we have called Hv.AGP.S.1. The plastidial AGPase activity in the mutant is unaffected. This shows that the cytosolic and plastidial small subunits of AGPase are encoded by separate genes. We purified the plastidial AGPase protein and, using amino acid sequence information, we identified the novel small subunit gene that encodes this protein. Studies of the Ris? 16 mutant revealed the following. First, the reduced starch content of the mutant showed that a cytosolic AGPase is required to achieve the normal rate of starch synthesis. Second, the mutant makes both A- and B-type starch granules, showing that the cytosolic AGPase is not necessary for the synthesis of these two granule types. Third, analysis of the phylogenetic relationships between the various small subunit proteins both within and between species, suggest that the cytosolic AGPase single small subunit gene probably evolved from a leaf single small subunit gene.     

Over-expression of AGPase genes enhances seed weight and starch content in transgenic maize

Cereal crops accumulate starch in the seed endosperm as an energy reserve. ADP-glucose pyrophosphorylase (AGPase) plays a key role in regulating starch biosynthesis in cereal seeds. The AGPase in the maize endosperm is a heterotetramer of two small subunits, encoded by Brittle2 (Bt2) gene, and two large subunits, encoded by the Shrunken2 (Sh2) gene. The two genes (Bt2, Sh2) from maize were introduced into two elite maize inbred lines, solely and in tandem, and under the control of endosperm-specific promoters for over-expression. PCR, Southern blotting, and real-time RT-PCR analysis indicated that the transgenes were integrated into the genome of transgenic plants and were over-expressed in their progeny. The over-expression of either gene enhanced AGPase activity, seed weight and starch content compared with the WT, but the amounts were lower than plants with over-expression of both Bt2 and Sh2. Developing seeds from co-expression transgenic maize plants had higher cytoplasmic AGPase activity: the 100-grain weight increased 15% over the wild type (WT), and the starch content increased to over 74% compared with the WT of 65%. These results indicate that over-expression of the genes in transgenic maize plants could improve kernel traits. This report provides a feasible approach for increasing starch content and seed weight in maize.     

Analysis of allosteric effector binding sites of potato ADP-glucose pyrophosphorylase through reverse genetics

ADP-glucose pyrophosphorylase (AGPase) is a key regulatory enzyme of bacterial glycogen and plant starch synthesis as it controls carbon flux via its allosteric regulatory behavior. Unlike the bacterial enzyme that is composed of a single subunit type, the plant AGPase is a heterotetrameric enzyme (alpha2beta2) with distinct roles for each subunit type. The large subunit (LS) is involved mainly in allosteric regulation through its interaction with the catalytic small subunit (SS). The LS modulates the catalytic activity of the SS by increasing the allosteric regulatory response of the hetero-oligomeric enzyme. To identify regions of the LS involved in binding of effector molecules, a reverse genetics approach was employed. A potato (Solanum tuberosum L.) AGPase LS down-regulatory mutant (E38A) was subjected to random mutagenesis using error-prone polymerase chain reaction and screened for the capacity to form an enzyme capable of restoring glycogen production in glgC(-) Escherichia coli. Dominant mutations were identified by their capacity to restore glycogen production when the LS containing only the second site mutations was co-expressed with the wild-type SS. Sequence analysis showed that most of the mutations were decidedly nonrandom and were clustered at conserved N- and C-terminal regions. Kinetic analysis of the dominant mutant enzymes indicated that the K(m) values for cofactor and substrates were comparable with the wild-type AGPase, whereas the affinities for activator and inhibitor were altered appreciably. These AGPase variants displayed increased resistance to P(i) inhibition and/or greater sensitivity toward 3-phosphoglyceric acid activation. Further studies of Lys-197, Pro-261, and Lys-420, residues conserved in AGPase sequences, by site-directed mutagenesis suggested that the effectors 3-phosphoglyceric acid and P(i) interact at two closely located binding sites.     

The different large subunit isoforms of Arabidopsis thaliana ADP-glucose pyrophosphorylase confer distinct kinetic and regulatory properties to the heterotetrameric enzyme

ADP-glucose pyrophosphorylase catalyzes the first and limiting step in starch biosynthesis and is allosterically regulated by the levels of 3-phosphoglycerate and phosphate in plants. ADP-glucose pyrophosphorylases from plants are heterotetramers composed of two types of subunits (small and large). In this study, the six Arabidopsis thaliana genes coding for ADP-glucose pyrophosphorylase isoforms (two small and four large subunits) have been cloned and expressed in an Escherichia coli mutant deficient in ADP-glucose pyrophosphorylase activity. The co-expression of the small subunit APS1 with the different Arabidopsis large subunits (APL1, APL2, APL3, and APL4) resulted in heterotetramers with different regulatory and kinetic properties. Heterotetramers composed of APS1 and APL1 showed the highest sensitivity to the allosteric effectors as well as the highest apparent affinity for the substrates (glucose-1-phosphate and ATP), whereas heterotetramers formed by APS1 and APL2 showed the lower response to allosteric effectors and the lower affinity for the substrates. No activity was detected for the second gene coding for a small subunit isoform (APS2) annotated in the Arabidopsis genome. This lack of activity is possibly due to the absence of essential amino acids involved in catalysis and/or in the binding of glucose-1-phosphate and 3-phosphoglycerate. Kinetic and regulatory properties of the different heterotetramers, together with sequence analysis has allowed us to make a distinction between sink and source enzymes, because the combination of different large subunits would provide a high plasticity to ADP-glucose pyrophosphorylase activity and regulation. This is the first experimental data concerning the role that all the ADP-glucose pyrophosphorylase isoforms play in a single plant species. This phenomenon could have an important role in vivo, because different large subunits would confer distinct regulatory properties to ADP-glucose pyrophosphorylase according to the necessities for starch synthesis in a given tissue.     

Modulation of Allosteric Regulation by E38K and G101N Mutations in the Potato Tuber ADP-glucose Pyrophosphorylase

The higher plant ADP-glucose (ADPG) pyrophosphorylase (AGPase), composed of two small subunits and two large subunits (LSs), produces ADPG, the sole substrate for starch biosynthesis from α-D-glucose 1-phosphate and ATP. This enzyme controls a key step in starch synthesis as its catalytic activity is activated by 3-phosphoglycerate (3-PGA) and inhibited by orthophosphate (Pi). Previously, two mutations in the LS of potato AGPase (PLS), PLS-E38K and PLS-G101N, were found to increase sensitivity to 3-PGA activation and tolerance to Pi inhibition. In the present study, the double mutated enzyme (PLS-E38K/G101N) was evaluated. In a complementation assay of ADPG synthesis in an Escherichia coli mutant defective in the synthesis of ADPG, expression of PLS-E38K/G101N mediated higher glycogen production than wild-type potato AGPase (PLS-WT) and the single mutant enzymes, PLS-E38K and PLS-G101N, individually. Purified PLS-E38K/G101N showed higher sensitivity to 3-PGA activation and tolerance to Pi inhibition than PLS-E38K or PLS-G101N. Moreover, the enzyme activities of PLS-E38K, PLS-G101N, and PLS-E38K/G101N were more readily stimulated by other major phosphate-ester metabolites, such as fructose 6-phosphate, fructose 2,6-bisphosphate, and ribose 5-phosphate, than was that of PLS-WT. Hence, although the specific enzyme activities of the LS mutants toward 3-PGA were impaired to some extent by the mutations, our results suggest that their enhanced allosteric regulatory properties and the broadened effector selectivity gained by the same mutations not only offset the lowered enzyme catalytic turnover rates but also increase the net performance of potato AGPase in vivo in view of increased glycogen production in bacterial cells.     

Starchless mutants of Chlamydomonas reinhardtii lack the small subunit of a heterotetrameric ADP-glucose pyrophosphorylase

ADP-glucose synthesis through ADP-glucose pyrophosphorylase defines the major rate-controlling step of storage polysaccharide synthesis in both bacteria and plants. We have isolated mutant strains defective in the STA6 locus of the monocellular green alga Chlamydomonas reinhardtii that fail to accumulate starch and lack ADP-glucose pyrophosphorylase activity. We show that this locus encodes a 514-amino-acid polypeptide corresponding to a mature 50-kDa protein with homology to vascular plant ADP-glucose pyrophosphorylase small-subunit sequences. This gene segregates independently from the previously characterized STA1 locus that encodes the large 53-kDa subunit of the same heterotetramer enzyme. Because STA1 locus mutants have retained an AGPase but exhibit lower sensitivity to 3-phosphoglyceric acid activation, we suggest that the small and large subunits of the enzyme define, respectively, the catalytic and regulatory subunits of AGPase in unicellular green algae. We provide preliminary evidence that both the small-subunit mRNA abundance and enzyme activity, and therefore also starch metabolism, may be controlled by the circadian clock.     

The two AGPase subunits evolve at different rates in angiosperms, yet they are equally sensitive to activity-altering amino acid changes when expressed in bacteria

The rate of protein evolution is generally thought to reflect, at least in part, the proportion of amino acids within the protein that are needed for proper function. In the case of ADP-glucose pyrophosphorylase (AGPase), this premise led to the hypothesis that, because the AGPase small subunit is more conserved compared with the large subunit, a higher proportion of the amino acids of the small subunit are required for enzyme activity compared with the large subunit. Evolutionary analysis indicates that the AGPase small subunit has been subject to more intense purifying selection than the large subunit in the angiosperms. However, random mutagenesis and expression of the maize (Zea mays) endosperm AGPase in bacteria show that the two AGPase subunits are equally predisposed to enzyme activity-altering amino acid changes when expressed in one environment with a single complementary subunit. As an alternative hypothesis, we suggest that the small subunit exhibits more evolutionary constraints in planta than does the large subunit because it is less tissue specific and thus must form functional enzyme complexes with different large subunits. Independent approaches provide data consistent with this alternative hypothesis.     

Enhancement of starch content by constitutive expression of <Emphasis Type="Italic">GmTrxF</Emphasis> in transgenic <Emphasis Type="Italic">Arabidopsis</Emphasis>

The plastidic thioredoxin F-type (TrxF) protein plays an important role in plant carbohydrate metabolism biosynthesis. In this study, a gene encoding the TrxF protein, named GmTrxF, was isolated from soybean. The open reading frame (ORF) contained 540 nucleotides encoding 179 amino acids. The coding region of GmTrxF was cloned into a binary vector under the control of 35S promoter and then transformed into Arabidopsis. The starch content in GmTrxF expressing plants was increased by 57–109% compared to that in wild-type (WT). Real-time quantitative PCR (qRT-PCR) analysis showed that constitutive expression of GmTrxF up-regulated the expression of phosphoglucomutase (AtPGM), ADP-glucose pyrophosphorylase (AGPase) small subunit (AtAGPase-S1 and AtAGPase-S2), AGPase large subunit (AtAGPase-L1 and AtAGPase-L2) and soluble starch synthases (AtSSS I, AtSSS II, AtSSS III and AtSSS IV) genes involved in starch biosynthesis in the transgenic Arabidopsis plants. Meanwhile, enzymatic analyses showed that the major enzymes (AGPase and SSS) involved in the starch biosynthesis exhibited higher activities in the transgenic plants compared to WT. These results suggest that GmTrxF may improve starch content of Arabidopsis by up-regulating the expression of the related genes and increasing the activities of the major enzymes invovled in starch biosynthesis. The manipulation of GmTrxF expression might be used for increasing starch accumulation of plants in the future.     

Starch biosynthesis in cereal endosperm

Stored starch generally consists of two d-glucose homopolymers, the linear polymer amylose and a highly branched glucan amylopectin that connects linear chains. Amylopectin structurally contributes to the crystalline organization of the starch granule in cereals. In the endosperm, amylopectin biosynthesis requires the proper execution of a coordinated series of enzymatic reactions involving ADP glucose pyrophosphorylase (AGPase), soluble starch synthase (SS), starch branching enzyme (BE), and starch debranching enzyme (DBE), whereas amylose is synthesized by AGPase and granule-bound starch synthase (GBSS). It is highly possible that plastidial starch phosphorylase (Pho1) plays an important role in the formation of primers for starch biosynthesis in the endosperm. Recent advances in our understanding of the functions of individual enzyme isoforms have provided new insights into how linear polymer chains and branch linkages are synthesized in cereals. In particular, genetic analyses of a suite of mutants have formed the basis of a new model outlining the role of various enzyme isoforms in cereal starch production. In our current review, we summarize the recent research findings related to starch biosynthesis in cereal endosperm, with a particular focus on rice.     

Characterization of an autonomously activated plant ADP-glucose pyrophosphorylase

ADP-glucose pyrophosphorylase (AGPase) catalyzes the rate-limiting step in starch biosynthesis in plants and changes in its catalytic and/or allosteric properties can lead to increased starch production. Recently, a maize (Zea mays)/potato (Solanum tuberosum) small subunit mosaic, MP [Mos(1–198)], containing the first 198 amino acids of the small subunit of the maize endosperm enzyme and the last 277 amino acids from the potato tuber enzyme, was expressed with the maize endosperm large subunit and was reported to have favorable kinetic and allosteric properties. Here, we show that this mosaic, in the absence of activator, performs like a wild-type AGPase that is partially activated with 3-phosphoglyceric acid (3-PGA). In the presence of 3-PGA, enzyme properties of Mos(1–198)/SH2 are quite similar to those of the wild-type maize enzyme. In the absence of 3-PGA, however, the mosaic enzyme exhibits greater activity, higher affinity for the substrates, and partial inactivation by inorganic phosphate. The Mos(1–198)/SH2 enzyme is also more stable to heat inactivation. The different properties of this protein were mapped using various mosaics containing smaller portions of the potato small subunit. Enhanced heat stability of Mos(1–198) was shown to originate from five potato-derived amino acids between 322 and 377. These amino acids were shown previously to be important in small subunit/large subunit interactions. These five potato-derived amino acids plus other potato-derived amino acids distributed throughout the carboxyl-terminal portion of the protein are required for the enhanced catalytic and allosteric properties exhibited by Mos(1–198)/SH2.     

Insights into subunit interactions in the heterotetrameric structure of potato ADP-glucose pyrophosphorylase

ADP-glucose pyrophosphorylase, a key allosteric enzyme involved in higher plant starch biosynthesis, is composed of pairs of large (LS) and small subunits (SS). Current evidence indicates that the two subunit types play distinct roles in enzyme function. The LS is involved in mainly allosteric regulation through its interaction with the catalytic SS. Recently the crystal structure of the SS homotetramer has been solved, but no crystal structure of the native heterotetrameric enzyme is currently available. In this study, we first modeled the three-dimensional structure of the LS to construct the heterotetrameric enzyme. Because the enzyme has a 2-fold symmetry, six different dimeric (either up-down or side-by-side) interactions were possible. Molecular dynamics simulations were carried out for each of these possible dimers. Trajectories obtained from molecular dynamics simulations of each dimer were then analyzed by the molecular mechanics/Poisson-Boltzmann surface area method to identify the most favorable dimers, one for up-down and the other for side-by-side. Computational results combined with site directed mutagenesis and yeast two hybrid experiments suggested that the most favorable heterotetramer is formed by LS-SS (side-by-side), and LS-SS (up-down). We further determined the order of assembly during the heterotetrameric structure formation. First, side-by-side LS-SS dimers form followed by the up-down tetramerization based on the relative binding free energies.