Differential inhibition of cytosolic PEPCK by substrate analogues. Kinetic and structural characterization of inhibitor recognition

Bidentate Substrate Binding Mode in Oxalate Decarboxylase

Oxalate decarboxylase is an Mn- and O2-dependent enzyme in the bicupin superfamily that catalyzes the redox-neutral disproportionation of the oxalate monoanion to form carbon dioxide and formate. Its best-studied isozyme is from Bacillus subtilis where it is stress-induced under low pH conditions. Current mechanistic schemes assume a monodentate binding mode of the substrate to the N-terminal active site Mn ion to make space for a presumed O2 molecule, despite the fact that oxalate generally prefers to bind bidentate to Mn. We report on X-band 13C-electron nuclear double resonance (ENDOR) experiments on 13C-labeled oxalate bound to the active-site Mn(II) in wild-type oxalate decarboxylase at high pH, the catalytically impaired W96F mutant enzyme at low pH, and Mn(II) in aqueous solution. The ENDOR spectra of these samples are practically identical, which shows that the substrate binds bidentate (κO, κO') to the active site Mn(II) ion. Domain-based local pair natural orbital coupled cluster singles and doubles (DLPNO-CCSD) calculations of the expected 13C hyperfine coupling constants for bidentate bound oxalate predict ENDOR spectra in good agreement with the experiment, supporting bidentate bound substrate. Geometry optimization of a substrate-bound minimal active site model by density functional theory shows two possible substrate coordination geometries, bidentate and monodentate. The bidentate structure is energetically preferred by ~4.7 kcal/mol. Our results revise a long-standing hypothesis regarding substrate binding in the enzyme and suggest that dioxygen does not bind to the active site Mn ion after substrate binds. The results are in agreement with our recent mechanistic hypothesis of substrate activation via a long-range electron transfer process involving the C-terminal Mn ion.

A structural perspective on the temperature-dependent activity of enzymes

Enzymes are biomolecular catalysts whose activity varies with temperature. Unlike for small-molecule catalysts, the structural ensembles of enzymes can vary substantially with temperature, and it is in general unclear how this modulates the temperature dependence of activity. Here multi-temperature X-ray crystallography was used to record structural changes from -20°C to 40°C for a mesophilic enzyme in complex with inhibitors mimicking substrate-, intermediate-, and product-bound states, representative of major complexes underlying the kinetic constantk c a t . Both inhibitors, substrates and catalytically relevant loop motifs increasingly populate catalytically competent conformations as temperature increases. These changes occur even in temperature ranges where kinetic measurements show roughly linear Arrhenius/Eyring behavior where parameters characterizing the system are assumed to be temperature independent. Simple analysis shows that linear Arrhenius/Eyring behavior can still be observed when the underlying activation energy / enthalpy values vary with temperature, e.g., due to structural changes, and that the underlying thermodynamic parameters can be far from values derived from Arrhenius/Eyring model fits. Our results indicate a critical role for temperature-dependent atomic-resolution structural data in interpreting temperature-dependent kinetic data from enzymatic systems.

Biochemical, structural, and kinetic characterization of PPi -dependent phosphoenolpyruvate carboxykinase from Propionibacterium freudenreichii

Phosphoenolpyruvate carboxykinases (PEPCK) are a well-studied family of enzymes responsible for the regulation of TCA cycle flux, where they catalyze the interconversion of oxaloacetic acid (OAA) and phosphoenolpyruvate (PEP) using a phosphoryl donor/acceptor. These enzymes have typically been divided into two nucleotide-dependent classes, those that use ATP and those that use GTP. In the 1960's and early 1970's, a group of papers detailed biochemical properties of an enzyme named phosphoenolpyruvate carboxytransphosphorylase (later identified as a third PEPCK) from Propionibacterium freudenreichii (PPi -PfPEPCK), which instead of using a nucleotide, utilized PPi to catalyze the same interconversion of OAA and PEP. The presented work expands upon the initial biochemical experiments for PPi -PfPEPCK and interprets these data considering both the current understanding of nucleotide-dependent PEPCKs and is supplemented with a new crystal structure of PPi -PfPEPCK in complex with malate at a putative allosteric site. Most interesting, the data are consistent with PPi -PfPEPCK being a Fe2+ activated enzyme in contrast with the Mn2+ activated nucleotide-dependent enzymes which in part results in some unique kinetic properties for the enzyme when compared to the more widely distributed GTP- and ATP-dependent enzymes.

Cost-effective selective deuteration of aromatic amino acid residues produces long-lived solution 1H NMR magnetization in proteins

Solution NMR studies of large proteins are hampered by rapid signal decay due to short-range dipolar 1H-1H and 1H-13C interactions. These are attenuated by rapid rotation in methyl groups and by deuteration (2H), so selective 1H,13C-isotope labelling of methyl groups in otherwise perdeuterated proteins, combined with methyl transverse relaxation optimized spectroscopy (methyl-TROSY), is now standard for solution NMR of large protein systems > 25 kDa. For non-methyl positions, long-lived magnetization can be introduced as isolated 1H-12C groups. We have developed a cost-effective chemical synthesis for producing selectively deuterated phenylpyruvate and hydroxyphenylpyruvate. Feeding these amino acid precursors to E. coli in D2O, along with selectively deuterated anthranilate and unlabeled histidine, results in isolated and long-lived 1H magnetization in the aromatic rings of Phe (HD, HZ), Tyr (HD), Trp (HH2, HE3) and His (HD2 and HE1). We are additionally able to obtain stereoselective deuteration of Asp, Asn, and Lys amino acid residues using unlabeled glucose and fumarate as carbon sources and oxalate and malonate as metabolic inhibitors. Combining these approaches produces isolated 1H-12C groups in Phe, Tyr, Trp, His, Asp, Asn, and Lys in a perdeuterated background, which is compatible with standard 1H-13C labeling of methyl groups in Ala, Ile, Leu, Val, Thr, Met. We show that isotope labeling of Ala is improved using the transaminase inhibitor L-cycloserine, and labeling of Thr is improved through addition of Cys and Met, which are known inhibitors of homoserine dehydrogenase. We demonstrate the creation of long-lived 1H NMR signals in most amino acid residues using our model system, the WW domain of human Pin1, as well as the bacterial outer membrane protein PagP.

Pharmacology and preclinical validation of a novel anticancer compound targeting PEPCK-M

BACKGROUND

Phosphoenolpyruvate carboxykinase (PEPCK) catalyzes the decarboxylation of oxaloacetate to phosphoenolpyruvate. The mitochondrial isozyme, PEPCK-M is highly expressed in cancer cells, where it plays a role in nutrient stress response. To date, pharmacological strategies to target this pathway have not been pursued.

METHODS

A compound embodying a 3-alkyl-1,8-dibenzylxanthine nucleus (iPEPCK-2), was synthesized and successfully probed in silico on a PEPCK-M structural model. Potency and target engagement in vitro and in vivo were evaluated by kinetic and cellular thermal shift assays (CETSA). The compound and its target were validated in tumor growth models in vitro and in murine xenografts.

RESULTS

Cross-inhibitory capacity and increased potency as compared to 3-MPA were confirmed in vitro and in vivo. Treatment with iPEPCK-2 inhibited cell growth and survival, especially in poor-nutrient environment, consistent with an impact on colony formation in soft agar. Finally, daily administration of the PEPCK-M inhibitor successfully inhibited tumor growth in two murine xenograft models as compared to vehicle, without weight loss, or any sign of apparent toxicity.

CONCLUSION

We conclude that iPEPCK-2 is a compelling anticancer drug targeting PEPCK-M, a hallmark gene product involved in metabolic adaptations of the tumor.

Characterization of 3-[(Carboxymethyl)thio]picolinic Acid: A Novel Inhibitor of Phosphoenolpyruvate Carboxykinase

Phosphoenolpyruvate carboxykinase (PEPCK) has traditionally been characterized for its role in the first committed step of gluconeogenesis. The current understanding of PEPCK's metabolic role has recently expanded to include it serving as a general mediator of tricarboxylic acid cycle flux. Selective inhibition of PEPCK in vivo and in vitro has been achieved with 3-mercaptopicolinic acid (MPA) (K_i ∼ 8 μM), whose mechanism of inhibition has been elucidated only recently. On the basis of crystallographic and mechanistic data of various inhibitors of PEPCK, MPA was used as the initial chemical scaffold to create a potentially more selective inhibitor, 3-[(carboxymethyl)thio]picolinic acid (CMP), which has been characterized both structurally and kinetically here. These data demonstrate that CMP acts as a competitive inhibitor at the OAA/PEP binding site, with its picolinic acid moiety coordinating directly with the M1 metal in the active site (K_i ∼ 29-55 μM). The extended carboxy tail occupies a secondary binding cleft that was previously shown could be occupied by sulfoacetate (K_i ∼ 82 μM) and for the first time demonstrates the simultaneous occupation of both OAA/PEP subsites by a single molecular structure. By occupying both the OAA/PEP binding subsites simultaneously, CMP and similar molecules can potentially be used as a starting point for the creation of additional selective inhibitors of PEPCK.

Computational evaluation of natural compounds as potential inhibitors of human PEPCK-M: an alternative for lung cancer therapy

Background

Lung cancer is the leading cause of cancer-related death worldwide. Among its subtypes, non-small cell lung cancer (NSCLC) is the most common. Recently, the mitochondrial isoform of the enzyme phosphoenolpyruvate carboxykinase (HsPEPCK-M) was identified as responsible for the metabolic adaptation in the NSCLC allowing tumor growth even under conditions of glucose deficiency. This adaptation is possible due to the role of HsPEPCK-M in gluconeogenesis, converting the oxaloacetate to phosphoenolpyruvate in the presence of GTP, which plays an important role in the energetic support of these tumors. In this context, it was shown that the inhibition or knockdown of this enzyme was able to induce apoptosis in NSCLC under low glucose conditions.

Purpose

In this study, novel putative inhibitors were proposed for the human PEPCK-M (HsPEPCK-M) based on a computer-aided approach.

Methods

Comparative modeling was used to generate 3D models for HsPEPCK-M. Subsequently, the set of natural compounds of the ZINC database was screened against HsPEPCK-M models using structure-based pharmacophore modeling and molecular docking approaches. The selected compounds were evaluated according to its chemical diversity and clustered based on chemical similarity.

Results

The pharmacophore hypotheses, generated based on known PEPCK inhibitors, were able to select 7,124 candidate compounds. These compounds were submitted to molecular docking studies using three conformations of HsPEPCK-M generated by comparative modeling. The aim was to select compounds with high predicted binding affinity for at least one of the conformations of HsPEPCK-M. After molecular docking, 612 molecules were selected as potential inhibitors of HsPEPCK-M. These compounds were clustered according to their structural similarity. Chemical profiling and binding mode analyses of these compounds allowed the proposal of four promising compounds: ZINC01656421, ZINC895296, ZINC00895535 and ZINC02571340.

Conclusion

These compounds may be considered as potential candidates for HsPEPCK-M inhibitors and may also be used as lead compounds for the development of novel HsPEPCK-M inhibitors.

Inhibition of Pig Phosphoenolpyruvate Carboxykinase Isoenzymes by 3-Mercaptopicolinic Acid and Novel Inhibitors

There exist two isoforms of cytosolic phosphoenolpyruvate carboxykinase (PEPCK-C) in pig populations that differ in a single amino acid (Met139Leu). The isoenzymes have different kinetic properties, affecting more strongly the K_m and V_max of nucleotides. They are associated to different phenotypes modifying traits of considerable economic interest. In this work we use inhibitors of phosphoenolpyruvate carboxykinase activity to search for further differences between these isoenzymes. On the one hand we have used the well-known inhibitor 3-mercaptopicolinic acid. Its inhibition patterns were the same for both isoenzymes: a three-fold decrease of the K_i values for GTP in 139Met and 139Leu (273 and 873 μM, respectively). On the other hand, through screening of a chemical library we have found two novel compounds with inhibitory effects of a similar magnitude to that of 3-mercaptopicolinic acid but with less solubility and specificity. One of these novel compounds, (N'1-({5-[1-methyl-5-(trifluoromethyl)-1H-pyrazol-3-yl]-2-thienyl}methylidene)-2,4-dichlorobenzene-1-carbohydrazide), exhibited significantly different inhibitory effects on either isoenzyme: it enhanced threefold the apparent K_m value for GTP in 139Met, whereas in 139Leu, it reduced it from 99 to 69 μM. The finding of those significant differences in the binding of GTP reinforces the hypothesis that the Met139Leu substitution affects strongly the nucleotide binding site of PEPCK-C.

Comparison of Newly Assembled Full Length HIV-1 Integrase With Prototype Foamy Virus Integrase: Structure-Function Prospective

Background

Drug design against human immunodeficiency virus type 1 (HIV-1) integrase through its mechanistic study is of great interest in the area in biological research. The main obstacle in this area is the absence of the full-length crystal structure for HIV-1 integrase to be used as a model. A complete structure, similar to HIV-1 of a prototype foamy virus integrase in complex with DNA, including all conservative residues, is available and has been extensively used in recent investigations.

Objectives

The aim of this study was to determine whether the above model is precisely representative of HIV-1 integrase. This would critically determine the success of any designed drug using the model in deactivation of integrase and AIDS treatment.

Materials and Methods

Primarily, a new structure for HIV-1 was constructed, using a crystal structure of prototype foamy virus as the starting structure. The constructed structure of HIV-1 integrase was simultaneously simulated with a prototype foamy virus integrase on a separate occasion.

Results

Our results indicate that the HIV-1 system behaves differently from the prototype foamy virus in terms of folding, hydration, hydrophobicity of binding site and stability.

Conclusions

Based on our findings, we can conclude that HIV-1 integrase is vastly different from the prototype foamy virus integrase and does not resemble it, and the modeling output of the prototype foamy virus simulations could not be simply generalized to HIV-1 integrase. Therefore, our HIV-1 model seems to be more representative and more useful for future research.

Mechanism of preferential packaging of negative sense genomic RNA by viral nucleoproteins in Crimean-Congo hemorrhagic Fever virus

The Crimean-Congo Hemorrhagic Fever (CCHF) is an infectious disease of high virulence and mortality caused by a negative sense RNA nairovirus. The genomic RNA of CCHFV is enwrapped by its nucleoprotein. Positively charged residues on CCHFV nucleoprotein provide multiple binding sites to facilitate genomic RNA encapsidation. In the present work, we investigated the mechanism underlying preferential packaging of the negative sense genomic RNA by CCHFV nucleoprotein in the presence of host cell RNAs during viral assembly. The work included genome sequence analyses for different families of negative and positive sense RNA viruses, using serial docking experiments and molecular dynamic simulations. Our results indicated that the main determinant parameter of the nucleoprotein binding affinity for negative sense RNA is the ratio of purine/pyrimidine in the RNA molecule. A negative sense RNA with a purine/pyrimidine ratio (>1) higher than that of a positive sense RNA (<1) exhibits higher affinity for the nucleoprotein. Our calculations revealed that a negative sense RNA expresses about 0.5 kJ/mol higher binding energy per nucleotide compared to a positive sense RNA. This energy difference produces a binding energy high enough to make the negative sense RNA, the preferred substrate for packaging by CCHFV nucleoprotein in the presence of cellular or complementary positive sense RNAs. The outcome of this study may contribute to ongoing researches on other viral diseases caused by negative sense RNA viruses such as Ebola virus which poses a security threat to all humanity.

Utilization of Substrate Intrinsic Binding Energy for Conformational Change and Catalytic Function in Phosphoenolpyruvate Carboxykinase

Phosphoenolpyruvate carboxykinase (PEPCK) is an essential metabolic enzyme operating in the gluconeogenesis and glyceroneogenesis pathways. Previous work has demonstrated that the enzyme cycles between a catalytically inactive open state and a catalytically active closed state. The transition of the enzyme between these states requires the transition of several active site loops to shift from mobile, disordered structural elements to stable ordered states. The mechanism by which these disorder-order transitions are coupled to the ligation state of the active site however is not fully understood. To further investigate the mechanisms by which the mobility of the active site loops is coupled to enzymatic function and the transitioning of the enzyme between the two conformational states, we have conducted structural and functional studies of point mutants of E89. E89 is a proposed key member of the interaction network of mobile elements as it resides in the R-loop region of the enzyme active site. These new data demonstrate the importance of the R-loop in coordinating interactions between substrates at the OAA/PEP binding site and the mobile R- and Ω-loop domains. In turn, the studies more generally demonstrate the mechanisms by which the intrinsic ligand binding energy can be utilized in catalysis to drive unfavorable conformational changes, changes that are subsequently required for both optimal catalytic activity and fidelity.

Inhibition and Allosteric Regulation of Monomeric Phosphoenolpyruvate Carboxykinase by 3-Mercaptopicolinic Acid

For almost 40 years, it has been known that tryptophan metabolites and picolinic acid analogues act as inhibitors of gluconeogenesis. Early studies observed that 3-mercaptopicolinic acid (MPA) was a potent hypoglycemic agent via inhibition of glucose synthesis through the specific inhibition of phosphoenolpyruvate carboxykinase (PEPCK) in the gluconeogenesis pathway. Despite prior kinetic investigation, the mechanism of the inhibition by MPA is unclear. To clarify the mechanism of inhibition exerted by MPA on PEPCK, we have undertaken structural and kinetic studies. The kinetic data in concert with crystallographic structures of PEPCK in complex with MPA and the substrates for the reaction illustrate that PEPCK is inhibited by the binding of MPA at two discrete binding sites: one acting in a competitive fashion with PEP/OAA (∼10 μM) and the other acting at a previously unidentified allosteric site (Ki ∼ 150 μM). The structural studies suggest that binding of MPA to the allosteric pocket stabilizes an altered conformation of the nucleotide-binding site that in turn reduces the affinity of the enzyme for the nucleotide.

The mitochondrial isoform of phosphoenolpyruvate carboxykinase (PEPCK-M) and glucose homeostasis: has it been overlooked?

BACKGROUND

Plasma glucose levels are tightly regulated within a narrow physiologic range. Insulin-mediated glucose uptake by tissues must be balanced by the appearance of glucose from nutritional sources, glycogen stores, or gluconeogenesis. In this regard, a common pathway regulating both glucose clearance and appearance has not been described. The metabolism of glucose to produce ATP is generally considered to be the primary stimulus for insulin release from beta-cells. Similarly, gluconeogenesis from phosphoenolpyruvate (PEP) is believed to be the primarily pathway via the cytosolic isoform of phosphoenolpyruvate carboxykinase (PEPCK-C). These models cannot adequately explain the regulation of insulin secretion or gluconeogenesis.

SCOPE OF REVIEW

A metabolic sensing pathway involving mitochondrial GTP (mtGTP) and PEP synthesis by the mitochondrial isoform of PEPCK (PEPCK-M) is associated with glucose-stimulated insulin secretion from pancreatic beta-cells. Here we examine whether there is evidence for a similar mtGTP-dependent pathway involved in gluconeogenesis. In both islets and the liver, mtGTP is produced at the substrate level by the enzyme succinyl CoA synthetase (SCS-GTP) with a rate proportional to the TCA cycle. In the beta-cell PEPCK-M then hydrolyzes mtGTP in the production of PEP that, unlike mtGTP, can escape the mitochondria to generate a signal for insulin release. Similarly, PEPCK-M and mtGTP might also provide a significant source of PEP in gluconeogenic tissues for the production of glucose. This review will focus on the possibility that PEPCK-M, as a sensor for TCA cycle flux, is a key mechanism to regulate both insulin secretion and gluconeogenesis suggesting conservation of this biochemical mechanism in regulating multiple aspects of glucose homeostasis. Moreover, we propose that this mechanism may be important for regulating insulin secretion and gluconeogenesis compared to canonical nutrient sensing pathways.

MAJOR CONCLUSIONS

PEPCK-M, initially believed to be absent in islets, carries a substantial metabolic flux in beta-cells. This flux is intimately involved with the coupling of glucose-stimulated insulin secretion. PEPCK-M activity may have been similarly underestimated in glucose producing tissues and could potentially be an unappreciated but important source of gluconeogenesis.

GENERAL SIGNIFICANCE

The generation of PEP via PEPCK-M may occur via a metabolic sensing pathway important for regulating both insulin secretion and gluconeogenesis. This article is part of a Special Issue entitled Frontiers of Mitochondrial Research.

Whiskers-less HIV-protease: a possible way for HIV-1 deactivation

BACKGROUND

Among viral enzymes, the human HIV-1 protease comprises the most interesting target for drug discovery. There are increasing efforts focused on designing more effective inhibitors for HIV-1 protease in order to prevent viral replication in AIDS patients. The frequent and continuous mutation of HIV-1 protease gene creates a formidable obstacle for enzyme inhibition which could not be overcome by the traditional single drug therapy. Nowadays, in vitro and in silico studies of protease inhibition constitute an advanced field in biological researches. In this article, we tried to simulate protease-substrate complexes in different states; a native state and states with whiskers deleted from one and two subunits. Molecular dynamic simulations were carried out in a cubic box filled with explicit water at 37°C and in 1atomsphere of pressure.

RESULTS

Our results showed that whisker truncation of protease subunits causes the dimer structure to decrease in compactness, disrupts substrate-binding site interactions and changes in flap status simultaneously.

CONCLUSIONS

Based on our findings we claim that whisker truncation even when applied to a single subunit, threats dimer association which probably leads to enzyme inactivation. We may postulate that inserting a gene to express truncated protease inside infected cells can interfere with protease dimerization. The resulted proteases would presumably have a combination of native and truncated subunits in their structures which exert no enzyme activities as evidenced by the present work. Our finding may create a new field of research in HIV gene therapy for protease inhibition, circumventing problems of drug resistance.

The Ω-loop lid domain of phosphoenolpyruvate carboxykinase is essential for catalytic function

Phosphoenolpyruvate carboxykinase (PEPCK) is an essential metabolic enzyme operating in the gluconeogenesis and glyceroneogenesis pathways. Recent studies have demonstrated that the enzyme contains a mobile active site lid domain that undergoes a transition between an open, disorded conformation and a closed, ordered conformation as the enzyme progresses through the catalytic cycle. The understanding of how this mobile domain functions in catalysis is incomplete. Previous studies showed that the closure of the lid domain stabilizes the reaction intermediate and protects the reactive intermediate from spurious protonation and thus contributes to the fidelity of the enzyme. To more fully investigate the roles of the lid domain in PEPCK function, we introduced three mutations that replaced the 11-residue lid domain with one, two, and three glycine residues. Kinetic analysis of the mutant enzymes demonstrates that none of the enzyme constructs exhibit any measurable kinetic activity, resulting in a decrease in the catalytic parameters of at least 10(6). Structural characterization of the mutants in complexes representing the catalytic cycle suggests that the inactivity is due to a role for the lid domain in the formation of the fully closed state of the enzyme that is required for catalytic function. In the absence of the lid domain, the enzyme is unable to achieve the fully closed state and is rendered inactive despite possessing all of the residues and substrates required for catalytic function. This work demonstrates how enzyme catalytic function can be abolished through the alteration of conformational equilibria despite all the elements required for chemical conversion of substrates to products remaining intact.

Nature of protein-CO2 interactions as elucidated via molecular dynamics

Rising global temperatures require innovative measures to reduce atmospheric concentrations of CO(2). The most successful carbon capture technology on Earth is the enzymatic capture of CO(2) and its sequestration in the form of glucose. Efforts to improve upon or mimic this naturally occurring process will require a rich understanding of protein-CO(2) interactions. Toward that end, extensive all-atom molecular dynamics (MD) simulations were performed on the CO(2)-utilizing enzyme phosphoenolpyruvate carboxykinase (PEPCK). Preliminary simulations were performed using implicit and explicit solvent models, which yielded similar results: arginine, lysine, tyrosine, and asparagine enhance the ability of a protein to bind carbon dioxide. Extensive explicit solvent simulations were performed for both wild-type PEPCK and five single-point PEPCK mutants, revealing three prevalent channels by which CO(2) enters (or exits) the active site cleft, as well as a fourth channel (observed only once), the existence of which can be rationalized in terms of the position of a key Arg residue. The strongest CO(2) binding sites in these simulations consist of appropriately positioned hydrogen bond donors and acceptors. Interactions between CO(2) and both Mn(2+) and Mg(2+) present in PEPCK are minimal due to the stable protein- and solvent-based coordination environments of these cations. His 232, suggested by X-ray crystallography as being a potential important CO(2) binding site, is indeed found to be particularly "CO(2)-philic" in these simulations. Finally, a recent mechanism, proposed on the basis of X-ray crystallography, for PEPCK active site lid closure is discussed in light of the MD trajectories. Overall, the results of this work will prove useful not only to scientists investigating PEPCK, but also to groups seeking to develop an environmentally benign, protein-based carbon capture, sequestration, and utilization system.

Biophysical limits of protein-ligand binding

In classic work, Kuntz et al. (Proc. Nat. Acad. Sci. USA 1999, 96, 9997-10002) introduced the concept of ligand efficiency. Though that study focused primarily on drug-like molecules, it also showed that metal binding led to the greatest ligand efficiencies. Here, the physical limits of binding are examined across the wide variety of small molecules in the Binding MOAD database. The complexes with the greatest ligand efficiencies share the trait of being small, charged ligands bound in highly charged, well buried binding sites. The limit of ligand efficiency is -1.75 kcal/mol·atom for the protein-ligand complexes within Binding MOAD, and 95% of the set have efficiencies below a "soft limit" of -0.83 kcal/mol·atom. On the basis of buried molecular surface area, the hard limit of ligand efficiency is -117 cal/mol·Å(2), which is in surprising agreement with the limit of macromolecule-protein binding. Close examination of the most efficient systems reveals their incredibly high efficiency is dictated by tight contacts between the charged groups of the ligand and the pocket. In fact, a misfit of 0.24 Å in the average contacts inherently decreases the maximum possible efficiency by at least 0.1 kcal/mol·atom.

Saccharomyces cerevisiae phosphoenolpyruvate carboxykinase: the relevance of Glu299 and Leu460 for nucleotide binding

A homology model of Saccharomyces cerevisiae phosphoenolpyruvate (PEP) carboxykinase (ATP + oxaloacetate right arrow over left arrow ADP + PEP + CO(2)) in complex with its substrates shows that the isobutyl group of Leu460 is in close proximity to the adenine ring of the nucleotide, while the carboxyl group of Glu299 is within hydrogen-bonding distance of the ribose 2'OH. The Leu460Ala mutation caused three-fold and seven-fold increases in the K (m) for ADPMn(-) and ATPMn(2-), respectively, while the Glu299Ala mutation had no effect. Binding studies showed losses of approximately 2 kcal mol(-1) in the nucleotide binding affinity due to the Leu460Ala mutation and no effect for the Glu299Ala mutation. PEP carboxykinase utilized 2'deoxyADP and 2'deoxyATP as substrates with kinetic and equilibrium dissociation constants very similar to those of ADP and ATP, respectively. These results show that the hydrophobic interaction between Leu460 and the adenine ring of the nucleotide significantly contributed to the nucleotide affinity of the enzyme. The 2'deoxy nucleotide studies and the lack of an effect of the Glu299Ala mutation in nucleotide binding suggest that the possible hydrogen bond contributed by Glu299 and the ribose 2'OH group may not be relevant for nucleotide binding.

Increasing the conformational entropy of the Omega-loop lid domain in phosphoenolpyruvate carboxykinase impairs catalysis and decreases catalytic fidelity

Many studies have shown that the dynamic motions of individual protein segments can play an important role in enzyme function. Recent structural studies of the gluconeogenic enzyme phosphoenolpyruvate carboxykinase (PEPCK) demonstrate that PEPCK contains a 10-residue Omega-loop domain that acts as an active site lid. On the basis of these structural studies, we have previously proposed a model for the mechanism of PEPCK catalysis in which the conformation of this mobile lid domain is energetically coupled to ligand binding, resulting in the closed conformation of the lid, necessary for correct substrate positioning, becoming more energetically favorable as ligands associate with the enzyme. Here we test this model by introducing a point mutation (A467G) into the center of the Omega-loop lid that is designed to increase the entropic penalty for lid closure. Structural and kinetic characterization of this mutant enzyme demonstrates that the mutation has decreased the favorability of the enzyme adapting the closed lid conformation. As a consequence of this shift in the equilibrium defining the conformation of the active site lid, the enzyme's ability to stabilize the reaction intermediate is weakened, resulting in catalytic defect. This stabilization is initially surprising, as the lid domain makes no direct contacts with the enolate intermediate formed during the reaction. Furthermore, during the conversion of OAA to PEP, the destabilization of the lid-closed conformation results in the reaction becoming decoupled as the enolate intermediate is protonated rather than phosphorylated, resulting in the formation of pyruvate. Taken together, the structural and kinetic characterization of A467G-PEPCK supports our model of the role of the active site lid in catalytic function and demonstrates that the shift in the lowest-energy conformation between open and closed lid states is a function of the free energy available to the enzyme through ligand binding and the entropic penalty for ordering of the 10-residue Omega-loop lid domain.

Electrostatic interactions play a significant role in the affinity of Saccharomyces cerevisiae phosphoenolpyruvate carboxykinase for Mn2+

Phosphoenolpyruvate (PEP) carboxykinases catalyse the reversible formation of oxaloacetate (OAA) and ATP (or GTP) from PEP, ADP (or GDP) and CO(2). They are activated by Mn(2+), a metal ion that coordinates to the protein through the epsilon-amino group of a lysine residue, the N(epsilon-2)-imidazole of a histidine residue, and the carboxylate from an aspartic acid residue. Neutrality in the epsilon-amino group of Lys213 of Saccharomyces cerevisiae PEP carboxykinase is expected to be favoured by the vicinity of ionised Lys212. Glu272 and Glu284, located close to Lys212, should, in turn, electrostatically stabilise its positive charge and hence assist in keeping the epsilon-amino group of Lys213 in a neutral state. The mutations Glu272Gln, Glu284Gln, and Lys212Met increased the activation constant for Mn(2+) in the main reaction of the enzyme up to seven-fold. The control mutation Lys213Gln increased this constant by ten-fold, as opposed to control mutation Lys212Arg, which did not affect the Mn(2+) affinity of the enzyme. These observations indicate a role for Glu272, Glu284, and Lys212 in assisting Lys213 to properly bind Mn(2+). In an unexpected result, the mutations Glu284Gln, Lys212Met and Lys213Gln changed the nucleotide-independent OAA decarboxylase activity of S. cerevisiae PEP carboxykinase into an ADP-requiring activity, implying an effect on the OAA binding characteristics of PEP carboxykinase.

Functional evaluation of serine 252 of Saccharomyces cerevisiae phosphoenolpyruvate carboxykinase

Saccharomyces cerevisiae phosphoenolpyruvate (PEP) carboxykinase mutant Ser252Ala, affecting the conserved Walker A serine residue, was characterized to elucidate the role of this serine residue. The substitution did not result in changes in the protein structure, as indicated by circular dichroism, intrinsic fluorescence spectroscopy, and gel-exclusion chromatography. Kinetic analysis of the mutated enzyme in both directions of the main reaction and in the two secondary reactions showed an approximately 50-fold increase in apparent K(m) for oxaloacetate with minor alterations in the other kinetic parameters. These results show that the hydroxyl group of serine 252 is required for proper oxaloacetate interaction.

Enzymes with lid-gated active sites must operate by an induced fit mechanism instead of conformational selection

The induced fit and conformational selection/population shift models are two extreme cases of a continuum aimed at understanding the mechanism by which the final key-lock or active enzyme conformation is achieved upon formation of the correctly ligated enzyme. Structures of complexes representing the Michaelis and enolate intermediate complexes of the reaction catalyzed by phosphoenolpyruvate carboxykinase provide direct structural evidence for the encounter complex that is intrinsic to the induced fit model and not required by the conformational selection model. In addition, the structural data demonstrate that the conformational selection model is not sufficient to explain the correlation between dynamics and catalysis in phosphoenolpyruvate carboxykinase and other enzymes in which the transition between the uninduced and the induced conformations occludes the active site from the solvent. The structural data are consistent with a model in that the energy input from substrate association results in changes in the free energy landscape for the protein, allowing for structural transitions along an induced fit pathway.

Stereochemistry of the carboxylation reaction catalyzed by the ATP-dependent phosphoenolpyruvate carboxykinases from Saccharomyces cerevisiae and Anaerobiospirillum succiniciproducens

The stereochemistry of CO(2) addition to phosphoenolpyruvate (PEP) to yield oxaloacetate catalyzed by ATP-dependent Saccharomyces cerevisiae and Anaerobiospirillum succiniciproducens PEP carboxykinases was determined using (Z)-3-fluorophosphoenolpyruvate ((Z)-F-PEP) as a substrate analog. A. succiniciproducens and S. cerevisiae PEP carboxykinases utilized (Z)-F-PEP with 1/14 and 1/47 the respective K(m) values for PEP. On the other hand, in the bacterial and yeast enzymes k(cat) was reduced to 1/67 and 1/48 the value with PEP, respectively. The binding affinity of pyridoxylphosphate-labeled S. cerevisiae and A. succiniciproducens PEP carboxykinases for PEP and (Z)-F-PEP was checked and found to be of similar magnitude for both substrates, suggesting that the lowered K(m) values for the fluorine-containing PEP analog are due to kinetic effects. The lowered k(cat) values when using (Z)-F-PEP as substrate suggest that the electron withdrawing effect of fluorine affects the nucleophilic attack of the double bond of (Z)-F-PEP to CO(2). For the stereochemical analyses, the carboxylation of (Z)-F-PEP was coupled to malate dehydrogenase to yield 3-fluoromalate, which was analyzed by (19)F NMR. The fluoromalate obtained was identified as (2R, 3R)-3-fluoromalate for both the A. succiniciproducens and S. cerevisiae PEP carboxykinases, thus indicating that CO(2) addition to (Z)-F-PEP, and hence PEP, takes place through the 2-si face of the double bond. These results, together with previously published data [Rose, I.A. et al. J. Biol. Chem. 244 (1969) 6130-6133; Hwang, S.H. and Nowak, T. Biochemistry 25 (1986) 5590-5595] indicate that PEP carboxykinases, no matter their nucleotide specificity, catalyze the carboxylation of PEP from the 2-si face of the double bond.