Structure-Based Network Analysis of Activation Mechanisms in the ErbB Family of Receptor Tyrosine Kinases: The Regulatory Spine Residues Are Global Mediators of Structural Stability and Allosteric Interactions

The ErbB protein tyrosine kinases are among the most important cell signaling families and mutation-induced modulation of their activity is associated with diverse functions in biological networks and human disease. We have combined molecular dynamics simulations of the ErbB kinases with the protein structure network modeling to characterize the reorganization of the residue interaction networks during conformational equilibrium changes in the normal and oncogenic forms. Structural stability and network analyses have identified local communities integrated around high centrality sites that correspond to the regulatory spine residues. This analysis has provided a quantitative insight to the mechanism of mutation-induced “superacceptor” activity in oncogenic EGFR dimers. We have found that kinase activation may be determined by allosteric interactions between modules of structurally stable residues that synchronize the dynamics in the nucleotide binding site and the αC-helix with the collective motions of the integrating αF-helix and the substrate binding site. The results of this study have pointed to a central role of the conserved His-Arg-Asp (HRD) motif in the catalytic loop and the Asp-Phe-Gly (DFG) motif as key mediators of structural stability and allosteric communications in the ErbB kinases. We have determined that residues that are indispensable for kinase regulation and catalysis often corresponded to the high centrality nodes within the protein structure network and could be distinguished by their unique network signatures. The optimal communication pathways are also controlled by these nodes and may ensure efficient allosteric signaling in the functional kinase state. Structure-based network analysis has quantified subtle effects of ATP binding on conformational dynamics and stability of the EGFR structures. Consistent with the NMR studies, we have found that nucleotide-induced modulation of the residue interaction networks is not limited to the ATP site, and may enhance allosteric cooperativity with the substrate binding region by increasing communication capabilities of mediating residues.     

Conformational regulation of the EGFR kinase core by the juxtamembrane and C-terminal tail: a molecular dynamics study

The catalytic domain of epidermal growth factor receptor (EGFR) is activated by dimerization, which requires allosteric coupling between distal dimerization and catalytic sites. Although crystal structures of EGFR kinases, solved in various conformational states, have provided important insights into EGFR activation by dimerization, the atomic details of how dimerization signals are dynamically coupled to catalytic regions of the kinase core are not fully understood. In this study, we have performed unrestrained and targeted molecular dynamics simulations on the active and inactive states of EGFR, followed by principal component analysis on the simulated trajectories, to identify correlated motions in the EGFR kinase domain upon dimerization. Our analysis reveals that the conformational changes associated with the catalytic functions of the kinase core are highly correlated with motions in the juxtamembrane (JM) and C-terminal tail, two flexible structural elements that play an active role in EGFR kinase activation and dimerization. In particular, the opening and closing of the ATP binding lobe relative to the substrate binding lobe is highly correlated with motions in the JM and C-terminal tail, suggesting that ATP and substrate binding can be coordinated with dimerization through conformational changes in the JM and C-terminal tail. Our study pinpoints key residues involved in this conformational coupling, and provides new insights into the role of the JM and C-terminal tail segments in EGFR kinase functions.     

Integration of Evolutionary Features for the Identification of Functionally Important Residues in Major Facilitator Superfamily Transporters

The identification of functionally important residues is an important challenge for understanding the molecular mechanisms of proteins. Membrane protein transporters operate two-state allosteric conformational changes using functionally important cooperative residues that mediate long-range communication from the substrate binding site to the translocation pathway. In this study, we identified functionally important cooperative residues of membrane protein transporters by integrating sequence conservation and co-evolutionary information. A newly derived evolutionary feature, the co-evolutionary coupling number, was introduced to measure the connectivity of co-evolving residue pairs and was integrated with the sequence conservation score. We tested this method on three Major Facilitator Superfamily (MFS) transporters, LacY, GlpT, and EmrD. MFS transporters are an important family of membrane protein transporters, which utilize diverse substrates, catalyze different modes of transport using unique combinations of functional residues, and have enough characterized functional residues to validate the performance of our method. We found that the conserved cores of evolutionarily coupled residues are involved in specific substrate recognition and translocation of MFS transporters. Furthermore, a subset of the residues forms an interaction network connecting functional sites in the protein structure. We also confirmed that our method is effective on other membrane protein transporters. Our results provide insight into the location of functional residues important for the molecular mechanisms of membrane protein transporters.     

Molecular Dynamics Simulation Studies of dTTP Binding and Catalysis Mediated by YhdE Dimerization

YhdE is a Maf-like (multicopy associated filamentation) protein that primarily acts as dTTPase to hydrolyze dTTP into dTMP and two phosphate molecules in cell metabolism pathway. Two crystal structures of YhdE have been previously determined, representing the open and closed active site conformations, respectively. Based on the structures, we have carried out molecular dynamics simulations and free energy calculations to investigate dTTP binding to and hydrolysis by YhdE. Our results suggest that YhdE closed state is structurally more compact than its open state at room temperature. YhdE open state is a favorable conformation for dTTP binding and closed state is a structurally favorable conformation for catalytic reaction. This observation is supported by the structure of YhdE homolog in complex with a nucleotide analog. Free energy calculations reveal that YhdE dimerization occurs preferentially in dTTP binding and is favorable for successive cooperative reaction. The key residues R11, R12 and K80, are found to contribute to the substrate stabilization. Further, YhdE dimerization and binding of dTTP induce the cooperative effect through a direct allosteric communication network in YhdE from the dTTP binding sites in the catalytic center to the intermolecular β-strand in YhdE dimer.     

Identification of regions of the tomato gamma-glutamyl kinase that are involved in allosteric regulation by proline

The first step of proline biosynthesis is catalyzed by gamma-glutamyl kinase (GK). To better understand the feedback inhibition properties of GK, we randomly mutagenized a plasmid carrying tomato tomPRO1 cDNA, which encodes proline-sensitive GK. A pool of mutagenized plasmids was transformed into an Escherichia coli GK mutant, and proline-overproducing derivatives were selected on minimal medium containing the toxic proline analog 3,4-dehydro-dl-proline. Thirty-two mutations that conferred 3,4-dehydro-dl-proline resistance were obtained. Thirteen different single amino acid substitutions were identified at nine different residues. The residues were distributed throughout the N-terminal two-thirds of the polypeptide, but 9 mutations affecting 6 residues were in a cluster of 16 residues. GK assays revealed that these amino acid substitutions caused varying degrees of diminished sensitivity to proline feedback inhibition and also resulted in a range of increased proline accumulation in vivo. GK belongs to a family of amino acid kinases, and a predicted three-dimensional model of this enzyme was constructed on the basis of the crystal structures of three related kinases. In the model, residues that were identified as important for allosteric control were located close to each other, suggesting that they may contribute to the structure of a proline binding site. The putative allosteric binding site partially overlaps the dimerization and substrate binding domains, suggesting that the allosteric regulation of GK may involve a direct structural interaction between the proline binding site and the dimerization and catalytic domains.     

Influence of Transmembrane Helix Mutations on Cytochrome P450-Membrane Interactions and Function

Human cytochrome P450 (CYP) enzymes play an important role in the metabolism of drugs, steroids, fatty acids, and xenobiotics. Microsomal CYPs are anchored in the endoplasmic reticulum membrane by an N-terminal transmembrane (TM) helix that is connected to the globular catalytic domain by a flexible linker sequence. However, the structural and functional importance of the TM-helix is unclear because it has been shown that CYPs can still associate with the membrane and have enzymatic activity in reconstituted systems after truncation or modification of the N-terminal sequence. Here, we investigated the effect of mutations in the N-terminal TM-helix residues of two human steroidogenic enzymes, CYP 17A1 and CYP 19A1, that are major drug targets for cancer therapy. These mutations were originally introduced to increase the expression of the proteins in Escherichia coli. To investigate the effect of the mutations on protein-membrane interactions and function, we carried out coarse-grained and all-atom molecular dynamics simulations of the CYPs in a phospholipid bilayer. We confirmed the orientations of the globular domain in the membrane observed in the simulations by linear dichroism measurements in a Nanodisc. Whereas the behavior of CYP 19A1 was rather insensitive to truncation of the TM-helix, mutations in the TM-helix of CYP 17A1, especially W2A and E3L, led to a gradual drifting of the TM-helix out of the hydrophobic core of the membrane. This instability of the TM-helix could affect interactions with the allosteric redox partner, cytochrome b5, required for CYP 17A1’s lyase activity. Furthermore, the simulations showed that the mutant TM-helix influenced the membrane interactions of the CYP 17A1 globular domain. In some simulations, the mutated TM-helix obstructed the substrate access tunnel from the membrane to the CYP active site, indicating a possible effect on enzyme function.     

A Combinatorial Approach to Detect Coevolved Amino Acid Networks in Protein Families of Variable Divergence

Communication between distant sites often defines the biological role of a protein: amino acid long-range interactions are as important in binding specificity, allosteric regulation and conformational change as residues directly contacting the substrate. The maintaining of functional and structural coupling of long-range interacting residues requires coevolution of these residues. Networks of interaction between coevolved residues can be reconstructed, and from the networks, one can possibly derive insights into functional mechanisms for the protein family. We propose a combinatorial method for mapping conserved networks of amino acid interactions in a protein which is based on the analysis of a set of aligned sequences, the associated distance tree and the combinatorics of its subtrees. The degree of coevolution of all pairs of coevolved residues is identified numerically, and networks are reconstructed with a dedicated clustering algorithm. The method drops the constraints on high sequence divergence limiting the range of applicability of the statistical approaches previously proposed. We apply the method to four protein families where we show an accurate detection of functional networks and the possibility to treat sets of protein sequences of variable divergence.     

Molecular Determinants Underlying Binding Specificities of the ABL Kinase Inhibitors: Combining Alanine Scanning of Binding Hot Spots with Network Analysis of Residue Interactions and Coevolution

Quantifying binding specificity and drug resistance of protein kinase inhibitors is of fundamental importance and remains highly challenging due to complex interplay of structural and thermodynamic factors. In this work, molecular simulations and computational alanine scanning are combined with the network-based approaches to characterize molecular determinants underlying binding specificities of the ABL kinase inhibitors. The proposed theoretical framework unveiled a relationship between ligand binding and inhibitor-mediated changes in the residue interaction networks. By using topological parameters, we have described the organization of the residue interaction networks and networks of coevolving residues in the ABL kinase structures. This analysis has shown that functionally critical regulatory residues can simultaneously embody strong coevolutionary signal and high network centrality with a propensity to be energetic hot spots for drug binding. We have found that selective (Nilotinib) and promiscuous (Bosutinib, Dasatinib) kinase inhibitors can use their energetic hot spots to differentially modulate stability of the residue interaction networks, thus inhibiting or promoting conformational equilibrium between inactive and active states. According to our results, Nilotinib binding may induce a significant network-bridging effect and enhance centrality of the hot spot residues that stabilize structural environment favored by the specific kinase form. In contrast, Bosutinib and Dasatinib can incur modest changes in the residue interaction network in which ligand binding is primarily coupled only with the identity of the gate-keeper residue. These factors may promote structural adaptability of the active kinase states in binding with these promiscuous inhibitors. Our results have related ligand-induced changes in the residue interaction networks with drug resistance effects, showing that network robustness may be compromised by targeted mutations of key mediating residues. This study has outlined mechanisms by which inhibitor binding could modulate resilience and efficiency of allosteric interactions in the kinase structures, while preserving structural topology required for catalytic activity and regulation.     

Prediction of amino acid residues participated in substrate recognition by cytochrome P450 subfamilies with broad substrate specificity

Cytochromes P450 comprise a large superfamily and several of their isoforms play a crucial role in metabolism of xenobiotics, including drugs. Although these enzymes demonstrate broad and cross‐substrate specificity, different cytochrome P450 subfamilies exhibit certain selectivity for some types of substrates. Analysis of amino acid residues of the active sites of six cytochrome subfamilies (CYP1А, CYP2А, CYP2С, CYP2D, CYP2E and CYP3А) enables to define subfamily‐specific patterns that consist of four residues. These residues are located on the periphery of the active sites of these cytochromes. We suggest that they can form a primary binding site at the entrance to the active site, defining cytochrome substrate recognition. Copyright © 2013 John Wiley & Sons, Ltd.     

NMR investigation of main-chain dynamics of the H80E mutant of bovine neurophysin-I: demonstration of dimerization-induced changes at the hormone-binding site

Neurophysins are hormone-binding proteins composed of two partially homologous domains. Ligand-binding (localized to the amino domain) and dimerization (involves both domains) are cooperatively linked by an as yet undefined allosteric mechanism. To help define this mechanism, we investigated the backbone dynamics of the unliganded monomeric state of the H80E mutant of bovine neurophysin-I by (15)N NMR. Model-free analysis of the NMR relaxation parameters indicated significantly greater flexibility in the carboxyl domain than in the amino domain, particularly at their dimerization interface segments. Amino domain residues critical to hormone binding were highly structured, constraining potential allosteric mechanisms. Model-free analysis additionally demonstrated chemical exchange effects, manifest as R(ex) terms, in 16 residues, 14 of which are located in the amino domain at, or immediately adjacent to, either the dimerization interface or the hormone-binding site. The chemical exchange process was further characterized using relaxation-compensated CPMG measurements, the results allowing assignment of the process to monomer-dimer exchange and calculation of the exchange kinetics, which were slow on the NMR time scale. An apparently different concentration-dependent process, distinguished from normal dimerization by its fast exchange behavior and pH-independence, also principally involved a subset of residues at and immediately adjacent to either the hormone-binding site or the amino domain dimerization interface. The data represent the first direct demonstration of an effect of dimerization in the unliganded state on neurophysin's hormone-binding site, the effect particularly involving residues that interact with hormone residue 2, and specifically identify Ser25 and Ile26 as likely intermediaries between the sites of dimerization and of hormone binding. Consistent with recent views of the role of anchor residues in protein interactions, we propose that dimerization proceeds by a fast pH-independent association of the well-structured amino domain interface that is rapidly communicated to the binding site for hormone residue 2, followed by a rate-determining pH-dependent interaction of the less structured carboxyl domain interface.     

Dimerization Interface of 3-Hydroxyacyl-CoA Dehydrogenase Tunes the Formation of Its Catalytic Intermediate

3-hydroxyacyl-CoA dehydrogenase (HAD, EC 1.1.1.35) is a homodimeric enzyme localized in the mitochondrial matrix, which catalyzes the third step in fatty acid β-oxidation. The crystal structures of human HAD and subsequent complexes with cofactor/substrate enabled better understanding of HAD catalytic mechanism. However, numerous human diseases were found related to mutations at HAD dimerization interface that is away from the catalytic pocket. The role of HAD dimerization in its catalytic activity needs to be elucidated. Here, we solved the crystal structure of Caenorhabditis elegans HAD (cHAD) that is highly conserved to human HAD. Even though the cHAD mutants (R204A, Y209A and R204A/Y209A) with attenuated interactions on the dimerization interface still maintain a dimerization form, their enzymatic activities significantly decrease compared to that of the wild type. Such reduced activities are in consistency with the reduced ratios of the catalytic intermediate formation. Further molecular dynamics simulations results reveal that the alteration of the dimerization interface will increase the fluctuation of a distal region (a.a. 60–80) that plays an important role in the substrate binding. The increased fluctuation decreases the stability of the catalytic intermediate formation, and therefore the enzymatic activity is attenuated. Our study reveals the molecular mechanism about the essential role of the HAD dimerization interface in its catalytic activity via allosteric effects.     

Allosteric regulation by membranes and hydrophobic subsites in phospholipase A2 enzymes determine their substrate specificity

Lipids play critical roles in several major chronic diseases of our times, including those that involve inflammatory sequelae such as metabolic syndrome including obesity, insulin sensitivity, and cardiovascular diseases. However, defining the substrate specificity of enzymes of lipid metabolism is a challenging task. For example, phospholipase A₂ (PLA₂) enzymes constitute a superfamily of degradative, biosynthetic, and signaling enzymes that all act stereospecifically to hydrolyze and release the fatty acids of membrane phospholipids. This review focuses on how membranes interact allosterically with enzymes to regulate cell signaling and metabolic pathways leading to inflammation and other diseases. Our group has developed “substrate lipidomics” to quantify the substrate phospholipid specificity of each PLA₂ and coupled this with molecular dynamics simulations to reveal that enzyme specificity is linked to specific hydrophobic binding subsites for membrane phospholipid substrates. We have also defined unexpected headgroup and acyl chain specificity for each of the major human PLA₂ enzymes, which explains the observed specificity at a structural level. Finally, we discovered that a unique hydrophobic binding site—and not each enzyme’s catalytic residues or polar headgroup binding site—predominantly determines enzyme specificity. We also discuss how PLA₂s release specific fatty acids after allosteric enzyme association with membranes and extraction of the phospholipid substrate, which can be blocked by stereospecific inhibitors. After decades of work, we can now correlate PLA₂ specificity and inhibition potency with molecular structure and physiological function.     

Evolutionarily conserved networks of residues mediate allosteric communication in proteins

A fundamental goal in cellular signaling is to understand allosteric communication, the process by which signals originating at one site in a protein propagate reliably to affect distant functional sites. The general principles of protein structure that underlie this process remain unknown. Here, we describe a sequence-based statistical method for quantitatively mapping the global network of amino acid interactions in a protein. Application of this method for three structurally and functionally distinct protein families (G protein-coupled receptors, the chymotrypsin class of serine proteases and hemoglobins) reveals a surprisingly simple architecture for amino acid interactions in each protein family: a small subset of residues forms physically connected networks that link distant functional sites in the tertiary structure. Although small in number, residues comprising the network show excellent correlation with the large body of mechanistic data available for each family. The data suggest that evolutionarily conserved sparse networks of amino acid interactions represent structural motifs for allosteric communication in proteins.     

Exploring Mechanisms of Allosteric Regulation and Communication Switching in the Multiprotein Regulatory Complexes of the Hsp90 Chaperone with Cochaperones and Client Proteins: Atomistic Insights from Integrative Biophysical Modeling and Network Analysis of Conformational Landscapes

Understanding molecular principles underlying Hsp90 chaperone functions and modulation of client activity is fundamental to dissect activation mechanisms of many proteins. In this work, we performed a computational investigation of the Hsp90-Hsp70-Hop-CR client complex to examine allosteric regulatory mechanisms underlying dynamic chaperone interactions and principles of chaperone-dependent client recognition and remodeling. Conformational dynamics analysis using high-resolution coarse-grained simulations and ensemble-based local frustration analysis suggest that the Hsp90 chaperone could recognize and recruit the GR client by invoking reciprocal dynamic exchanges near the intermolecular interfaces with the client. Using mutational scanning of the intermolecular residues in the Hsp90-Hsp70-Hop-GR complex, we identified binding energy hotspots in the regulatory complex. Perturbation-based network analysis and dynamic fluctuations-based modeling of allosteric residue potentials are employed for a detailed analysis of allosteric interaction networks and identification of conformational communication switches. We found that allosteric interactions between the Hsp90, the client-bound Hsp70 and Hop cochaperone can define two allosteric residue clusters that control client recruitment in which the intrinsic Hsp70 allostery is exploited to mediate integration of the Hsp70-bound client into the Hsp90 chaperone system. The results suggest a model of dynamics-driven allostery that enables efficient client recruitment and loading through allosteric couplings between intermolecular interfaces and communication switch centers. This study showed that the Hsp90 interactions with client proteins may operate under dynamic-based allostery in which ensembles of preexisting conformational states and intrinsic allosteric pathways present in the Hsp90 and Hsp70 chaperones can be exploited for recognition and integration of substrate proteins.     

Evolutionary variation and adaptation in a conserved protein kinase allosteric network: Implications for inhibitor design

The activation of protein kinases involves conformational changes in key functional regions of the kinase domain, a detailed understanding of which is essential for the design of selective protein kinase inhibitors. Through statistical analysis of protein kinase sequences and crystal structures from diverse organisms, we recently proposed that the activation of protein kinases involves a hidden strain switch in the catalytic loop. Specifically, we demonstrated that the backbone torsion-angles of residues in the catalytic loop switch from a “relaxed” to “strained” conformation upon kinase activation and the strained geometry results in a network of hydrogen bonds involving conserved non-catalytic residues in the ATP and substrate binding lobes. Here, we further explore this activation mechanism by analyzing families that lack the canonical hydrogen bonding interactions with the strained backbone. We find that alternative mechanisms have evolved to maintain catalytic loop strain. In PIM kinase, for example, two water molecules account for the lack of a conserved aspartate in the substrate binding by hydrogen bonds to the strained backbone. We discuss the relevance of these findings in the design of family-specific allosteric inhibitors, and in predicting the structural and functional impact of cancer mutations that alter the strain associated hydrogen bonding network. This article is part of a Special Issue entitled: Inhibitors of Protein Kinases (2012).     

Exploring the effect of N308D mutation on protein tyrosine phosphatase-2 cause gain-of-function activity by a molecular dynamics study

One of the most common protein tyrosine phosphatase-2 (SHP2) mutations in Noonan syndrome is the N308D mutation, and it increases the activity of the protein. However, the molecular basis of the activation of N308D mutation on SHP2 conformations is poorly understood. Here, molecular dynamic simulations were performed on SHP2 and SHP2-N308D to explore the effect of N308D mutation on SHP2 cause gain of function activity, respectively. The principal component analysis, dynamic cross-correlation map, secondary structure analysis, residue interaction networks, and solvent accessible surface area analysis suggested that the N308D mutation distorted the residues interactions network between the allosteric site (residue Gly244-Gly246) and C-SH2 domain, including the hydrogen bond formation and the binding energy. Meanwhile, the activity of catalytic site (residue Gly503-Val505) located in the Q-loop in mutant increased due to this region's high fluctuations. Therefore, the substrate had more chances to access to the catalytic activity site of the precision time protocol domain of SHP2-N308D, which was easy to be exposed. In addition, we had speculated that the Lys244 located in the allosteric site was the key residue which lead to the protein conformation changes. Consequently, overall calculations presented in this study ultimately provide a useful understanding of the increased activity of SHP2 caused by the N308D mutation.     

DNA interaction and dimerization of eukaryotic SMC hinge domains

The eukaryotic SMC1/SMC3 heterodimer is essential for sister chromatid cohesion and acts in DNA repair and recombination. Dimerization depends on the central hinge domain present in all SMC proteins, which is flanked at each side by extended coiled-coil regions that terminate in specific globular domains. Here we report on DNA interactions of the eukaryotic, heterodimeric SMC1/SMC3 hinge regions, using the two known isoforms, SMC1alpha/SMC3 and the meiotic SMC1beta/SMC3. Both dimers bind DNA with a preference for double-stranded DNA and DNA rich in potential secondary structures. Both dimers form large protein-DNA networks and promote reannealing of complementary DNA strands. DNA binding but not dimerization depends on approximately 20 amino acids of transitional sequence into the coiled-coil region. Replacement of three highly conserved glycine residues, thought to be required for dimerization, in one of the two hinge domains still allows formation of a stable dimer, but if two hinge domains are mutated dimerization fails. Single-mutant dimers bind DNA, but hinge monomers do not. Together, we show that eukaryotic hinge dimerization does not require conserved glycines in both hinge domains, that only the transition into the coiled-coil region rather than the entire coiled-coil region is necessary for DNA binding, and that dimerization is required but not sufficient for DNA binding of the eukaryotic hinge heterodimer.     

Explicit water near the catalytic I helix Thr in the predicted solution structure of CYP2A4

The solution structure of mouse cytochrome P450 2A4 (CYP2A4), a monooxygenase of deoxysteroids, was obtained using homology modeling and molecular dynamics. The solvent-equilibrated CYP2A4 preserves the essential features of CYP450s. A comparison of the models CYP2A4 and CYP2A4 with testosterone bound CYP2A4/T illustrates the changes induced by the binding of the substrate. Experimental evidence links four amino acid residues to the catalytic activity, substrate specificity, and regioselectivity of this enzyme. Three of the four amino acids are found within contact distance of the testosterone substrate, and therefore may control the binding of the substrate through direct interaction. Remarkably, a water complex previously observed in x-ray crystal structure forms near the bulge in the central I helix that contains a conserved Thr. The properties of the I helix are computed in the context of the presence or absence of ligand.     

Rational engineering of enzyme allosteric regulation through sequence evolution analysis

Control of enzyme allosteric regulation is required to drive metabolic flux toward desired levels. Although the three-dimensional (3D) structures of many enzyme-ligand complexes are available, it is still difficult to rationally engineer an allosterically regulatable enzyme without decreasing its catalytic activity. Here, we describe an effective strategy to deregulate the allosteric inhibition of enzymes based on the molecular evolution and physicochemical characteristics of allosteric ligand-binding sites. We found that allosteric sites are evolutionarily variable and comprised of more hydrophobic residues than catalytic sites. We applied our findings to design mutations in selected target residues that deregulate the allosteric activity of fructose-1,6-bisphosphatase (FBPase). Specifically, charged amino acids at less conserved positions were substituted with hydrophobic or neutral amino acids with similar sizes. The engineered proteins successfully diminished the allosteric inhibition of E. coli FBPase without affecting its catalytic efficiency. We expect that our method will aid the rational design of enzyme allosteric regulation strategies and facilitate the control of metabolic flux.