Enzymatic transition states: thermodynamics, dynamics and analogue design

Kinetic isotope effects and computational chemistry have defined the transition state structures for several members of the N-ribosyltransferase family. Transition state analogues designed to mimic their cognate transition state structures are among the most powerful enzyme inhibitors. In complexes of N-ribosyltransferases with their transition state analogues, the dynamic nature of the transition state is converted to an ordered, thermodynamic structure closely related to the transition state. This phenomenon is documented by peptide bond H/D exchange, crystallography and computational chemistry. Complexes with substrate, transition state and product analogues reveal reaction coordinate motion and catalytic interactions. Isotope-edited spectroscopic analysis and binding specificity of these complexes provides information about specific enzyme-transition state contacts. In combination with protein dynamic QM/MM models, it is proposed that the transition state is reached by stochastic dynamic excursions of the protein groups near the substrates in the closed conformation. Examples from fully dissociated (D(N) *A(N)), hybrid (D(N)A(N)) and symmetric nucleophilic displacement (A(N)D(N)) transition states are found in the N-ribosyltransferases. The success of transition state analogue inhibitor design based on kinetic isotope effects validates this approach to understanding enzymatic transition states.     

Structural and biochemical evidence for a boat-like transition state in beta-mannosidases

Enzyme inhibition through mimicry of the transition state is a major area for the design of new therapeutic agents. Emerging evidence suggests that many retaining glycosidases that are active on alpha- or beta-mannosides harness unusual B2,5 (boat) transition states. Here we present the analysis of 25 putative beta-mannosidase inhibitors, whose Ki values range from nanomolar to millimolar, on the Bacteroides thetaiotaomicron beta-mannosidase BtMan2A. B2,5 or closely related conformations were observed for all tightly binding compounds. Subsequent linear free energy relationships that correlate log Ki with log Km/kcat for a series of active center variants highlight aryl-substituted mannoimidazoles as powerful transition state mimics in which the binding energy of the aryl group enhances both binding and the degree of transition state mimicry. Support for a B2,5 transition state during enzymatic beta-mannosidase hydrolysis should also facilitate the design and exploitation of transition state mimics for the inhibition of retaining alpha-mannosidases--an area that is emerging for anticancer therapeutics.     

Exploring the role of large conformational changes in the fidelity of DNA polymerase beta

The relationships between the conformational landscape, nucleotide insertion catalysis and fidelity of DNA polymerase beta are explored by means of computational simulations. The simulations indicate that the transition states for incorporation of right (R) and wrong (W) nucleotides reside in substantially different protein conformations. The protein conformational changes that reproduce the experimentally observed fidelity are significantly larger than the small rearrangements that usually accompany motions from the reactant state to the transition state in common enzymatic reactions. Once substrate binding has occurred, different constraints imposed on the transition states for insertion of R and W nucleotides render it highly unlikely that both transition states can occur in the same closed structure, because the predicted fidelity would then be many orders of magnitude too large. Since the conformational changes reduce the transition state energy of W incorporation drastically they decrease fidelity rather than increase it. Overall, a better agreement with experimental data is attained when the R is incorporated through a transition state in a closed conformation and W is incorporated through a transition state in one or perhaps several partially open conformations. The generation of free energy surfaces for R and W also allow us to analyze proposals about the relationship between induced fit and fidelity.     

Effects of high pressure on enzymatic activity

Effects of high pressure on enzymatic reactions are poised to revolutionize enzyme kinetics. The reason for this is that experimental designs are at hand to separate effects on equilibria between reactant states from effects on catalytic transition states and both yield new information. The first of the former runs contrary to Pauling's hypothesis that substrates are bound more tightly in the transition state, while the latter penetrates the 'black box' of catalysis, the stabilized transition state itself, and returns a precise measure of a physical parameter, deltaV. This in turn opens the door to new forms of structure-activity relationships. The first of these has been described, the effect of pressure on isotope effects, with the surprising finding that the entire isotope effect comes from a transition state phenomenon such as quantum mechanical hydrogen tunneling.     

Over-the-barrier transition state analogues and crystal structure with Mycobacterium tuberculosis purine nucleoside phosphorylase

Stable chemical analogues of enzymatic transition states are imperfect mimics since they lack the partial bond character of the transition state. We synthesized structural variants of the Immucillins as transition state analogues for purine nucleoside phosphorylase and characterized them with the enzyme from Mycobacterium tuberculosis (MtPNP). PNPs form transition states with ribooxacarbenium ion character and catalyze nucleophilic displacement reactions by migration of the cationic ribooxacarbenium carbon between the enzymatically immobilized purine and phosphate nucleophiles. As bond-breaking progresses, carbocation character builds on the ribosyl group, the distance between the purine and the carbocation increases, and the distance between carbocation and phosphate anion decreases. Transition state analogues were produced with carbocation character and increased distance between the ribooxacarbenium ion and the purine mimics by incorporating a methylene bridge between these groups. Immucillin-H (ImmH), DADMe-ImmH, and DADMe-ImmG mimic the transition state of MtPNP and are slow-onset, tight-binding inhibitors of MtPNP with equilibrium dissociation constants of 650, 42, and 24 pM. Crystal structures of MtPNP complexes with ImmH and DADMe-ImmH reveal an ion-pair between the inhibitor cation and the nucleophilic phosphoryl anion. The stronger ion-pair (2.7 A) is found with DADMe-ImmH. The position of bound ImmH resembles the substrate side of the transition state barrier, and DADMe-ImmH more closely resembles the product side of the barrier. The ability to probe both substrate and product sides of the transition state barrier provides expanded opportunities to explore transition state analogue design in N-ribosyltransferases. This approach has resulted in the highest affinity transition state analogues known for MtPNP.     

Ionic states of substrates and transition state analogues at the catalytic sites of N-ribosyltransferases

Purine nucleoside phosphorylase (PNP) and hypoxanthine-guanine phosphoribosyltransferase (HGPRTase) catalyze N-ribosidic bond cleavage in purine nucleosides and nucleotides, with addition of phosphate or pyrophosphate to form phosphorylated alpha-D-ribose products. The transition states have oxacarbenium ion character with a positive charge near 1'-C and ionic stabilization from nearby phosphoryl anions. Immucillin-H (ImmH) and Immucillin-H 5'-PO(4) (ImmHP) resemble the transition state charge when protonated at 4'-N and bind tightly to these enzymes with K(d) values of 20 pM to 1 nM. It has been proposed that Immucillins bind as the 4'-N neutral form and are protonated in the slow-onset step. Solution and solid-state NMR spectra of ImmH, ImmHP, guanosine, and GMP in complexes with two PNPs and a HGPRTase have been used to characterize their ionization states. Results with PNP*ImmH*PO(4) and HGPRTase*ImmHP*MgPP(i) indicate protonation at N-4' for the tightly bound inhibitors. The 1'-(13)C and 1'-(1)H resonances of bound Immucillins showed large downfield shifts as compared to Michaelis complexes, suggesting distortion of 1'-C toward sp(2) geometry. The Immucillins act as transition state mimics by binding with neutral iminoribitol groups followed by 4'-N protonation during slow-onset inhibition to form carbocationic mimics of the transition states. The ability of the Immucillins to mimic both substrate and transition state features contributes to their capture of transition state binding energy. Enzyme-activated phosphoryl nucleophiles bound to PNP and HGPRTase suggest enhanced electrostatic stabilization of the cationic transition states. Distortion of the oxacarbenium ion mimic toward transition state geometry is a common feature of the three distinct enzymatic complexes analyzed here. Substrate complexes, even in catalytically cycling equilibrium mixtures, do not reveal similar distortions.     

The dynamic role of distal side residues in heme hydroperoxidase catalysis. Interplay between X-ray crystallography and ab initio MD simulations

The enzymatic cycle of hydroperoxidases involves the resting Fe(III) state of the enzyme and the high-valent iron intermediates Compound I and Compound II. These states might be characterized by X-ray crystallography and the transition pathways between each state can be investigated using atomistic simulations. Here we review our recent work in the modeling of two key steps of the enzymatic reaction of hydroperoxidases: the formation of Cpd I in peroxidase and the reduction of Cpd I in catalase. It will be shown that small conformational motions of distal side residues (His in peroxidases and His/Asn in catalases), not,or only partially, revealed by the available X-ray structures, play an important role in the catalytic processes examined.     

Transition-state structures for N-glycoside hydrolysis of AMP by acid and by AMP nucleosidase in the presence and absence of allosteric activator

The mechanism of acid and enzymatic hydrolysis of the N-glycosidic bond of AMP has been investigated by fitting experimentally observed kinetic isotope effects [Parkin, D. W., & Schramm, V. L. (1987) Biochemistry (preceding paper in this issue)] to calculated kinetic isotope effects for proposed transition-state structures. The sensitivity of the transition-state calculations was tested by "arying the transition-state structure and comparing changes in the calculated kinetic isotope effects with the experimental values of the isotope effect measurements. The kinetic isotope effects for the acid-catalyzed hydrolysis of AMP are best explained by a transition state with considerable oxycarbonium character in the ribose ring, significant bonding remaining to the departing adenine ring, participation of a water nucleophile, and protonation of the adenine ring. A transition-state structure without preassociation of the water nucleophile cannot be eliminated by the data. Enzymatic hydrolysis of the N-glycosidic bond of AMP by AMP nucleosidase from Azotobacter vinelandii was analyzed in the absence and presence of MgATP, the allosteric activator that increases Vmax approximately 200-fold. The transition states for enzyme-catalyzed hydrolysis that best explain the kinetic isotope effects involve early SN1 transition states with significant bond order in the glycosidic bond and protonation of the adenine base. The enzyme enforces participation of an enzyme-bound water molecule, which has weak bonding to C1' in the transition state. Activation of AMP nucleosidase by MgATP causes the bond order of the glycosidic bond in the transition state to increase significantly. Hyperconjugation in the ribosyl group is altered by enzymatic stabilization of the oxycarbonium ion. This change is consistent with the interaction of an amino acid on the enzyme. Together, these changes stabilize a carboxonium-like transition-state complex that occurs earlier in the reaction pathway than in the absence of allosteric activator. In addition to the allosteric changes that alter transition-state structure, the presence of other inductive effects that are unobserved by kinetic isotope measurements is also likely to increase the catalytic rate.     

Mutation of Arg-166 of alkaline phosphatase alters the thio effect but not the transition state for phosphoryl transfer. Implications for the interpretation of thio effects in reactions of phosphatases

It has been suggested that the mechanism of alkaline phosphatase (AP) is associative, or triester-like, because phosphorothioate monoesters are hydrolyzed by AP approximately 10(2)-fold slower than phosphate monoesters. This "thio effect" is similar to that observed for the nonenzymatic hydrolysis of phosphate triesters, and is the inverse of that observed for the nonenzymatic hydrolysis of phosphate monoesters. The latter reactions proceed by loose, dissociative transition states, in contrast to reactions of triesters, which have tight, associative transition states. Wild-type alkaline phosphatase catalyzes the hydrolysis of p-nitrophenyl phosphate approximately 70 times faster than p-nitrophenyl phosphorothioate. In contrast, the R166A mutant alkaline phosphatase enzyme, in which the active site arginine at position 166 is replaced with an alanine, hydrolyzes p-nitrophenyl phosphate only about 3 times faster than p-nitrophenyl phosphorothioate. Despite this approximately 23-fold change in the magnitude of the thio effects, the magnitudes of Bronsted beta(lg) for the native AP (-0.77 +/- 0.09) and the R166A mutant (-0.78 +/- 0. 06) are the same. The identical values for the beta(lg) indicate that the transition states are similar for the reactions catalyzed by the wild-type and the R166A mutant enzymes. The fact that a significant change in the thio effect is not accompanied by a change in the beta(lg) indicates that the thio effect is not a reliable reporter for the transition state of the enzymatic phosphoryl transfer reaction. This result has important implications for the interpretation of thio effects in enzymatic reactions.     

Role of protein conformational dynamics in the catalysis by 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase

Enzymatic catalysis has conflicting structural requirements of the enzyme. In order for the enzyme to form a Michaelis complex, the enzyme must be in an open conformation so that the substrate can get into its active center. On the other hand, in order to maximize the stabilization of the transition state of the enzymatic reaction, the enzyme must be in a closed conformation to maximize its interactions with the transition state. The conflicting structural requirements can be resolved by a flexible active center that can sample both open and closed conformational states. For a bisubstrate enzyme, the Michaelis complex consists of two substrates in addition to the enzyme. The enzyme must remain flexible upon the binding of the first substrate so that the second substrate can get into the active center. The active center is fully assembled and stabilized only when both substrates bind to the enzyme. However, the side-chain positions of the catalytic residues in the Michaelis complex are still not optimally aligned for the stabilization of the transition state, which lasts only approximately 10(-13) s. The instantaneous and optimal alignment of catalytic groups for the transition state stabilization requires a dynamic enzyme, not an enzyme which undergoes a large scale of movements but an enzyme which permits at least a small scale of adjustment of catalytic group positions. This review will summarize the structure, catalytic mechanism, and dynamic properties of 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase and examine the role of protein conformational dynamics in the catalysis of a bisubstrate enzymatic reaction.     

Mathematical analysis of multienzyme systems. II. Steady state and transient control.

The control theory of steady states, previously presented for linear enzymatic systems (Heinrich and Rapoport, 1974) is extended to nonlinear systems. On the basis of three theorems a new procedure for the calculation of the control strength and of the control matrix is developed. The theory is applied to the extended model of glycolysis of erythrocytes, which includes also ATP-consuming processes. Also in this model the glycolytic flux is mainly controlled by the hexokinase-phosphofructokinase-system. The control strengths of the pyruvate kinase and of the enzymes of the 2.3 P2G-bypass are negligibly small. The control strength of the ATPase is negative, i.e. an activation of this enzyme leads to a decrease of the flux. For transition states of multienzyme systems definitions are given for the mean time required for the transition of the metabolites and for the "transient control" of enzymes. Enzymes with a pronounced influence on the transition time are called time-limiting enzymes. Enzymes which excert strong control on the time-dependent processes may have little influence under steady state conditions and vice versa. The transition times of ATP have been calculated for transient states of glycolysis.     

Secondary 18O isotope effects as a tool for studying reactions of phosphate mono-, di-, and triesters

Secondary 18O isotope effects have been developed as a tool for determining transition state structures in enzymatic and nonenzymatic phosphoryl transfer reactions. 18O substitution in the nonbridge oxygens of a phosphoryl group makes the reaction go faster when the bond order is higher to these oxygens in the transition state than in the reactant, whereas the reaction goes slower if the bond order is less. The isotope effects are measured by the remote label method, using an isotope ratio mass spectrometer for analysis. The bond order to p-nitrophenolate ion when it is the leaving group is indicated by the secondary 15N isotope effect in the nitro group, with a value of 1.0028 representing nearly complete bond cleavage. It appears that the transition states for phosphoryl transfer have no more than one negative charge on the nonbridge oxygens, so that reactions of monoesters are dissociative, reactions of triesters are associative, and reactions of diesters are SN2 with half bond order to entering and leaving groups.     

Action of cobra venom phospholipase A2 on the gel and liquid crystalline states of dimyristoyl and dipalmitoyl phosphatidylcholine vesicles

The activity of phospholipase A2 from cobra venom toward phospholipid in single-walled, sonicated vesicles was analyzed, particularly with respect to its activity toward the saturated phosphatidylcholines in the gel and liquid crystalline states. When egg phosphatidylcholine vesicles are used as substrate, the phospholipase has an apparent Km of 4.4 mM, an apparent Vmax of 100 mumol min-1 mg-1 of protein, and a pH optimum of 5.0 at 40 degrees C. The phospholipase hydrolyzed the gel state of dimyristoyl phosphatidylcholine vesicles and dipalmitoyl phosphatidylcholine vesicles at a rate 2 to 3 times greater than the liquid crystalline state, taking into account temperature effects on the enzymatic reaction itself. The results suggest that, toward sonicated vesicles, there is no specific enhancement of the rate when the both liquid crystalline and gel states are present together, as has been suggested to occur for multibilayers studied with other phospholipases. An apparent stimulation of activity as the reaction proceeded was observed above the phase transition temperature. This might be attributed to an increase in the phase transition temperature caused by free fatty acids so that, in the presence of reaction products, the enzyme is actually hydrolyzing gel state phospholipid which was found to be the preferred lipid state for phospholipase activity.     

Multiparameter study of denaturation of 20 beta-hydroxysteroid dehydrogenase by urea in the presence of stabilizing agents

We investigated the denaturation of tetrameric 20 beta-hydroxysteroid dehydrogenase (20R)-17 beta,20 beta,21-trihydroxysteroid:NAD+ oxidoreductase, EC 1.1.1.53) to find out whether intermediate states are formed during the process. The denaturation process was studied in the presence and absence of stabilizers, both specific, such as NADH, and non-specific, such as the salting-out anion phosphate. Changes in enzymatic activity, intrinsic protein fluorescence and far-ultraviolet circular dichroism were monitored. When NADH was present, denaturation of the enzyme by urea was a one-step transition between the native and the completely denatured state. In dilute phosphate, and even more so in concentrated phosphate, the existence of intermediate states with different stability is evidenced by the noncoincidence of the transition curves that probe for different functional and conformational aspects of the enzyme. Therefore, for 20 beta-hydroxysteroid dehydrogenase the formation of intermediates can be prevented by adding NADH, or enhanced by adding concentrated phosphate.     

Atomic motion in enzymatic reaction coordinates

Atomic excursions of reactants in enzymatic catalytic sites can be estimated from high-resolution crystal structures of enzyme complexes with substrates, transition state analog inhibitors and products. Transition state structures, defined from kinetic isotope effect studies, are compared to crystallographic structures to validate the properties of the transition state analog. Atomic excursions in enzymatic catalytic sites can differ from those in solution and define the role of the enzymatic catalyst in directing atomic motion.     

Population density analysis for determining the protonation state of the catalytic dyad in BACE1-tertiary carbinamine-based inhibitor complex

BACE1 is an aspartyl protease with a very relevant role in medicinal chemistry related to Alzheimer Disease since it has demonstrated to be a promising therapeutic target for inhibition and possible control for the progress of the peptide accumulation characteristic of this pathology. The enzymatic activity of this protein is given by the aspartic dyad, Asp93 and Asp289, which can adopt several protonation states depending on the chemical nature of its inhibitors, this is, monoprotonated, diprotonated and di-deprotonated states. In the present study, the analysis of the population density, for a series of protein-inhibitor molecular dynamics simulations, was carried out to identify the most feasible protonation state adopted by the catalytic dyad in the presence of tertiary carbinamine (TC) transition state analog inhibitors. The results revealed that the monoprotonated Asp289i state, in which the Asp93 and Asp289 residue side chains are deprotonated and protonated on the inner oxygen, respectively, is the most preferred in the presence of TC family inhibitors. This result was obtained after evaluating, for all 9 possible protonation state configurations, the individual and combined population densities of a set of parameters sensitive to protonation state of the Aspartic dyad, using an X-ray experimental BACE1/TC crystallographic structure as reference. This case study demonstrates again the usefulness of the concept of population density as a quantitative tool to establish the most stable system settings, among all possible, by measuring the level of occurrence of simultaneous events obtained from a sampling over time. These results will help to clear the phenomena related to the TCs inhibitory pathway, as well as assist in the design of better TC inhibitors against Alzheimer’s protease.     

Molecular modeling of the mechanochemical triggering mechanism for catalysis of carbon-cobalt bond homolysis in coenzyme B12

The possible contributions of the mechanochemical triggering effect to the enzymatic activation of the carbon-cobalt bond of coenzyme B12 (5'-deoxyadenosylcobalamin, AdoCbl) for homolytic cleavage have been studied by molecular modeling and semiempirical molecular orbital calculations. Classically, this effect has envisioned enzymatic compression of the axial Co-N bond in the ground state to cause upward folding of the corrin ring and subsequent sterically induced distortion of the Co-C bond leading to its destabilization. The models of this process show that in both methylcobalamin (CH3Cbl) and AdoCbl, compression of the axial Co-N bond does engender upward folding of the corrin ring, and that the extent of such upward folding is smaller in an analog in which the normal 5,6-dimethylbenzimidazole axial ligand is replaced by the sterically smaller ligand, imidazole (CH3(lm)Cbl and Ado(lm)Cbl). Furthermore, in AdoCbl, this upward folding of the corrin is accompanied by increases in the carbon-cobalt bond length and in the Co-C-C bond angle (which are also less pronounced in Ado(Im)Cbl), and which indicate that the Co-C bond is indeed destabilized by this mechanism. However, these effects on the Co-C bond are small, and destabilization of this bond by this mechanism is unlikely to contribute more than ca. 3 kcal mol(-1) towards the enzymatic catalysis of Co-C bond homolysis, far short of the observed ca. 14 kcal mol(-1). A second version of mechanochemical triggering, in which compression of the axial Co-N bond in the transition state for Co-C bond homolysis stabilizes the transition state by increased Co-N orbital overlap, has also been investigated. Stretching the Co-C bond to simulate the approach to the transition state was found to result in an upward folding of the corrin ring, a slight decrease in the axial Co-N bond length, a slight displacement of the metal atom from the plane of the equatorial nitrogens towards the "lower" axial ligand, and a decrease in strain energy amounting to about 8 kcal mol(-1) for both AdoCbl and Ado(Im)Cbl. In such modeled transition states, compression of the axial Co-N bond to just below 2.0 A (the distance subsequently found to provide maximal stabilization of the transition state by increased orbital overlap) required about 4 kcal mol(-1) for AdoCbl, and about 2.5 kcal mol(-1) for Ado(Im)Cbl. ZINDO/1 calculations on slightly simplified structures showed that maximal electronic stabilization of the transition state by about 10 kcal mol(-1) occurred at an axial Co-N bond distance of 1.96 A for both AdoCbl and Ado(Im)Cbl. The net result is that this type of transition state mechanochemical triggering can provide 14 kcal mol(-1) of transition state stabilization for AdoCbl, and about 15.5 kcal mol(-1) for the Ado(Im)Cbl, enough to completely explain the observed enzymatic catalysis. These results are discussed in the light of current knowledge about class I AdoCbl-dependent enzymes, in which the coenzyme is bound in its "base-off" conformation, with the lower axial ligand position occupied by the imidazole moiety of an active site histidine residue, and the class II enzymes, in which AdoCbl binds to the enzyme in its "base-on" conformation, and the pendent 5,6-dimethylbenzimidazole base remains coordinated to the metal during Co-C bond activation.     

Circularization changes the folding transition state of the src SH3 domain

Native state topology has been implicated as a major determinant of protein-folding mechanisms. Here, we test experimentally the robustness of the src SH3-domain folding transition state to changes in topology by covalently constraining regions of the protein with disulfide crosslinks and then performing kinetic analysis on point mutations in the context of these modified proteins. Circularization (crosslinking the N and C termini) of the src SH3 domain makes the protein topologically symmetric and causes delocalization of structure in the transition state ensemble suggesting a change in the folding mechanism. In contrast, crosslinking a single structural element (the distal beta-hairpin) which is an essential part of the transition state, results in a protein that folds 30 times faster, but does not change the distribution of structure in the transition state. As the transition states of distantly related SH3 domains were previously found to be very similar, we conclude that the free energy landscape of this protein family contains deep features which are relatively insensitive to sequence variations but can be altered by changes in topology.     

A coarse-grained model for force-induced protein deformation and kinetics

Force-induced changes in protein conformation are thought to be responsible for certain cellular responses to mechanical force. Changes in conformation subsequently initiate a biochemical response by alterations in, for example, binding affinity to another protein or enzymatic activity. Here, a model of protein extension under external forcing is created inspired by Kramers' theory for reaction rate kinetics in liquids. The protein is assumed to have two distinct conformational states: a relaxed state, C(1), preferred in the absence of external force, and an extended state, C(2), favored under force application. In the context of mechanotransduction, the extended state is a conformation from which the protein can initiate signaling. Appearance and persistence of C(2) are assumed to lead to transduction of the mechanical signal into a chemical one. The protein energy landscape is represented by two harmonic wells of stiffness kappa(1) and kappa(2), whose minima correspond to conformations C(1) and C(2). First passage time t(f) from C(1) to C(2) is determined from the Fokker-Plank equation employing several different approaches found in the literature. These various approaches exhibit significant differences in behavior as force increases. Although the level of applied force and the energy difference between states largely determine equilibrium, the dominant influence on t(f) is the height of the transition state. Distortions in the energy landscape due to force can also have a significant influence, however, exhibiting a weaker force dependence than exponential as previously reported, approaching a nearly constant value at a level of force that depends on the ratio kappa(1)/kappa(2). Two model systems are used to demonstrate the utility of this approach: a short alpha-helix undergoing a transition between two well-defined states and a simple molecular motor.