Temperature-Dependent Kinetic Isotope Effects in R67 Dihydrofolate Reductase from Path-Integral Simulations

Heavy water inhibits DNA double-strand break repairs and disturbs cellular transcription, presumably via quantum-level mechanisms of kinetic isotope effects on hydrolytic enzyme reactions

Heavy water, containing the heavy hydrogen isotope, is toxic to cells, although the underlying mechanism remains incompletely understood. In addition, certain enzymatic proton transfer reactions exhibit kinetic isotope effects attributed to hydrogen isotopes and their temperature dependencies, indicative of quantum tunneling phenomena. However, the correlation between the biological effects of heavy water and the kinetic isotope effects mediated by hydrogen isotopes remains elusive. In this study, we elucidated the kinetic isotope effects arising from hydrogen isotopes of water and their temperature dependencies in vitro, focusing on deacetylation, DNA cleavage, and protein cleavage, which are crucial enzymatic reactions mediated by hydrolysis. Intriguingly, the intracellular isotope effects of heavy water, related to the in vitro kinetic isotope effects, significantly impeded multiple DNA double-strand break repair mechanisms crucial for cell survival. Additionally, heavy water exposure enhanced histone acetylation and associated transcriptional activation in cells, consistent with the in vitro kinetic isotope effects observed in histone deacetylation reactions. Moreover, as observed for the in vitro kinetic isotope effects, the cytotoxic effect on cell proliferation induced by heavy water exhibited temperature-dependency. These findings reveal the substantial impact of heavy water-induced isotope effects on cellular functions governed by hydrolytic enzymatic reactions, potentially mediated by quantum-level mechanisms underlying kinetic isotope effects.

From a binding module to essential catalytic activity: how nature stumbled on a good thing

Enzymes are complex macromolecules capable of catalyzing a wide variety of chemical reactions with high efficiency. Nonetheless, biological catalysis can be rudimentary. Here, we describe an enzyme that is built from a simple protein fold. This short protein sequence - almost a peptide - belongs to the ancient SH3 family of binding modules. Surprisingly, this binding module catalyzes the specific reduction of dihydrofolate using NADPH as a reducing cofactor, making this a dihydrofolate reductase. Too small to provide all the required binding and catalytic machinery on its own, it homotetramerizes, thus creating a large, central active site environment. Remarkably, none of the active site residues is essential to the catalytic function. Instead, backbone interactions juxtapose the reducing cofactor proximal to the target imine of the folate substrate, and a specific motion of the substrate promotes formation of the transition state. In this feature article, we describe the features that make this small protein a functional enzyme capable of catalyzing a metabolically essential reaction, highlighting the characteristics that make it a model for the evolution of primitive enzymes from binding modules.

Reaction within the coulomb-cage; science in retrospect

The Coulomb-cage is defined as the space where the electrostatic interaction between two bodies is more intensive than the thermal energy (k_BT). For small molecule, the Coulomb-cage is a small sphere, extending only few water molecules towards the bulk and its radius is sensitive to the ionic strength of the solution. For charged proteins or membranal structures, the Coulomb-cage can engulf large fraction of the surface and provides a preferred pathway for ion propagation along the surface. Similarly, electrostatic potential at the inner space of a channel can form preferential trajectories passage for ions. The dynamics of ions inside the Coulomb-cage of ions was formulated by the studies of proton-anion recombination of excited photoacids. In the present article, we recount the study of intra- Coulomb-cage reaction taking place on the surface of macro-molecular bodies like micelles, membranes, proteins and intra-protein cavities. The study progressed stepwise, tracing the dynamics of a proton ejected from a photo-acid molecule located at defined sites (on membrane, inter-membrane space, active site of enzyme, inside Large Pore Channels etc.). Accumulation of experimental observations encouraged us to study of the reaction mechanism by molecular dynamics simulations of ions within the Coulomb-cage of proteins surface or inside large pores. The intra-Coulomb-cage proton transfer events follows closely the fine structure of the electrostatic field inside the cage and reflects the shape of nearby dielectric boundaries, the temporal ordering of the solvent molecules and the structural fluctuations of the charged side chains. The article sums some 40 years of research, which in retrospect clarifies the intra-Coulomb-cage reaction mechanism.

Solvent Effects on the Temperature Dependence of Hydride Kinetic Isotope Effects: Correlation to the Donor-Acceptor Distances

Protein structural effects on the temperature (T) dependence of kinetic isotope effects (KIEs) in H-tunneling reactions have recently been used to discuss about the role of enzyme thermal motions in catalysis. Frequently observed nearly T-independent KIEs in the wild-type enzymes and T-dependent KIEs in variants suggest that H-tunneling in the former is assisted by the naturally evolved protein constructive vibrations that help sample short donor-acceptor distances (DADs) needed. This explanation that correlates the T-dependence of KIEs with DAD sampling has been highly debated as simulations following other H-tunneling models sometimes gave alternative explanations. In this paper, solvent effects on the T-dependence of KIEs of two hydride tunneling reactions of NADH/NAD⁺ analogues (represented by ΔE_a = E_aD - E_aH) were determined in attempts to replicate the observations in enzymes and test the protein vibration-assisted DAD sampling concept. Effects of selected aprotic solvents on the DAD_PRC's of the productive reactant complexes (PRCs) and the DAD_TRS's of the activated tunneling ready states (TRSs) were obtained through computations and analyses of the kinetic data, including 2° KIEs, respectively. A weaker T-dependence of KIEs (i.e., smaller ΔE_a) was found in a more polar aprotic solvent in which the system has a shorter average DAD_PRC and DAD_TRS. Further results show that a charge-transfer (CT) complexation made of a stronger donor/acceptor gives rise to a smaller ΔE_a. Overall, the shorter and less broadly distributed DADs resulting from the stronger CT complexation vibrations give rise to a smaller ΔE_a. Our results appear to support the explanation that links the T-dependence of KIEs to the donor-acceptor rigidity in enzymes.

The catalytic mechanism of the mitochondrial methylenetetrahydrofolate dehydrogenase/cyclohydrolase (MTHFD2)

Methylenetetrahydrofolate dehydrogenase/cyclohydrolase (MTHFD2) is a new drug target that is expressed in cancer cells but not in normal adult cells, which provides an Achilles heel to selectively kill cancer cells. Despite the availability of crystal structures of MTHFD2 in the inhibitor- and cofactor-bound forms, key information is missing due to technical limitations, including (a) the location of absolutely required Mg2+ ion, and (b) the substrate-bound form of MTHFD2. Using computational modeling and simulations, we propose that two magnesium ions are present at the active site whereby (i) Arg233, Asp225, and two water molecules coordinate MgA2+, while MgA2+ together with Arg233 stabilize the inorganic phosphate (Pi); (ii) Asp168 and three water molecules coordinate MgB2+, and MgB2+ further stabilizes Pi by forming a hydrogen bond with two oxygens of Pi; (iii) Arg201 directly coordinates the Pi; and (iv) through three water-mediated interactions, Asp168 contributes to the positioning and stabilization of MgA2+, MgB2+ and Pi. Our computational study at the empirical valence bond level allowed us also to elucidate the detailed reaction mechanisms. We found that the dehydrogenase activity features a proton-coupled electron transfer with charge redistribution connected to the reorganization of the surrounding water molecules which further facilitates the subsequent cyclohydrolase activity. The cyclohydrolase activity then drives the hydration of the imidazoline ring and the ring opening in a concerted way. Furthermore, we have uncovered that two key residues, Ser197/Arg233, are important factors in determining the cofactor (NADP+/NAD+) preference of the dehydrogenase activity. Our work sheds new light on the structural and kinetic framework of MTHFD2, which will be helpful to design small molecule inhibitors that can be used for cancer treatment.

Pairing structural reconstruction with catalytic competence to evaluate the mechanisms of key enzymes in the folate-mediated one-carbon pathway

Mammalian metabolism comprises a series of interlinking pathways that include two major cycles: the folate and methionine cycles. The folate-mediated metabolic cycle uses several oxidation states of tetrahydrofolate to carry activated one-carbon units to be readily used and interconverted within the cell. They are required for nucleotide synthesis, methylation and metabolism, and particularly for proliferation of cancer cells. Based on the latest progress in genome-wide CRISPR loss-of-function viability screening of 789 cell lines, we focus on the most cancer-dependent enzymes in this pathway, especially those that are hyperactivated in cancer, to provide new insight into the chemical basis for cancer drug development. Since the complete 3D structure of several of these enzymes of the one-carbon pathway in their active form are not available, we used homology modelling integrated with the interpretation of the reaction mechanism. In addition, have reconstructed the most likely scenario for the reactions taking place paired with their catalytic competence that provides a testable framework for this pathway.

Direct Correlation between Donor-Acceptor Distance and Temperature Dependence of Kinetic Isotope Effects in Hydride-Tunneling Reactions of NADH/NAD+ Analogues

Recent study of structural effects on primary kinetic isotope effects (1° KIEs) of H-transfer reactions in enzymes and solution revealed that a more rigid reaction system gave rise to a weaker temperature dependence of 1° KIEs, i.e., a smaller isotopic activation energy difference (ΔE_a = E_aD - E_aH). This has been explained within the contemporary vibrationally assisted activated H-tunneling (VA-AHT) model in which rigidity is defined according to the density of donor-acceptor distance (DAD_TRS) populations at the tunneling ready state (TRS) sampled by heavy atom motions. To test the relationship between DAD_TRS and ΔE_a in the model, we developed a computational method to obtain the TRS structures for H-transfer reactions. The method was applied to three hydride transfer reactions of NADH/NAD⁺ analogues for which the ΔE_a's as well as secondary (2°) KIEs have been reported. The 2° KIEs computed from each TRS structure were fitted to the observed values to obtain the optimal TRSs/DAD_TRS's. It was found that a shorter DAD_TRS does correspond with a smaller ΔE_a. This appears to support the VA-AHT model. Moreover, an analysis of hybridizations at the bent TRS structures shows that rehybridizations at the donor-acceptor centers are much more advanced than predicted from the classical mechanism. This implies that more orbital preparations are required for the nonclassical H-tunneling to take place.

Large Anharmonic Effects on Tunneling and Kinetics: Reaction of Propane with Muonium

Calculations of kinetic isotope effects (KIEs) provide challenging tests of quantal mass effects on reaction rates, and muonium KIEs are the most challenging. Here, we show that it can be very important to include reaction-coordinate-dependent vibrational anharmonicity along the whole reaction path to calculate tunneling probabilities and KIEs. For the reaction of propane with Mu, this decreases both the height and width of the vibrationally adiabatic potential barrier, with both effects increasing the rate constants. Our results agree well with the experimental observations.