Involvement of Val 315 located in the C-terminal region of thermolysin in its expression in Escherichia coli and its thermal stability

Thermolysin is a thermophilic and halophilic zinc metalloproteinase that consists of β-rich N-terminal (residues 1–157) and α-rich C-terminal (residues 158–316) domains. Expression of thermolysin variants truncated from the C-terminus was examined in E. coli culture. The C-terminal Lys316 residue was not significant in the expression, but Val315 was critical. Variants in which Val315 was substituted with fourteen amino acids were prepared. The variants substituted with hydrophobic amino acids such as Leu and Ile were almost the same as wild-type thermolysin (WT) in the expression amount, α-helix content, and stability. Variants with charged (Asp, Glu, Lys, and Arg), bulky (Trp), or small (Gly) amino acids were lower in these characteristics than WT. All variants exhibited considerably high activities (50–100% of WT) in hydrolyzing protein and peptide substrates. The expression amount, helix content, and stability of variants showed good correlation with hydropathy indexes of the amino acids substituted for Val315. Crystallographic study of thermolysin has indicated that V315 is a member of the C-terminal hydrophobic cluster. The results obtained in the present study indicate that stabilization of the cluster increases thermolysin stability and that the variants with higher stability are expressed more in the culture. Although thermolysin activity was not severely affected by the variation at position 315, the stability and specificity were modified significantly, suggesting the long-range interaction between the C-terminal region and active site.     

Engineering the substrate specificity of Escherichia coli asparaginase. II. Selective reduction of glutaminase activity by amino acid replacements at position 248

The use of Escherichia coli asparaginase II as a drug for the treatment of acute lymphoblastic leukemia is complicated by the significant glutaminase side activity of the enzyme. To develop enzyme forms with reduced glutaminase activity, a number of variants with amino acid replacements in the vicinity of the substrate binding site were constructed and assayed for their kinetic and stability properties. We found that replacements of Asp248 affected glutamine turnover much more strongly than asparagine hydrolysis. In the wild-type enzyme, N248 modulates substrate binding to a neighboring subunit by hydrogen bonding to side chains that directly interact with the substrate. In variant N248A, the loss of transition state stabilization caused by the mutation was 15 kJ mol(-1) for L-glutamine compared to 4 kJ mol(-1) for L-aspartic beta-hydroxamate and 7 kJ mol(-1) for L-asparagine. Smaller differences were seen with other N248 variants. Modeling studies suggested that the selective reduction of glutaminase activity is the result of small conformational changes that affect active-site residues and catalytically relevant water molecules.     

Evaluation at atomic resolution of the role of strain in destabilizing the temperature‐sensitive T4 lysozyme mutant Arg 96 → His

Mutant R96H is a classic temperature‐sensitive mutant of bacteriophage T4 lysozyme. It was in fact the first variant of the protein to be characterized structurally. Subsequently, it has been studied extensively by a variety of experimental and computational techniques, but the reasons for the loss of stability of the mutant protein remain controversial. In the crystallographic refinement of the mutant structure at 1.9 Å resolution one of the bond angles at the site of substitution appeared to be distorted by about 11^°, and it was suggested that this steric strain was one of the major factors in destabilizing the mutant. Different computationally‐derived models of the mutant structure, however, did not show such distortion. To determine the geometry at the site of mutation more reliably, we have extended the resolution of the data and refined the wildtype (WT) and mutant structures to be better than 1.1 Å resolution. The high‐resolution refinement of the structure of R96H does not support the bond angle distortion seen in the 1.9 Å structure determination. At the same time, it does confirm other manifestations of strain seen previously including an unusual rotameric state for His96 with distorted hydrogen bonding. The rotamer strain has been estimated as about 0.8 kcal/mol, which is about 25% of the overall reduction in stability of the mutant. Because of concern that contacts from a neighboring molecule in the crystal might influence the geometry at the site of mutation we also constructed and analyzed supplemental mutant structures in which this crystal contact was eliminated. High‐resolution refinement shows that the crystal contacts have essentially no effect on the conformation of Arg96 in WT or on His96 in the R96H mutant.     

Evidence of the coevolution of antigenicity and host cell tropism of foot-and-mouth disease virus in vivo

In this work we analyze the antigenic properties and the stability in cell culture of virus mutants recovered upon challenge of peptide-vaccinated cattle with foot-and-mouth disease virus (FMDV) C3 Arg85. Previously, we showed that a significant proportion of 29 lesions analyzed (41%) contained viruses with single amino acid replacements (R141G, L144P, or L147P) within a major antigenic site located at the G-H loop of VP1, known to participate also in interactions with integrin receptors. Here we document that no replacements at this site were found in viruses from 12 lesions developed in six control animals upon challenge with FMDV C3 Arg85. Sera from unprotected, vaccinated animals exhibited poor neutralization titers against mutants recovered from them. Sequence analyses of the viruses recovered upon 10 serial passages in BHK-21 and FBK-2 cells in the presence of preimmune (nonneutralizing) sera revealed that mutants reverted to the parental sequence, suggesting an effect of the amino acid replacements in the interaction of the viruses with cells. Parallel passages in the presence of subneutralizing concentrations of immune homologous sera resulted in the maintenance of mutations R141G and L147P, while mutation L144P reverted to the C3 Arg85 sequence. Reactivity with a panel of FMDV type C-specific monoclonal antibodies indicated that mutant viruses showed altered antigenicity. These results suggest that the selective pressure exerted by host humoral immune response can play a role in both the selection and stability of antigenic FMDV variants and that such variants can manifest alterations in cell tropism.     

Maturation efficiency, trypsin sensitivity, and optical properties of Arg96, Glu222, and Gly67 variants of green fluorescent protein

Spontaneous chromophore biosynthesis in green fluorescent protein (GFP) is initiated by a main-chain cyclization reaction catalyzed by the protein fold. To investigate the structural prerequisites for chromophore formation, we have substituted the conserved residues Arg96, Glu222, and Gly67. Upon purification, the variants can be ordered based on their decreasing extent of chromophore maturation according to the series EGFP, E222Q, R96K, G67A, and R96M. Arg96 and Glu222 appear to play catalytic roles, whereas Gly67 is likely important in interior packing to enforce correct hydrogen bonding to Arg96. The effect of Arg96 can be partially compensated for by a lysine, but not by a methionine residue, confirming its electrophilic role. Limited trypsinolysis data suggest that protein stability is largely unaffected by the presence of the chromophore, inconsistent with the mechanical compression hypothesis. Trends in optical properties may be related to the degree of chromophore charge delocalization, which is modulated by residue 96.     

Characterization of an L-phosphinothricin resistant glutamine synthetase from Exiguobacterium sp. and its improvement

A glutamine synthetase (GS; 1341 bp) gene with potent L-phosphinothricin (PPT) resistance was isolated and characterized from a marine bacterium Exiguobacterium sp. Molecular docking analysis indicated that the substitution of residues Glu60 and Arg64 may lead to significant changes in binding pocket. To enhance the enzymatic property of GS, variants E60A and R64G were obtained by site-directed mutagenesis. The results revealed a noteworthy change in the thermostability and activity in comparison to the wild type (WT). WT exhibited optimum activity at 35 °C, while E60A and R64G exhibited optimum activity at 45 and 40 °C, respectively. The mutant R64G was 4.3 times more stable at 70 °C in comparison to WT, while E60A was 5.7 times more stable. Kinetic analysis revealed that the k _cat value of R64G mutant was 8.10-, 7.25- and 7.63-fold that of WT for ADP, glutamine and hydroxylamine, respectively. The kinetic inhibition (K _i, 4.91 ± 0.42 mM) of R64G was 2.02-fold that of WT (2.43 ± 0.14 mM) for L-phosphinothricin. The analysis of structure and function relationship showed that the binding pocket underwent dramatic changes when Arg site of 64 was substituted by Gly, thus promoting the rapid capture of substrates and leading to increase in activity and PPT-resistance of mutant R64G. The rearrangements of the residues at the molecular level formed new hydrogen bonds around the active site, which contributed to the increase of thermostability of enzymes. This study provides new insights into substrate binding mechanism of glutamine synthetase and the improved GS gene also has a potential for application in transgenic crops with L-phosphinothricin tolerance.

     

Functional analysis of three amino acid residues of purR repressor, Trpl47, Gln-218 and Gln-292 inSalmonella typhimurium

The amber mutation sites of 6 purR(am) mutants were determined by cloning and DNA sequencing. The results showed that the mutations were distributed at three different sites in PurR coding region, G⁷²¹(→A), C⁹³³(→T) and C¹¹⁵⁵(→T), which respectively turn Trp-147,Gln-218 and Gln-292 of PurR into TAG terminal codon. To determine the effect of the three amino acid residues on regulatory function of PurR protein 5 different kinds of tRNA suppressor genes, Su3, Su4, Su6, Su7 and Su9 were used for creating the PurR protein variants with single amino acid substitution. The results indicated that Cys, Glu, Gly, His and Arg which substituted Trp-147 respectively all could not recover the regulation function of PurR. It confirmed that Trp-147 is a critical amino acid for the PurR function. Gln-292 substituted respectively by the same amino acids also could not recover the PurR function, demonstrating that Gln-292 is also an important amino acid residue in PurR.     

Functional analysis of three amino acid residues of purR repressor, Trpl47, Gln-218 and Gln-292 in Salmonella typhimurium

The amber mutation sites of 6 purR(am) mutants were determined by cloning and DNA sequencing. The results showed that the mutations were distributed at three different sites in PurR coding region, G⁷²¹(→A), C⁹³³(→T) and C¹¹⁵⁵(→T), which respectively turn Trp-147, Gln-218 and Gln-292 of PurR into TAG terminal codon. To determine the effect of the three amino acid residues on regulatory function of PurR protein 5 different kinds of tRNA suppressor genes, Su3, Su4, Su6, Su7 and Su9 were used for creating the PurR protein variants with single amino acid substitution. The results indicated that Cys, Glu, Gly, His and Arg which substituted Trp-147 respectively all could not recover the regulation function of PurR. It confirmed that Trp-147 is a critical amino acid for the PurR function. Gln-292 substituted respectively by the same amino acids also could not recover the PurR function, demonstrating that Gln-292 is also an important amino acid residue in PurR.     

Sequences flanking Arg336 in factor VIIIa modulate factor Xa-catalyzed cleavage rates at this site and cofactor function

Factor (F)VIII can be activated to FVIIIa by FXa following cleavages at Arg(372), Arg(740), and Arg(1689). FXa also cleaves FVIII/FVIIIa at Arg(336) and Arg(562) resulting in inactivation of the cofactor. These inactivating cleavages occur on a slower time scale than the activating ones. We assessed the contributions to cleavage rate and cofactor function of residues flanking Arg(336), the primary site yielding FVIII(a) inactivation, following replacement of these residues with those flanking the faster-reacting Arg(740) and Arg(372) sites and the slower-reacting Arg(562) site. Replacing P4-P3' residues flanking Arg(336) with those from Arg(372) or Arg(740) resulted in ～4-6-fold increases in rates of FXa-catalyzed inactivation of FVIIIa, which paralleled the rates of proteolysis at Arg(336). Examination of partial sequence replacements showed a predominant contribution of prime residues flanking the scissile bonds to the enhanced rates. Conversely, replacement of this sequence with residues flanking the slow-reacting Arg(562) site yielded inactivation and cleavage rates that were ～40% that of the WT values. The capacity for FXa to activate FVIII variants where cleavage at Arg(336) was accelerated due to flanking sequence replacement showed marked reductions in peak activity, whereas reducing the cleavage rate at this site enhanced peak activity. Furthermore, plasma-based thrombin generation assays employing the variants revealed significant reductions in multiple parameter values with acceleration of Arg(336) cleavage suggesting increased down-regulation of FXase. Overall, these results are consistent with a model of competition for activating and inactivating cleavages catalyzed by FXa that is modulated in large part by sequences flanking the scissile bonds.     

Substrate specificity at the P1' site of Escherichia coli OmpT under denaturing conditions

Though OmpT has been reported to mainly cleave the peptide bond between consecutive basic amino acids, we identified more precise substrate specificity by using a series of modified substrates, termed PRX fusion proteins, consisting of 184 residues. The cleavage site of the substrate PRR was Arg140-Arg141 and the modified substrates PRX substituted all 19 natural amino acids at the P1' site instead of Arg141. OmpT under denaturing conditions (in the presence of 4 M urea) cleaved not only between two consecutive basic amino acids but also at the carboxyl side of Arg140 except for the Arg140-Asp141, -Glu141, and -Pro141 pairs. In addition to Arg140 at the P1 site, similar results were obtained when Lys140 was substituted into the P1 site. In the absence of urea, an aspartic acid residue at the P1' site was unfavorable for OmpT cleavage of synthetic decapeptides, the enzyme showed a preference for a dibasic site.     

Dissection of helix capping in T4 lysozyme by structural and thermodynamic analysis of six amino acid substitutions at Thr 59.

Threonine 59, a helix-capping residue at the amino terminus of the longest helix in T4 phage lysozyme, was substituted with valine, alanine, glycine, serine, asparagine, and aspartic acid. The valine, alanine, and glycine replacements were observed to be somewhat more destabilizing than serine, asparagine, and aspartic acid. The crystal structures of the different variants showed that changes in conformation occurred at the site of substitution, including Asp 61, which is nearby, as well as displacement of a solvent molecule that is hydrogen-bonded to the gamma-oxygen of Thr 59 in wild-type lysozyme. Neither the structures nor the stabilities of the mutant proteins support the hypothesis of Serrano and Fersht (1989) that glycine and alanine are better helix-capping residues than valine because a smaller-sized residue allows better hydration at the end of the helix. In the aspartic acid and asparagine replacements the substituted side chains form hydrogen bonds with the end of the helix, as does threonine and serine at this position. In contrast, however, the Asp and Asn side chains also make unusually close contacts with carbon atoms in Asp 61. This suggests a structural basis for the heretofore puzzling observations that asparagine is more frequently observed as a helix-capping residue than threonine [Richardson, J. S., & Richardson, D. C. (1988) Science 240, 1648-1652] yet Thr----Asn replacements at N-cap positions in barnase were found to be destabilizing [Serrano, L., & Fersht, A. R. (1989) Nature 342, 296-299].(ABSTRACT TRUNCATED AT 250 WORDS)     

Stabilization of protein structure through π–π interaction in the second coordination sphere of pseudoazurin

Noncovalent, weak interactions in the second coordination sphere of the copper active site of Pseudoazurin (PAz) from Achromobacter cycloclastes were examined using a series of Met16X variants. In this study, the differences in protein stability due to the changes in the nature of the 16th amino acid (Met, Phe, Val, Ile) were investigated by electrospray ionization mass spectrometry (ESI‐MS) and far‐UV circular dichroism (CD) as a result of acid denaturation. The percentage of native states (folded holo forms) of Met16Phe variants was estimated to be 75% at pH 2.9 although the wild‐type (WT), Met16Val and Met16Ile PAz, became completely unfolded. The high stability under acidic conditions is correlated with the result of the active site being stabilized by the aromatic substitution of the Met16 residue. The π–π interaction in the second coordination sphere makes a significant contribution to the stability of active site and the protein matrix.     

ATP as effector of inorganic pyrophosphatase of <Emphasis Type="Italic">Escherichia coli</Emphasis>. The role of residue Lys112 in binding effectors

It has been shown that PP_i, methylenediphosphonate, and ATP act as effectors of Escherichia coli inorganic pyrophosphatase (E-PPase), and that they compete for binding at the allosteric regulatory site. On the basis of chemical modification and computer modeling of a structure of the enzyme-ATP complex, a number of amino acid residues presumably involved in binding effectors has been revealed. Mutant variants Lys112Gln, Lys112Gln/Lys148Gln, and Lys112Gln/Lys115Ala of E-PPase have been obtained, as well as a modified variant of wild type E-PPase (^Adwt PPase) with a derivative of ATP chemically attached to the amino group of Lys146. Kinetic properties of these variants have been investigated and compared to the earlier described variants Lys115Ala, Arg43Gln, and Lys148Gln. Analysis of the data confirms the proposed location of an effector binding site in a cluster of positively charged amino acid residues including the side chains of Arg43, Lys146 (subunit A), Lys112, and Lys115 (subunit B). Lys112 is supposed to play a key role in forming contacts with the phosphate groups of the three studied effectors. Published in Russian in Biokhimiya, 2007, Vol. 72, No. 1, pp. 118–127.     

Cofactor Activity in Factor VIIIa of the Blood Clotting Pathway Is Stabilized by an Interdomain Bond between His281 and Ser524 Formed in Factor VIII

The factor VIII (FVIII) crystal structure suggests a possible bonding interaction of His²⁸¹ (A1 domain) with Ser⁵²⁴ (A2 domain), although the resolution of the structure (∼4 Å) does not firmly establish this bonding. To establish that side chains of these residues participate in an interdomain bond, we prepared and examined the functional properties of a residue swap variant (H281S/S524H) where His²⁸¹ and Ser⁵²⁴ residues were exchanged with one another and a disulfide-bridged variant (H281C/S524C) where the two residues were replaced with Cys. The latter variant showed efficient disulfide bonding of the A1 and A2 domains. The swap variant showed WT-like FVIII and FVIIIa stability, which were markedly reduced for H281A and S524A variants in an earlier study. The disulfide-bridged variant showed ∼20% increased FVIII stability, and FVIIIa did not decay during the time course measured. This variant also yielded 35% increased thrombin peak values compared with WT in a plasma-based thrombin generation assay. Binding analyses of H281S-A1/A3C1C2 dimer with S524H-A2 subunit yielded a near WT-like affinity value, whereas combining the variant dimer or A2 subunit with the WT complement yielded ∼5- and ∼10-fold reductions, respectively, in affinity. Other functional properties including thrombin generation potential, FIXa binding affinity, K_m for FX of FXase complexes, thrombin activation efficiency, and down-regulation by activated protein C showed similar results for the two variants compared with WT FVIII. These results indicate that the side chains of His²⁸¹ and Ser⁵²⁴ are in close proximity and contribute to a bonding interaction in FVIII that is retained in FVIIIa.     

Aromatic interactions with naphthylalanine in a β‐hairpin peptide

Stable peptides have been explored as epitope mimics for protein–protein and protein–nucleic acid interactions; however, presentation of a regular structure is critical. Aromatic interactions are ubiquitous and are competent at stabilizing a β‐hairpin fold. The greatest stabilization has been reported from pairs of tryptophan side chains. Naphthylalanine residues are often used as tryptophan replacements, but it is not clear if 1‐naphthylalanine or 2‐naphthylalanine is adequate at replicating the geometry and stability observed with tryptophan aromatic interactions. Herein, a 12‐residue peptide has been constructed with laterally disposed aromatic amino acids. A direct comparison is made between tryptophan and other bicyclic, unnatural amino acids. Significant stabilization is gained from all bicyclic amino acids; however, geometric analysis shows that only 1‐naphthylalanine adopts a similar edge to face geometry as tryptophan, whereas the 2‐naphthylalanine appears most similar to a substituted phenylalanine. Copyright © 2013 European Peptide Society and John Wiley & Sons, Ltd.     

An Unorthodox Sensory Adaptation Site in the Escherichia coli Serine Chemoreceptor

The serine chemoreceptor of Escherichia coli contains four canonical methylation sites for sensory adaptation that lie near intersubunit helix interfaces of the Tsr homodimer. An unexplored fifth methylation site, E502, lies at an intrasubunit helix interface closest to the HAMP domain that controls input-output signaling in methyl-accepting chemotaxis proteins. We analyzed, with in vivo Förster resonance energy transfer (FRET) kinase assays, the serine thresholds and response cooperativities of Tsr receptors with different mutationally imposed modifications at sites 1 to 4 and/or at site 5. Tsr variants carrying E or Q at residue 502, in combination with unmodifiable D and N replacements at adaptation sites 1 to 4, underwent both methylation and demethylation/deamidation, although detection of the latter modifications required elevated intracellular levels of CheB. These Tsr variants could not mediate a chemotactic response to serine spatial gradients, demonstrating that adaptational modifications at E502 alone are not sufficient for Tsr function. Moreover, E502 is not critical for Tsr function, because only two amino acid replacements at this residue abrogated serine chemotaxis: Tsr-E502P had extreme kinase-off output and Tsr-E502I had extreme kinase-on output. These large threshold shifts are probably due to the unique HAMP-proximal location of methylation site 5. However, a methylation-mimicking glutamine at any Tsr modification site raised the serine response threshold, suggesting that all sites influence signaling by the same general mechanism, presumably through changes in packing stability of the methylation helix bundle. These findings are consistent with control of input-output signaling in Tsr through dynamic interplay of the structural stabilities of the HAMP and methylation bundles.     

QM/MM studies on the calcium‐assisted β‐elimination mechanism of pectate lyase from bacillus subtilis

Pectate lyase utilizes the anti‐β‐elimination chemistry to catalyze the cleavage of α‐1,4 glycosidic bond between _D‐galacturonate regions during the degradation of plant polysaccharide pectin. We report here detailed mechanistic studies of the Bacillus subtilis pectate lyase (BsPel) using QM/MM calculations. It was found that the residue Arg279 serves as the catalytic base to abstract the α‐proton from C₅₂ atom of substrate Ada2 subsite, forming an unstable carbanion intermediate. The glycosidic bond of this intermediate is scissile to generate the 4,5‐unsaturated digalacturonate product and a negatively charged β‐leaving group. Two active site residues (Lys247 and Arg279) and two Ca²⁺ ions (Ca2 and Ca3) form hydrogen‐bonding and coordination interactions with C₅₂? COO^? of Ada2, respectively, which facilitate the proton abstraction and stabilize the generated carbanion intermediates. Arg284 is not the potential proton donor to saturate the leaving group. Actually, the proton source of leaving group is the solvent water molecule rather than any active site acidic residues. In addition, the calculation results suggest that careful selections of QM‐ and Active‐regions are essential to accurately explore the enzymatic reactions. Proteins 2016; 84:1606–1615. © 2016 Wiley Periodicals, Inc.     

Constitutively active glutaminase variants provide insights into the activation mechanism of anthranilate synthase

The glutamine amidotransferase (GATase) family comprises enzyme complexes which consist of glutaminase and synthase subunits that catalyze in a concerted reaction the incorporation of nitrogen within various metabolic pathways. An important feature of GATases is the strong stimulation of glutaminase activity by the associated synthase. To understand the mechanism of this tight activity regulation, we probed by site-directed mutagenesis four residues of the glutaminase subunit TrpG from anthranilate synthase that are located between the catalytic Cys-His-Glu triad and the synthase subunit TrpE. In order to minimize structural perturbations induced by the introduced exchanges, the amino acids from TrpG were substituted with the corresponding residues of the closely related glutaminase HisH from imidazole glycerol phosphate synthase. Steady-state kinetic characterization showed that, in contrast to wild-type TrpG, two TrpG variants with single exchanges constitutively hydrolyzed glutamine in the absence of TrpE. A reaction assay performed with hydroxylamine as a stronger nucleophile replacing water and a filter assay with radiolabeled glutamine indicated that the formation of the thioester intermediate is the rate-limiting step of constitutive glutamine hydrolysis. Molecular dynamics simulations with wild-type TrpG and constitutively active TrpG variants suggest that the introduced amino acid exchanges result in a distance reduction between the active site Cys-His pair, which facilitates the deprotonation of the sulfhydryl group of the catalytic cysteine and thus enables its nucleophilic attack onto the carboxamide group of the glutamine side chain. We propose that native TrpG in the anthranilate synthase complex is activated by a similar mechanism.     

Shifting the optimum pH of Bacillus circulans xylanase towards acidic side by introducing arginine

Electrostatic interactions are important in protein folding, binding, flexibility, stability and function. The pH at which the enzyme is maximally active is determined by the pK_as of the active site residues, which are modulated by several factors including the change in electrostatics in its vicinity. As the acidic xylanases are important in food and animal feed industries, electrostatic interactions are introduced in Bacillus circulans xylanase to shift their pH optima towards the acidic side. Arg substitutions are made to modulate the pKas of the active site residues. Neutral residues are substituted by Arg in such a way that the substituted residue can make direct interaction with the catalytic residues. However, the mutations with other titratable residues (Asp, Arg, Lys, His, Tyr, and Ser) present in between the catalytic sites and the substituted sites are avoided. Site directed mutagenesis was conducted to confirm the strategy. The results show the shift in pH optima of the mutants towards the acidic side by 0.5–1.5 unit. Molecular dynamics simulation of the mutant V37R reveals that the decrease in activity is due to the increase in distance between the substrate oxygen atoms and catalytic glutamates.     

Structural basis for the complete loss of GSK3β catalytic activity due to R96 mutation investigated by molecular dynamics study

Many Ser/Thr protein kinases, to be fully activated, are obligated to introduce a phospho‐Ser/Thr in their activation loop. Presently, the similarity of activation loop between two crystal complexes, i.e. glycogen synthase kinase 3β (GSK3β)‐AMPNP and GSK3β‐sulfate ion complex, indicates that the activation segment of GSK3β is preformed requiring neither a phosphorylation event nor conformational changes. GSK3β, when participated in glycogen synthesis and Wnt signaling pathways, possesses a unique feature with the preference of such substrate with a priming phosphate. Experimental mutagenesis proved that the residue arginine at amino acid 96 mutations to lysine (R96K) or alanine (R96A) selectively abolish activity on the substrates involved in glycogen synthesis signaling pathway. Based on two solved crystal structures, wild type (WT) and two mutants (R96K and R96A) GSK3β‐ATP‐phospho‐Serine (pSer) complexes were modeled. Molecular dynamics simulations and energy analysis were employed to investigate the effect of pSer involvement on the GSK3β structure in WT, and the mechanisms of GSK3β deactivation due to R96K and R96A mutations. The results indicate that the introduction of pSer to WT GSK3β generates a slight lobe closure on GSK3β without any remarkable changes, which may illuminate the experimental conclusion, whereas the conformations of GSK3β and ATP undergo significant changes in two mutants. As to GSK3β, the affected positions distribute over activation loop, α‐helix, and glycine‐rich loop. Based on coupling among the mentioned positions, the allosteric mechanisms for distorted ATP were proposed. Energy decomposition on the residues of activation loop identified the important residues Arg96 and Arg180 in anchoring the phosphate group. Proteins 2009. © 2008 Wiley‐Liss, Inc.