Open Access Article
Xinjia Yu†
a,
Yutong Li†a,
Fengxi Li
a,
Shenhan Xiea,
Liang Lic,
Hong Zhang*b,
Zhi Wang
*a and
Lei Wang
*a
aKey Laboratory of Molecular Enzymology and Engineering of Ministry of Education, School of Life Sciences, Jilin University, Changchun 130023, P. R. China. E-mail: wangzhi@jlu.edu.cn; w_lei@jlu.edu.cn
bInstitute for Interdisciplinary Biomass Functional Materials Studies, Jilin Engineering Normal University, Kaixuan Road 3050, Changchun, 130052, P. R. China. E-mail: zhanghong0825@jlenu.edu.cn
cDepartment of Industry and Information Technology of Jilin Province, 199 Jianshe Street, Kuancheng District, Changchun 130000, P. R. China. E-mail: 15810006169@163.com
First published on 12th January 2026
In this study, we report the development of an environmentally friendly artificial Vitreoscilla hemoglobin (VHb) for the synthesis of 3-amino-[1,2,4]-triazole and [4,3-a]pyridine. We employed a strategy that combines porphyrin substitution with axial ligand mutations to create a highly active VHb oxidase containing cobalt protoporphyrin IX (Co(ppIX)), while simultaneously introducing double mutations (H85Y, P54C). This artificial enzyme catalyzes the cyclization desulfurization reaction of the corresponding 2-hydrazinopyridine and isothiocyanate in PBS containing 10% DMSO (v/v) under aerobic conditions at room temperature. This method addresses the limitations of catalytic activity observed in natural hemoglobin and provides a novel pathway for the green synthesis of nitrogen-containing heterocycles. Furthermore, the porphyrin ligand substitution strategy broadens the application scope of artificial metalloenzymes in non-natural reactions.
Green foundation1. A rationally engineered artificial metalloenzyme (VHb-Co(H85Y, P54C)) enables biocatalytic green synthesis of pharmaceutical 3-amino-[1,2,4]-triazolo[4,3-a]pyridines via aqueous cyclodesulfurization, overcoming traditional route drawbacks through biocatalysis's mild, specific, low-toxicity and optimized hemoglobin active site cofactors.2. Engineered VHb-Co (H85Y, P54C) delivers a high yield for target model heterocycle substrates with a low loading of catalyst and an aqueous solvent system. The room-temperature, O2-driven process eschews toxic reagents and metals, and harnesses biocatalysis’ merits alongside optimized hemoglobin active sites and cofactors to overcome traditional synthesis flaws. 3. Subsequent studies could explore cascade reactions using the same artificial metalloenzyme to synthesize target heterocycles directly from simple, renewable precursors eliminating intermediate isolation steps |
Biocatalysis, characterized by its environmental friendliness, high efficiency, and specificity, is fundamentally revolutionizing the paradigms of traditional organic synthesis and emerging as an increasingly favored synthetic strategy in organic chemistry.6 Among these, enzymes as green biocatalysts exhibit exceptional catalytic performance, enabling a series of non-natural transformations under mild conditions and overcoming the dependencies of conventional chemical synthesis on metal catalysts, toxic reagents, and harsh reaction conditions.7 Notably, hemoglobins, a class of multifunctional biomacromolecules containing heme cofactors, have recently emerged as a research hotspot in the field of green biocatalysis due to their unique catalytic potential. They efficiently mediate various reaction types, such as epoxidation, hydroxylation, and carbene transfer, through their iron-porphyrin active sites while leveraging their distinctive protein cavities to precisely control the stereoselectivity and regioselectivity of reactions.8 However, the catalytic activity of natural hemoglobin is constrained by the inherent reactivity of its heme iron center. Through directed substitution with artificial metal cofactors, one can precisely regulate the redox potential, ligand field effect, and steric hindrance of the metal center, thereby endowing the protein with novel catalytic functions.9
For instance, in 2014, Fasan's group replaced the heme cofactor of myoglobin (Mb) with Mn(ppIX) or Co(ppIX) to enable efficient catalysis of intramolecular C(sp3)–H amination of aryl sulfonyl azides (Fig. 1c-1). In 2017, they further confirmed that engineered Mb conjugated with chlorin iron e6 could catalyze stereoselective cyclopropanation of olefins (Fig. 1c-2).10 In light of recent advancements in hemoglobin protein engineering, our team has conducted related research. We replaced the cofactor of Vitreoscilla hemoglobin with various metal porphyrins and, utilizing directed evolution techniques, developed a novel artificial metal oxidase. These artificial metal enzymes were employed to efficiently synthesize a range of nitrogen-containing heterocycles, including 2-amino benzo[d]oxazole/thiazole and thiazole.11 Building on our group's prior research findings, VHb and its artificial metal enzymes have exhibited significant potential in catalyzing novel non-natural reactions. Consequently, in this study, we further explored the application of VHb artificial metal oxidase in the cyclization desulfurization reaction, successfully achieving the efficient synthesis of 3-amino-[1,2,4]-triazolo[4,3-a]pyridine (Fig. 1d). Notably, this reaction was performed at room temperature using PBS as the reaction medium and O2 as the oxidant. Compared to previously reported methods, our approach presents substantial advantages in terms of low toxicity and high process safety, while minimizing the reliance on organic solvents.
To further enhance catalytic activity, we focused on engineering modifications of the active site of VHb. By integrating data from previous studies on VHb mutants, we constructed an artificial mutant enzyme library. These mutations may alter the coordination environment with metal porphyrin cofactors and the binding modes with reactants by changing the amino acid residues at the active site of VHb. Among these residues, the proximally conserved histidine (H85) is crucial for VHb function, as it not only directly coordinates axially with the metal ligands but also fine-tunes the electronic density on the metal through the electronic effects of its side chain, thereby further influencing the redox properties of VHb.13 Consequently, axial ligands represent a key factor that dictates the catalytic potential of the active site. Based on this structural analysis, we first selected VHb variants with mutations at the H85 residue to recombine with Co(ppIX) and evaluated the catalytic efficiency of the mutants and recombinant proteins against substrates 1a and 2a, respectively (Fig. 2c-1).
The results indicated that the vast majority of H85 mutants were unable to establish effective axial coordination with Co(ppIX), resulting in catalytic product yields ranging from 6% to 25%. Notably, the catalytic product yields of these H85 mutant enzymes were comparable to those of the wild-type VHb. In contrast, the VHb (H85Y) variant exhibited significantly enhanced catalytic performance, with a yield of 50%, and an even higher yield of 61% when bound to Co(ppIX). This substantial difference can be attributed to the fact that the phenolic hydroxyl side chain of tyrosine forms more stable axial coordination bonds with the cobalt at the center of Co(ppIX) through its oxygen atoms, compared to histidine. Additionally, the moderate electron-donating ability of tyrosine adjusts the electron density at the center of Co(ppIX) to an optimal level for substrate activation and electron transfer. Consequently, we performed further screening based on the H85Y variant (Fig. 2c-2). The results revealed that, except for the double mutants VHb-Co (H85Y, P54L) and VHb-Co (H85Y, P54C), the catalytic performance of all other double mutants was inferior to that of the single mutant VHb (H85Y). These findings suggest that mutations at non-axial sites can only fine-tune the microenvironment of the active cavity and cannot substitute for the essential role of functional axial ligands. Notably, the VHb-Co (H85Y, P54C) variant achieved a catalytic product yield of 84% and a turnover number (TON) of 1680 (Table S2, entry 13), demonstrating substantially enhanced catalytic performance that significantly outperforms all other engineered variants. Sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) exhibited a single, distinct band at the expected molecular weight (17 kDa), which matched the theoretical molecular weight of the VHb variant. The electrophoretic patterns obtained confirmed the high homogeneity of the protein preparations. Collectively, these results confirm that the mutant enzyme achieved proper folding and maintained structural integrity throughout expression and purification, thereby providing a critical structural foundation for its superior catalytic performance. Detailed data for screening VHb mutants, including comprehensive catalytic performance data for both single and double mutants within the constructed enzyme library, are provided in Table S2 of the SI.
To further understand the reaction pathway and provide a theoretical basis for exploring substrate applicability, we conducted systematic mechanistic investigations of the cyclic desulfurization reaction catalyzed by VHb-Co (H85Y, P54C). Subsequently, a series of control experiments were performed, and the results are presented in Fig. 3a. First, a radical trapping experiment was conducted in which three equivalents of the radical scavenger TEMPO were added to the reaction system while keeping other parameters unchanged (Fig. 3a-1). The results indicated that only a trace amount of the target product 3a was detected, even after continuous stirring for 12 hours under standard conditions. This finding suggests that the reaction may involve radical intermediates. When nitrogen (N2) was used to replace oxygen (O2) in the system, no target product 3a was detected (Fig. 3a-2), indicating that O2 is essential for maintaining the catalytic activity of VHb-Co (H85Y, P54C). Furthermore, under catalyst-free and oxygen-free conditions, intermediate I was obtained with a yield of 72% (Fig. 3a-3); however, after adding VHb-Co (H85Y, P54C) and bubbling air, intermediate I was ultimately converted to 3a with a yield of 86% (Fig. 3a-4). In summary, intermediate II is a key intermediate in this reaction pathway. Subsequently, utilizing the property of hydrogen peroxide (H2O2) to rapidly oxidize ferrous iron (Fe2+) to ferric iron (Fe3+), Mohr's salt solution was added to detect the generation of hydrogen peroxide during the reaction14 (Fig. 3b). Under standard reaction conditions, with a total reaction duration of 12 hours, the reaction system was centrifuged after 8 hours. Following centrifugation, Mohr's salt solution was added to the supernatant, and the rapid formation of ferric hydroxide flocs was immediately observed. After further centrifugation, distinct precipitates were noted. This result indicates that H2O2 is produced during this reaction process.
To further investigate substrate binding to proteins and analyze the stability and interactions of substrate-protein complexes, we docked intermediate II into the active cavities of WT-VHb-Co and VHb-Co (H85Y, P54C) respectively using the AutoDock Vina tool within Chimera (the structure of VHb-Co (H85Y, P54C) was constructed by homology modeling strategy using SWISS-MODEL).15 Fig. 3d illustrates the interactions between intermediate II and both WT-VHb-Co and VHb-Co (H85Y, P54C). As shown in Fig. 3d-1, within the active pocket of WT-VHb-Co, intermediate II is only embedded in the marginal region of the hydrophobic cavity, exhibiting limited interactions with surrounding residues. Additionally, the sulfur atom on intermediate II is at a relatively large distance from the metal center of Co(ppIX) (5.7 Å), and this steric hindrance significantly restricts the efficiency of the oxidation step. In contrast, as a key axial amino acid residue in the active cavity that is crucial for VHb-substrate binding and mediating the catalytic reaction, the H85Y mutation exerts a synergistic effect together with the P54C mutation, which plays an auxiliary optimization role, collectively inducing significant conformational rearrangement and optimization of the active cavity in VHb-Co (H85Y, P54C) and thereby forming a more compact and highly specific interaction network with intermediate II. Specifically, this structural adjustment stabilizes the conformation of the intermediate through hydrogen bonds, pi-sulfur interactions, and pi-alkyl interactions among the surrounding residues. Notably, the H85Y axial mutation precisely optimizes the spatial arrangement of the active cavity along the axial dimension, directly reducing the steric hindrance of intermediate II and effectively shortening the distance between intermediate II and the Co(ppIX) metal center. These conformational changes decrease the steric hindrance of intermediate II, bringing it closer to the metal center of Co(ppIX) (3.8 Å), which may further enhance catalytic efficiency (Fig. 3d-2).
Based on literature reports and the aforementioned experimental results, a plausible reaction mechanism for the VHb-Co (H85Y, P54C)-catalyzed synthesis of 3-amino-[1,2,4]-triazolo[4,3-a]pyridine compounds is proposed (Fig. 3c). Initially, 2-hydrazinopyridine (1a) undergoes a spontaneous nucleophilic addition reaction with phenyl isothiocyanate (2a), resulting in the formation of thiourea intermediate I. Under aerobic conditions, the Co(ppIX) cofactor interacts with O2 to generate a Co3+–superoxide intermediate, which subsequently acidifies intermediate I. The sulfur radical II is formed via an electron transfer process. Intermediate II is then converted to intermediate III through proton transfer, followed by intramolecular cyclization of intermediate III to form intermediate IV, accompanied by the generation of superoxide anions. Finally, intermediate IV undergoes proton exchange with the solvent, ultimately completing the cyclodesulfurization reaction to yield the target product 3a.
Following the above screening and mechanistic clarification, the substrate scope of the reaction was explored using the VHb-Co (H85Y, P54C) mutant (Fig. 4). The applicability of substituents on the N-atom of isothiocyanates was first evaluated. Experiments showed that under standard reaction conditions, 2-hydrazinopyridine (1a) and phenyl isothiocyanate derivatives (2) could be efficiently converted to a series of 1,2,4-triazolopyridine derivatives (3) catalyzed by the biocatalyst. For monosubstituted phenyl isothiocyanates, ortho-, meta-, and para-substituents had minimal effect on the yield (3b–3d, 63%–72%). The influence of electron-donating (EDG) and electron-withdrawing (EWG) groups on the reaction yield was negligible (3e–3g), and even with a strong EWG like CF3, the corresponding product yield still reached 68%. The yield decrease for 3h and 3i was likely attributed to the steric hindrance of 2,6-dimethyl and 2,4,6-trimethyl substituents, which hindered the substrate from binding to the active cavity of VHb and is consistent with the molecular docking results that steric hindrance restricts substrate-enzyme interaction efficiency. Additionally, the lower catalytic product yields with aliphatic isothiocyanates (3j–3l, 51%–62%) might be due to the poor reactivity of aliphatic substrates. Furthermore, other 2-hydrazinopyridine substrates were evaluated in this study. The results indicated that the nature of substituents on the pyridine ring exerted a significant effect on the reaction efficiency, 2-hydrazinopyridines bearing electron-donating groups or halogen substituents afforded good catalytic product yields (3m–3p, 71%–79%). In contrast, substrates with strong electron-withdrawing groups (such as –CN and –CF3) showed a marked decrease in product yields. This discrepancy was attributed to the electron-withdrawing groups significantly reducing the electron density of the pyridine ring, thereby decreasing substrate reactivity and inhibiting the cyclodesulfurization reaction (3q–3r, 34%–45%).
These results are consistent with the proposed mechanism, in which electron transfer and the formation of intermediates are highly dependent on the electron density of the substrate. Furthermore, we systematically evaluated the industrial application potential of the cyclodesulfurization reaction catalyzed by VHb-Co (H85Y, P54C), and this study investigated the applicability of the VHb-catalyzed reaction at the gram-scale. Specifically, a gram-scale experiment was performed using 5.0 mmol of 1a and 2a, with both the catalyst and solvent scaled up proportionally from the small-scale reaction. The mixture was stirred at 800 rpm in a 500 mL glass vessel. The target product 3a was successfully synthesized with a yield of 79%.
Footnote |
| † These authors contributed equally to this work. |
| This journal is © The Royal Society of Chemistry 2026 |