On the mechanism of action of ribonuclease A: relevance of enzymatic studies with a p-nitrophenylphosphate ester and a thiophosphate ester.

It has been reported that His-119 of ribonuclease A plays a major role as an imidazolium ion acid catalyst in the cyclization/cleavage of normal dinucleotides but that it is not needed for the cyclization/cleavage of 3'-uridyl p-nitrophenyl phosphate. We see that this is also true for simple buffer catalysis, where imidazole (as in His-12 of the enzyme), but not imidazolium ion, plays a significant catalytic role with the nitrophenyl substrate, but both are catalytic for normal dinucleotides such as uridyluridine. Rate studies show that the enzyme catalyzes the cyclization of the nitrophenylphosphate derivative 47,000,000 times less effectively (kcat/kuncat) than it does uridyladenosine, indicating that approximately 50% of the catalytic free energy change is lost with this substrate. This suggests that the nitrophenyl substrate is not correctly bound to take full advantage of the catalytic groups of the enzyme and is thus not a good guide to the mechanism used by normal nucleotides. The published data on kinetic effects with ribonuclease A of substituting thiophosphate groups for the phosphate groups of normal substrates has been discussed elsewhere, and it was argued that these effects are suggestive of the classical mechanism for ribonuclease action, not the novel mechanism we have recently proposed. The details of these rate effects, including stereochemical preferences in the thiophosphate series, can be invoked as support for our newer mechanism.

Progress toward chemical accuracy in the computer simulation of condensed phase reactions.

We describe a procedure for the generation of chemically accurate computer-simulation models to study chemical reactions in the condensed phase. The process involves (i) the use of a coupled semiempirical quantum and classical molecular mechanics method to represent solutes and solvent, respectively; (ii) the optimization of semiempirical quantum mechanics (QM) parameters to produce a computationally efficient and chemically accurate QM model; (iii) the calibration of a quantum/classical microsolvation model using ab initio quantum theory; and (iv) the use of statistical mechanical principles and methods to simulate, on massively parallel computers, the thermodynamic properties of chemical reactions in aqueous solution. The utility of this process is demonstrated by the calculation of the enthalpy of reaction in vacuum and free energy change in aqueous solution for a proton transfer involving methanol, methoxide, imidazole, and imidazolium, which are functional groups involved with proton transfers in many biochemical systems. An optimized semiempirical QM model is produced, which results in the calculation of heats of formation of the above chemical species to within 1.0 kcal/mol (1 kcal = 4.18 kJ) of experimental values. The use of the calibrated QM and microsolvation QM/MM (molecular mechanics) models for the simulation of a proton transfer in aqueous solution gives a calculated free energy that is within 1.0 kcal/mol (12.2 calculated vs. 12.8 experimental) of a value estimated from experimental pKa values of the reacting species.