Energetics of RNA cleavage: implications for the mechanism of action of ribozymes.

A new class of ribozymes produce 2',3'-cyclic phosphate upon self-catalyzed cleavage of RNA molecules, similar to those observed during enzymatic (RNase-catalyzed) as well as non-enzymatic hydrolyses of RNAs. This product suggests that the reaction intermediate/transition state is a pentacoordinated oxyphosphorane. In order to elucidate the energetics of these RNA cleaving reactions, the reaction coordinate has been simulated and a pentacoordinated intermediate has been characterized via ab initio molecular orbital calculations utilizing the dianionic hydrolysis-intermediate of methyl ethylene phosphate as a model compound. The calculated reaction coordinate indicates that the transition state for the P-O(2') bond cleavage is lower in energy than that for the P-O(5') bond cleavage under uncatalyzed conditions. Thus, the dianionic pentacoordinated phosphorus intermediate tends to revert back to the starting RNA by cleaving the P-O(2') bond rather than productively cleaving the P-O(5') bond. In order for ribozymes to effectively cleave RNA molecules, it is therefore mandatory to stabilize the leaving 5'-oxygen, e.g. by means of a divalent magnesium ion.     

RNA hydrolysis via an oxyphosphorane intermediate.

From calculations of a model reaction scheme for base-catalyzed RNA hydrolysis, a pentacoodinate dianionic intermediate 2a (Storer, et al., J. Am. Chem. Soc., 1991, 113, 5216-5219) as well as two transition states, TS1 and TS2, to the intermediate have been located by ab initio calculations at the 3-21G* level. Although the intermediate, which has the well depth on the order of kBT, is unlikely to be kinetically significant, the overall rate-limiting transition state structure TS2 obtained at 3-21G* level is very close to the corresponding structure at the STO-3G level; it has an extended P-O(5') bond breaking character. These gas-phase calculation results are used to qualitatively interpret mutagenesis results of Barnase and RNase T1 where water molecules are absent from the active site.     

P-O bond destabilization accelerates phosphoenzyme hydrolysis of sarcoplasmic reticulum Ca2+ -ATPase

The phosphate group of the ADP-insensitive phosphoenzyme (E2-P) of sarcoplasmic reticulum Ca2+ -ATPase (SERCA1a) was studied with infrared spectroscopy to understand the high hydrolysis rate of E2-P. By monitoring an autocatalyzed isotope exchange reaction, three stretching vibrations of the transiently bound phosphate group were selectively observed against a background of 50,000 protein vibrations. They were found at 1194, 1137, and 1115 cm(-1). This information was evaluated using the bond valence model and empirical correlations. Compared with the model compound acetyl phosphate, structure and charge distribution of the E2-P aspartyl phosphate resemble somewhat the transition state in a dissociative phosphate transfer reaction; the aspartyl phosphate of E2-P has 0.02 A shorter terminal P-O bonds and a 0.09 A longer bridging P-O bond that is approximately 20% weaker, the angle between the terminal P-O bonds is wider, and -0.2 formal charges are shifted from the phosphate group to the aspartyl moiety. The weaker bridging P-O bond of E2-P accounts for a 10(11)-10(15)-fold hydrolysis rate enhancement, implying that P-O bond destabilization facilitates phosphoenzyme hydrolysis. P-O bond destabilization is caused by a shift of noncovalent interactions from the phosphate oxygens to the aspartyl oxygens. We suggest that the relative positioning of Mg2+ and Lys684 between phosphate and aspartyl oxygens controls the hydrolysis rate of the ATPase phosphoenzymes and related phosphoproteins.     

Mechanisms of guanosine triphosphate hydrolysis by Ras and Ras-GAP proteins as rationalized by ab initio QM/MM simulations

The hydrolysis reaction of guanosine triphosphate (GTP) by p21(ras) (Ras) has been modeled by using the ab initio type quantum mechanical-molecular mechanical simulations. Initial geometry configurations have been prompted by atomic coordinates of the crystal structure (PDBID: 1QRA) corresponding to the prehydrolysis state of Ras in complex with GTP. Multiple searches of minimum energy geometry configurations consistent with the hydrogen bond networks have been performed, resulting in a series of stationary points on the potential energy surface for reaction intermediates and transition states. It is shown that the minimum energy reaction path is consistent with an assumption of a two-step mechanism of GTP hydrolysis. At the first stage, a unified action of the nearest residues of Ras and the nearest water molecules results in a substantial spatial separation of the gamma-phosphate group of GTP from the rest of the molecule (GDP). This phase of hydrolysis process proceeds through the low barrier (16.7 kcal/mol) transition state TS1. At the second stage, the inorganic phosphate is formed in consequence of proton transfers mediated by two water molecules and assisted by the Gln61 residue from Ras. The highest transition state at this segment, TS3, is estimated to have an energy 7.5 kcal/mol above the enzyme-substrate complex. The results of simulations are compared to the previous findings for the GTP hydrolysis in the Ras-GAP (p21(ras)-p120(GAP)) protein complex. Conclusions of the modeling lead to a better understanding of the anticatalytic effect of cancer causing mutation of Gln61 from Ras, which has been debated in recent years.     

Co-operative interactions during protein folding.

The theory for measuring co-operativity between interactions in proteins by protein engineering experiments is developed by introducing a procedure for analysing increasing orders of synergy in a protein with increasing numbers of residues. The (pairwise) interaction energy (delta 2Gint) between two side-chains may be measured experimentally by a double-mutant cycle consisting of the wild-type protein, the two single mutants and the double mutant. This procedure may be extended to three residues to give a value for delta 3Gint for a triple-mutant cube, and to higher orders using multi-dimensional mutant space. We now show that delta 3Gint is the excess energy of adding all three chains compared with the sum of all the pairwise values of delta 2Gint for each of the constituent double-mutant cycles and the sum of all the single addition energies. This physical interpretation extends to higher orders of mutation. delta nGint (i.e. the interaction energy for n residues), thus, reveals the layers of synergy in interactions as a protein is built up. This procedure is applied to measuring changes in synergy during the refolding of barnase for the triad of salt-linked residues Asp8, Asp12 and Arg110, which are mutated to alanine residues. The value of delta 3Gint in the folded structure is 0.77(+/- 0.06) kcal mol-1 (i.e. the triad is 0.77 kcal mol-1 more stable than expected from the sum of the individual pairwise interactions and single contributions). The value of delta 3Gint is still significant in the transition state for unfolding (0.60(+/- 0.07) kcal mol-1) and in the folding intermediate (0.60(+/- 0.13 kcal mol-1)). These results show that synergistic interactions exist in barnase, in its transition state for unfolding and in a refolding intermediate. A direct measurement of the change of co-operativity between the folded state and the transition state for unfolding shows a decrease of 0.17(+/- 0.04) kcal mol-1, suggesting that the initial stages of protein unfolding may be accompanied by some loosening of structure in parts that still interact. The similar extent of co-operativity in the transition state for unfolding and the intermediate in refolding suggests that the intermediate is homogeneous, at least in the region of the salt-linked triad, as heterogeneity would lower the co-operativity.     

Probing the transition-state structure of dual-specificity protein phosphatases using a physiological substrate mimic

Dual-specificity phosphatases (DSPs) belong to the large family of protein tyrosine phosphatases that contain the active-site motif (H/V)CxxGxxR(S/T), but unlike the tyrosine-specific enzymes, DSPs are able to catalyze the efficient hydrolysis of both phosphotyrosine and phosphoserine/threonine found on signaling proteins, as well as a variety of small-molecule aryl and alkyl phosphates. It is unclear how DSPs accomplish similar reaction rates for phosphoesters, whose reactivity (i.e., pK(a) of the leaving group) can vary by more than 10(8). Here, we utilize the alkyl phosphate m-nitrobenzyl phosphate (mNBP), leaving-group pK(a) = 14.9, as a physiological substrate mimic to probe the mechanism and transition state of the DSP, Vaccinia H1-related (VHR). Detailed pH and kinetic isotope effects of the V/K value for mNBP indicates that VHR reacts with the phosphate dianion of mNBP and that the nonbridge phosphate oxygen atoms are unprotonated in the transition state. (18)O and solvent isotope effects indicate differences in the respective timing of the proton transfer to the leaving group and P-O fission; with the alkyl ester substrate, protonation is ahead of P-O fission, while with the aryl substrate, the two processes are more synchronous. Kinetic analysis of the general-acid mutant D92N with mNBP was consistent with the requirement of Asp-92 in protonating the ester oxygen, either in a step prior to significant P-O bond cleavage or in a concerted but asynchronous mechanism in which protonation is ahead of P-O bond fission. Collectively, the data indicate that VHR and likely all DSPs can match leaving-group potential with the timing of the proton transfer to the ester oxygen, such that diverse aryl and alkyl phosphoesters are turned over with similar catalytic efficiency.     

Kinetic Isotope Effects and Stereochemical Studies on a Ribonuclease Model: Hydrolysis Reactions of Uridine 3'-Nitrophenyl Phosphate

The reactions of a ribonuclease model substrate, the compound uridine-3'-p-nitrophenyl phosphate, have been examined using heavy-atom isotope effects and stereochemical analysis. The cyclization of this compound is subject to catalysis by general base (by imidazole buffer), specific base (by carbonate buffer), and by acid. All three reactions proceed by the same mechanistic sequence, via cyclization to cUMP, which is stable under basic conditions but which is rapidly hydrolyzed to a mixture of 2'- and 3'-UMP under acid conditions. The isotope effects indicate that the specific base-catalyzed reaction exhibits an earlier transition state with respect to bond cleavage to the leaving group compared to the general base-catalyzed reaction. Stereochemical analysis indicates that both of the base-catalyzed reactions proceed with the same stereochemical outcome. It is concluded that the difference in the nucleophile in the two base-catalyzed reactions results in a difference in the transition state structure but both reactions are most likely concerted, with no phosphorane intermediate. The (15)N isotope effects were also measured for the reaction of the substrate with ribonuclease A. The results indicate that considerably less negative charge develops on the leaving group in the transition state than for the general base-catalyzed reaction in solution. Copyright 2000 Academic Press.     

Recruitment of both uniform and differential binding energy in enzymatic catalysis: xylanases from families 10 and 11

The contributions of enzyme-substrate hydrogen-binding interactions to catalysis by two different families of xylanases were evaluated through kinetic studies with two representative wild-type enzymes, Cellulomonas fimi xylanase (Cex) and Bacillus circulans xylanase (Bcx), on a series of monodeoxygenated and monodeoxyfluorinated p-nitrophenyl xylobioside substrates. Effects of substitution in the distal (-2 subsite) sugar on kcat/Km for Cex were moderately large (up to 2.9 kcal mol-1), with no effect seen on kcat. By contrast, substantial effects upon both kcat and kcat/Km were seen for substrates modified in the proximal (-1 subsite) sugar. Very similar results were obtained with Bcx. Kinetic analyses with a series of eight mutants of Cex in which active site residues interacting with the substrate were mutated yielded complementary insights. Again, interactions with the distal (-2) sugar were seen to contribute substantially to kcat/Km (up to 3.7 kcal mol-1), thus to the formation of the glycosyl-enzyme intermediate, but not to kcat, thus to the hydrolysis of the glycosyl-enzyme. Interactions with the proximal (-1) sugar are much more significant, contributing up to 6.7 kcal mol-1 to both kcat/Km and kcat. These results together indicate that interactions with the distal sugar maintain similar magnitudes in the transition states for glycosylation and deglycosylation as well as in the glycosyl-enzyme intermediate and can be referred to as "uniform binding interactions" in the parlance of Albery and Knowles (Albery, W. J., and Knowles, J. R. (1976) Biochemistry 15, 5631-5640). Interactions with the proximal sugar are considerably stronger at the deglycosylation transition state than in the intermediate, and fall into the category of differential binding interactions. This behavior likely has its origins in the changes in ring conformation of the proximal sugar but not of the distal sugar between the ground state and the reaction transition state. Correlation of these individual interaction energies with the hydrogen-bonding pattern seen in the glycosyl-enzyme intermediate allows for the assignment of hydrogen-bond strengths to each interaction, with good correlation between the two approaches. These findings are relevant to the discussion of remote binding effects upon enzymatic catalysis.     

The unfolding pathway for Apo Escherichia coli aspartate aminotransferase is dependent on the choice of denaturant

The guanidine hydrochloride (GdnHCl) mediated denaturation pathway for the apo form of homodimeric Escherichia coli aspartate aminotransferase (eAATase) (molecular mass = 43.5 kDa/monomer) includes a partially folded monomeric intermediate, M* [Herold, M., and Kirschner, K. (1990) Biochemistry 29, 1907-1913; Birolo, L., Dal Piaz, F., Pucci, P., and Marino, G. (2002) J. Biol. Chem. 277, 17428-17437]. The present investigation of the urea-mediated denaturation of eAATase finds no evidence for an M* species but uncovers a partially denatured dimeric form, D*, that is unpopulated in GdnHCl. Thus, the unfolding process is a function of the employed denaturant. D* retains less than 50% of the native secondary structure (circular dichroism), conserves significant quaternary and tertiary interactions, and unfolds cooperatively (mD*<==>U = 3.4 +/- 0.3 kcal mol-1 M-1). Therefore, the following equilibria obtain in the denaturation of apo-eAATase: D <==> 2M 2M* <==> 2U in GdnHCl and D <==> D* <==> 2U in urea (D = native dimer, M = folded monomer, and U = unfolded state). The free energy of unfolding of apo-eAATase (D <==> 2U) is 36 +/- 3 kcal mol-1, while that for the D* 2U transition is 24 +/- 2 kcal mol-1, both at 1 M standard state and pH 7.5.     

A phosphate-binding histidine of binuclear metallophosphodiesterase enzymes is a determinant of 2',3'-cyclic nucleotide phosphodiesterase activity

Binuclear metallophosphoesterases are an enzyme superfamily defined by a shared fold and a conserved active site. Although many family members have been characterized biochemically or structurally, the physiological substrates are rarely known, and the features that determine monoesterase versus diesterase activity are obscure. In the case of the dual phosphomonoesterase/diesterase enzyme CthPnkp, a phosphate-binding histidine was implicated as a determinant of 2',3'-cyclic nucleotide phosphodiesterase activity. Here we tested this model by comparing the catalytic repertoires of Mycobacterium tuberculosis Rv0805, which has this histidine in its active site (His(98)), and Escherichia coli YfcE, which has a cysteine at the equivalent position (Cys(74)). We find that Rv0805 has a previously unappreciated 2',3'-cyclic nucleotide phosphodiesterase function. Indeed, Rv0805 was 150-fold more active in hydrolyzing 2',3'-cAMP than 3',5'-cAMP. Changing His(98) to alanine or asparagine suppressed the 2',3'-cAMP phosphodiesterase activity of Rv0805 without adversely affecting hydrolysis of bis-p-nitrophenyl phosphate. Further evidence for a defining role of the histidine derives from our ability to convert the inactive YfcE protein to a vigorous and specific 2',3'-cNMP phosphodiesterase by introducing histidine in lieu of Cys(74). YfcE-C74H cleaved the P-O2' bond of 2',3'-cAMP to yield 3'-AMP as the sole product. Rv0805, on the other hand, hydrolyzed either P-O2' or P-O3' to yield a mixture of 3'-AMP and 2'-AMP products, with a bias toward 3'-AMP. These reaction outcomes contrast with that of CthPnkp, which cleaves the P-O3' bond of 2',3'-cAMP to generate 2'-AMP exclusively. It appears that enzymic features other than the phosphate-binding histidine can influence the orientation of the cyclic nucleotide and thereby dictate the choice of the leaving group.     

The merits of bipartite transition-state mimics for inhibition of uracil DNA glycosylase

The glycosidic bond hydrolysis reaction of the enzyme uracil DNA glycosylase (UDG) occurs by a two-step mechanism involving complete bond breakage to the uracil anion leaving group in the first step, formation of a discrete glycosyl cation-uracil anion intermediate, followed by water attack in a second transition-state leading to the enzyme-bound products of uracil and abasic DNA. We have synthesized and determined the binding affinities of unimolecular mimics of the substrate and first transition-state (TS1) in which the uracil base is covalently attached to the sugar, and in addition, bimolecular mimics of the second addition transition state (TS2) in which the base and sugar are detached. We find that the bipartite mimics of TS2 are superior to the TS1 mimics. These results indicate that bipartite TS2 inhibitors could be useful for inhibition of glycosylases that proceed by stepwise reaction mechanisms.     

The association of actin and myosin in the presence of gamma-amido-ATP proceeds mainly via a complex with myosin in the closed conformation

The interaction of gamma-amido-ATP (ATPN) and its 2'(3')-O-methylanthraniloyl derivative (mantATPN) with skeletal myosin subfragment 1 (S1) and actomyosin (actoS1) was studied in stopped-flow experiments. Tryptophan fluorescence and fluorescence of the mant label or light scattering were measured simultaneously. Information about the binding of mant nucleotides was obtained from the quenching of tryptophan fluorescence by the mant label. The parameters of various kinetic models were fitted to the experimental traces. The high-fluorescence state of S1 forms with ATPN at a rate of 95 s-1 ("open-closed" transition); the transition is only slowly reversible, in contrast to the very fast equilibrium seen with its better known isomer AMPPNP [Urbanke, C., and Wray, J. (2001) Biochem. J. 358, 165-173]. The stabilization of the closed state of myosin by ATPN may be due to the formation of a complex with a pentacoordinated amido-gamma-phosphate, from which ATPN can dissociate at a rate of 0.005 s-1 or be hydrolyzed by cleavage of the beta-gamma bond at a rate of 2.5 x 10(-4) s-1. A corresponding actoS1-ATPN complex with myosin in the "closed" conformation is the first detectable intermediate in the association of actin and S1-ATPN, giving an experimental access to a state analogous to a key intermediate in the cross-bridge cycle.     

Mechanistic studies of protein tyrosine phosphatases YopH and Cdc25A with m-nitrobenzyl phosphate

Protein tyrosine phosphatases (PTPs) constitute a large family of signaling enzymes that include both tyrosine specific and dual-specificity phosphatases that hydrolyze pSer/Thr in addition to pTyr. Previous mechanistic studies of PTPs have relied on the highly activated substrate p-nitrophenyl phosphate (pNPP), an aryl phosphate with a leaving group pK(a) of 7. In the study presented here, we employ m-nitrobenzyl phosphate (mNBP), an alkyl phosphate with a leaving group pK(a) of 14.9, which mimics the physiological substrates of the PTPs. We have carried out pH dependence and kinetic isotope effect measurements to characterize the mechanism of two important members of the PTP superfamily: Yersinia PTP (YopH) and Cdc25A. Both YopH and Cdc25A exhibit bell-shaped pH-rate profiles for the hydrolysis of mNBP, consistent with general acid catalysis. The slightly inverse (18)(V/K)(nonbridge) isotope effects (0.9999 for YopH and 0.9983 for Cdc25A) indicate a loose transition state with little nucleophilic participation for both enzymes. The smaller (18)(V/K)(bridge) primary isotope effects (0.9995 for YopH and 1.0012 for Cdc25A) relative to the corresponding isotope effects for pNPP hydrolysis suggest that protonation of the leaving group oxygen at the transition state by the general acid is ahead of P-O bond fission with the alkyl substrate, while general acid catalysis of pNPP by YopH is more synchronous with P-O bond fission. The isotope effect data also confirm findings from previous studies that Cdc25A utilizes general acid catalysis for substrates with a leaving group pK(a) of >8, but not for pNPP. Interestingly, the difference in the kinetic isotope effects for the reactions of aryl phosphate pNPP and alkyl phosphate mNBP by the PTPs parallels what is observed in the uncatalyzed reactions of their monoanions. In these reactions, the leaving group is protonated in the transition state, as is the case in PTP-catalyzed reactions. Also, the phosphoryl group in the transition states of the enzymatic reactions does not differ substantially from those of the uncatalyzed reactions. These results provide further evidence that these enzymes do not change the transition state but simply stabilize it.     

Pentacoordinated phosphorus revisited by high‐level QM/MM calculations

Enzymes catalyzing phosphoryl transfer reactions are extremely efficient and are involved in crucial biochemical processes. The mechanisms of these enzymes are complex due to the diversity of substrates that are involved. The reaction can proceed through a pentacoordinated phosphorus species that is either a stable intermediate or a transition state (TS). Because of this, the first X‐ray structure of a pentacoordinated phosphorus intermediate in the β‐phosphoglucomutase enzyme aroused great interest but also much controversy. To provide new insights into the nature of that structure, we have determined the reaction path of the phosphorylation step using high‐level QM/MM calculations, and have also calculated the geometry of a complex with a transition state analogue (TSA) that has been suggested to be the actual species in the crystal. The protein crystalline environment has been modeled so as to mimic the experimental conditions. We conclude that the pentacoordinated phosphorus formed in this enzyme is not a stable species but a TS, which gives an activation energy for phosphorylation in agreement with kinetic results. We also show that the TSA is a good mimic of the true TS. We have performed a new crystallographic refinement of the original diffraction map of the pentacoordinated phosphorus structure with the MgF TSA. The new fit improves significantly with respect to the original one, which strongly supports that Allen and coworkers wrongly assigned the X‐ray structure to a pentavalent phosphorane. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

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.     

Stereochemical effects of non-ionic methylphosphonates on nucleotide conformations.

The preferred conformations of deoxyribo and ribonucleoside 3'-methylphosphonates are analysed by minimizing the conformational energy as a function of all the major parameters including the sugar ring for both the S- and R-isomers. The results show that neither the substitution nor the nature of the diastereomer affects significantly the preferred conformations compared to the naturally occurring nucleoside 3'-phosphates. The preferred range of C3'-O3' bond torsions or the phase angles of pseudorotation (P) of the sugar are unaffected. The chiral substitution on the phosphate always adopts a conformation distal to the secondary C3' carbon atom in the minimum energy conformational state. Further, it introduces certain restrictions on the preferred range of P-O3' torsions depending on the methylphosphonate configuration. Methylphosphonate, especially the S-isomer, renders the normal gauche- range of P-O3' bond torsions responsible for the stacked helical duplexes to be energetically unfavourable besides introducing a high energy barrier between trans and gauche conformations. Therefore it is suggested that duplexes with S-methylphosphonate may favour extended phosphodiester conformations. These factors explain the observed lower melting temperature as well as the downfield shifts in the 31P signals in duplexes containing the S-isomer.     

Theoretical investigation of the enzymatic phosphoryl transfer of β-phosphoglucomutase: revisiting both steps of the catalytic cycle

Enzyme catalyzed phosphate transfer is a part of almost all metabolic processes. Such reactions are of central importance for the energy balance in all organisms and play important roles in cellular control at all levels. Mutases transfer a phosphoryl group while nucleases cleave the phosphodiester linkages between two nucleotides. The subject of our present study is the Lactococcus lactis β-phosphoglucomutase (β-PGM), which effectively catalyzes the interconversion of β-D-glucose-1-phosphate (β-G1P) to β-D-glucose-6-phosphate (β-G6P) and vice versa via stabile intermediate β-D-glucose-1,6-(bis)phosphate (β-G1,6diP) in the presence of Mg(2+). In this paper we revisited the reaction mechanism of the phosphoryl transfer starting from the bisphosphate β-G1,6diP in both directions (toward β-G1P and β-G6P) combining docking techniques and QM/MM theoretical method at the DFT/PBE0 level of theory. In addition we performed NEB (nudged elastic band) and free energy calculations to optimize the path and to identify the transition states and the energies involved in the catalytic cycle. Our calculations reveal that both steps proceed via dissociative pentacoordinated phosphorane, which is not a stabile intermediate but rather a transition state. In addition to the Mg(2+) ion, Ser114 and Lys145 also play important roles in stabilizing the large negative charge on the phosphate through strong coordination with the phosphate oxygens and guiding the phosphate group throughout the catalytic process. The calculated energy barrier of the reaction for the β-G1P to β-G1,6diP step is only slightly higher than for the β-G1,6diP to β-G6P step (16.10 kcal mol(-1) versus 15.10 kcal mol(-1)) and is in excellent agreement with experimental findings (14.65 kcal mol(-1)).     

Molecular mechanics studies of dinucleoside methylphosphonates: influence of methylphosphonate and its chirality on the phosphodiester conformation.

Molecular mechanics studies are performed on single stranded as well as base paired forms of dinucleoside methylphosphonates comprising different base sequences for both the S- and R-isomers of methylphosphonate (MP). S-MP produces noticeable distortions in the geometry, locally at the phosphate center, and this enables the stereochemical feasibility of compact g- g- phosphodiester. Besides, it tends to perturb the conformations around the P-O3' and glycosyl bonds, causing minor variations in stacking interactions. In single stranded dinucleosides, the gain in adjacent base stacking interaction energies seems to be sufficient to overcome the barrier to P-O3' bond rotation arising due to S-MP...sugar interaction, and this results in transition to a compact phosphodiester (BI-type) from an initial extended phosphodiester (BII-type) conformation. Such a thing seems rather difficult in base pair constrained duplexes. Dinucleosides with R-MP behave analogous to normal phosphate duplexes as the methyl group is away from the sugar. It is found that dinucleoside methylphosphonates are energetically less favoured than the corresponding dinucleoside phosphates mainly due to the depletion of contributions from electrostatic attractive interactions involving the base and sugar with the methylphosphonate consequent to the nonionic nature of the latter. Neither S-MP nor R-MP seem to significantly alter the stereochemistry of duplex structure.     

Energetics of the calcium-transporting ATPase

A thermodynamic cycle for catalysis of calcium transport by the sarcoplasmic reticulum ATPase is described, based on equilibrium constants for the microscopic steps of the reaction shown in Equation 1 under a single set of experimental (formula; see text) conditions (pH 7.0, 25 degrees C, 100 mM KCl, 5 mM MgSO4): KCa = 5.9 X 10(-12) M2, K alpha ATP = 15 microM, Kint = 0.47, K alpha ADP = 0.73 mM, K'int = 1.7, K"Ca = 2.2 X 10(-6) M2, and Kp = 37 mM. The value of K"Ca was calculated by difference, from the free energy of hydrolysis of ATP. The spontaneous formation of an acylphosphate from Pi and E is made possible by the expression of 12.5 kcal mol-1 of noncovalent binding energy in E-P. Only 1.9 kcal mol-1 of binding energy is expressed in E X Pi. There is a mutual destabilization of bound phosphate and calcium in E-P X Ca2, with delta GD = 7.6 kcal mol-1, that permits transfer of phosphate to ADP and transfer of calcium to a concentrated calcium pool inside the vesicle. It is suggested that the ordered kinetic mechanism for the dissociation of E-P X Ca2, with phosphate transfer to ADP before calcium dissociation outside and phosphate transfer to water after calcium dissociation inside, preserves the Gibbs energies of these ligands and makes a major contribution to the coupling in the transport process. A lag (approximately 5 ms) before the appearance of E-P after mixing E and Pi at pH 6 is diminished by ATP and by increased [Pi]. This suggests that ATP accelerates the binding of Pi. The weak inhibition by ATP of E-P formation at equilibrium also suggests that ATP and phosphate can bind simultaneously to the enzyme at pH 6. Rate constants are greater than or equal to 115 s-1 for all the steps in the reaction sequence to form E-32P X Ca2 from E-P, Ca2+ and [32P]ATP at pH 7. E-P X Ca2 decomposes with kappa = 17 s-1, which shows that it is a kinetically competent intermediate. The value of kappa decreases to 4 s-1 if the intermediate is formed in the presence of 2 mM Ca2+. This decrease and inhibition of turnover by greater than 0.1 mM Ca2+ may result from slow decomposition of E-P X Ca3.     

Arginine 127 stabilizes the transition state in carboxypeptidase

Crystallographic studies suggest that Arg-127 is a key amino acid in the hydrolysis of peptides and esters by carboxypeptidase A. The guanidinium group of Arg-127 is hypothesized to stabilize the oxyanion of the tetrahedral intermediate formed by the attack of water on the scissile carbonyl bond. We have replaced this amino acid in rat carboxypeptidase A1 with lysine (R127K), methionine (R127M), and alanine (R127A), in order to define the role of Arg-127 in carboxypeptidase catalyzed hydrolysis. The wild-type and mutant enzymes were expressed in yeast and purified. Kinetic studies show that Arg-127 substitution decreases kcat for both ester and amide substrates, whereas Km is relatively unchanged; for R127M and R127A this corresponds to a 6 kcal/mol decrease in transition state stabilization of the rate-limiting step. The binding affinity for the phosphonate transition state analog, Cbz-Phe-Ala(P)-OAla, was decreased by 5.4 kcal/mol, whereas binding affinity for the ground state inhibitor, DL-benzylsuccinic acid, was decreased by only 1.7 kcal/mol for R127M. Electrostatic calculations employing a finite difference solution to the Poisson-Boltzmann equation predict that the positive charge of Arg-127 should stabilize the transition state by 6-8 kcal/mol. Therefore, the experimental and theoretical data suggest that the primary role of Arg-127 is stabilization of the transition state through electrostatic interaction with the oxyanion.