Substrate Orientation in 4-Oxalocrotonate Tautomerase and Its Effect on QM/MM Energy Profiles

Combining Evolutionary Conservation and Quantum Topological Analyses To Determine Quantum Mechanics Subsystems for Biomolecular Quantum Mechanics/Molecular Mechanics Simulations

Selection of residues and other molecular fragments for inclusion in the quantum mechanics (QM) region for QM/molecular mechanics (MM) simulations is an important step for these calculations. Here, we present an approach that combines protein sequence/structure evolution and electron localization function (ELF) analyses. The combination of these two analyses allows the determination of whether a residue needs to be included in the QM subsystem or can be represented by the MM environment. We have applied this approach on two systems previously investigated by QM/MM simulations, 4-oxalocrotonate tautomerase (4OT) and ten-eleven translocation-2 (TET2), that provide examples where fragments may or may not need to be included in the QM subsystem. Subsequently, we present the use of this approach to determine the appropriate QM subsystem to calculate the minimum energy path (MEP) for the reaction catalyzed by human DNA polymerase λ (Polλ) with a third cation in the active site. Our results suggest that the combination of protein evolutionary and ELF analyses provides insights into residue/molecular fragment selection for QM/MM simulations.

Using Atomic Confining Potentials for Geometry Optimization and Vibrational Frequency Calculations in Quantum-Chemical Models of Enzyme Active Sites

Quantum-chemical studies of enzymatic reaction mechanisms sometimes use truncated active-site models as simplified alternatives to mixed quantum mechanics molecular mechanics (QM/MM) procedures. Eliminating the MM degrees of freedom reduces the complexity of the sampling problem, but the trade-off is the need to introduce geometric constraints in order to prevent structural collapse of the model system during geometry optimizations that do not contain a full protein backbone. These constraints may impair the efficiency of the optimization, and care must be taken to avoid artifacts such as imaginary vibrational frequencies. We introduce a simple alternative in which terminal atoms of the model system are placed in soft harmonic confining potentials rather than being rigidly constrained. This modification is simple to implement and straightforward to use in vibrational frequency calculations, unlike iterative constraint-satisfaction algorithms, and allows the optimization to proceed without constraint even though the practical result is to fix the anchor atoms in space. The new approach is more efficient for optimizing minima and transition states, as compared to the use of fixed-atom constraints, and also more robust against unwanted imaginary frequencies. We illustrate the method by application to several enzymatic reaction pathways where entropy makes a significant contribution to the relevant reaction barriers. The use of confining potentials correctly describes reaction paths and facilitates calculation of both vibrational zero-point and finite-temperature entropic corrections to barrier heights.

Catalytic mechanism of 4-oxalocrotonate tautomerase: significances of protein-protein interactions on proton transfer pathways

4-Oxalocrotonate tautomerase (4-OT), a member of tautomerase superfamily, is an essential enzyme in the degradative metabolism pathway occurring in the Krebs cycle. The proton transfer process catalyzed by 4-OT has been explored previously using both experimental and theoretical methods; however, the elaborate catalytic mechanism of 4-OT still remains unsettled. By combining classical molecular mechanics with quantum mechanics, our results demonstrate that the native hexametric 4-OT enzyme, including six protein monomers, must be employed to simulate the proton transfer process in 4-OT due to protein-protein steric and electrostatic interactions. As a consequence, only three out of the six active sites in the 4-OT hexamer are observed to be occupied by three 2-oxo-4-hexenedioates (2o4hex), i.e., half-of-the-sites occupation. This agrees with experimental observations on negative cooperative effect between two adjacent substrates. Two sequential proton transfers occur: one proton from the C3 position of 2o4hex is initially transferred to the nitrogen atom of the general base, Pro1. Subsequently, the same proton is shuttled back to the position C5 of 2o4hex to complete the proton transfer process in 4-OT. During the catalytic reaction, conformational changes (i.e., 1-carboxyl group rotation) of 2o4hex may occur in the 4-OT dimer model but cannot proceed in the hexametric structure. We further explained that the docking process of 2o4hex can influence the specific reactant conformations and an alternative substrate (2-hydroxymuconate) may serve as reactant under a different reaction mechanism than 2o4hex.

Catalytic mechanism and roles of Arg197 and Thr183 in the Staphylococcus aureus sortase A enzyme

The sortase A enzyme, which catalyzes the peptidoglycan cell wall anchoring reaction of LPXTG surface proteins, has been proposed to be a universal target for therapeutic agents against Gram-positive bacteria. The catalytic mechanism of the Staphylococcus aureus sortase A enzyme has been systematically studied using molecular dynamics simulations, ONIOM(DFT:MM) calculations, and QM/MM charge deletion analysis. The catalytic roles of Arg197 and Thr183 were analyzed. Our calculations show that Arg197 has several important roles in the mechanism. It is crucial for substrate binding, and is capable of reversible shift of its hydrogen bonds between the LP and TG carbonyls of the LPXTG substrate motif, depending on the protonation state of the catalytic Cys184-His120 dyad. Arg197 stabilizes the catalytic dyad in the active ion pair form but at the same time raises the barrier to acylation by approximately 8 kcal/mol. Thr183 is also essential for the catalytic reaction in that it correspondingly lowers the barrier by the same amount via electrostatic interactions. The catalytic mechanism proceeds via proton transfer from His120, followed by nucleophilic attack from the thiolate anion of Cys184. The data thus supports the proposed reverse protonation mechanism, and disproves the hypothesis of the Arg197 generating an oxyanion hole to stabilize the tetrahedral intermediate of the reaction.

QM/MM methods for biomolecular systems

Combined quantum-mechanics/molecular-mechanics (QM/MM) approaches have become the method of choice for modeling reactions in biomolecular systems. Quantum-mechanical (QM) methods are required for describing chemical reactions and other electronic processes, such as charge transfer or electronic excitation. However, QM methods are restricted to systems of up to a few hundred atoms. However, the size and conformational complexity of biopolymers calls for methods capable of treating up to several 100,000 atoms and allowing for simulations over time scales of tens of nanoseconds. This is achieved by highly efficient, force-field-based molecular mechanics (MM) methods. Thus to model large biomolecules the logical approach is to combine the two techniques and to use a QM method for the chemically active region (e.g., substrates and co-factors in an enzymatic reaction) and an MM treatment for the surroundings (e.g., protein and solvent). The resulting schemes are commonly referred to as combined or hybrid QM/MM methods. They enable the modeling of reactive biomolecular systems at a reasonable computational effort while providing the necessary accuracy.

Critical role of substrate conformational change in the proton transfer process catalyzed by 4-oxalocrotonate tautomerase

4-Oxalocrotonate tautomerase enzyme (4-OT) catalyzes the isomerization of 2-oxo-4-hexenedioate to 2-oxo-3-hexenedioate. The chemical process involves two proton transfers, one from a carbon of the substrate to the nitrogen of Pro1 and another from this nitrogen atom to a different carbon of the substrate. In this paper the isomerization has been studied using the combined quantum mechanical and molecular mechanical method with a dual-level treatment of the quantum subsystem employing the MPW1BK density functional as the higher level. Exploration of the potential energy surface shows that the process is stepwise, with a stable intermediate state corresponding to the deprotonated substrate and a protonated proline. The rate constant of the overall process has been evaluated using ensemble-averaged variational transition state theory, including the quantized vibrational motion of a primary zone of active-site atoms and a transmission coefficient based on an ensemble of optimized reaction coordinates to account for recrossing trajectories and optimized multidimensional tunneling. The two proton-transfer steps have similar free energy barriers, but the transition state associated with the first proton transfer is found to be higher in energy. The calculations show that reaction progress is coupled to a conformational change of the substrate, so it is important that the simulation allows this flexibility. The coupled conformational change is promoted by changes in the electron distribution of the substrate that take place as the proton transfers occur.

Quantum chemical modeling of enzymatic reactions: the case of 4-oxalocrotonate tautomerase

The reaction mechanism of 4-oxalocrotonate tautomerase (4-OT) is studied using the density functional theory method B3LYP. This enzyme catalyzes the isomerisation of unconjugated alpha-keto acids to their conjugated isomers. Two different quantum chemical models of the active site are devised and the potential energy curves for the reaction are computed. The calculations support the proposed reaction mechanism in which Pro-1 acts as a base to shuttle a proton from the C3 to the C5 position of the substrate. The first step (proton transfer from C3 to proline) is shown to be the rate-limiting step. The energy of the charge-separated intermediate (protonated proline-deprotonated substrate) is calculated to be quite low, in accordance with measured pKa values. The results of the two models are used to evaluate the methodology employed in modeling enzyme active sites using quantum chemical cluster models.