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.

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.

Most enzymes present a catalytic mechanism where one or more proton transfer events occur coupled tightly together with the enzymatic chemical reaction. We show here that inactivating mutations decouple this proton transfer step from the phosphate cleavage reaction in dUTPase. Homotrimeric dUTPase enzymes catalyze the hydrolysis of dUTP to dUMP and pyrophosphate, using largely similar structural and functional groups as most AAA+ enzymes. dUTPases typically use a single Mg2+ ion as a cofactor in the active site that is formed by direct protein–protein contacts including all three protomers. Here we focus on the C-terminal arm structural motif, which has sequence and functional similarities to P-loop motifs and is required for catalysis. In this work, we have studied the functional roles of the C-terminal arm in ligand binding and catalysis by using QM/MM (quantum mechanics/molecular mechanics) calculations in conjunction with site-directed mutagenesis experiments. We also present a new method to assess the metal ion coordination symmetry during the catalytic reaction. Using this new implementation, we identified that the coordination symmetry follows a consistent pattern in the three systems studied, reaching the most symmetrical state near the transition states. We found that the phosphate cleavage proceeds with a concerted bimolecular (ANDN) mechanism with a loose dissociative transition state and that it is coupled with a proton transfer step involving a unanimously conserved Asp residue. We show that the main mechanistic effect of the lack of the C-terminal arm is to decouple the phosphate cleavage from the subsequent proton transfer step, resulting in a high-barrier altered reaction pathway.Keywords: ab initio; continuous symmetry measure; coupled proton transfer; DFT (density functional theory); dUTPase; one-metal ion catalytic mechanism; QM/MM (quantum mechanics/molecular mechanics)

The weighted histogram analysis method (WHAM) is widely used to obtain accurate free energies from biased molecular simulations. However, WHAM free energies can exhibit significant errors if some of the biasing windows are not fully equilibrated. To account for the lack of full equilibration, we develop the dynamic histogram analysis method (DHAM). DHAM uses a global Markov state model to obtain the free energy along the reaction coordinate. A maximum likelihood estimate of the Markov transition matrix is constructed by joint unbiasing of the transition counts from multiple umbrella-sampling simulations along discretized reaction coordinates. The free energy profile is the stationary distribution of the resulting Markov matrix. For this matrix, we derive an explicit approximation that does not require the usual iterative solution of WHAM. We apply DHAM to model systems, a chemical reaction in water treated using quantum-mechanics/molecular-mechanics (QM/MM) simulations, and the Na+ ion passage through the membrane-embedded ion channel GLIC. We find that DHAM gives accurate free energies even in cases where WHAM fails. In addition, DHAM provides kinetic information, which we here use to assess the extent of convergence in each of the simulation windows. DHAM may also prove useful in the construction of Markov state models from biased simulations in phase-space regions with otherwise low population.

Most phosphate-processing enzymes require Mg2+ as a cofactor to catalyze nucleotide cleavage and transfer reactions. Ca2+ ions inhibit many of these enzymatic activities, despite Ca2+ and Mg2+ having comparable binding affinities and overall biological abundances. Here we study the molecular details of the calcium inhibition mechanism for phosphodiester cleavage, an essential reaction in the metabolism of nucleic acids and nucleotides, by comparing Ca2+- and Mg2+ catalyzed reactions. We study the functional roles of the specific metal ion sites A and B in enabling the catalytic cleavage of an RNA/DNA hybrid substrate by B. halodurans ribonuclease (RNase) H1 using hybrid quantum-mechanics/molecular mechanics (QM/MM) free energy calculations. We find that Ca2+ substitution of either of the two active-site Mg2+ ions substantially increases the height of the reaction barrier and thereby abolishes the catalytic activity. Remarkably, Ca2+ at the A site is inactive also in Mg2+-optimized active-site structures along the reaction path, whereas Mg2+ substitution recovers activity in Ca2+-optimized structures. Geometric changes resulting from Ca2+ substitution at metal ion site A may thus be a secondary factor in the loss of catalytic activity. By contrast, at metal ion site B geometry plays a more important role, with only a partial recovery of activity after Mg2+ substitution in Ca2+-optimized structures. Ca2+-substitution also leads to a change in mechanism, with deprotonation of the water nucleophile requiring a closer approach to the scissile phosphate, which in turn increases the barrier. As a result, Ca2+ is less efficient in activating the water. As a likely cause for the different reactivities of Mg2+ and Ca2+ ions in site A, we identify differences in charge transfer to the ions and the associated decrease in the pKa of the oxygen nucleophile attacking the phosphate group.

Understanding the relationship between the adiabatic free energy profiles of chemical reactions and the underlining diabatic states is central to the description of chemical reactivity. The diabatic states form the theoretical basis of linear free energy relationships (LFERs) and thus play a major role in physical organic chemistry and related fields. However, the theoretical justification for some of the implicit LFER assumptions has not been fully established by quantum mechanical studies. This study follows our earlier works(1, 2) and uses the ab initio frozen density functional theory (FDFT) method(3) to evaluate both the diabatic and the adiabatic free energy surfaces and to determine the corresponding off-diagonal coupling matrix elements for a series of SN2 reactions. It is found that the off-diagonal coupling matrix elements are almost the same regardless of the nucleophile and the leaving group but change upon changing the central group. Furthermore, it is also found that the off-diagonal elements are basically the same in gas phase and in solution, even when the solvent is explicitly included in the ab initio calculations. Furthermore, our study establishes that the FDFT diabatic profiles are parabolic to a good approximation, thus providing a first-principles support to the origin of LFER. These findings further support the basic approximation of the empirical valence bond treatment.