The regioselective activation of C–H bonds in complex molecules containing several of them is still an exciting challenge in chemistry. However, many enzymes catalyze these processes, and much can be learned from the way they do it. For example, each mammalian lipoxygenase isoenzyme abstracts a hydrogen atom from essentially a unique carbon position. This paper presents a comprehensive study at the B3LYP­(QM)/CHARMM­(MM) level of the hydrogen abstraction reaction from arachidonic acid (AA) catalyzed by rabbit 15-lipoxygenase (15-rLO). Most of the products of this reaction arise from the initial hydrogen abstraction from the carbon C13 of AA. Nevertheless, we have shown that 15-rLO seems able to catalyze not only the abstraction of H13 but also that of H10. After having studied 20 of these hydrogen transfers initiated from different snapshots of an extensive MD sampling of the Michaelis complex, we have even concluded that the reaction mechanisms for both abstractions are identical (proton-coupled electron transfer processes), with transition state structures matching their geometries around the shifting hydrogen. Despite that similarity, the average potential energy barrier for the H13 abstractions is 4.0 kcal/mol lower than for the H10 abstractions, in good agreement with the experimental C13:C10 ratio of 97:3. We have found that a subtle steric hindrance by Leu597 and Ile663 is the main cause for that difference. Driving the strict regiospecificity exhibited by 15-rLO appears to be the essential function of the bulky side chains of those conserved residues, in this way making possible the vital physiological role of 15-rLO and, probably, of all the mammalian lipoxygenase isoenzymes. The understanding of how nature uses residues with the bulkiest aliphatic side chains to achieve the selective activation of C–H bonds can stimulate the design of efficient biocatalysts to that aim.

The basis set and the functional dependence of one-bond carbon-carbon NMR spin-spin coupling constants (SSCC) have been analyzed using density functional theory. Four basis sets (6-311G**, TZVP, EPR-III, and aug-cc-pVTZ-J) and four functionals (PBE, PW91, B3LYP, and B3P86) are tested by comparison with 70 experimental values corresponding to 49 molecules that represent multiple types of hybridization of the carbon atoms. The two hybrid functionals B3P86 and B3LYP combined either EPR-III or aug-cc-pVTZ-J basis sets lead to the best accuracy of calculated SSCC. However, a simple linear regression allows for the obtaining of scaled coupling constants that fit much better with the experimental data and where the differences between the different basis sets and/or functional results are significantly reduced. For large molecules the TZVP basis set can be an appropriate election presenting a good compromise between quality of results and computational cost.

Lipoxygenases (LOXs) are a family of enzymes involved in the biosynthesis of several lipid mediators. In the case of human 15-LOX, the 15-LOX-1 and 15-LOX-2 isoforms show slightly different reaction regiospecificity and substrate specificity, indicating that substrate binding and recognition may be different, a fact that could be related to their different biological role. Here, we have used long molecular dynamics simulations, QM(DFT)/MM potential energy and free energy calculations (using the newly developed DHAM method), to investigate the binding mode of the arachidonic acid (AA) substrate into 15-LOX-2 and the rate-limiting hydrogen-abstraction reaction 15-LOX-2 catalyzes. Our results strongly indicate that hydrogen abstraction from C13 in 15-LOX-2 is only consistent with the "tail-first" orientation of AA, with its carboxylate group interacting with Arg429, and that only the pro-S H13 hydrogen will be abstracted (being the pro-R H13 and H10 too far from the acceptor oxygen atom). At the B3LYP/6-31G(d) level the potential and free energy barriers for the pro-S H13 abstraction of AA by 15-LOX-2 are 18.0 and 18.6 kcal/mol, respectively. To analyze the kinetics of the hydrogen abstraction process, we determined a Markov model corresponding to the unbiased simulations along the state-discretized reaction coordinate. The calculated rates based on the second largest eigenvalue of the Markov matrices agree well with experimental measurements, and also provide the means to directly determine the pre-exponential factor for the reaction by comparing with the free energy barrier height. Our calculated pre-exponential factor is close to the value of kBT/h. On the other hand, our results suggest that the spin inversion of the complete system (including the O2 molecule) that is required to happen at some point along the full process to lead to the final hydroperoxide product, is likely to take place during the hydrogen transfer, which is a proton coupled electron transfer. Overall, a different binding mode from the one accepted for 15-LOX-1 is proposed, which provides a molecular basis for 15-LOX-2 exclusive 15-HPETE production in front of the double (although highly 15-) 12/15 regiospecificity of 15-LOX-1. Understanding how these different isoenzymes achieve their regiospecificity is expected to help in specific inhibitor design.

Arginine finger is a highly conserved and essential residue in many GTPase and AAA+ ATPase enzymes that completes the active site from a distinct protomer, forming contacts with the γ-phosphate of the nucleotide. To date, no pyrophosphatase has been identified that employs an arginine finger fulfilling all of the above properties; all essential arginine fingers are used to catalyze the cleavage of the γ-phosphate. Here, we identify and unveil the role of a conserved arginine residue in trimeric dUTPases that meets all the criteria established for arginine fingers. We found that the conserved arginine adjacent to the P-loop-like motif enables structural organization of the active site for efficient catalysis via its nucleotide coordination, while its direct electrostatic role in transition state stabilization is secondary. An exhaustive structure-based comparison of analogous, conserved arginines from nucleotide hydrolases and transferases revealed a consensus amino acid location and orientation for contacting the γ-phosphate of the substrate nucleotide. Despite the structurally equivalent position, functional differences between arginine fingers of dUTPases and NTPases are explained on the basis of the unique chemistry performed by the pyrophosphatase dUTPases.

Optimized shifting and/or scaling factors for calculating one-bond carbon-hydrogen spin-spin coupling constants have been determined for 35 combinations of representative functionals (PBE, B3LYP, B3P86, B97-2 and M06-L) and basis sets (TZVP, HIII-su3, EPR-III, aug-cc-pVTZ-J, ccJ-pVDZ, ccJ-pVTZ, ccJ-pVQZ, pcJ-2 and pcJ-3) using 68 organic molecular systems with 88 (1)JCH couplings including different types of hybridized carbon atoms. Density functional theory assessment for the determination of (1)JCH coupling constants is examined, comparing the computed and experimental values. The use of shifting constants for obtaining the calculated coupling improves substantially the results, and most models become qualitatively similar. Thus, for the whole set of couplings and for all approaches excluding those using the M06 functional, the root-mean-square deviations lie between 4.7 and 16.4 Hz and are reduced to 4-6.5 Hz when shifting constants are considered. Alternatively, when a specific rovibrational contribution of 5 Hz is subtracted from the experimental values, good results are obtained with PBE, B3P86 and B97-2 functionals in combination with HIII-su3, aug-cc-pVTZ-J and pcJ-2 basis sets.