Association between molecular variation in the phosphoglucose isomerase (Pgi) locus and migration propensity in the Glanville fritillary butterfly

Habitat fragmentation produces patches of suitable habitat surrounded by unfavourable matrix habitat. A species may persist in such a fragmented landscape in an equilibrium between the extinctions and recolonizations of local populations, thus forming a metapopulation. Migration between local populations is necessary for the long-term persistence of a metapopulation. The Glanville fritillary butterfly (Melitaea cinxia) forms a metapopulation in the Åland islands in Finland. There is migration between the populations, the extent of which is affected by several environmental factors and variation in the phenotype of individual butterflies. Different allelic forms of the glycolytic enzyme phosphoglucose isomerase (Pgi) has been identified as a possible genetic factor influencing flight performance and migration rate in this species. The frequency of a certain Pgi allele, Pgi-f, follows the same pattern in relation to population age and connectivity as migration propensity. Furthermore, variation in flight metabolic performance, which is likely to affect migration propensity, has been linked to genetic variation in Pgi or a closely linked locus. The aim of this study was to investigate the association between Pgi genotype and the migration propensity in the Glanville fritillary both at the individual and population levels using a statistical modelling approach. A mark-release-recapture (MRR) study was conducted in a habitat patch network of M. cinxia in Åland to collect data on the movements of individual butterflies. Larval samples from the study area were also collected for population level examinations. Each butterfly and larva was genotyped at the Pgi locus. The MRR data was parameterised with two mathematical models of migration: the Virtual Migration Model (VM) and the spatially explicit diffusion model. VM model predicted and observed numbers of emigrants from populations with high and low frequencies of Pgi-f were compared. Posterior predictive data sets were simulated based on the parameters of the diffusion model. Lack-of-fit of observed values to the model predicted values of several descriptors of movements were detected, and the effect of Pgi genotype on the deviations was assessed by randomizations including the genotype information. This study revealed a possible difference in the effect of Pgi genotype on migration propensity between the two sexes in the Glanville fritillary. The females with and males without the Pgi-f allele moved more between habitat patches, which is probably related to differences in the function of flight in the two sexes. Females may use their high flight capacity to migrate between habitat patches to find suitable oviposition sites, whereas males may use it to acquire mates by keeping a territory and fighting off other intruding males, possibly causing them to emigrate. The results were consistent across different movement descriptors and at the individual and population levels. The effect of Pgi is likely to be dependent on the structure of the landscape and the prevailing environmental conditions.

Unfolding of Plasmodium falciparum triosephosphate isomerase in urea and guanidinium chloride solutions. Evidence for a novel disulfide exchange reaction in a covalently crosslinked mutant

The conformational stability of Plasmodium falciparum triosephosphate isomerase (TIMWT) enzyme has been investigated in urea and guanidinium chloride (GdmCl) solutions using circular dichroism, fluorescence, and size-exclusion chromatography. The dimeric enzyme is remarkably stable in urea solutions. It retains considerable secondary, tertiary, and quaternary structure even in 8 M urea. In contrast, the unfolding transition is complete by 2.4 M GdmCl. Although the secondary as well as the tertiary interactions melt before the perturbation of the quaternary structure, these studies imply that the dissociation of the dimer into monomers ultimately leads to the collapse of the structure, suggesting that the interfacial interactions play a major role in determining multimeric protein stability. The Cm(urea)/Cm(GdmCl) ratio (where Cm is the concentration of the denaturant required at the transition midpoint) is unusually high for triosephosphate isomerase as compared to other monomeric and dimeric proteins. A disulfide cross-linked mutant protein (Y74C) engineered to form two disulfide cross-links across the interface (13-74‘) and (13‘-74) is dramatically destablized in urea. The unfolding transition is complete by 6 M urea and involves a novel mechanism of dimer dissociation through intramolecular thiol−disulfide exchange.

Engineered dimer interface mutants of triosephosphate isomerase: the role of inter-subunit interactions in enzyme function and stability

The role of inter-subunit interactions in maintaining optimal catalytic activity in triosephosphate isomerase (TIM) has been probed, using the Plasmodium falciparum enzyme as a model. Examination of subunit interface contacts in the crystal structures suggests that residue 75 (Thr, conserved) and residue 13 (Cys, variable) make the largest number of inter-subunit contacts. The mutants Cys13Asp (C13D) and Cys13Glu (C13E) have been constructed and display significant reduction in catalytic activity when compared with wild-type (WT) enzyme (similar to 7.4-fold decrease in k(cat) for the C13D and similar to 3.3-fold for the C13E mutants). Analytical gel filtration demonstrates that the C13D mutant dissociates at concentrations < 1.25 mu M, whereas the WT and the C13E enzymes retain the dimeric structure. The order of stability of the mutants in the presence of chemical denaturants, like urea and guanidium chloride, is WT > Cys13Glu > Cys13Asp. Irreversible thermal precipitation temperatures follow the same order as well. Modeling studies establish that the Cys13Asp mutation is likely to cause a significantly greater structural perturbation than Cys13Glu. Analysis of sequence and structural data for TIMs from diverse sources suggests that residues 13 and 82 form a pair of proximal sites, in which a limited number of residue pairs may be accommodated.

Probing the role of the fully conserved Cys126 in triosephosphate isomerase by site-specific mutagenesis - distal effects on dimer stability

Cys126 is a completely conserved residue in triosephosphate isomerase that is proximal to the active site but has been ascribed no specific role in catalysis. A previous study of the C126S and C126A mutants of yeast TIM reported substantial catalytic activity for the mutant enzymes, leading to the suggestion that this residue is implicated in folding and stability [Gonzalez-Mondragon E et al. (2004) Biochemistry43, 3255–3263]. We re-examined the role of Cys126 with the Plasmodium falciparum enzyme as a model. Five mutants, C126S, C126A, C126V, C126M, and C126T, were characterized. Crystal structures of the 3-phosphoglycolate-bound C126S mutant and the unliganded forms of the C126S and C126A mutants were determined at a resolution of 1.7–2.1 Å. Kinetic studies revealed an approximately five-fold drop in kcat for the C126S and C126A mutants, whereas an approximately 10-fold drop was observed for the other three mutants. At ambient temperature, the wild-type enzyme and all five mutants showed no concentration dependence of activity. At higher temperatures (> 40 °C), the mutants showed a significant concentration dependence, with a dramatic loss in activity below 15 μm. The mutants also had diminished thermal stability at low concentration, as monitored by far-UV CD. These results suggest that Cys126 contributes to the stability of the dimer interface through a network of interactions involving His95, Glu97, and Arg98, which form direct contacts across the dimer interface.

Probing the role of the fully conserved Cys126 in triosephosphate isomerase by site-specific mutagenesis - distal effects on dimer stability

DatabaseStructural data are available in the Protein Data Bank under the accession numbers 3PVF, 3PY2, and 3PWA. Structured digital abstract Tim binds to Tim by x-ray crystallography (View interaction).

Revisiting the Mechanism of the Triosephosphate Isomerase Reaction: The Role of the Fully Conserved Glutamic Acid 97 Residue

An analysis of 503 available triosephosphate isomerase sequences revealed nine fully conserved residues. Of these, four residues-K12, H95, E97 and E165-are capable of proton transfer and are all arrayed around the dihydroxyacetone phosphate substrate in the three-dimensional structure. Specific roles have been assigned to the residues K12, H95 and E165, but the nature of the involvement of E97 has not been established. Kinetic and structural characterization is reported for the E97Q and E97D mutants of Plasmodium falciparum triosephosphate isomerase (Pf TIM). A 4000-fold reduction in k(cat) is observed for E97Q, whereas the E97D mutant shows a 100-fold reduction. The control mutant, E165A, which lacks the key catalytic base, shows an approximately 9000-fold drop in activity. The integrity of the overall fold and stability of the dimeric structure have been demonstrated by biophysical studies. Crystal structures of E97Q and E97D mutants have been determined at 2.0 angstrom resolution. In the case of the isosteric replacement of glutamic acid by glutamine in the E97Q mutant a large conformational change for the critical K12 side chain is observed, corresponding to a trans-to-gauche transition about the C gamma-C delta (chi(3)) bond. In the E97D mutant, the K12 side chain maintains the wild-type orientation, but the hydrogen bond between K12 and D97 is lost. The results are interpreted as a direct role for E97 in the catalytic proton transfer cycle. The proposed mechanism eliminates the need to invoke the formation of the energetically unfavourable imidazolate anion at H95, a key feature of the classical mechanism.

The Structure of the Thioredoxin-Triosephosphate Isomerase Complex Provides Insights into the Reversible Glutathione-Mediated Regulation of Triosephosphate Isomerase

Protein-protein interactions are crucial for many biological functions. The redox interactome encompasses numerous weak transient interactions in which thioredoxin plays a central role. Proteomic studies have shown that thioredoxin binds to numerous proteins belonging to various cellular processes, including energy metabolism. Thioredoxin has cross talk with other redox mechanisms involving glutathionylation and has functional overlap with glutaredoxin in deglutathionylation reactions. In this study, we have explored the structural and biochemical interactions of thioredoxin with the glycolytic enzyme, triosephosphate isomerase. Nuclear magnetic resonance chemical shift mapping methods and molecular dynamics-based docking have been applied in deriving a structural model of the thioredoxin-triosephosphate isomerase complex. The spatial proximity of active site cysteine residues of thioredoxin to reactive thiol groups on triosephosphate isomerase provides a direct link to the observed deglutathionylation of cysteine 217 in triosephosphate isomerase, thereby reversing the inhibitory effect of S-glutathionylation of triosephosphate isomerase.

The Structure of the Thioredoxin−Triosephosphate Isomerase Complex Provides Insights into the Reversible Glutathione-Mediate Regulation of Triosephosphate Isomerase

Protein−protein interactions are crucial for many biological functions. The redox interactome encompasses numerous weak transient interactions in which thioredoxin plays a central role. Proteomic studies have shown that thioredoxin binds to numerous proteins belonging to various cellular processes, including energy metabolism. Thioredoxin has cross talk with other redox mechanisms involving glutathionylation and has functional overlap with glutaredoxin in deglutathionylation reactions. In this study, we have explored the structural and biochemical interactions of thioredoxin with the glycolytic enzyme, triosephosphate isomerase. Nuclear magnetic resonance chemical shift mapping methods and molecular dynamics-based docking have been applied in deriving a structural model of the thioredoxin−triosephosphate isomerase complex. The spatial proximity of active site cysteine residues of thioredoxin to reactive thiol groups on triosephosphate isomerase provides a direct link to the observed deglutathionylation of cysteine 217 in triosephosphate isomerase, thereby reversing the inhibitory effect of S-glutathionylation of triosephosphate isomerase.

Structural insights into N-terminal to C-terminal interactions and implications for thermostability of a (beta/alpha)(8)-triosephosphate isomerase barrel enzyme

Although several factors have been suggested to contribute to thermostability, the stabilization strategies used by proteins are still enigmatic. Studies on a recombinant xylanase from Bacilllus sp. NG-27 (RBSX), which has the ubiquitous (beta/alpha)(8)-triosephosphate isomerase barrel fold, showed that just a single mutation, V1L, although not located in any secondary structural element, markedly enhanced the stability from 70 degrees C to 75 degrees C without loss of catalytic activity. Conversely, the V1A mutation at the same position decreased the stability of the enzyme from 70 degrees C to 68 degrees C. To gain structural insights into how a single extreme N-terminus mutation can markedly influence the thermostability of the enzyme, we determined the crystal structure of RBSX and the two mutants. On the basis of computational analysis of their crystal structures, including residue interaction networks, we established a link between N-terminal to C-terminal contacts and RBSX thermostability. Our study reveals that augmenting N-terminal to C-terminal noncovalent interactions is associated with enhancement of the stability of the enzyme. In addition, we discuss several lines of evidence supporting a connection between N-terminal to C-terminal noncovalent interactions and protein stability in different proteins. We propose that the strategy of mutations at the termini could be exploited with a view to modulate stability without compromising enzymatic activity, or in general, protein function in diverse folds where N and C termini are in close proximity. Database The coordinates of RBSX, V1A and V1L have been deposited in the PDB database under the accession numbers 4QCE, 4QCF, and 4QDM, respectively

Probing the role of highly conserved residues in triosephosphate isomerase-analysis of site specific mutants at positions 64 and 75 in the Plasmodial enzyme

Highly conserved residues in enzymes are often found to be clustered close to active sites, suggesting that functional constraints dictate the nature of amino acid residues accommodated at these sites. Using the Plasmodiumfalciparum triosephosphate isomerase (PfTIM) enzyme () as a template, we have examined the effects of mutations at positions 64 and 75, which are not directly involved in the proton transfer cycle. Thr (T) occurring at position 75 is completely conserved, whereas only Gln (Q) and Glu (E) are accommodated at position 64. Biophysical and kinetic data are reported for four T75 (T75S/V/C/N) and two Q64 (Q64N/E) mutants. The dimeric structure is weakened in the Q64E and Q64N mutants, whereas dimer integrity is unimpaired in all four T75 mutants. Measurement of the concentration dependence of enzyme activity permits an estimate of K-d values for dimer dissociation (Q64N=73.79.2nm and Q64E=44.6 +/- 8.4nm). The T75S/V/C mutants have activities comparable to the wild-type enzyme, whereas a fourfold drop is observed for T75N. All four T75 mutants show a dramatic fall in activity between 35 degrees C and 45 degrees C. Crystal structure determination of the T75S/V/N mutants provides insights into the variations in local interactions, with the T75N mutant showing the largest changes. Hydrogen-bond interactions determine dimer stability restricting the choice of residues at position 64 to Gln (Q) and Glu (E). At position 75, the overwhelming preference for Thr (T) may be dictated by the imperative of maintaining temperature stability of enzyme activity.

Connecting Active-Site Loop Conformations and Catalysis in Triosephosphate Isomerase: Insights from a Rare Variation at Residue96 in the Plasmodial Enzyme

Despite extensive research into triosephosphate isomerases (TIMs), there exists a gap in understanding of the remarkable conjunction between catalytic loop-6 (residues 166-176) movement and the conformational flip of Glu165 (catalytic base) upon substrate binding that primes the active site for efficient catalysis. The overwhelming occurrence of serine at position96 (98% of the 6277 unique TIM sequences), spatially proximal to E165 and the loop-6 residues, raises questions about its role in catalysis. Notably, Plasmodium falciparum TIM has an extremely rare residuephenylalanineat this position whereas, curiously, the mutant F96S was catalytically defective. We have obtained insights into the influence of residue96 on the loop-6 conformational flip and E165 positioning by combining kinetic and structural studies on the PfTIM F96 mutants F96Y, F96A, F96S/S73A, and F96S/L167V with sequence conservation analysis and comparative analysis of the available apo and holo structures of the enzyme from diverse organisms.

Mass spectrometric analysis of dimer-disrupting mutations in Plasmodium triosephosphate isomerase

Electrospray ionization mass spectrometry (ESI MS) under nanospray conditions has been used to examine the effects of mutation at two key dimer interface residues, Gln (Q) 64 and Thr (T) 75, in Plasmodium falciparum triosephosphate isomerase. Both residues participate in an intricate network of intra- and intersubunit hydrogen bonds. The gas phase distributions of dimeric and monomeric protein species have been examined for the wild type enzyme (TWT) and three mutants, Q64N, Q64E, and 175S, under a wide range of collision energies (40-160 eV). The results established the order of dimer stability as TWT > T75S > Q64E similar to Q64N. The mutational effects on dimer stability are in good agreement with the previously reported estimates, based on the concentration dependence of enzyme activity. Additional experiments in solution, using inhibition of activity by a synthetic dimer interface peptide, further support the broad agreement between gas phase and solution studies. (C) 2016 Elsevier Inc. All rights reserved.

Plastidic Phosphoglucose Isomerase Is an Important Determinant of Starch Accumulation in Mesophyll Cells, Growth, Photosynthetic Capacity, and Biosynthesis of Plastidic Cytokinins in Arabidopsis

Phosphoglucose isomerase (PGI) catalyzes the reversible isomerization of glucose-6-phosphate and fructose-6-phosphate. It is involved in glycolysis and in the regeneration of glucose-6-P molecules in the oxidative pentose phosphate pathway (OPPP). In chloroplasts of illuminated mesophyll cells PGI also connects the Calvin-Benson cycle with the starch biosynthetic pathway. In this work we isolated pgi1-3, a mutant totally lacking pPGI activity as a consequence of aberrant intron splicing of the pPGI encoding gene, PGI1. Starch content in pgi1-3 source leaves was ca. 10-15% of that of wild type (WT) leaves, which was similar to that of leaves of pgi1-2, a T-DNA insertion pPGI null mutant. Starch deficiency of pgi1 leaves could be reverted by the introduction of a sex1 null mutation impeding beta-amylolytic starch breakdown. Although previous studies showed that starch granules of pgi1-2 leaves are restricted to both bundle sheath cells adjacent to the mesophyll and stomata guard cells, microscopy analyses carried out in this work revealed the presence of starch granules in the chloroplasts of pgi1-2 and pgi1-3 mesophyll cells. RT-PCR analyses showed high expression levels of plastidic and extra-plastidic beta-amylase encoding genes in pgi1 leaves, which was accompanied by increased beta-amylase activity. Both pgi1-2 and pgi1-3 mutants displayed slow growth and reduced photosynthetic capacity phenotypes even under continuous light conditions. Metabolic analyses revealed that the adenylate energy charge and the NAD(P) H/NAD(P) ratios in pgi1 leaves were lower than those of WT leaves. These analyses also revealed that the content of plastidic 2-C-methyl-D-erythritol 4-phosphate (MEP)-pathway derived cytokinins (CKs) in pgi1 leaves were exceedingly lower than in WT leaves. Noteworthy, exogenous application of CKs largely reverted the low starch content phenotype of pgi1 leaves. The overall data show that pPGI is an important determinant of photosynthesis, energy status, growth and starch accumulation in mesophyll cells likely as a consequence of its involvement in the production of OPPP/glycolysis intermediates necessary for the synthesis of plastidic MEP-pathway derived hormones such as CKs.