Kinetic mechanism of human DNA ligase I reveals magnesium-dependent changes in the rate-limiting step that compromise ligation efficiency

DNA ligase I (LIG1) catalyzes the ligation of single-strand breaks to complete DNA replication and repair. The energy of ATP is used to form a new phosphodiester bond in DNA via a reaction mechanism that involves three distinct chemical steps: enzyme adenylylation, adenylyl transfer to DNA, and nick sealing. We used steady state and pre-steady state kinetics to characterize the minimal mechanism for DNA ligation catalyzed by human LIG1. The ATP dependence of the reaction indicates that LIG1 requires multiple Mg(2+) ions for catalysis and that an essential Mg(2+) ion binds more tightly to ATP than to the enzyme. Further dissection of the magnesium ion dependence of individual reaction steps revealed that the affinity for Mg(2+) changes along the reaction coordinate. At saturating concentrations of ATP and Mg(2+) ions, the three chemical steps occur at similar rates, and the efficiency of ligation is high. However, under conditions of limiting Mg(2+), the nick-sealing step becomes rate-limiting, and the adenylylated DNA intermediate is prematurely released into solution. Subsequent adenylylation of enzyme prevents rebinding to the adenylylated DNA intermediate comprising an Achilles' heel of LIG1. These ligase-generated 5'-adenylylated nicks constitute persistent breaks that are a threat to genomic stability if they are not repaired. The kinetic and thermodynamic framework that we have determined for LIG1 provides a starting point for understanding the mechanism and specificity of mammalian DNA ligases.     

Cytochrome P450 7A1 cholesterol 7alpha-hydroxylation: individual reaction steps in the catalytic cycle and rate-limiting ferric iron reduction

Cytochrome P450 (P450) 7A1 is well known as the cholesterol 7α-hydroxylase, the first enzyme involved in bile acid synthesis from cholesterol. The human enzyme has been reported to have the highest catalytic activity of any mammalian P450. Analyses of individual steps of cholesterol 7α-hydroxylation reaction revealed several characteristics of this reaction: (i) two-step binding of cholesterol to ferric P450, with an apparent K(d) of 0.51 μM, (ii) a rapid reduction rate in the presence of cholesterol (～10 s(-1) for the fast phase), (iii) rapid formation of a ferrous P450-cholesterol-O(2) complex (29 s(-1)), (iv) the lack of a non-competitive kinetic deuterium isotope effect, (v) the lack of a kinetic burst, and (vi) the lack of a deuterium isotope effect when the reaction was initiated with the ferrous P450-cholesterol complex. A minimum kinetic model was developed and is consistent with all of the observed phenomena and the rates of cholesterol 7α-hydroxylation and H(2)O and H(2)O(2) formation. The results indicate that the first electron transfer step, although rapid, becomes rate-limiting in the overall P450 7A1 reaction. This is a different phenomenon compared with other P450s that have much lower rates of catalysis, attributed to the much more efficient substrate oxidation steps in this reaction.     

Mechanical modulation of catalytic power on F1-ATPase

The conformational fluctuation of enzymes has a crucial role in reaction acceleration. However, the contribution to catalysis enhancement of individual substates with conformations far from the average conformation remains unclear. We studied the catalytic power of the rotary molecular motor F(1)-ATPase from thermophilic Bacillus PS3 as it was stalled in transient conformations far from a stable pausing angle. The rate constants of ATP binding and hydrolysis were determined as functions of the rotary angle. Both rates exponentially increase with rotation, revealing the molecular basis of positive cooperativity among three catalytic sites: elementary reaction steps are accelerated via the mechanical rotation driven by other reactions on neighboring catalytic sites. The rate enhancement induced by ATP binding upon rotation was greater than that brought about by hydrolysis, suggesting that the ATP binding step contributes more to torque generation than does the hydrolysis step. Additionally, 9% of the ATP-driven rotary step was supported by thermal diffusion, suggesting that acceleration of the ATP docking process occurs via thermally agitated conformational fluctuations.     

Affinity purification and catalytic properties of a soluble, Ca2+-independent, diacylglycerol kinase

We used a new procedure that involved selective enzyme binding to lipid vesicles to partially purify a soluble diacylglycerol kinase, then studied the relation between enzyme-vesicle binding and activity in vesicle-based assays. The vesicle-binding procedure required about 2 h, increased the enzyme's specific activity 50-fold with a 50% yield of activity, and combined well with additional purification steps. Studies of the activity of the partially purified diacylglycerol kinase toward vesicle-associated diacylglycerols revealed linear reaction kinetics that reflected enzyme binding to the vesicles; factors known to influence enzyme binding to the vesicles affected enzyme activity only indirectly, not by influencing the diacylglycerol kinase reaction itself. On the other hand, special incubation experiments that caused both substrate depletion in vesicles and enzyme stalling provided evidence that the diacylglycerol kinase could desorb from these vesicles, adsorb to freshly added, substrate-containing vesicles, and resume catalysis of phosphorylation reactions. The molecular basis for this enzyme-vesicle "hopping" behavior remains to be clarified. But enzyme-catalyzed conversion of diacylglycerol to phosphatidic acid may not have been a contributing factor because separate, enzyme-vesicle binding experiments showed that the enzyme had only a marginally higher affinity for diacylglycerol-containing vesicles than it did for vesicles that contained comparable amounts of phosphatidic acid. The combined results of our experiments suggest that the linear rates of diacylglycerol phosphorylation observed in standard assays with diacylglycerol-containing vesicles may have been combined functions of both the rate of enzyme hopping among vesicles and the rate of diacylglycerol phosphorylation by enzyme that was bound transiently on substrate-containing vesicles.     

Chemical transformation is not rate-limiting in the reaction catalyzed by Escherichia coli 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) catalyzes the transfer of pyrophosphate from ATP to 6-hydroxymethyl-7,8-dihydropterin (HMDP). Because HPPK is essential for microorganisms but is absent from human and animals, the enzyme is an excellent target for developing antimicrobial agent. Thermodynamic analysis shows that Mg(2+) is important not only for the binding of nucleotides but also for the binding of HMDP. Transient kinetic analysis shows that a step or steps after the chemical transformation are rate-limiting in the reaction catalyzed by HPPK. The pre-steady-state kinetics is composed of a burst phase and a steady-state phase. The rate constant for the burst phase is approximately 50 times larger than that for the steady-state phase. The latter is very similar to the k(cat) value measured by steady-state kinetics. A set of rate constants for the individual steps of the HPPK-catalyzed reaction has been determined by a combination of stopped-flow and quench-flow analyses. These results form a thermodynamic and kinetic framework for dissecting the roles of active site residues in the substrate binding and catalysis by HPPK.     

Arginine coordination in enzymatic phosphoryl transfer: evaluation of the effect of Arg166 mutations in Escherichia coli alkaline phosphatase

Arginine residues are commonly found in the active sites of enzymes catalyzing phosphoryl transfer reactions. Numerous site-directed mutagenesis experiments establish the importance of these residues for efficient catalysis, but their role in catalysis is not clear. To examine the role of arginine residues in the phosphoryl transfer reaction, we have measured the consequences of mutations to arginine 166 in Escherichia coli alkaline phosphatase on hydrolysis of ethyl phosphate, on individual reaction steps in the hydrolysis of the covalent enzyme-phosphoryl intermediate, and on thio substitution effects. The results show that the role of the arginine side chain extends beyond its positive charge, as the Arg166Lys mutant is as compromised in activity as Arg166Ser. Through measurement of individual reaction steps, we construct a free energy profile for the hydrolysis of the enzyme-phosphate intermediate. This analysis indicates that the arginine side chain strengthens binding by approximately 3 kcal/mol and provides an additional 1-2 kcal/mol stabilization of the chemical transition state. A 2.1 A X-ray diffraction structure of Arg166Ser AP is presented, which shows little difference in enzyme structure compared to the wild-type enzyme but shows a significant reorientation of the bound phosphate. Altogether, these results support a model in which the arginine contributes to catalysis through binding interactions and through additional transition state stabilization that may arise from complementarity of the guanidinum group to the geometry of the trigonal bipyramidal transition state.     

Kinetics of ATP-stimulated nuclease activity of the Escherichia coli RecBCD enzyme

The RecBCD enzyme is an ATP-dependent nuclease on both single-stranded and double-stranded DNA substrates. We have investigated the kinetics of the RecBCD-catalyzed reaction with small, single-stranded oligodeoxyribonucleotide substrates under single-turnover conditions using rapid-quench flow techniques. RecBCD-DNA complexes were allowed to form in pre-incubation mixtures. The nuclease reactions were initiated by mixing with ATP. The reaction time-courses were fit to several possible reaction mechanisms and quantitative estimates were obtained for rate constants for individual reaction steps. The relative rates of forward reaction versus dissociation from the DNA, and the fact that inclusion of excess non-radiolabeled single-stranded DNA to trap free RecBCD has no effect on the nuclease reaction, indicates that the reaction is processive. The reaction products show that the reaction begins near the 3'-end of the [5'-32P]DNA substrates and the major cleavage sites are two to four phosphodiester bonds apart. The product distribution is unchanged as the ATP concentration varies from 10 microM to 100 microM ATP, while the overall reaction rate varies by about tenfold. These observations suggest that DNA cleavage is tightly coordinated with movement of the enzyme along the DNA. The reaction time-courses at low concentrations of ATP (10 microM and 25 microM) have a significant lag before cleavage products appear. We propose that the lag represents ATP-dependent movement of the DNA from an initial binding site in the helicase domain of the RecB subunit to the nuclease active site in a separate domain of RecB. The extent of reaction of the substrate is limited (approximately 50%) under all conditions. This may indicate the formation of a non-productive RecBCD-DNA complex that does not dissociate in the 1-2 s time-scale of our experiments.     

Pre-steady-state studies of phosphite dehydrogenase demonstrate that hydride transfer is fully rate limiting

Phosphite dehydrogenase (PTDH) is a unique NAD-dependent enzyme that catalyzes the oxidation of inorganic phosphite to phosphate. The enzyme has great potential for cofactor regeneration, and mechanistic studies have provided some insight into the residues that are important for catalysis. In this investigation, pre-steady-state studies were performed on the His6-tagged wild-type (WT) enzyme, several active site mutants, a thermostable mutant (12X-PTDH), and a thermostable mutant with dual cofactor specificity (NADP-12X-PTDH). Stopped-flow kinetic experiments indicate that slow steps after hydride transfer do not significantly limit the rate of reaction for the WT enzyme, the active site mutants, or the thermostable mutant. Pre-steady-state kinetic isotope effects (KIEs) and single-turnover experiments further confirm that slow steps after the chemical step do not significantly limit the rate of reaction for any of these proteins. Collectively, these results suggest that the hydride transfer step is fully rate determining in PTDH and that the observed KIE on kcat is the intrinsic effect in WT PTDH and the mutants examined. In contrast, a slow step after catalysis may partially limit the rate of phosphite oxidation by NADP-12X-PTDH with NADP as the cofactor. Finally, site-directed mutagenesis of Asp79 indicates that this residue is important in orienting Arg237 for proper interaction with phosphite.     

Dehydration is catalyzed by glutamate-136 and aspartic acid-135 active site residues in Escherichia coli dTDP-glucose 4,6-dehydratase

The dTDP-glucose 4,6-dehydratase catalyzed conversion of dTDP-glucose to dTDP-4-keto-6-deoxyglucose occurs in three sequential chemical steps: dehydrogenation, dehydration, and rereduction. The enzyme contains the tightly bound coenzyme NAD(+), which mediates the dehydrogenation and rereduction steps of the reaction mechanism. In this study, we have determined that Asp135 and Glu136 are the acid and base catalysts, respectively, of the dehydration step. Identification of the acid catalyst was performed using an alternative substrate, dTDP-6-fluoro-6-deoxyglucose (dTDP-6FGlc), which undergoes fluoride ion elimination instead of dehydration, and thus does not require protonation of the leaving group. The steady-state rate of conversion of dTDP-6FGlc to dTDP-4-keto-6-deoxyglucose by each Asp135 variant was identical to that of wt, in contrast to turnover using dTDP-glucose where differences in rates of up to 2 orders of magnitude were observed. These results demonstrate Asp135's role in protonating the glucosyl-C6(OH) during dehydration. The base catalyst was identified using a previously uncharacterized, enzyme-catalyzed glucosyl-C5 hydrogen-solvent exchange reaction of product, dTDP-4-keto-6-deoxyglucose. Base catalysis of this exchange reaction is analogous to that occurring at C5 during the dehydration step of net catalysis. Thus, the decrease in the rate of catalysis ( approximately 2 orders of magnitude) of the exchange reaction observed with Glu136 variants demonstrates this residue's importance in base catalysis of dehydration.     

Mechanism of formation of bound alpha-iminoglutarate from alpha-ketoglutarate in the glutamate dehydrogenase reaction. A chemical basis for ammonia recognition

Two mechanisms have been postulated for the formation of bound alpha-iminoglutarate intermediate during the glutamate dehydrogenase-catalyzed reductive amination of alpha-ketoglutarate; one involves the nucleophilic attack of ammonia on a covalently bound Schiff base in the enzyme-NADPH-alpha-ketoglutarate complex, and the other involves the reaction of ammonia with the carbonyl group of alpha-ketoglutarate in the ternary complex. We have measured the rates of carbonyl oxygen exchange in the complex to unambiguously distinguish between these two mechanisms. We find that the loss of label in the carbonyl oxygen-labeled ternary complex is at least 10(5) times slower than the rate of the reductive amination reaction. Therefore, the former mechanism cannot be operative. We also find that (i) the carbonyl oxygen exchange in free alpha-ketoglutarate proceeds without any significant catalysis by its gamma-carboxylate group; (ii) this exchange reaction has energy parameters which are comparable to those observed for the hydration of simple aliphatic ketones; and (iii) the carbonyl oxygen exchange in bound alpha-ketoglutarate is slower than that in the free keto acid over a wide pH range. We conclude that the oxygen exchange in the free and bound alpha-ketoglutarate must occur via a gem-diol intermediate. The observation that the enzyme inhibits the reaction of water with alpha-ketoglutarate while it catalyzes the reaction of ammonia with the same keto acid points to an extraordinary recognition of ammonia by the enzyme. We interpret this observation by assuming that the enzyme-NADPH-alpha-ketoglutarate complex exists in two forms, a predominant form which is produced rapidly upon mixing the components together and an unstable form which is produced in trace amounts from the predominant form via a gem-diol intermediate. These two forms are presumed to differ in the spatial relationship of the carbonyl group to the enzyme functional groups. The carbonyl group in the unstable form is assumed to be surrounded by the same enzyme groups as the iminium ion is in the bound iminoglutarate complex. We ascribe the remarkable catalysis of the ammonia reaction and the inhibition of the water reaction by the enzyme to the opposing interactions of the iminium and carbonyl groups with these surrounding enzyme groups.     

Kinetic mechanisms of cholesterol oxidase from Streptomyces hygroscopicus and Brevibacterium sterolicum.

The kinetic properties of two cholesterol oxidases, one from Brevibacterium sterolicum (BCO) the other from Streptomyces hygroscopicus (SCO) were investigated. BCO works via a ping-pong mechanism, whereas the catalytic pathway of SCO is sequential. The turnover numbers at infinite cholesterol and oxygen concentrations are 202 s-1 and 105 s-1 for SCO and BCO, respectively. The rates of flavin reduction extrapolated to saturating substrate concentration, under anaerobic conditions, are 235 s-1 for BCO and 232 s-1 for SCO (in the presence of 1% Thesit and 10% 2-propanol). With reduced SCO the rate of Delta5-6-->Delta4-5 isomerization of the intermediate 5-cholesten-3-one to final product is slow (0.3 s-1). With oxidized SCO and BCO the rate of isomerization is much faster ( approximately 300 s-1), thus it is not rate-limiting for catalysis. The kinetic behaviour of both reduced COs towards oxygen is unusual in that they exhibit apparent saturation with increasing oxygen concentrations (extrapolated rates approximately 250 s-1 and 1.3 s-1, for BCO and SCO, respectively): too slow to account for catalysis. For BCO the kinetic data are compatible with a step preceding the reaction with oxygen, involving interconversion of reactive and nonreactive forms of the enzyme. We suggest that the presence of micelles in the reaction medium, due to the necessary presence of detergents to solubilize the substrate, influence the availability or reactivity of oxygen towards the enzyme. The rate of re-oxidation of SCO in the presence of product is also too slow to account for catalysis, probably due to the impossibility of producing quantitatively the reduced enzyme-product complexes.     

Inhibitory potency and specificity of subtilase-like pro-protein convertase (SPC) prodomains.

The SPCs (subtilisin-like pro-protein convertases) are a family of enzymes responsible for the proteolytic processing of numerous precursor proteins of the constitutive and regulated secretory pathways. SPCs are themselves synthesized as inactive zymogens. Activation of SPCs occurs via the intramolecular autocatalytic removal of the prodomain. SPC prodomains have been proposed as templates in the development of potent and specific SPC inhibitors. In this study, we investigated the specificity and potency of complete prodomains and short C-terminal prodomain peptides of each SPC on highly purified, soluble enzyme preparations of human SPC1, SPC6, and SPC7. Progress curve kinetic analysis of prodomain peptides and complete prodomains showed competitive inhibitory profiles in the low nanomolar range. Complete prodomains were 5-100 times more potent than C-terminal prodomain peptides, suggesting that N-terminal determinants are involved in the recognition process. However, complete prodomains and prodomain peptides exhibit only a partial specificity toward their cognate enzyme. Ala-scan structure activity studies indicated the importance of basic residues in the P(4), P(5), and P(6) positions for inhibition of SPC1. In contrast, hydrophobic residues in P(6) and P(7), as well as basic residues in P(4) and P(5), were critical for inhibition of SPC7. Our data demonstrated that the use of prodomains as specific inhibitors acting in trans would be of limited usefulness, unless modified into more specific compounds.     

Insight into microwave irradiation and enzyme catalysis in enantioselective resolution of dl-(±)-3-phenyllactic acid

Lipase catalyzed kinetic resolution of DL-(±)-3-phenyllactic acid (DL-(±)-3-PLA) was investigated to study the synergistic effect of microwave irradiation and enzyme catalysis. Lipases from different sources were employed for the transesterification of DL-(±)-3-PLA under otherwise similar conditions, among which Novozyme 435 efficiently catalyzed the resolution of DL-(±)-3-PLA to L-(-)-O-acetyl-3-PLA using vinyl acetate as the acyl donor, showing excellent conversion (49?%) and enantiomeric excess (>99?%). The effect of various parameters affecting the initial rate, conversion and enantiomeric excess of the reaction were studied to establish kinetics and mechanism. There is a synergism between enzyme catalysis and microwave irradiation; an increase in initial rates up to 1.8-fold was observed under microwave irradiation than that under conventional heating. The analysis of initial rate data showed that reaction obeys ternary complex (ordered bi-bi) mechanism with inhibition by DL-(±)-3-PLA. The calculated and simulated rates match very well showing the validity of the proposed kinetic model.     

Interaction of NADP(H) with oxidized and reduced P450 reductase during catalysis. Studies with nucleotide analogues

Previous studies have shown that the interaction of P450 reductase with bound NADP(H) is essential to ensure fast electron transfer through the two flavin cofactors. In this study we investigated in detail the interaction of the house fly flavoprotein with NADP(H) and a number of nucleotide analogues. 1,4,5,6-Tetrahydro-NADP, an analogue of NADPH, was used to characterize the interaction of P450 reductase with the reduced nucleotide. This analogue is inactive as electron donor, but its binding affinity and rate constant of release are very close to those for NADPH. The 2'-phosphate contributes about 5 kcal/mol of the binding energy of NADP(H). Oxidized nicotinamide does not interact with the oxidized flavoprotein, while reduced nicotinamide contributes 1.3 kcal/mol of the binding energy. Oxidized P450 reductase binds NADPH with a K(d) of 0.3 microM, while the affinity of the reduced enzyme is considerably lower, K(d) = 1.9 microM. P450 reductase catalyzes a transhydrogenase reaction between NADPH and oxidized nucleotides, such as thionicotinamide-NADP(+), acetylpyridine-NADP(+), or [(3)H]NADP(+). The reverse reaction, reduction of [(3)H]NADP(+) by the reduced analogues, is also catalyzed by P450 reductase. We define the mechanism of the transhydrogenase reaction as follows: NADPH binding, hydride ion transfer, and release of the NADP(+) formed. An NADP(+) or its analogue binds to the two-electron-reduced flavoprotein, and the electron-transfer steps reverse to transfer hydride ion to the oxidized nucleotide, which is released. Measurements of the flavin semiquinone content, rate constant for NADPH release, and transhydrogenase turnover rates allowed us to estimate the steady-state distribution of P450 reductase species during catalysis, and to calculate equilibrium constants for the interconversion of catalytic intermediates. Our results demonstrate that equilibrium redox potentials of the flavin cofactors are not the sole factor governing rapid electron transfer during catalysis, but conformational changes must be considered to understand P450 reductase catalysis.     

Enzymatic Activity versus Structural Dynamics: The Case of Acetylcholinesterase Tetramer

The function of many proteins, such as enzymes, is modulated by structural fluctuations. This is especially the case in gated diffusion-controlled reactions (where the rates of the initial diffusional encounter and of structural fluctuations determine the overall rate of the reaction) and in oligomeric proteins (where function often requires a coordinated movement of individual subunits). A classic example of a diffusion-controlled biological reaction catalyzed by an oligomeric enzyme is the hydrolysis of synaptic acetylcholine (ACh) by tetrameric acetylcholinesterase (AChEt). Despite decades of efforts, the extent to which enzymatic efficiency of AChEt (or any other enzyme) is modulated by flexibility is not fully determined. This article attempts to determine the correlation between the dynamics of AChEt and the rate of reaction between AChEt and ACh. We employed equilibrium and nonequilibrium electro-diffusion models to compute rate coefficients for an ensemble of structures generated by molecular dynamics simulation. We found that, for the static initial model, the average reaction rate per active site is ∼22-30% slower in the tetramer than in the monomer. However, this effect of tetramerization is modulated by the intersubunit motions in the tetramer such that a complex interplay of steric and electrostatic effects either guides or blocks the substrate into or from each of the four active sites. As a result, the rate per active site calculated for some of the tetramer structures is only ∼15% smaller than the rate in the monomer. We conclude that structural dynamics minimizes the adverse effect of tetramerization, allowing the enzyme to maintain similar enzymatic efficiency in different oligomerization states.     

Kinetic studies on beta-site amyloid precursor protein-cleaving enzyme (BACE). Confirmation of an iso mechanism

The steady-state kinetic mechanism of beta-amyloid precursor protein-cleaving enzyme (BACE)-catalyzed proteolytic cleavage was evaluated using product and statine- (Stat(V)) or hydroxyethylene-containing (OM99-2) peptide inhibition data, solvent kinetic isotope effects, and proton NMR spectroscopy. The noncompetitive inhibition pattern observed for both cleavage products, together with the independence of Stat(V) inhibition on substrate concentration, suggests a uni-bi-iso kinetic mechanism. According to this mechanism, the enzyme undergoes multiple conformation changes during the catalytic cycle. If any of these steps are rate-limiting to turnover, an enzyme form preceding the rate-limiting conformational change should accumulate. An insignificant solvent kinetic isotope effect (SKIE) on k(cat)/K(m), a large inverse solvent kinetic isotope effect on k(cat), and the absence of any SKIE on the inhibition onset by Stat(V) during catalysis together indicate that the rate-limiting iso-step occurs after formation of a tetrahedral intermediate. A moderately short and strong hydrogen bond (at delta 13.0 ppm and phi of 0.6) has been observed by NMR spectroscopy in the enzyme-hydroxyethylene peptide (OM99-2) complex that presumably mimics the tetrahedral intermediate of catalysis. Collapse of this intermediate, involving multiple steps and interconversion of enzyme forms, has been suggested to impose a rate limitation, which is manifested in a significant SKIE on k(cat). Multiple enzyme forms and their distribution during catalysis were evaluated by measuring the SKIE on the noncompetitive (mixed) inhibition constants for the C-terminal reaction product. Large, normal SKIE values were observed for these inhibition constants, suggesting that both kinetic and thermodynamic components contribute to the K(ii) and K(is) expressions, as has been suggested for other iso-mechanism featuring enzymes. We propose that a conformational change related to the reprotonation of aspartates during or after the bond-breaking event is the rate-limiting segment in the catalytic reaction of beta-amyloid precursor protein-cleaving enzyme, and ligands binding to other than the ground-state forms of the enzyme might provide inhibitors of greater pharmacological relevance.     

Correlated conformational fluctuations during enzymatic catalysis: Implications for catalytic rate enhancement

Correlated enzymatic conformational fluctuations are shown to contribute to the rate of enhancement achieved during catalysis. Cytidine deaminase serves as a model system. Crystallographic temperature factor data for this enzyme complexed with substrate analog, transition-state analog, and product are available, thereby establishing a measure of atomic scale conformational fluctuations along the (approximate) reaction coordinate. First, a neural network-based algorithm is used to visualize the decreased conformational fluctuations at the transition state. Second, a dynamic diffusion equation along the reaction coordinate is solved and shows that the flux velocity through the associated enzymatic conformation space is greatest at the transition state. These results suggest (1) that there are both dynamic and energetic restrictions to conformational fluctuations at the transition state, (2) that enzymatic catalysis occurs on a fluctuating potential energy surface, and (3) a form for the potential energy. The Michaelis-Menten equations are modified to describe catalysis on this fluctuating potential energy profile, leading to enhanced catalytic rates when fluctuations along the reaction coordinate are appropriately correlated. This represents a dynamic tuning of the enzyme for maximally effective transformation of the ES complex into EP.     

On the relationship between thermal stability and catalytic power of enzymes

The possible relationship between the thermal stability and the catalytic power of enzymes is of great current interest. In particular, it has been suggested that thermophilic or hyperthermophilic (Tm) enzymes have lower catalytic power at a given temperature than the corresponding mesophilic (Ms) enzymes, because the thermophilic enzymes are less flexible (assuming that flexibility and catalysis are directly correlated). These suggestions presume that the reduced dynamics of the thermophilic enzymes is the reason for their reduced catalytic power. The present paper takes the specific case of dihydrofolate reductase (DHFR) and explores the validity of the above argument by simulation approaches. It is found that the Tm enzymes have restricted motions in the direction of the folding coordinate, but this is not relevant to the chemical process, since the motions along the reaction coordinate are perpendicular to the folding motions. Moreover, it is shown that the rate of the chemical reaction is determined by the activation barrier and the corresponding reorganization energy, rather than by dynamics or flexibility in the ground state. In fact, as far as flexibility is concerned, we conclude that the displacement along the reaction coordinate is larger in the Tm enzyme than in the Ms enzyme and that the general trend in enzyme catalysis is that the best catalyst involves less motion during the reaction than the less optimal catalyst. The relationship between thermal stability and catalysis appears to reflect the fact that to obtain small electrostatic reorganization energy it is necessary to invest some folding energy in the overall preorganization process. Thus, the optimized catalysts are less stable. This trend is clearly observed in the DHFR case.     

The kinetics of conformation change as determinant of Rubisco's specificity

The molecular basis of Rubisco's specificity is investigated in terms of the structure and kinetics of the enzyme. We propose that the rates of the conformational changes (closing/opening) of the binding niche exert a crucial influence on apparent binding rates and the enzyme's specificity. An extended reaction scheme for binding and conformational kinetics is presented and expressed in a mathematical model. The closed conformation, known from X-ray structures, is assumed to be necessary for binding of the gaseous substrates (carbon dioxide and oxygen) and for catalysis. Opening the niche interrupts catalysis and enables a fast exchange of those molecules between the internal cavity and the surrounding solvent. Our model predicts that specificity of Rubisco for CO₂ increases with the rate by which the niche opens. This is due to the fact that binding of the carbon dioxide is faster than oxygen binding, which is hampered by spin inversion. The apparent rate of carbon dioxide binding correlates with the repetition rate of the conformational change, and the rate of oxygen binding with the probability of the closed state. This revised version was published online in June 2006 with corrections to the Cover Date.     

DNA cytosine C5 methyltransferase Dnmt1: catalysis-dependent release of allosteric inhibition

We followed the cytosine C(5) exchange reaction with Dnmt1 to characterize its preference for different DNA substrates, its allosteric regulation, and to provide a basis for comparison with the bacterial enzymes. We determined that the methyl transfer is rate-limiting, and steps up to and including the cysteine-cytosine covalent intermediate are in rapid equilibrium. Changes in these rapid equilibrium steps account for many of the previously described features of Dnmt1 catalysis and specificity including faster reactions with premethylated DNA versus unmethylated DNA, faster reactions with DNA in which guanine is replaced with inosine [poly(dC-dG) vs poly(dI-dC)], and 10-100-fold slower catalytic rates with Dnmt1 relative to the bacterial enzyme M.HhaI. Dnmt1 interactions with the guanine within the CpG recognition site can prevent the premature release of the target base and solvent access to the active site that could lead to mutagenic deamination. Our results suggest that the beta-elimination step following methyl transfer is not mediated by free solvent. Dnmt1 shows a kinetic lag in product formation and allosteric inhibition with unmethylated DNA that is not observed with premethylated DNA. Thus, we suggest the enzyme undergoes a slow relief from allosteric inhibition upon initiation of catalysis on unmethylated DNA. Notably, this relief from allosteric inhibition is not caused by self-activation through the initial methylation reaction, as the same effect is observed during the cytosine C(5) exchange reaction in the absence of AdoMet. We describe limitations in the Michaelis-Menten kinetic analysis of Dnmt1 and suggest alternative approaches.