Exploring nucleoside hydrolase catalysis in silico: molecular dynamics study of enzyme-bound substrate and transition state

Structure, Oligomerization and Activity Modulation in N-Ribohydrolases

Enzymes catalyzing the hydrolysis of the N-glycosidic bond in nucleosides and other ribosides (N-ribohydrolases, NHs) with diverse substrate specificities are found in all kingdoms of life. While the overall NH fold is highly conserved, limited substitutions and insertions can account for differences in substrate selection, catalytic efficiency, and distinct structural features. The NH structural module is also employed in monomeric proteins devoid of enzymatic activity with different physiological roles. The homo-oligomeric quaternary structure of active NHs parallels the different catalytic strategies used by each isozyme, while providing a buttressing effect to maintain the active site geometry and allow the conformational changes required for catalysis. The unique features of the NH catalytic strategy and structure make these proteins attractive targets for diverse therapeutic goals in different diseases.

Structural explanation for the tunable substrate specificity of an E. coli nucleoside hydrolase: insights from molecular dynamics simulations

Parasitic protozoa rely on nucleoside hydrolases that play key roles in the purine salvage pathway by catalyzing the hydrolytic cleavage of the N-glycosidic bond that connects nucleobases to ribose sugars. Cytidine-uridine nucleoside hydrolase (CU-NH) is generally specific toward pyrimidine nucleosides; however, previous work has shown that replacing two active site residues with Tyr, specifically the Thr223Tyr and Gln227Tyr mutations, allows CU-NH to process inosine. The current study uses molecular dynamics (MD) simulations to gain atomic-level insight into the activity of wild-type and mutant E. coli CU-NH toward inosine. By examining systems that differ in the identity and protonation states of active site catalytic residues, key enzyme-substrate interactions that dictate the substrate specificity of CU-NH are identified. Regardless of the wild-type or mutant CU-NH considered, our calculations suggest that inosine binding is facilitated by interactions of the ribose moiety with active site residues and Ca²⁺, and π-interactions between two His residues (His82 and His239) and the nucleobase. However, the lack of observed activity toward inosine for wild-type CU-NH is explained by no residue being correctly aligned to stabilize the departing nucleobase. In contrast, a hydrogen-bonding network between hypoxanthine and a newly identified general acid (Asp15) is present when the two Tyr mutations are engineered into the active site. Investigation of the single CU-NH mutants reveals that this hydrogen-bonding network is only maintained when both Tyr mutations are present due to a π-interaction between the residues. These results rationalize previous experiments that show the single Tyr mutants are unable to efficiently hydrolyze inosine and explain how the Tyr residues work synergistically in the double mutant to stabilize the nucleobase leaving group during hydrolysis. Overall, our simulations provide a structural explanation for the substrate specificity of nucleoside hydrolases, which may be used to rationally develop new treatments for kinetoplastid diseases.

QM/MM and MM MD Simulations on the Pyrimidine-Specific Nucleoside Hydrolase: A Comprehensive Understanding of Enzymatic Hydrolysis of Uridine

The pyrimidine-specific nucleoside hydrolase Yeik (CU-NH) from Escherichia coli cleaves the N-glycosidic bond of uridine and cytidine with a 10²-10⁴-fold faster rate than that of purine nucleoside substrates, such as inosine. Such a remarkable substrate specificity and the plausible hydrolytic mechanisms of uridine have been explored by using QM/MM and MM MD simulations. The present calculations show that the relatively stronger hydrogen-bond interactions between uridine and the active-site residues Gln227 and Tyr231 in CU-NH play an important role in enhancing the substrate binding and thus promoting the N-glycosidic bond cleavage, in comparison with inosine. The estimated energy barrier of 30 kcal/mol for the hydrolysis of inosine is much higher than 22 kcal/mol for uridine. Extensive MM MD simulations on the transportation of substrates to the active site of CU-NH indicate that the uridine binding is exothermic by ∼23 kcal/mol, more remarkable than inosine (∼12 kcal/mol). All of these arise from the noncovalent interactions between the substrate and the active site featured in CU-NH, which account for the substrate specificity. Quite differing from other nucleoside hydrolases, here the enzymatic N-glycosidic bond cleavage of uridine is less influenced by its protonation.

Active site plasticity revealed from the structure of the enterobacterial N-ribohydrolase RihA bound to a competitive inhibitor

Background

Pyrimidine-preferring N-ribohydrolases (CU-NHs) are a class of Ca²⁺-dependent enzymes that catalyze the hydrolytic cleavage of the N-glycosidic bond in pyrimidine nucleosides. With the exception of few selected organisms, their physiological relevance in prokaryotes and eukaryotes is yet under investigation.

Results

Here, we report the first crystal structure of a CU-NH bound to a competitive inhibitor, the complex between the Escherichia coli enzyme RihA bound to 3, 4-diaminophenyl-iminoribitol (DAPIR) to a resolution of 2.1 Å. The ligand can bind at the active site in two distinct orientations, and the stabilization of two flexible active site regions is pivotal to establish the interactions required for substrate discrimination and catalysis.

Conclusions

A comparison with the product-bound RihA structure allows a rationalization of the structural rearrangements required for an enzymatic catalytic cycle, highlighting a substrate-assisted cooperative motion, and suggesting a yet overlooked role of the conserved His82 residue in modulating product release. Differences in the structural features of the active sites in the two homologous CU-NHs RihA and RihB from E. coli provide a rationale for their fine differences in substrate specificity. These new findings hint at a possible role of CU-NHs in the breakdown of modified nucleosides derived from RNA molecules.

Effects of nucleophile, oxidative damage, and nucleobase orientation on the glycosidic bond cleavage in deoxyguanosine

Deglycosylation of nucleotides occurs during many essential biological processes, including DNA repair, and is initiated by a variety of nucleophiles. In the present work, density functional theory (B3LYP) was used to investigate the thermodynamics and kinetics of the glycosidic bond cleavage reaction in the model nucleoside forms of guanine and its major oxidation product, 8-oxoguanine. Base excision facilitated by four different nucleophiles (hydroxyl anion (fully activated water), formate-water complex (partially activated water), lysine, and proline) was considered, which spans nucleophiles involved in a collection of spontaneous and enzyme-catalyzed processes. Because some enzymes that catalyze deglycosylation can accommodate more than one orientation of the base with respect to the sugar moiety, the effects of the (anti/syn) base orientation on the barrier height were also considered. We find that the nucleophile has a very large effect on the overall (gas-phase) reaction energetics. Although this effect decreases in different (polar) environments, the nucleophile has the greatest influence on the overall reaction as compared to whether the base is damaged or to the base orientation. Furthermore, the effects are significant in environments that most closely resemble (nonpolar) enzymatic active sites. Our results provide a greater understanding of the relative effects of the nucleophile, damage to the nucleobase, and the nucleobase orientation with respect to the sugar moiety on the deglycosylation pathway, which provide qualitative explanations for relative base excision rates observed in some biological systems.

Modeling the dissociative hydrolysis of the natural DNA nucleosides

Two-dimensional PCM-B3LYP/6-31+G(d) potential energy surfaces for the hydrolysis of the four natural 2'-deoxyribonucleosides (2'-deoxyadenosine, 2'-deoxyguanosine, 2'-deoxycytidine, and thymidine) are characterized using a model that includes both implicit (bulk) solvent effects and (three or four) explicit water molecules in the optimization routine. For the first time, the experimentally predicted dissociative (S(N)1) mechanism is found to be favored over the synchronous (S(N)2) pathway for all nucleosides studied. Due to the success of our model in stabilizing the charge-separated intermediates along the S(N)1 pathway, it is proposed that the new model presented here is the smallest system capable of generating the experimentally predicted oxacarbenium cation intermediate. We therefore stress that dissociative mechanisms should be studied with methodologies that account for the (bulk) environment in the optimization routine, where these effects are often only included as a correction to the energy in the current literature. In addition to accounting for charge stabilization through implicit solvation, nucleophile activation and leaving group stabilization should also be explicitly introduced into the model to further stabilize the system. Our work also emphasizes the importance of studying the Gibbs surface, which in some cases provides a better description of chemically important regions of the reaction surface or changes the calculated trend in the magnitude of dissociative barriers. In addition, it is proposed that the methodology presented in this study can be used to calculate uncatalyzed deglycosylation barriers for a range of DNA nucleosides, which when compared to the corresponding enzyme-catalyzed reactions, will allow the prediction of the rate enhancement (barrier reduction) due to the enzyme.

Glycosidic bond cleavage in deoxynucleotides — A density functional study

Density functional theory was used to study the glycosidic bond cleavage in deoxynucleotides with the main goal to determine the effects of the nucleobase, hydrogen bonding with the nucleobase, and the (bulk) environment on the reaction energetics. Since direct glycosidic bond cleavage is a high-energy process, two nucleophile models were considered (HCOO–···H2O and HO–), which represent different stages of activation of a water nucleophile. The glycosidic bond cleavage barriers were found to decrease, while the reaction exothermicity increases, with an increase in the nucleobase acidity. The gas-phase barriers and reaction energies for bond cleavage in all deoxynucleotides were found to be significantly affected by hydrogen-bonding interactions with the nucleobase (by up to 30 kJ mol–1 depending on the nucleophile). Although the barriers increase and reaction energies become less exothermic in enzymatic and aqueous environments, the effects of the bulk environment are similar in the presence and absence of small molecules bound to the nucleobase. Therefore, the effects of hydrogen bonding with the bases are approximately the same in all environments. Our results suggest that hydrogen bonding with the nucleobase may play an important role in the glycosidic bond cleavage in both pyrimidine and purine nucleotides in a variety of environments.

Transition-state Complex of the Purine-specific Nucleoside Hydrolase of T.vivax: Enzyme Conformational Changes and Implications for Catalysis

Nucleoside hydrolases cleave the N-glycosidic bond of ribonucleosides. Crystal structures of the purine-specific nucleoside hydrolase from Trypanosoma vivax have previously been solved in complex with inhibitors or a substrate. All these structures show the dimeric T. vivax nucleoside hydrolase with an "open" active site with a highly flexible loop (loop 2) in its vicinity. Here, we present the crystal structures of the T. vivax nucleoside hydrolase with both soaked (TvNH-ImmH(soak)) and co-crystallised (TvNH-ImmH(co)) transition-state inhibitor immucillin H (ImmH or (1S)-1-(9-deazahypoxanthin-9-yl)-1,4-dideoxy-1,4-imino-D-ribitol) to 2.1 A and 2.2 A resolution, respectively. In the co-crystallised structure, loop 2 is ordered and folds over the active site, establishing previously unobserved enzyme-inhibitor interactions. As such this structure presents the first complete picture of a purine-specific NH active site, including leaving group interactions. In the closed active site, a water channel of highly ordered water molecules leads out from the N7 of the nucleoside toward bulk solvent, while Trp260 approaches the nucleobase in a tight parallel stacking interaction. Together with mutagenesis results, this structure rules out a mechanism of leaving group activation by general acid catalysis, as proposed for base-aspecific nucleoside hydrolases. Instead, the structure is consistent with the previously proposed mechanism of leaving group protonation in the T. vivax nucleoside hydrolase where aromatic stacking with Trp260 and an intramolecular O5'-H8C hydrogen bond increase the pKa of the N7 sufficiently to allow protonation by solvent. A mechanism that couples loop closure to the positioning of active site residues is proposed based on a comparison of the soaked structure with the co-crystallized structure. Interestingly, the dimer interface area increases by 40% upon closure of loop 2, with loop 1 of one subunit interacting with loop 2 of the other subunit, suggesting a relationship between the dimeric form of the enzyme and its catalytic activity.

Catalysis by nucleoside hydrolases

Nucleoside hydrolases cleave the N-glycosidic bond of ribonucleosides. Because of their vital role in the protozoan purine salvage pathway, nucleoside hydrolases from parasitic protozoa in particular have been studied extensively by X-ray crystallography, kinetic methods and site-directed mutagenesis. An elaborate network of conserved interactions between the metalloenzyme and the ribose enables steric and electrostatic stabilisation of the oxocarbenium-ion-like transition state. Activation of the leaving group by protonation before the formation of the transition state is a recurring catalytic strategy of enzymes that cleave N-glycosidic bonds. However, the mechanisms underlying leaving group activation are still the subject of debate for the nucleoside hydrolases.