A very short hydrogen bond provides only moderate stabilization of an enzyme-inhibitor complex of citrate synthase

Association of the malate dehydrogenase-citrate synthase metabolon is modulated by intermediates of the Krebs tricarboxylic acid cycle

Mitochondrial malate dehydrogenase (MDH)-citrate synthase (CS) multi-enzyme complex is a part of the Krebs tricarboxylic acid (TCA) cycle ‘metabolon’ which is enzyme machinery catalyzing sequential reactions without diffusion of reaction intermediates into a bulk matrix. This complex is assumed to be a dynamic structure involved in the regulation of the cycle by enhancing metabolic flux. Microscale Thermophoresis analysis of the porcine heart MDH-CS complex revealed that substrates of the MDH and CS reactions, NAD⁺ and acetyl-CoA, enhance complex association while products of the reactions, NADH and citrate, weaken the affinity of the complex. Oxaloacetate enhanced the interaction only when it was present together with acetyl-CoA. Structural modeling using published CS structures suggested that the binding of these substrates can stabilize the closed format of CS which favors the MDH-CS association. Two other TCA cycle intermediates, ATP, and low pH also enhanced the association of the complex. These results suggest that dynamic formation of the MDH-CS multi-enzyme complex is modulated by metabolic factors responding to respiratory metabolism, and it may function in the feedback regulation of the cycle and adjacent metabolic pathways.

Optimal pH shift of the NADH oxidase from Lactobacillus rhamnosus with a single mutation

OBJECTIVE

To improve the activity of a water-forming NADH oxidase from Lactobacillus rhamnosus under neutral or alkaline pH for coupling NAD⁺-dependent dehydrogenases with an alkaline optimal pH.

RESULTS

The water-forming NADH oxidase from Lactobacillus rhamnosus was engineered by replacing the aspartic acid or glutamic acid with arginine on the surface. The mutant D251R improved the activity with a 112%, 111%, and 244% relative activity to the wild-type at pH 6.5, pH 7.0, and pH 7.5, respectively. Docking substrate into the D251R mutant reveals that the NADH is access to the substrate-binding site with a larger substrate loop due to the enhanced electrostatic repulsion between ARG-251 and ARG-243. In the D251R-NADH complex, the carboxyl of NADH additionally forms two hydrogen bonds (2.6 and 2.9 Å) with G154 due to the changed interaction of substrate and the residues in the catalytic sites, and the hydrogen bond with the oxygen of carbonyl in P295 is shortened from 2.9 to 2.0 Å, which could account for the enhanced specific activity.

CONCLUSIONS

The D251R mutant displayed higher catalytic activity than the wild-type in the pH range 6.5-7.5, and further insight into those shorter and newly formed hydrogen bonds in substrate docking analysis could account for the higher bind affinity and catalytic efficiency of D251R mutant.

On the Case of the Misplaced Hydrogens

Few other elements play a more central role in biology than hydrogen. The interactions, bonding and movement of hydrogen atoms are central to biological catalysis, structure and function. Yet owing to the elusive nature of a single hydrogen atom few experimental and computational techniques can precisely determine its location. This is exemplified in short hydrogen bonds (SHBs) where the location of the hydrogen atom is indicative of the underlying strength of the bonds, which can vary from 1-5 kcal/mol in canonical hydrogen bonds, to an almost covalent nature in single-well hydrogen bonds. Owing to the often-times inferred position of hydrogen, the role of SHBs in biology has remained highly contested and debated. This has also led to discrepancies in computational, biochemical and structural studies of proteins thought to use SHBs in performing chemistry and stabilizing interactions. Herein, we discuss in detail two distinct examples, namely the conserved catalytic triad and the photoreceptor, photoactive yellow protein, where studies of these SHB-containing systems have permitted contextualization of the role these unique hydrogen bonds play in biology.

Low barrier hydrogen bonds in protein structure and function

Low-barrier hydrogen bonds (LBHBs) are a special type of short hydrogen bond (HB) that is characterized by the equal sharing of a hydrogen atom. The existence and catalytic role of LBHBs in proteins has been intensely contested. Advancements in X-ray and neutron diffraction methods has revealed delocalized hydrogen atoms involved in potential LBHBs in a number of proteins, while also demonstrating that short HBs are not necessarily LBHBs. More importantly, a series of experiments on ketosteroid isomerase (KSI) have suggested that LBHBs are significantly stronger than standard HBs in the protein microenvironment in terms of enthalpy, but not free energy. The discrepancy between the enthalpy and free energy of LBHBs offers clues to the challenges, and potential solutions, of the LBHB debate, where the unique strength of LBHBs plays a special role in the kinetic processes of enzyme function and structure, together with other molecular forces in a pre-organized environment.

Identification of the active site residues in ATP-citrate lyase's carboxy-terminal portion

ATP-citrate lyase (ACLY) catalyzes production of acetyl-CoA and oxaloacetate from CoA and citrate using ATP. In humans, this cytoplasmic enzyme connects energy metabolism from carbohydrates to the production of lipids. In certain bacteria, ACLY is used to fix carbon in the reductive tricarboxylic acid cycle. The carboxy(C)-terminal portion of ACLY shows sequence similarity to citrate synthase of the tricarboxylic acid cycle. To investigate the roles of residues of ACLY equivalent to active site residues of citrate synthase, these residues in ACLY from Chlorobium limicola were mutated, and the proteins were investigated using kinetics assays and biophysical techniques. To obtain the crystal structure of the C-terminal portion of ACLY, full-length C. limicola ACLY was cleaved, first non-specifically with chymotrypsin and subsequently with Tobacco Etch Virus protease. Crystals of the C-terminal portion diffracted to high resolution, providing structures that show the positions of active site residues and how ACLY tetramerizes.

A Combined N-terminomics and Shotgun Proteomics Approach to Investigate the Responses of Human Cells to Rapamycin and Zinc at the Mitochondrial Level

All but thirteen mammalian mitochondrial proteins are encoded by the nuclear genome, translated in the cytosol and then imported into the mitochondria. For a significant proportion of the mitochondrial proteins, import is coupled with the cleavage of a presequence called the transit peptide, and the formation of a new N-terminus. Determination of the neo N-termini has been investigated by proteomic approaches in several systems, but generally in a static way to compile as many N-termini as possible. In the present study, we have investigated how the mitochondrial proteome and N-terminome react to chemical stimuli that alter mitochondrial metabolism, namely zinc ions and rapamycin. To this end, we have used a strategy that analyzes both internal and N-terminal peptides in a single run, the dN-TOP approach. We used these two very different stressors to sort out what could be a generic response to stress and what is specific to each of these stressors. Rapamycin and zinc induced different changes in the mitochondrial proteome. However, convergent changes to key mitochondrial enzymatic activities such as pyruvate dehydrogenase, succinate dehydrogenase and citrate synthase were observed for both treatments. Other convergent changes were seen in components of the N-terminal processing system and mitochondrial proteases. Investigations into the generation of neo-N-termini in mitochondria showed that the processing system is robust, as indicated by the lack of change in neo N-termini under the conditions tested. Detailed analysis of the data revealed that zinc caused a slight reduction in the efficiency of the N-terminal trimming system and that both treatments increased the degradation of mitochondrial proteins. In conclusion, the use of this combined strategy allowed a detailed analysis of the dynamics of the mitochondrial N-terminome in response to treatments which impact the mitochondria.

An active site-tail interaction in the structure of hexahistidine-tagged Thermoplasma acidophilum citrate synthase

Citrate synthase (CS) plays a central metabolic role in aerobes and many other organisms. The CS reaction comprises two half-reactions: a Claisen aldol condensation of acetyl-CoA (AcCoA) and oxaloacetate (OAA) that forms citryl-CoA (CitCoA), and CitCoA hydrolysis. Protein conformational changes that `close' the active site play an important role in the assembly of a catalytically competent condensation active site. CS from the thermoacidophile Thermoplasma acidophilum (TpCS) possesses an endogenous Trp fluorophore that can be used to monitor the condensation reaction. The 2.2 Å resolution crystal structure of TpCS fused to a C-terminal hexahistidine tag (TpCSH6) reported here is an `open' structure that, when compared with several liganded TpCS structures, helps to define a complete path for active-site closure. One active site in each dimer binds a neighboring His tag, the first nonsubstrate ligand known to occupy both the AcCoA and OAA binding sites. Solution data collectively suggest that this fortuitous interaction is stabilized by the crystalline lattice. As a polar but almost neutral ligand, the active site-tail interaction provides a new starting point for the design of bisubstrate-analog inhibitors of CS.

Contribution of a low-barrier hydrogen bond to catalysis is not significant in ketosteroid isomerase

Low-barrier hydrogen bonds (LBHBs) have been proposed to have important influences on the enormous reaction rate increases achieved by many enzymes. Δ(5)-3-ketosteroid isomerase (KSI) catalyzes the allylic isomerization of Δ(5)-3-ketosteroid to its conjugated Δ(4)-isomers at a rate that approaches the diffusion limit. Tyr14, a catalytic residue of KSI, has been hypothesized to form an LBHB with the oxyanion of a dienolate steroid intermediate generated during the catalysis. The unusual chemical shift of a proton at 16.8 ppm in the nuclear magnetic resonance spectrum has been attributed to an LBHB between Tyr14 Oη and C3-O of equilenin, an intermediate analogue, in the active site of D38N KSI. This shift in the spectrum was not observed in Y30F/Y55F/D38N and Y30F/Y55F/Y115F/D38N mutant KSIs when each mutant was complexed with equilenin, suggesting that Tyr14 could not form LBHB with the intermediate analogue in these mutant KSIs. The crystal structure of Y30F/Y55F/Y115F/D38N-equilenin complex revealed that the distance between Tyr14 Oη and C3-O of the bound steroid was within a direct hydrogen bond. The conversion of LBHB to an ordinary hydrogen bond in the mutant KSI reduced the binding affinity for the steroid inhibitors by a factor of 8.1-11. In addition, the absence of LBHB reduced the catalytic activity by only a factor of 1.7-2. These results suggest that the amount of stabilization energy of the reaction intermediate provided by LBHB is small compared with that provided by an ordinary hydrogen bond in KSI.

Using unnatural amino acids to probe the energetics of oxyanion hole hydrogen bonds in the ketosteroid isomerase active site

Hydrogen bonds are ubiquitous in enzyme active sites, providing binding interactions and stabilizing charge rearrangements on substrate groups over the course of a reaction. But understanding the origin and magnitude of their catalytic contributions relative to hydrogen bonds made in aqueous solution remains difficult, in part because of complexities encountered in energetic interpretation of traditional site-directed mutagenesis experiments. It has been proposed for ketosteroid isomerase and other enzymes that active site hydrogen bonding groups provide energetic stabilization via "short, strong" or "low-barrier" hydrogen bonds that are formed due to matching of their pKa or proton affinity to that of the transition state. It has also been proposed that the ketosteroid isomerase and other enzyme active sites provide electrostatic environments that result in larger energetic responses (i.e., greater "sensitivity") to ground-state to transition-state charge rearrangement, relative to aqueous solution, thereby providing catalysis relative to the corresponding reaction in water. To test these models, we substituted tyrosine with fluorotyrosines (F-Tyr's) in the ketosteroid isomerase (KSI) oxyanion hole to systematically vary the proton affinity of an active site hydrogen bond donor while minimizing steric or structural effects. We found that a 40-fold increase in intrinsic F-Tyr acidity caused no significant change in activity for reactions with three different substrates. F-Tyr substitution did not change the solvent or primary kinetic isotope effect for proton abstraction, consistent with no change in mechanism arising from these substitutions. The observed shallow dependence of activity on the pKa of the substituted Tyr residues suggests that the KSI oxyanion hole does not provide catalysis by forming an energetically exceptional pKa-matched hydrogen bond. In addition, the shallow dependence provides no indication of an active site electrostatic environment that greatly enhances the energetic response to charge accumulation, consistent with prior experimental results.

The mechanism of citryl-coenzyme A formation catalyzed by citrate synthase

The enzyme citrate synthase is used by all living cells to catalyze the first step of the citric acid cycle. In this work, we have investigated the enolization and condensation steps catalyzed by citrate synthase, using ab initio (B3LYP/def2-TZVP and MP2/aug-cc-pVDZ) quantum chemical/molecular mechanical hybrid potentials in conjunction with reaction-path-location algorithms and molecular dynamics free energy simulations. The results of the latter indicate that the catalytic His238 residue is in its neutral form, and also argue strongly for the presence of a water molecule in the enzyme's catalytic center. Such a water is observed in some, but not all, of the experimentally resolved structures of the protein. The mechanism itself starts with an enolization that proceeds via an enolate intermediate rather than the enol form, which is much more unstable. This is in agreement with the results of other workers. For the condensation step, we investigated two mechanisms in which there is a direct nucleophilic attack of the enolate intermediate on the oxaloacetate carbonyl carbon, and found the one in which there is no proton transfer from the neighboring arginine to be preferred. Although this residue, Arg329, is not implicated directly in the reaction, it helps to stabilize the negative citryl-CoA formed during the condensation step.

Enzyme-substrate complexes of allosteric citrate synthase: evidence for a novel intermediate in substrate binding

The citrate synthase (CS) of Escherichia coli is an allosteric hexameric enzyme specifically inhibited by NADH. The crystal structure of wild type (WT) E. coli CS, determined by us previously, has no substrates bound, and part of the active site is in a highly mobile region that is shifted from the position needed for catalysis. The CS of Acetobacter aceti has a similar structure, but has been successfully crystallized with bound substrates: both oxaloacetic acid (OAA) and an analog of acetyl coenzyme A (AcCoA). We engineered a variant of E. coli CS wherein five amino acids in the mobile region have been replaced by those in the A. aceti sequence. The purified enzyme shows unusual kinetics with a low affinity for both substrates. Although the crystal structure without ligands is very similar to that of the WT enzyme (except in the mutated region), complexes are formed with both substrates and the allosteric inhibitor NADH. The complex with OAA in the active site identifies a novel OAA-binding residue, Arg306, which has no functional counterpart in other known CS-OAA complexes. This structure may represent an intermediate in a multi-step substrate binding process where Arg306 changes roles from OAA binding to AcCoA binding. The second complex has the substrate analog, S-carboxymethyl-coenzyme A, in the allosteric NADH-binding site and the AcCoA site is not formed. Additional CS variants unable to bind adenylates at the allosteric site show that this second complex is not a factor in positive allosteric activation of AcCoA binding.

Inhibition of human α-methylacyl CoA racemase (AMACR): a target for prostate cancer

The enzyme α-methylacyl CoA racemase (AMACR) is involved in the metabolism of branched-chain fatty acids and has been identified as a promising therapeutic target for prostate cancer. By using the recently available human AMACR from HEK293 kidney cell cultures, we tested a series of new rationally designed inhibitors to determine the structural requirements in the acyl component. An N-methylthiocarbamate (Ki=98 nM), designed to mimic the proposed enzyme-bound enolate, was found to be the most potent AMACR inhibitor reported to date.

Evolution of the genetic code by incorporation of amino acids that improved or changed protein function

Fifty years have passed since the genetic code was deciphered, but how the genetic code came into being has not been satisfactorily addressed. It is now widely accepted that the earliest genetic code did not encode all 20 amino acids found in the universal genetic code as some amino acids have complex biosynthetic pathways and likely were not available from the environment. Therefore, the genetic code evolved as pathways for synthesis of new amino acids became available. One hypothesis proposes that early in the evolution of the genetic code four amino acids-valine, alanine, aspartic acid, and glycine-were coded by GNC codons (N = any base) with the remaining codons being nonsense codons. The other sixteen amino acids were subsequently added to the genetic code by changing nonsense codons into sense codons for these amino acids. Improvement in protein function is presumed to be the driving force behind the evolution of the code, but how improved function was achieved by adding amino acids has not been examined. Based on an analysis of amino acid function in proteins, an evolutionary mechanism for expansion of the genetic code is described in which individual coded amino acids were replaced by new amino acids that used nonsense codons differing by one base change from the sense codons previously used. The improved or altered protein function afforded by the changes in amino acid function provided the selective advantage underlying the expansion of the genetic code. Analysis of amino acid properties and functions explains why amino acids are found in their respective positions in the genetic code.

A hybrid elastic band string algorithm for studies of enzymatic reactions

A common challenge in theoretical biophysics is the identification of a minimum energy path (MEP) for the rearrangement of a group of atoms from one stable configuration to another. The structure with maximum energy along the MEP approximates the transition state for the process and the energy profile itself permits estimation of the transition rates. In this work we describe a computationally efficient algorithm for the identification of minimum energy paths in complicated biosystems. The algorithm is a hybrid of the nudged elastic band (NEB) and string methods. It has been implemented in the pDynamo simulation program and tested by examining elementary steps in the reaction mechanisms of three enzymes: citrate synthase, RasGAP, and lactate dehydrogenase. Good agreement is found for the energies and geometries of the species along the reaction profiles calculated using the new algorithm and previous versions of the NEB and string techniques, and also those obtained by the common method of adiabatic exploration of the potential energy surface as a function of predefined reaction coordinates. Precisely refined structures of the saddle points along the paths may be subsequently obtained with the climbing image variant of the NEB algorithm. Directions in which the utility of the methods that we have implemented can be further improved are discussed.

Direct carbon-carbon bond formation via soft enolization: a biomimetic asymmetric Mannich reaction of phenylacetate thioesters

An asymmetric Mannich reaction of phenylacetate thioesters and sulfonylimines using cinchona alkaloid-based amino (thio)urea catalysts is reported that employs proximity-assisted soft enolization. This approach to enolization is based on the cooperative action of a carbonyl-activating hydrogen bonding (thio)urea moiety and an amine base contained within a single catalytic entity to facilitate intracomplex deprotonation. Significantly, this allows thioesters over a range of acidity to react efficiently, thereby opening the door to the development of a general mode of enolization-based organocatalysis of monocarboxylic acid derivatives.

The partial substrate dethiaacetyl-coenzyme A mimics all critical carbon acid reactions in the condensation half-reaction catalyzed by Thermoplasma acidophilum citrate synthase

Citrate synthase (CS) performs two half-reactions: the mechanistically intriguing condensation of acetyl-CoA with oxaloacetate (OAA) to form citryl-CoA and the subsequent, slower hydrolysis of citryl-CoA that generally dominates steady-state kinetics. The condensation reaction requires the abstraction of a proton from the methyl carbon of acetyl-CoA to generate a reactive enolate intermediate. The carbanion of that intermediate then attacks the OAA carbonyl to furnish citryl-CoA, the initial product. Using stopped-flow and steady-state fluorescence methods, kinetic substrate isotope effects, and mutagenesis of active site residues, we show that all of the processes that occur in the condensation half-reaction performed by Thermoplasma acidophilum citrate synthase (TpCS) with the natural thioester substrate, acetyl-CoA, also occur with the ketone inhibitor dethiaacetyl-CoA. Free energy profiles demonstrate that the nonhydrolyzable product of the condensation reaction, dethiacitryl-CoA, forms a particularly stable complex with TpCS but not pig heart CS.

Substrate polarization in enzyme catalysis: QM/MM analysis of the effect of oxaloacetate polarization on acetyl-CoA enolization in citrate synthase

Citrate synthase is an archetypal carbon-carbon bond forming enzyme. It promotes the conversion of oxaloacetate (OAA) to citrate by catalyzing the deprotonation (enolization) of acetyl-CoA, followed by nucleophilic attack of the enolate form of this substrate on OAA to form a citryl-CoA intermediate and subsequent hydrolysis. OAA is strongly bound to the active site and its alpha-carbonyl group is polarized. This polarization has been demonstrated spectroscopically, [(Kurz et al., Biochemistry 1985;24:452-457; Kurz and Drysdale, Biochemistry 1987;26:2623-2627)] and has been suggested to be an important catalytic strategy. Substrate polarization is believed to be important in many enzymes. The first step, formation of the acetyl-CoA enolate intermediate, is thought to be rate-limiting in the mesophilic (pig/chicken) enzyme. We have examined the effects of substrate polarization on this key step using quantum mechanical/molecular mechanical (QM/MM) methods. Free energy profiles have been calculated by AM1/CHARMM27 umbrella sampling molecular dynamics (MD) simulations, together with potential energy profiles. To study the influence of OAA polarization, profiles were calculated with different polarization of the OAA alpha-carbonyl group. The results indicate that OAA polarization influences catalysis only marginally but has a larger effect on intermediate stabilization. Different levels of treatment of OAA are compared (MM or QM), and its polarization in the protein and in water analyzed at the B3LYP/6-31+G(d)/CHARMM27 level. Analysis of stabilization by individual residues shows that the enzyme mainly stabilizes the enolate intermediate (not the transition state) through electrostatic (including hydrogen bond) interactions: these contribute much more than polarization of OAA.

How accurate and statistically robust are catalytic site predictions based on closeness centrality?

Background

We examine the accuracy of enzyme catalytic residue predictions from a network representation of protein structure. In this model, amino acid α-carbons specify vertices within a graph and edges connect vertices that are proximal in structure. Closeness centrality, which has shown promise in previous investigations, is used to identify important positions within the network. Closeness centrality, a global measure of network centrality, is calculated as the reciprocal of the average distance between vertex i and all other vertices.

Results

We benchmark the approach against 283 structurally unique proteins within the Catalytic Site Atlas. Our results, which are inline with previous investigations of smaller datasets, indicate closeness centrality predictions are statistically significant. However, unlike previous approaches, we specifically focus on residues with the very best scores. Over the top five closeness centrality scores, we observe an average true to false positive rate ratio of 6.8 to 1. As demonstrated previously, adding a solvent accessibility filter significantly improves predictive power; the average ratio is increased to 15.3 to 1. We also demonstrate (for the first time) that filtering the predictions by residue identity improves the results even more than accessibility filtering. Here, we simply eliminate residues with physiochemical properties unlikely to be compatible with catalytic requirements from consideration. Residue identity filtering improves the average true to false positive rate ratio to 26.3 to 1. Combining the two filters together has little affect on the results. Calculated p-values for the three prediction schemes range from 2.7E-9 to less than 8.8E-134. Finally, the sensitivity of the predictions to structure choice and slight perturbations is examined.

Conclusion

Our results resolutely confirm that closeness centrality is a viable prediction scheme whose predictions are statistically significant. Simple filtering schemes substantially improve the method's predicted power. Moreover, no clear effect on performance is observed when comparing ligated and unligated structures. Similarly, the CC prediction results are robust to slight structural perturbations from molecular dynamics simulation.

Ab initio QM/MM modelling of acetyl-CoA deprotonation in the enzyme citrate synthase

The first step of the reaction catalysed by the enzyme citrate synthase is studied here with high level combined quantum mechanical/molecular mechanical (QM/MM) methods (up to the MP2/6-31+G(d)//6-31G(d)/CHARMM level). In the first step of the reaction, acetyl-CoA is deprotonated by Asp375, producing an intermediate, which is the nucleophile for attack on the second substrate, oxaloacetate, prior to hydrolysis of the thioester bond of acetyl-CoA and release of the products. A central question has been whether the nucleophilic intermediate is the enolate of acetyl-CoA, the enol, or an 'enolic' intermediate stabilized by a 'low-barrier' hydrogen bond with His274 at the active site. The imidazole sidechain of His274 is neutral, and donates a hydrogen bond to the carbonyl oxygen of acetyl-CoA in substrate complexes. We have investigated the identity of the nucleophilic intermediate by QM/MM calculations on the substrate (keto), enolate, enol and enolic forms of acetyl-CoA at the active site of citrate synthase. The transition states for proton abstraction from acetyl-CoA by Asp375, and for transfer of the hydrogen bonded proton between His274 and acetyl-CoA have been modelled approximately. The effects of electron correlation are included by MP2/6-31G(d) and MP2/6-31+G(d) calculations on active site geometries produced by QM/MM energy minimization. The results do not support the hypothesis that a low-barrier hydrogen bond is involved in catalysis in citrate synthase, in agreement with earlier calculations. The acetyl-CoA enolate is identified as the only intermediate consistent with the experimental barrier for condensation, stabilized by conventional hydrogen bonds from His274 and a water molecule.

Transition state modeling and catalyst design for hydrogen bond-stabilized enolate formation

[reaction: see text] A catalyst for enolate formation was designed that incorporates an amine base along with a thiourea to bind to the oxygen atom of the substrate and enolate through hydrogen bonding. A computational model of the transition state was developed in which the thiourea (modeled initially as a urea) and amine were separate molecules. This model and models incorporating one or two methanol molecules in place of the urea showed an out-of-plane hydrogen bond, apparently to the carbonyl pi-bond, in addition to an in-plane hydrogen bond to an unshared electron pair. In contrast, optimized complexes of the ketone and the fully formed enolate showed only in-plane hydrogen bonding. The transition state model with the urea and amine was used to define a database search with the computer program CAVEAT to identify structures suitable for linking the amine and urea/thiourea moieties in the transition state. On the basis of a group of structures identified from this search, a flexible but conformationally biased linker was designed to connect the two catalytic moieties. The molecule having the amine and thiourea moieties connected by this linker was synthesized and was shown to catalyze proton exchange between methanol and deuterated acetone. The catalyst was about 5-fold more efficient than the amine and thiourea as separate molecules and relative to a similar but less conformationally biased catalyst.

An intramolecular ionic hydrogen bond stabilizes a cis amide bond rotamer of a ring-opened rapamycin-degradation product

Rapamycin (1), a macrolide immunosuppressant, undergoes degradation into ring-opened acid products 2 and 3 under physiologically relevant conditions. The unsaturated product (3) was isolated and studied in this work. Unlike 1, which has its amide primarily in a trans conformation in solution, 3 has both cis and trans conformations in approximately a 1:1 ratio in dimethyl sulfoxide (DMSO). The amount of cis rotamer was increased dramatically in the presence of an organic base such as triethylamine. The detailed NMR results indicate that the cis rotamer is stabilized through an intramolecular ionic hydrogen bond of the carboxylate anion with the tertiary alcohol as part of a nine-membered ring system. This hydrogen bond was characterized further in organic media and the trans-cis rotamer equilibria were used to estimate the relative bond strengths in several solvents. The additional stabilization arising from this ionic hydrogen bond in the cis rotamer was determined to be 1.4 kcal mol(-1) in DMSO-d6, 2.0 kcal mol(-1) in CD3CN and 1.1 kcal mol(-1) in CD3OD.

The importance of the amide bond nearest the thiol group in enzymatic reactions of coenzyme A

Analogues of coenzyme A (CoA) and of CoA thioesters have been prepared in which the amide bond nearest the thiol group has been modified. An analogue of acetyl-CoA in which this amide bond is replaced with an ester linkage was a good substrate for the enzymes carnitine acetyltransferase, chloramphenicol acetyltransferase, and citrate synthase, with K(m) values 2- to 8-fold higher than those of acetyl-CoA and V(max) values from 14 to >80% those of the natural substrate. An analogue in which an extra methylene group was inserted between the amide bond and the thiol group showed less than 4-fold diminished binding to the three enzymes but exhibited less than 1% activity relative to acetyl-CoA with carnitine acetyltransferase and no measurable activity with the other two enzymes. Analogues of several CoA thioesters in which the amide bond was replaced with a hemithioacetal linkage exhibited no measurable activity with the appropriate enzymes. The results indicate that some aspects of the amide bond and proper distance between this amide and the thiol/thioester moiety are critical for activity of CoA ester-utilizing enzymes.

Different roles of electrostatics in heat and in cold: adaptation by citrate synthase

Electrostatics plays a major role in heat adaptation by thermophilic proteins. Here we ask whether electrostatics similarly contributes to cold adaptation in psychrophilic proteins. We compare the sequences and structures of citrate synthases from the psychrophile Arthobacter Ds2-3R, from chicken, and from the hyperthermophile Pyrococcus furiosus. The three enzymes share similar packing, burial of nonpolar surface area, and main-chain hydrogen bonding. However, both psychrophilic and hyperthermophilic citrate synthases contain more charged residues, salt bridges, and salt-bridge networks than the mesophile. The electrostatic free-energy contributions toward protein stability by individual charged residues show greater variabilities in the psychrophilic citrate synthase than in the hyperthermophilic enzyme. The charged residues in the active-site regions of the psychrophile are more destabilizing than those in the active-site regions of the hyperthermophile. In the hyperthermophilic enzyme, salt bridges and their networks largely cluster in the active-site regions and at the dimer interface. In contrast, in the psychrophile, they are more dispersed throughout the structure. On average, salt bridges and their networks provide greater electrostatic stabilization to the thermophilic citrate synthase at 100 degrees C than to the psychrophilic enzyme at 0 degrees C. Electrostatics appears to play an important role in both heat and cold adaptation of citrate synthase. However, remarkably, the role may be different in the two types of enzyme: In the hyperthermophile, it may contribute to the integrity of both the protein dimer and the active site by possibly countering conformational disorder at high temperatures. On the other hand, in the psychrophile at low temperatures, electrostatics may contribute to enhance protein solvation and to ensure active-site flexibility.

Structure of the Escherichia coli malate synthase G:pyruvate:acetyl-coenzyme A abortive ternary complex at 1.95 A resolution

Malate synthase, an enzyme of the glyoxylate pathway, catalyzes the condensation and subsequent hydrolysis of acetyl-coenzyme A (acetyl-CoA) and glyoxylate to form malate and CoA. In the present study, we present the 1.95 A-resolution crystal structure of Escherichia coli malate synthase isoform G in complex with magnesium, pyruvate, and acetyl-CoA, and we compare it with previously determined structures of substrate and product complexes. The results reveal how the enzyme recognizes and activates the substrate acetyl-CoA, as well as conformational changes associated with substrate binding, which may be important for catalysis. On the basis of these results and mutagenesis of active site residues, Asp 631 and Arg 338 are proposed to act in concert to form the enolate anion of acetyl-CoA in the rate-limiting step. The highly conserved Cys 617, which is immediately adjacent to the presumed catalytic base Asp 631, appears to be oxidized to cysteine-sulfenic acid. This can explain earlier observations of the susceptibility of the enzyme to inactivation and aggregation upon X-ray irradiation and indicates that cysteine oxidation may play a role in redox regulation of malate synthase activity in vivo. There is mounting evidence that enzymes of the glyoxylate pathway are virulence factors in several pathogenic organisms, notably Mycobacterium tuberculosis and Candida albicans. The results described in this study add insight into the mechanism of catalysis and may be useful for the design of inhibitory compounds as possible antimicrobial agents.

A low-barrier hydrogen bond between histidine of secreted phospholipase A2 and a transition state analog inhibitor

This work describes in-depth NMR characterization of a unique low-barrier hydrogen bond (LBHB) between an active site residue from the enzyme and a bound inhibitor: the complex between secreted phospholipase A(2) (sPLA(2), from bee venom and bovine pancreas) and a transition-state analog inhibitor HK32. A downfield proton NMR resonance, at 17-18 ppm, was observed in the complex but not in the free enzyme. On the basis of site-specific mutagenesis and specific 15N-decoupling, this downfield resonance was assigned to the active site H48, which is part of the catalytic dyad D99-H48. These results led to a hypothesis that the downfield resonance represents the proton (H(epsilon 2) of H48) involved in the H-bonding between D99 and H48, in analogy with serine proteases. However, this was shown not to be the case by use of the bovine enzyme labeled with specific [15N(epsilon 2)]His. Instead, the downfield resonance arises from H(delta1) of H48, which forms a hydrogen bond with a non-bridging phosphonate oxygen of the inhibitor. Further studies showed that this proton displays a fractionation factor of 0.62(+/-0.06), and an exchange rate protection factor of >100 at 285 K and >40 at 298 K, which are characteristic of a LBHB. The pK(a) of the imidazole ring of H48 was shown to be shifted from 5.7 for the free enzyme to an apparent value of 9.0 in the presence of the inhibitor. These properties are very similar to those of the Asp em leader His LBHBs in serine proteases. Possible structural bases and functional consequences for the different locations of the LBHB between these two types of enzymes are discussed. The results also underscore the importance of using specific isotope labeling, rather than extrapolation of NMR results from other enzyme systems, to assign the downfield proton resonance to a specific hydrogen bond. Although our studies did not permit the strength of the LBHB to be accurately measured, the data do not provide support for an unusually strong hydrogen bond strength (i.e. >10 kcal/mol).

Elaborate manifold of short hydrogen bond arrays mediating binding of active site-directed serine protease inhibitors

An extensive structural manifold of short hydrogen bond-mediated, active site-directed, serine protease inhibition motifs is revealed in a set of over 300 crystal structures involving a large suite of small molecule inhibitors (2-(2-phenol)-indoles and 2-(2-phenol)-benzimidazoles) determined over a wide range of pH (3.5-11.4). The active site hydrogen-bonding mode was found to vary markedly with pH, with the steric and electronic properties of the inhibitor, and with the type of protease (trypsin, thrombin or urokinase type plasminogen activator (uPA)). The pH dependence of the active site hydrogen-bonding motif is often intricate, constituting a distinct fingerprint of each complex. Isosteric replacements or minor substitutions within the inhibitor that modulate the pK(a) of the phenol hydroxyl involved in short hydrogen bonding, or that affect steric interactions distal to the active site, can significantly shift the pH-dependent structural profile characteristic of the parent scaffold, or produce active site-binding motifs unique to the bound analog. Ionization equilibria at the active site associated with inhibitor binding are probed in a series of the protease-inhibitor complexes through analysis of the pH dependence of the structure and environment of the active site-binding groups involved in short hydrogen bond arrays. Structures determined at high pH (>11), suggest that the pK(a) of His57 is dramatically elevated, to a value as high as approximately 11 in certain complexes. K(i) values involving uPA and trypsin determined as a function of pH for a set of inhibitors show pronounced parabolic pH dependence, the pH for optimal inhibition governed by the pK(a) of the inhibitor phenol involved in short hydrogen bonds. Comparison of structures of trypsin, thrombin and uPA, each bound by the same inhibitor, highlights important structural variations in the S1 and active sites accessible for engineering notable selectivity into remarkably small molecules with low nanomolar K(i) values.

Molecular dynamics and quantum chemical studies on the catalytic mechanism of Delta5-3-ketosteroid isomerase: the catalytic diad versus the cooperative hydrogen bond mechanism

To further understand Delta(5)-3-ketosteroid isomerase (KSI) catalysis, we carried out molecular dynamics (MD) simulations of the KSI dimer ligated with a substrate and reaction intermediate analogue and high level ab initio calculations on relevant enzymatic reaction models. Simulation of the enzyme-substrate complex dimer systems showed asymmetric dynamics between the two monomers, in which the hydrogen bond pattern between the substrate and active site residues in the first and the second subunits supported the cooperative hydrogen bond (CH) and the catalytic diad (CD) mechanisms, respectively. On the other hand, only the CH mechanism was supported in the MD simulation of the enzyme-intermediate complex dimer. From MP2/6-31+G**//RHF/6-31G** calculations, we found the kinetic barriers for the two reaction mechanisms were similar. The CH route afforded a greater stabilization to the enolate intermediate than did the CD counterpart. Thus, the present computational studies indicate that the CH mechanism would be favored over the CD one in the catalytic action of KSI. However, the latter could not be ruled out conclusively because of the explicit appearance of a CD configuration in the MD trajectories of the enzyme-substrate complex and because of the similar intrinsic activation barrier for the CH and CD mechanisms. The appearance of configurations that favor the CD pathway is rationalized in terms of a model in which the KSI-substrate complex does not have a strong preference for one hydrogen bonding pattern over another, while the KSI-intermediate complex favors a cooperative hydrogen bond pattern in order to stabilize the reaction intermediate. This hypothesis is supported by the ab initio calculations which indicate that the CH intermediate is more stable than the CD one by approximately 6.3 kcal/mol.

Contribution of Multicentered Short Hydrogen Bond Arrays to Potency of Active Site-Directed Serine Protease Inhibitors

We describe and compare the pH dependencies of the potencies and of the bound structures of two inhibitor isosteres that form multicentered short hydrogen bond arrays at the active sites of trypsin, thrombin, and urokinase type plasminogen activator (urokinase or uPA) over certain ranges of pH. Depending on the pH, short hydrogen bond arrays at the active site are mediated by two waters, one in the oxyanion hole (H(2)O(oxy)) and one on the other (S2) side of the inhibitor (H(2)O(S2)), by one water (H(2)O(oxy)), or by no water. The dramatic variation in the length of the active site hydrogen bonds as a function of pH, of inhibitor, and of enzyme, along with the involvement or absence of ordered water, produces a large structural manifold of active site hydrogen bond motifs. Diverse examples of multicentered and two-centered short hydrogen bond arrays, both at and away from the active site, recently discovered in several protein crystal systems, suggest that short hydrogen bonds in proteins may be more common than has been recognized. The short hydrogen bond arrays resemble one another with respect to ionic nature, highly polar environment, multitude of associated ordinary hydrogen bonds, and disparate pK(a) values of participating groups. Comparison of structures and K(i) values of trypsin complexes at pH values where the multicentered short hydrogen bond arrays mediating inhibitor binding are present or absent indicate that these arrays have a minor effect on inhibitor potency. These features suggest little covalent nature within the short hydrogen bonds, despite their extraordinary shortness (as short as 2.0 A).

The biosynthesis of coenzyme A in bacteria

Coenzyme A (I) and enzyme-bound phosphopantetheine (II) function as acyl carriers and as carbonyl activating groups for Claisen reactions as well as for amide-, ester-, and thioester-forming reactions in the cell. In so doing, these cofactors play a key role in the biosynthesis and breakdown of fatty acids and in the biosynthesis of polyketides and nonribosomal peptides. Coenzyme A is biosynthesized in bacteria in nine steps. The biosynthesis begins with the decarboxylation of aspartate to give beta-alanine. Pantoic acid is formed by the hydroxymethylation of alpha-ketoisovalerate followed by reduction. These intermediates are then condensed to give pantothenic acid. Phosphorylation of pantothenic acid followed by condensation with cysteine and decarboxylation gives 4'-phosphopantetheine. Adenylation and phosphorylation of 4'-phosphopantetheine completes the biosynthesis of coenzyme A. This review will focus on the mechanistic enzymology of coenzyme A biosynthesis in bacteria.

A novel serine protease inhibition motif involving a multi-centered short hydrogen bonding network at the active site

We describe a new serine protease inhibition motif in which binding is mediated by a cluster of very short hydrogen bonds (<2.3 A) at the active site. This protease-inhibitor binding paradigm is observed at high resolution in a large set of crystal structures of trypsin, thrombin, and urokinase-type plasminogen activator (uPA) bound with a series of small molecule inhibitors (2-(2-phenol)indoles and 2-(2-phenol)benzimidazoles). In each complex there are eight enzyme-inhibitor or enzyme-water-inhibitor hydrogen bonds at the active site, three of which are very short. These short hydrogen bonds connect a triangle of oxygen atoms comprising O(gamma)(Ser195), a water molecule co-bound in the oxyanion hole (H(2)O(oxy)), and the phenolate oxygen atom of the inhibitor (O6'). Two of the other hydrogen bonds between the inhibitor and active site of the trypsin and uPA complexes become short in the thrombin counterparts, extending the three-centered short hydrogen-bonding array into a tetrahedral array of atoms (three oxygen and one nitrogen) involved in short hydrogen bonds. In the uPA complexes, the extensive hydrogen-bonding interactions at the active site prevent the inhibitor S1 amidine from forming direct hydrogen bonds with Asp189 because the S1 site is deeper in uPA than in trypsin or thrombin. Ionization equilibria at the active site associated with inhibitor binding are probed through determination and comparison of structures over a wide range of pH (3.5 to 11.4) of thrombin complexes and of trypsin complexes in three different crystal forms. The high-pH trypsin-inhibitor structures suggest that His57 is protonated at pH values as high as 9.5. The pH-dependent inhibition of trypsin, thrombin, uPA and factor Xa by 2-(2-phenol)benzimidazole analogs in which the pK(a) of the phenol group is modulated is shown to be consistent with a binding process involving ionization of both the inhibitor and the enzyme. These data further suggest that the pK(a) of His57 of each protease in the unbound state in solution is about the same, approximately 6.8. By comparing inhibition constants (K(i) values), inhibitor solubilities, inhibitor conformational energies and corresponding structures of short and normal hydrogen bond-mediated complexes, we have estimated the contribution of the short hydrogen bond networks to inhibitor affinity ( approximately 1.7 kcal/mol). The structures and K(i) values associated with the short hydrogen-bonding motif are compared with those corresponding to an alternate, Zn(2+)-mediated inhibition motif at the active site. Structural differences among apo-enzymes, enzyme-inhibitor and enzyme-inhibitor-Zn(2+) complexes are discussed in the context of affinity determinants, selectivity development, and structure-based inhibitor design.

A QM/MM study of the racemization of vinylglycolate catalyzed by mandelate racemase enzyme

The experimentally postulated mechanism for the interconversion between (S)-vinylglycolate and (R)-vinylglycolate catalyzed by mandelate racemase enzyme consists of a two-step quite symmetric process through a dianionic enolic intermediate that is formed after the abstraction of the alpha-proton of vinylglycolate by a basic enzymatic residue and is then reprotonated by another residue. The challenging problem behind this reaction is how the enzyme manages to stabilize such an intermediate, that is, how it lowers enough the high pK(a) of the alpha-proton for the reaction to take place. The QM/MM simulations performed in this paper indicate that catalysis is based on the stabilization of the negative charge developed on the substrate along the reaction. We have identified three different reaction mechanisms starting from different quasi-degenerate structures of the substrate-enzyme complex. In two of them the stabilizing role is done by means of a catalytic proton transfer that avoids the formation of a dianionic intermediate, and they involve six steps instead of the two experimentally proposed. On the contrary, the third mechanism passes through a dianionic species stabilized by the concerted approach of a protonated enzymatic residue during the proton abstraction. The potential energy barriers theoretically found along these mechanisms are qualitatively in good agreement with the experimental free energy barriers determined for racemization of vinylglycolate and mandelate. The theoretical study of the effect of the mutation of Glu317 by Gln317 in the kinetics of the reaction reveals the important role in the catalysis of the hydrogen bond formed by Glu317 in the native enzyme, as only one of the mechanisms, the slower one, is able to produce the racemization in the active site of the mutant. However, we have found that this hydrogen bond is not an LBHB within our model.

Alternate modes of binding in two crystal structures of alkaline phosphatase-inhibitor complexes

Two high resolution crystal structures of Escherichia coli alkaline phosphatase (AP) in the presence of phosphonate inhibitors are reported. The phosphonate compounds, phosphonoacetic acid (PAA) and mercaptomethylphosphonic acid (MMP), bind competitively to AP with dissociation constants of 5.5 and 0.6 mM, respectively. The structures of the complexes of AP with PAA and MMP were refined at high resolution to crystallographic R-values of 19.0 and 17.5%, respectively. Refinement of the AP-inhibitor complexes was carried out using X-PLOR. The final round of refinement was done using SHELXL-97. Crystallographic analyses of the inhibitor complexes reveal different binding modes for the two phosphonate compounds. The significant difference in binding constants can be attributed to these alternative binding modes observed in the high resolution X-ray structures. The phosphinyl group of PAA coordinates to the active site zinc ions in a manner similar to the competitive inhibitor and product inorganic phosphate. In contrast, MMP binds with its phosphonate moiety directed toward solvent. Both enzyme-inhibitor complexes exhibit close contacts, one of which has the chemical and geometrical potential to be considered an unconventional hydrogen bond of the type C-H...X.

Crystallographic analysis of the reaction pathway of Zoogloea ramigera biosynthetic thiolase

Biosynthetic thiolases catalyze the biological Claisen condensation of two acetyl-CoA molecules to form acetoacetyl-CoA. This is one of the fundamental categories of carbon skeletal assembly patterns in biological systems and is the first step in many biosynthetic pathways including those which generate cholesterol, steroid hormones and ketone body energy storage molecules. High resolution crystal structures of the tetrameric biosynthetic thiolase from Zoogloea ramigera were determined (i) in the absence of active site ligands, (ii) in the presence of CoA, and (iii) from protein crystals which were flash frozen after a short soak with acetyl-CoA, the enzyme's substrate in the biosynthetic reaction. In the latter structure, a reaction intermediate was trapped: the enzyme was found to be acetylated at Cys89 and a molecule of acetyl-CoA was bound in the active site pocket. A comparison of the three new structures and the two previously published thiolase structures reveals that small adjustments in the conformation of the acetylated Cys89 side-chain allow CoA and acetyl-CoA to adopt identical modes of binding. The proximity of the acetyl moiety of acetyl-CoA to the sulfur atom of Cys378 supports the hypothesis that Cys378 is important for proton exchange in both steps of the reaction. The thioester oxygen atom of the acetylated enzyme points into an oxyanion hole formed by the nitrogen atoms of Cys89 and Gly380, thus facilitating the condensation reaction. The interaction between the thioester oxygen atom of acetyl-CoA and His348 assists the condensation step of catalysis by stabilizing a negative charge on the thioester oxygen atom. Our structure of acetyl-CoA bound to thiolase also highlights the importance in catalysis of a hydrogen bonding network between Cys89 and Cys378, which includes the thioester oxygen atom of acetyl-CoA, and extends from the catalytic site through the enzyme to the opposite molecular surface. This hydrogen bonding network is different in yeast degradative thiolase, indicating that the catalytic properties of each enzyme may be modulated by differences in their hydrogen bonding networks.

Hydrogen bonding in enzymatic catalysis: analysis of energetic contributions

Enzymes can provide catalysis by increasing the strengthening of hydrogen bonds to groups undergoing charge rearrangement in the course of reaction relative to the strengthening of the hydrogen bonds in the corresponding solution reactions. This can be accomplished by using hydrogen bond donors and acceptors that are stronger than water and by lowering the effective dielectric relative to that in aqueous solution. We suggest that these electrostatic effects are of general significance in enzymatic catalysis. The effective dielectric is lowered by the overall "rigidity" of the folded enzyme, which facilitates the formation of active site interactions, and by the fixation of active site functional groups within the enzyme x substrate complex. This underscores the fundamental interconnection of catalytic mechanisms in enzymatic catalysis.

Development of a Second Generation Coenzyme A Analogue Synthon

We have previously reported a general synthetic approach to analogues of coenzyme A (CoA) which involves enzymatic synthesis of a general CoA analogue synthon having a thioester linkage in place of the amide bond nearest the thiol group (Martin et al. J. Am. Chem. Soc. 1994, 116, 4660). We report here the synthesis of a second CoA analogue synthon 1c which has the amide bond more distant from the thiol group replaced with a thioester. This analogue was prepared by nonenzymatic synthesis of a racemic phosphopantetheine analogue followed by enzymatic conversion to the corresponding CoA analogue. Stereochemical analysis showed that the natural enantiomer of the phosphopantetheine analogue was selectively converted to product by the enzyme phosphopantetheine adenylyltansferase, yielding a product that possessed the desired stereoconfiguration. Reaction of the new synthon 1c with a primary amine results in amide bond formation to form the CoA analogue of interest. This new methodology provides access to an even broader array of CoA analogues modified in the beta-alanylcysteamine moiety. This has been demonstrated in the synthesis of an analogue having an extra methylene group in the beta-alanine moiety and two analogues in which the amide bond nearest the thiol group is replaced with a pair of methylene groups.

Amino acid substitutions in the subunit interface enhancing thermostability of Thermoplasma acidophilum citrate synthase

We have used citrate synthase from Thermoplasma (Tp.) acidophilum as a thermostable model system to investigate the role of hydrophobic interactions in dimer interface for maintaining high temperature stability. Three mutant enzymes were constructed by single amino acid substitutions in the interface helices: Ala97-->Ser, Ala104-->Thr, and Gly209-->Ala. All of the mutations enhanced the thermostability of Tp. citrate synthase, while improving its catalytic properties (Km, Vmax, and specific activity). The highest thermostability was achieved by the Gly209-->Ala substitution. The half-life of irreversible inactivation of the G209A mutant enzyme at 85 degreesC was about 57 min, and the midpoint of guanidinium chloride (GdmCl) induced irreversible denaturation was at 2.0 M GdmCl. Our results showed that amino acid substitutions increasing or decreasing interface hydrophobicity could further increase the thermostability of the Tp. citrate synthase. Thus, interface substitutions affecting the entropy of the unfolded state did not prove to be so critical in protein thermostabilization at higher temperatures.