Calculations enzymatic reactions

A final important area is the calculation of free energies with quantum mechanical models [72] or hybrid quanmm mechanics/molecular mechanics models (QM/MM) [9]. Such models are being used to simulate enzymatic reactions and calculate activation free energies, providing unique insights into the catalytic efficiency of enzymes. They are reviewed elsewhere in this volume (see Chapter 11). 

A non-linear regression analysis is employed using die Solver in Microsoft Excel spreadsheet to determine die values of and in die following examples. Example 1-5 (Chapter 1) involves the enzymatic reaction in the conversion of urea to ammonia and carbon dioxide and Example 11-1 deals with the interconversion of D-glyceraldehyde 3-Phosphate and dihydroxyacetone phosphate. The Solver (EXAMPLEll-l.xls and EXAMPLEll-3.xls) uses the Michaehs-Menten (MM) formula to compute v i- The residual sums of squares between Vg(,j, and v j is then calculated. Using guessed values of and the Solver uses a search optimization technique to determine MM parameters. The values of and in Example 11-1 are ... 

The values determined from Figure 5.23 agree well with the values calculated from the equations (Table 5.5), with an error of 3.81% for the slope and 4.65% for the intersect, respectively. The obtained experimental data were consistent with the proposed enzymatic reaction and the reaction mechanisms with uncompetitive substrate inhibition and the noncompetitive product inhibition model. 

The simplest method of coupling enzymatic reactions to electrochemical detection is to monitor an off-line reaction using FIAEC or LCEC. The enzymatic reaction is carried out in a test tube under controlled conditions with aliquots being taken at timed intervals. These aliquots are then analyzed for the electroactive product and the enzyme activity in the sample calculated from the generated kinetic information. 

Another method that has been applied by our group to the study of enzymatic reactions is the Effective Fragment Potential (EFP) method [19]. The EFP method (developed at Mark Gordon s group at Iowa State University) allows the explicit inclusion of environment effects in quantum chemical calculations. The solvent, which may consist of discrete solvent molecules, protein fragments or other material, is treated explicitly using a model potential that incorporates electrostatics, polarization, and exchange repulsion effects. The solute, which can include some... 

We have also presented the development of a method that enables the determination of free energy profiles associated with calculated MEPs for enzymatic reaction... 

Carloni et al.91 applied the DFT(PZ) calculations to investigate the electronic structure of various models of oxydized and reduced Cu, Zn superoxide dismutase. The first stage of the enzymatic reaction involves the electron transfer from Cu" ion to superoxide. The theoretical investigations provided a detailed description of the electronic structure of the molecules involved in the electron transfer process. The effect of charged groups, present in the active center, on the electron transfer process were analyzed and the Argl41 residue was shown to play a crucial role. 

A more complete list of early applications of QM/MM methods to enzymatic reactions can be found elsewhere [18, 35, 83, 84], Gao [85] has reviewed QM/MM studies of a variety of solution phenomena. QM/MM methods have also been used to study the spectra of small molecules in different solvents [86] and electrochemical properties of photosynthetic reaction centers within a protein environment [87-89], An approach has also been developed for calculation of NMR shielding tensors by use of a QM/ MM method [90]. 

Kinetic complexity definition, 43 Klinman s approach, 46 Kinetic isotope effects, 28 for 2,4,6-collidine, 31 a-secondary, 35 and coupled motion, 35, 40 in enzyme-catalyzed reactions, 35 as indicators of quantum tunneling, 70 in multistep enzymatic reactions, 44-45 normal temperature dependence, 37 Northrop notation, 45 Northrop s method of calculation, 55 rule of geometric mean, 36 secondary effects and transition state, 37 semiclassical treatment for hydrogen transfer,... 

From the second-order rate constant for imidazole-catalysed cyclization of the ethyl ester (34) (8 x 10 M s" ) and the rate constant for acylation of a-chy mo trypsin by N-acetyltyrosine ethyl ester (1600 s ), it can be calculated that, in order to attain a rate constant of the magnitude seen in the a-chymotrypsin reaction, a neighbouring imidazole would have to possess an effective molarity of 200,000 M. An effective concentration of this magnitude is not unreasonable, but it is probable that other factors are also important in the enzymatic reaction. 

Oxidation by direct H transfer from the a-carbon of alcohols to the pyrroloquinoline quinone (PQQ) cofactor of alcohol dehydrogenases was studied using ab initio quantum mechanical methods <2001JCC1732>. Energies and geometries were calculated at the 6-31G(d,p) level of theory, results were compared to available structural and spectroscopic data, and the role of calcium in the enzymatic reaction was explored. Transition state searches at the semi-empirical and STO-3G(d) level of theory provided evidence that direct transfer from the alcohol to C-5 of PQQ is energetically feasible. 

A requirement for determination of optical purity (P) is that the specific rotation of the pure enantiomer, [a]max, is known with certainty. This maximum rotation can be established by calculation, e.g., via competitive reaction methods41,42, or by direct determination employing an enantiomerically pure sample. Enantiomerically pure samples are generally believed to arise from enzymatic reactions performed on biogenic substrates, an assumption which in some instances is incorrect31, or are obtained in most cases by crystallization2. As many direct, nonchiroptical methods are available for determining enantiomeric purities, maximum rotations can also be extrapolated from the specific rotation of a sample of known ee. 

As the next enzymes of the heme pathway are active, uroporphyrinogen is converted to heptacarboxyporphyrinogen during the enzymatic reaction. This can be overcome if the incubation temperature is raised to 45°C, which inactivates uroporphyrinogen decarboxylase. Yet, 37°C is the standard temperature for enzyme activity determinations. Despite our neglect of the additionally formed porphyrins by uroporphyrinogen decarboxylase in the activity calculations, we found a highly reliable test performance. 

Using these data, calculate the protein-normalized F max and ANP MM for this enzymatic reaction. 

The mechanism of an enzymatic reaction is ultimately defined when all the intermediates, complexes, and conformational states of the enzyme are characterized and the rate constants for their interconversion are determined. The task of the kineticist in this elucidation is to detect the number and sequence of these intermediates and processes, define their approximate nature (that is, whether covalent intermediates are formed or conformational changes occur), measure the rate constants, and, from studying pH dependence, search for the participation of acidic and basic groups. The chemist seeks to identify the chemical nature of the intermediates, by what chemical paths they form and decay, and the types of catalysis that are involved. These results can then be combined with those from x-ray diffraction and NMR studies and calculations by theoretical chemists to give a complete description of the mechanism. 

---