Reaction path of phosphofructo-1-kinase is altered by mutagenesis and alternative substrates

Tunnel-like region observed as a potential allosteric site in Staphylococcus aureus Glyceraldehyde-3-phosphate dehydrogenase

Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) catalyzing the sixth step of glycolysis has been investigated for allosteric features that might be used as potential target for specific inhibition of Staphylococcus aureus (S.aureus). X-ray structure of bacterial enzyme for which a tunnel-like opening passing through the center previously proposed as an allosteric site has been subjected to six independent 500 ns long Molecular Dynamics simulations. Harmonic bond restraints were employed at key residues to underline the allosteric feature of this region. A noticeable reduction was observed in the mobility of NAD+ binding domains when restrictions were applied. Also, a substantial decrease in cross-correlations between distant Cα fluctuations was detected throughout the structure. Mutual information (MI) analysis revealed a similar decrease in the degree of correspondence in positional fluctuations in all directions everywhere in the receptor. MI between backbone and side-chain torsional variations changed its distribution profile and decreased considerably around the catalytic sites when restraints were employed. Principal component analysis clearly showed that the restrained state sampled a narrower range of conformations than apo state, especially in the first principal mode due to restriction in the conformational flexibility of NAD+ binding domain. Clustering the trajectory based on catalytic site residues displayed a smaller repertoire of conformations for restrained state compared to apo. Representative snapshots subjected to k-shortest pathway analysis revealed the impact of bond restraints on the allosteric communication which displayed distinct optimal and suboptimal pathways for two states, where observed frequencies of critical residues Gln51 and Val283 at the proposed site changed considerably.

Divalent cations in human liver pyruvate kinase exemplify the combined effects of complex-equilibrium and allosteric regulation

There is growing recognition that the functional outcome of binding of an allosteric regulator to a protein/enzyme is influenced by the presence of other ligands. Here, this complexity is exemplified in the allosteric regulation of human liver pyruvate kinase (hLPYK) that is influenced by the presence of a range of divalent cation types and concentrations. For this system, fructose-1,6-bisphosphate (activator) and alanine (inhibitor) both influence the protein's affinity for the substrate, phosphoenolpyruvate (PEP). Mg2+, Mn2+, Ni2+, and Co2+ were the primary divalent cations evaluated, although Zn2+, Cd2+, V2+, Pb2+, Fe2+, and Cu2+also supported activity. Allosteric coupling between Fru-1,6-BP and PEP and between Ala and PEP varied depending on divalent cation type and concentration. Due to complicating interactions among small molecules, we did not attempt the fitting of response trends and instead we discuss a range of potential mechanisms that may explain those observed trends. Specifically, observed "substrate inhibition" may result from substrate A in one active site acting as an allosteric regulator for the affinity for substrate B in a second active site of a multimer. We also discuss apparent changes in allosteric coupling that can result from a sub-saturating concentration of a third allosteric ligand.

Structural analysis and functional study of phosphofructokinase B (PfkB) from Mycobacterium marinum

Phosphofructokinase B (PfkB) belongs to the ribokinase family, which uses the phosphorylated sugar as substrate, and catalyzes fructose-6-phosphate into fructose-1,6-diphosphate. However, the structural basis of Mycobacterium marinum PfkB is not clear. Here, we found that the PfkB protein was monomeric in solution, which was different from most enzymes in this family. The crystal structure of PfkB protein from M. marinum was solved at a resolution of 2.21 Å. The PfkB structure consists of two domains, a major three-layered α/β/α sandwich-like domain characteristic of the ribokinase-like superfamily, and a second domain composed of four-stranded β sheets. Structural comparison analysis suggested that residues G236, A237, G238, and D239 could be critical for ATP catalysis and substrate binding of PfkB. Our current work provides new insights into understanding the mechanism of the glycolysis in M. marinum.

Potential allosteric sites captured in glycolytic enzymes via residue-based network models: Phosphofructokinase, glyceraldehyde-3-phosphate dehydrogenase and pyruvate kinase

Likelihood of new allosteric sites for glycolytic enzymes, phosphofructokinase (PFK), glyceraldehyde-3-phosphate dehydrogenase (GADPH) and pyruvate kinase (PK) was evaluated for bacterial, parasitic and human species. Allosteric effect of a ligand binding at a site was revealed on the basis of low-frequency normal modes via C_α-harmonic residue network model. In bacterial PFK, perturbation of the proposed allosteric site outperformed the known allosteric one, producing a high amount of stabilization or reduced dynamics, on all catalytic regions. Another proposed allosteric spot at the dimer interface in parasitic PFK exhibited major stabilization effect on catalytic regions. In parasitic GADPH, the most desired allosteric response was observed upon perturbation of its tunnel region which incorporated key residues for functional regulation. Proposed allosteric site in bacterial PK produced a satisfactory allosteric response on all catalytic regions, whereas in human and parasitic PKs, a partial inhibition was observed. Residue network model based solely on contact topology identified the 'hub residues' with high betweenness tracing plausible allosteric communication pathways between distant functional sites. For both bacterial PFK and PK, proposed sites accommodated hub residues twice as much as the known allosteric site. Tunnel region in parasitic GADPH with the strongest allosteric effect among species, incorporated the highest number of hub residues. These results clearly suggest a one-to-one correspondence between the degree of allosteric effect and the number of hub residues in that perturbation site, which increases the likelihood of its allosteric nature.

Identification of Alternative Allosteric Sites in Glycolytic Enzymes for Potential Use as Species-Specific Drug Targets

Three allosteric glycolytic enzymes, phosphofructokinase, glyceraldehyde-3 phosphate dehydrogenase and pyruvate kinase, associated with bacterial, parasitic and human species, were explored to identify potential allosteric sites that would be used as prime targets for species-specific drug design purposes using a newly developed approach which incorporates solvent mapping, elastic network modeling, sequence and structural alignments. The majority of binding sites detected by solvent mapping overlapped with the interface regions connecting the subunits, thus appeared as promising target sites for allosteric regulation. Each binding site was then evaluated by its ability to alter the global dynamics of the receptor defined by the percentage change in the frequencies of the lowest-frequency modes most significantly and as anticipated, the most effective ones were detected in the vicinity of the well-reported catalytic and allosteric sites. Furthermore, some of our proposed regions intersected with experimentally resolved sites which are known to be critical for activity regulation, which further validated our approach. Despite the high degree of structural conservation encountered between bacterial/parasitic and human glycolytic enzymes, the majority of the newly presented allosteric sites exhibited a low degree of sequence conservation which further increased their likelihood to be used as species-specific target regions for drug design studies.

Looking for a sequence based allostery definition: a statistical journey at different resolution scales

The aim of this work was to detect allosteric hotspots signatures characterizing protein regions acting as the 'key drivers' of global allosteric conformational change. We computationally estimated the relative strength of intra-molecular interaction in allosteric proteins between two putative allostery-susceptible sites using a co-evolution model based upon the optimization of the cross-correlation in terms of free-energy-transfer hydrophobicity scale (Tanford scale) distribution along the chain. Cross-Recurrence Quantification Analysis (Cross-RQA) applied on the sequences of allostery susceptible sites showed evidence of strong interaction amongst allosteric susceptible sites. This could be due to transient weak molecular bonds between allostery susceptible patches enabling regions far-apart to come together. Further, using a large protein dataset, by comparing allosteric protein set with a randomly generated sequence population as well as a generic protein set, we reconfirmed our earlier findings that hydrophobicity patterning (as formalized by Recurrence Quantification Analysis (RQA) descriptors) may serve as determinant of allostery and its relevance in the transmission of allosteric conformational change. We applied RQA to free-energy-transfer hydrophobicity-transformed amino acid sequences of the allostery dataset to extract allostery specific global sequence features. These free-energy-transfer hydrophobicity-based RQA markers proved to be representative of allosteric signatures and not related to the differences between randomly generated and real proteins. These free-energy-transfer hydrophobicity-based RQA markers when evaluated by pattern recognition tools could distinguish allosteric proteins with 92% accuracy.

The crystal complex of phosphofructokinase-2 of Escherichia coli with fructose-6-phosphate: kinetic and structural analysis of the allosteric ATP inhibition

Substrate inhibition by ATP is a regulatory feature of the phosphofructokinases isoenzymes from Escherichia coli (Pfk-1 and Pfk-2). Under gluconeogenic conditions, the loss of this regulation in Pfk-2 causes substrate cycling of fructose-6-phosphate (fructose-6-P) and futile consumption of ATP delaying growth. In the present work, we have broached the mechanism of ATP-induced inhibition of Pfk-2 from both structural and kinetic perspectives. The crystal structure of Pfk-2 in complex with fructose-6-P is reported to a resolution of 2 Å. The comparison of this structure with the previously reported inhibited form of the enzyme suggests a negative interplay between fructose-6-P binding and allosteric binding of MgATP. Initial velocity experiments show a linear increase of the apparent K(0.5) for fructose-6-P and a decrease in the apparent k(cat) as a function of MgATP concentration. These effects occur simultaneously with the induction of a sigmoidal kinetic behavior (n(H) of approximately 2). Differences and resemblances in the patterns of fructose-6-P binding and the mechanism of inhibition are discussed for Pfk-1 and Pfk-2, as an example of evolutionary convergence, because these enzymes do not share a common ancestor.

On the chemical mechanism of succinic semialdehyde dehydrogenase (GabD1) from Mycobacterium tuberculosis

Succinic semialdehyde dehydrogenases (SSADHs) are ubiquitous enzymes that catalyze the NAD(P)+-coupled oxidation of succinic semialdehyde (SSA) to succinate, the last step of the γ-aminobutyrate shunt. Mycobacterium tuberculosis encodes two paralogous SSADHs (gabD1 and gabD2). Here, we describe the first mechanistic characterization of GabD1, using steady-state kinetics, pH-rate profiles, ¹H NMR, and kinetic isotope effects. Our results confirmed SSA and NADP+ as substrates and demonstrated that a divalent metal, such as Mg²+, linearizes the time course. pH-rate studies failed to identify any ionizable groups with pK(a) between 5.5 and 10 involved in substrate binding or rate-limiting chemistry. Primary deuterium, solvent and multiple kinetic isotope effects revealed that nucleophilic addition to SSA is very fast, followed by a modestly rate-limiting hydride transfer and fast thioester hydrolysis. Proton inventory studies revealed that a single proton is associated with the solvent-sensitive rate-limiting step. Together, these results suggest that product dissociation and/or conformational changes linked to it are rate-limiting. Using structural information for the human homolog enzyme and ¹H NMR, we further established that nucleophilic attack takes place at the Si face of SSA, generating a thiohemiacetal with S stereochemistry. Deuteride transfer to the Pro-R position in NADP+ generates the thioester intermediate and [4A-²H, 4B-¹H] NADPH. A chemical mechanism based on these data and the structural information available is proposed.

Crystallographic structure of phosphofructokinase-2 from Escherichia coli in complex with two ATP molecules. Implications for substrate inhibition

Phosphofructokinase-1 and -2 (Pfk-1 and Pfk-2, respectively) from Escherichia coli belong to different homologous superfamilies. However, in spite of the lack of a common ancestor, they share the ability to catalyze the same reaction and are inhibited by the substrate MgATP. Pfk-2, an ATP-dependent 6-phosphofructokinase member of the ribokinase-like superfamily, is a homodimer of 66 kDa subunits whose oligomerization state is necessary for catalysis and stability. The presence of MgATP favors the tetrameric form of the enzyme. In this work, we describe the structure of Pfk-2 in its inhibited tetrameric form, with each subunit bound to two ATP molecules and two Mg ions. The present structure indicates that substrate inhibition occurs due to the sequential binding of two MgATP molecules per subunit, the first at the usual site occupied by the nucleotide in homologous enzymes and the second at the allosteric site, making a number of direct and Mg-mediated interactions with the first. Two configurations are observed for the second MgATP, one of which involves interactions with Tyr23 from the adjacent subunit in the dimer and the other making an unusual non-Watson-Crick base pairing with the adenine in the substrate ATP. The oligomeric state observed in the crystal is tetrameric, and some of the structural elements involved in the binding of the substrate and allosteric ATPs are also participating in the dimer-dimer interface. This structure also provides the grounds to compare analogous features of the nonhomologous phosphofructokinases from E. coli.

Is allostery an intrinsic property of all dynamic proteins?

Allostery involves coupling of conformational changes between two widely separated binding sites. The common view holds that allosteric proteins are symmetric oligomers, with each subunit existing in "at least" two conformational states with a different affinity for ligands. Recent observations such as the allosteric behavior of myoglobin, a classical example of a nonallosteric protein, call into question the existing allosteric dogma. Here we argue that all (nonfibrous) proteins are potentially allosteric. Allostery is a consequence of re-distributions of protein conformational ensembles. In a nonallosteric protein, the binding site shape may not show a concerted second-site change and enzyme kinetics may not reflect an allosteric transition. Nevertheless, appropriate ligands, point mutations, or external conditions may facilitate a population shift, leading a presumably nonallosteric protein to behave allosterically. In principle, practically any potential drug binding to the protein surface can alter the conformational redistribution. The question is its effectiveness in the redistribution of the ensemble, affecting the protein binding sites and its function. Here, we review experimental observations validating this view of protein allostery.

Structural and functional roles of Cys-238 and Cys-295 in Escherichia coli phosphofructokinase-2

Modification of Escherichia coli phosphofructokinase-2 (Pfk-2) with pyrene maleimide (PM) results in a rapid inactivation of the enzyme. The loss of enzyme activity correlates with the incorporation of 2 mol of PM/mol of subunit and the concomitant dissociation of the dimeric enzyme. The two modified residues were identified as Cys-238 and Cys-295. In the presence of the negative allosteric effector, MgATP, Cys-238 was the only modified cysteine residue. Kinetic characterization of the Cys-238-labelled Pfk-2 indicates that the enzyme is fully active, with the kinetic constants ( K(m), kcat) being almost identical to the ones obtained for the native enzyme. The modified enzyme is a monomer in the absence of ligands and, like the native enzyme, behaves as a tetramer in the presence of the nucleotide. However, in the presence of fructose-6-phosphate (fru-6-P) and ATP(-4), the enzyme behaves as a dimer, suggesting that the monomers undergo re-association in the presence of the substrates and that the active species is a dimer. Modification of Pfk-2 with eosin-5-maleimide (EM) results in the labelling of Cys-295. This modified enzyme is inactive and is not able to bind to the allosteric effector, remaining as a dimer in its presence. Nonetheless, Cys-295-labelled Pfk-2 is able to bind to the substrate fru-6-P in an hyperbolic fashion with a K(d) value that is 6-fold higher than the one determined for the native enzyme. These are the first residues to be implicated in the activity and/or structure of the Pfk-2.