Improving stereoselectivity of phosphotriesterase (PTE) for kinetic resolution of chiral phosphates

Specific stereoisomer is paramount as it is vital for optimizing drug efficacy and safety. The quest for the isolation of desired stereoisomer of active pharmaceutical ingredients or key intermediates drives innovation in drug synthetic and biocatalytic methods. Chiral phosphoramidate is an important building block for the synthesis of antiviral drugs such as remdesivir and sofosbuvir. Given the clinical potency of the (Sp)-diastereomer of the drugs, an enzyme capable of completely hydrolyzing the (Rp)-diastereomer is needed to achieve the purified diastereomers via biocatalytic reaction. In this study, protein engineering of phosphotriesterase (PTE) was aimed to improve the specificity. Employing rational design and site-directed mutagenesis, we generated a small library comprising 24 variants for activity screening. Notably, W131M and I106A/W131M variants demonstrated successful preparation of pure (Sp)-diastereomer of remdesivir and sofosbuvir precursors within a remarkably short hydrolysis time (<20 min). Our work unveils a promising methodology for producing pure stereoisomeric compounds, utilizing novel biocatalysts to enable the chemoenzymatic synthesis of phosphoramidate nucleoside prodrugs.


Introduction
Biocatalysis has emerged as a highly valuable tool, revolutionizing the production of intermediates, desired pure diastereomer building blocks, and active pharmaceutical ingredients (APIs) in the pharmaceutical industry (Alcántara et al., 2022;Rossino et al., 2022).Enzymes, as biocatalysts, offer several advantages over traditional chemical synthesis.Due to their specificity, efficiency, and environmental friendliness, enzymes promote greener and more efficient manufacturing processes (Adams et al., 2019;Sheldon et al., 2020;Lewis et al., 2023).Moreover, the U.S. Food and Drug Administration enforces strict rules on chiral drugs due to the significant impact of chirality on pharmaceutical properties (FDA, 1992).Different enantiomers can have varying efficacy and safety profiles, making it crucial to ensure the purity and proper identification of the active form.The regulatory environment compels manufacturers to promote enantiopure drugs whenever feasible.The conventional methods used for chiral separation, such as crystallization and chiral chromatography, are challenging and inefficient in terms of both speed and scalability (Rougeot and Hein, 2015;Sui et al., 2023).This further advocates for utilization of enzymes as enantioselective catalysts.
Stereoisomerically pure phosphoramidate nucleoside prodrugs (ProTide) represent an advancement in synthetic intermediates for nucleoside analogues (Figure 1) used to treat viral infection and cancer (Dousson, 2018;Slusarczyk et al., 2018;Sabat et al., 2022).Designing ProTide compounds is to enhance the pharmacological properties of nucleoside analogs, such as improved bioavailability, cell permeability, and resistance to enzymatic degradation within cells (Mehellou et al., 2018).The FDA approved ProTide are such as remdesivir, and sofosbuvir.Many ProTide have a stereogenic phosphorus center thus the therapeutic effect of the ProTide is predominantly attributed to one stereoisomer.For examples, the clinical prodrug form of remdesivir predominantly consists of the (Sp)-diastereomer because of its higher selectivity and broader therapeutic range (Liang et al., 2020).Similarly, (Sp)-diastereomer of sofosbuvir is more potent than the (Rp)-diastereomer.(Sp)-diastereomer of sofosbuvir demonstrates activity 18 times higher than the (Rp)diastereomer against hepatitis C (Sofia et al., 2010).As such, efficient methods for the chemical synthesis and isolation of stereochemically pure ProTides are imperative.
Recently, engineered phosphotriesterase (PTE also known as parathion hydrolase) from Brevundimonas diminuta (formerly known as Pseudomonas diminuta) was reported toward purification of precursors for ProTide (Xiang et al., 2019;Bigley et al., 2020).PTE catalyzes the hydrolysis of p-nitrophenyl phosphotriesters with substituents at the phosphorus center (Figure 2A).Variants of PTE were screened to specifically hydrolyze either (Rp)-or (Sp)-diastereomeric precursor of remdesivir and sofosbuvir, which are hereinafter referred to as (Rp, Sp)-rem and (Rp, Sp)-sof, respectively.Variant G60A showed 165-fold preference for hydrolyzing the (Rp)-sof.While a variant In1W (F132L/H254S/H257W/257-258insSAIGLDPIPN) impressively exhibited 1,400-fold for (Sp)-sof (Xiang et al., 2019).In1W also exhibited >200-fold preference for hydrolyzing (Sp)-rem.On the other hand, G60A did not demonstrate any observable stereoselective hydrolysis of (Rp)rem (Bigley et al., 2020).Stereoselective PTE toward the counterpart of the desired stereoisomer could enhance the production of highly pure stereoisomer precursors for ProTide synthesis.In case of preparative stereoisomer isolation of remdesivir and sofosbuvir, specific hydrolysis of (Rp)diastereomer would provide the pure (Sp)-diastereomer which is the potent isoform (Figure 2B).Although previous studies have demonstrated remarkably success in isolating purified diastereomers, challenges remain in improving hydrolysis of (Rp)-diastereomer of these prodrugs.
In this study, a small numbers of PTE variants were rationally designed and prepared to screen stereoselectivity toward (Rp)diastereomer of remdesivir and sofosbuvir precursors ((Rp)-rem and (Rp)-sof).The mutants with the improved specific activity were subjected to biochemical characterization in detail.We successfully obtained PTE variants that showed improved enantioselectivity toward both (Rp)-rem and (Rp)-sof.These findings underscore the potential of the variants for producing other stereoisomerically pure ProTides and organic compounds.

Structural analysis and molecular docking
Structural analysis and molecular docking were employed to investigate the interactions between PTE and substrates.The complex structure of PTE and paraoxon analog was retrieved from the protein databank (PDB ID: 1DPM) (Vanhooke et al., 1996;Berman et al., 2000).Protein and ligand molecules were prepared with default structural preparation and protonation procedure in MOE (Molecular Operating Environment) program (Chemical Computing Group ULC, 2023).Molecular docking was performed using Triangle Matcher and GBVI/WSA rescoring.Visualization of the docked complexes facilitated the identification of interactions between the protein and ligand molecules and key residues involved in binding.Molecular dynamic (MD) simulation was used to validate the binding interactions between the docked poses and proteins (Supplementary data).

Plasmid construction and PTE variants preparation
The metalloenzyme PTE derived from B. diminuta demonstrates remarkable efficiency in hydrolyzing diverse compounds.However, the investigation of PTE has faced challenges due to the inefficient expression of its recombinant form, despite its significant potential.This study utilizes an evolved variant of PTE, known for its successful heterologous expression, as the standard reference, denoted as wild type (WT), to address this issue (Roodveldt and Tawfik, 2005).The nucleotide sequence was obtained from Genbank (accession number: KU746636.1)and synthesized by GenScript without the leader peptide-coding fragment (99 bp) at N-terminal and tagged with six consecutive histidine residues (6×His) at C-terminal.The target gene was inserted into pMAL.c5xvector at NdeI and HindIII size (pMAL.c5x-WT-PTE).
In1W-PTE (F132L, H254S, H257W, and insertion of SAIGLDPIPN between amino acid residues at 257 and 258) was synthesized by GenScript and transformed into Escherichia coli BL21DE3 cells.Other PTE variants were constructed using the overlap extension PCR.Recombinant plasmid pMAL.c5x-WT-PTEwas used as a PCR template.To introduce a single-site mutation via substitution, two pairs of primers were designed: (i) a forward-NdeI and a reverse primer tailored for the desired mutation, and (ii) a reverse-HindIII primer and a forward primer tailored for the intended mutation (Supplementary Table S1).The obtained PCR products were utilized as the template for the second PCR procedure, employing the forward-NdeI and reverse-HindIII primer pair.The PCR reactions were performed in the following order: initial denaturation for 5 min at 95 °C, followed by 30 cycles of 95 °C for 30 s, 55 °C for 30 s, 72 °C for 1 min, and final extension 72 °C for 7 min.Each PTE mutant was inserted into the pMAL.c5xvector at the NdeI and HindIII site and subsequently transformed into E. coli DH5α.Randomly selected colonies were sequenced to ensure the presence of amino acid substitutions.The resulting recombinant plasmids were then used to transform into E. coli BL21DE3 cells, allowing for protein expression upon isopropyl thio-ß-D-galactoside (IPTG) induction.

Screening for diastereoselective hydrolytic PTE variants
The expression of host cells harboring pMAL.c5xplasmid encoding the mutants PTE were cultured in 50 mL of Luria broth (LB) containing ampicillin 50 μg/mL at 37 °C and 200 rpm.The target protein expression  was induced by 0.2 mM IPTG when the optical density at 600 nm (OD 600 ) reached 1.0-1.3.After that, the culture was incubated at 16 °C for 18-20 h.The bacteria cells were harvested by centrifugation at 8,000 rpm for 10 min and resuspended in 5 mL of 50 mM Tris-HCl buffer (pH 8.0).The crude enzyme was extracted by sonication and removed cell debris by centrifuging at 17,000 rpm for 15 min at 4 °C.The supernatant was then screened for an enantioselective variant.
To screen the enantioselectivity for hydrolysis of sofosbuvir and remdesivir precursor, the samples were performed and investigated by colorimetric assay.The reaction mixture (total volume 1 mL) contained 200 µg of crude enzyme extract, 50 mM CHES (pH 9.0), 0.1 mM CoCl 2 , and 60 µM substrate (two diastereomer of two precursors including (Sp)-sof, (Sp)-rem, (Rp)-sof or (Rp)-rem dissolved in 100% DMSO).The reaction was incubated at room temperature (25 °C) for 60 min, then detected p-nitrophenol at 400 nm.The substrate complete hydrolysis with KOH was defined as 100% hydrolysis and was the control of the experiment.Each analysis was conducted in triplicate.

Purification of recombinant PTE
The expression system for studying the enzyme kinetics of PTE was changed from pMAL.c5x plasmid, which contained MBP, to pET28a.While the fusion of PTE with MBP can enhance protein solubility, the relatively large size of MBP (42 kDa) (Raran-Kurussi and Waugh, 2012) may interfere with kinetic activity of the PTE.PTE has a molecular weight of ~36 kDa (Bigley and Raushel, 2019).To facilitate biochemical characterization of PTE, MBP cleavage was achieved using the protease Factor Xa (Thermo Fisher Scientific Inc., United States) over a period of several days.Given the inherent instability of PTE, mutant PTE genes were subsequently integrated into the pET28a plasmid at the identical restriction site positioned between NdeI and HindIII sites.The constructed plasmids were transformed into E. coli BL21DE3 cells for recombinant PTE expression.The host cells were cultured in 1 L of LB containing kanamycin 50 μg/mL at 37 °C and 200 rpm.The cultured cells reached OD 600 0.4-0.6,then added IPTG 0.1 mM, decreased the temperature to 16 °C, and cultured for 20 h.The cultured cells were harvested and sonicated.The clear supernatant (100 mL) was added imidazole to final concentration of 20 mM and filtered through 0.2 µm before being applied to a HiPrep FF 16/10 column (GE Healthcare) which was also equilibrated by the equilibration buffer (20 mM Tris-HCl, pH 8.0 contained 0.1 M NaCl and 20 mM imidazole).The column was washed with the same buffer at a flow rate of 5 mL/min and the bound protein was eluted with a stepwise gradient imidazole (20-500 mM).The protein pattern of the eluted protein fractions was analyzed by using 12% sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE).All purified PTE fractions were pooled and dialyzed against 50 mM Tris-HCl buffer (pH 8.0) at 4 °C with gentle stirring for desalting.The dialyzed enzyme was concentrated by ultrafiltration and immediately kept the enzyme in 20% glycerol at 4 °C.

Determination of enzyme kinetic parameters
Preliminary assays with varying enzyme concentrations were conducted, and the initial rates were measured to determine the appropriate concentration of enzyme in kinetic studies (Supplementary Table S2).The optimal concentration is selected from the linear range of the initial rate versus enzyme concentration curve, ensuring reliable measurements without saturation.This process balances sensitivity and linearity for accurate kinetic analysis.
The purified PTEs were subjected to characterization of steadystate kinetic constants using Rp-or Sp-diastereomer of sofosbuvir precursor ((Rp, Sp)-sof) and remdesivir precursors ((Rp, Sp)-rem).The reaction mixtures (1 mL) consisted of 50 mM CHES (pH 9.0), 0.1 mM CoCl 2 , and substrate concentrations ranging from 10 μM to 250 µM (dissolved in 100% DMSO).Incubation was carried out at 30 °C for 10 min.Reactions were initiated by adding the appropriate enzyme concentration (Supplementary Table S2) and monitoring the release of p-NP at 400 nm with a spectrophotometer.All experiments were conducted in triplicate, with a control reaction performed in the absence of an enzyme.The Michaelis-Menten constant (k m ) and maximum velocity (V max ) were determined from the Lineweaver-Burk plot using a computer program provided using GraFit Version 7 Software.Thus, the catalytic efficiency of the enzyme (k cat /K m ) was calculated using the equation: where k cat is the turnover number, [E t ] is the total enzyme concentration in moles (with the molecular weight of PTE being 38,128.55g/mol), and V max is maximum velocity or the reaction rate when the enzyme is fully saturated by the substrate.To facilitate the comparison of substrate stereoselectivity, the ratio of the kinetic parameters was calculated as follows.
A K m ratio less than one indicates that Rp substrate exhibits higher affinity for the enzyme compared to Sp substrate.Conversely, a ratio greater than one suggests that Sp substrate has higher affinity.A ration close to one indicates similar affinities for both isoforms.
A k cat ratio or k cat /K m ratio greater than one indicates that the enzyme catalyzes the Rp substrate more efficient than the Sp substrate, and vice versa.
The samples were filtered through 0.2 µm before the determination of diastereomer by HPLC using CHIRALPAK IG-U column.The mobile phase A and B were 95% (v/v) methanol and 100% acetonitrile, respectively.The column was performed at 40 °C before injecting the sample volume of 50 µL.The elution was performed with 100% of mobile phase A for 3 min and then programmed gradient solvent systems with 0%-35% of mobile phase B for 3-10 min.The HPLC peaks were detected with the wavelength at 270 nm.Hydrolytic activity was determined based on the area under the curve (AUC) formed by the decline in enantiomeric substrate over time (0-20 min) using the formula: AUC of isomer R or S ( )from the enzymatic reaction at X min AUC of isomer R or S ( )at 0 min × 100 _

Synthesis of racemic remdesivir and racemic sofosbuvir precursors
The synthesis of (Rp, Sp)-sof was slightly modified from previous reports (Li and Sha, 2008;Cho et al., 2014), whereas the synthesis of (Rp, Sp)-rem was adapted from (Warren et al., 2016) (Supplementary data).Briefly, both racemic compounds were prepared through two-step synthesis.For the racemic sofosbuvir precursor, the first step involved synthesizing L-alanine isopropyl ester hydrochloride, while for racemic remdesivir precursor, L-alanine 2-ethylbutyl ester hydrochloride was used as starting material.These diastereomeric precursors were produced in the second step under controlled conditions of solvent, temperature, stirrer, and reaction time.The resulting product appeared as a lightyellow syrup, comprising a mixture of diastereomeric products (Sp)and (Rp)-in an approximate 1:1 ratio, as confirmed by HPLC and NMR analysis (Supplementary Figures S1-S8).

Separation of Sp-diastereomer
The 0.5 g mixture of diastereomeric precursor for sofosbuvir ((Rp, Sp)-sof) was dissolved in diisopropyl ether and stirred at 5 °C in an ice bath.While stirring, 0.5 mL hexane was added to the solution.The mixture was allowed to stand in a freezer (5 °C) for 12 h.The solid product was collected by filtration, washed with a precooled 1:1 mixture of diisopropyl ether and hexane, and dried under vacuum.The obtained compound was recrystallized again under the same condition.Likewise, the 1 g mixture of diastereomeric precursor for remdesivir ((Rp, Sp)-rem) was dissolved in 4 mL diisopropyl ether.The solution was gently stirred at room temperature for 5 h.The diastereomers thus formed were filtrated off, washed with a precooled diisopropyl ether, and dried under vacuum (Supplementary data).The pure Sp-sof and Sp-rem were confirmed by HPLC and NMR analysis (Supplementary Figures S9-S16).

Computational design of mutations
Enhancing enzyme specificity is to fine-tune enzyme-substrate interactions.In this study, computational structural analysis was employed to modify selectivity of PTE for ProTide synthesis.Synthetic pathways to produce remdesivir and sofosbuvir exist in multiple routes, each involving a complex series of chemical reactions and the use of various solvents (Peifer et al., 2014;Barth et al., 2016;Hu et al., 2022;Kumar Palli et al., 2022).The selected synthetic route used in this study is presented in Figure 2.During the synthesis, racemic mixtures are produced, and subsequent steps are required for isolating and purifying the desired isomer.One of the key chemical intermediates that determine the isomer of these ProTides are 2-ethylbutyl 2-[[(4-nitrophenoxy) (phenoxy) phosphoryl]amino] propanoate ((Rp, Sp)-rem, Figure 2B), and propan-2-yl 2-[[(4nitrophenoxy) (phenoxy)phosphoryl] amino] propanoate ((Rp, Sp)sof, Figure 2B).PTE variants were employed in the precursor isolation step to specifically hydrolyze (Rp)-isomer therefore the desired (Sp)isomer will remain.
PTE is promiscuous on a broad spectrum of diverse substrates.The substrate promiscuity makes it attractive for various applications.PTE enzymes are renowned for their capability to detoxicate organophosphate agents such as paraoxon and parathion (Briseño-Roa et al., 2011).Several studies have shown that it possible to alter PTE specificity for organophosphorus insecticide (Hong and Raushel, 1999;Chen-Goodspeed et al., 2001b;Naqvi et al., 2014;Kronenberg et al., 2024).The PTE have three binding pockets (small, large, and leaving group, Figure 3) (Chen-Goodspeed et al., 2001a;Bigley and Raushel, 2013).The small pocket is composed of G60, I106, L303, and S308.The large pocket and the leaving pocket are highly exposed to solvent.The leaving pocket is composed of hydrophobic side chain of W131, F132, F306 and Y309 whereas the large pocket is defined by H254, H257, L271 and M317 (Figure 3).The hydrolysis of PTE requires two divalent metal ions, commonly zinc ions, for enzyme activity (Hong and Raushel, 1996).These metal ions activate a water molecule, converting it into a hydroxide ion, which then attacks the phosphorus atom of the substrate, leading to the cleavage of the ester bond (Aubert et al., 2004;Bigley and Raushel, 2013;2019).A product such as p-nitrophenol is released (Figure 2A).
Instead of investigating all possible mutations of amino acids aligned at the binding site, structure-based analysis and design was employed to improve specificity.Molecular docking of each stereoisomer of remdesivir and sofosbuvir precursors were used to specify possible binding poses.The docked poses were selected based on the lowest docking scores with preferred interactions where phosphoryl oxygen should interact with metal ions (Aubert et al., 2004;Bigley and Raushel, 2013;2019).MD simulations were performed to validate the docked poses.Binding free energies were also calculated to evaluate whether the mutation favored Rp-rem or Rp-sof or both (Supplementary data, Supplementary Figures S25, 26; Supplementary Table S3).
To improve selectivity toward (Rp)-isomer, the mutations were designed either to block the binding of (Sp)-isomer or to improve the accommodation of (Rp)-isomer.A small set of PTE variants were generated.Enlarging the size of the small pocket was carried out by substituent I106 with alanine (A) and valine (V) which have smaller side chain but possess hydrophobic side chain akin to isoleucine (I).Similarly, replacement of L303 or S308 with alanine (A) was conducted.To enhance potential pi-interactions, phenylalanine (F) and tyrosine (Y) having benzene ring were introduced at L303 and S308, respectively.W131, an amino acid at leaving group pocket, was substituted with methionine (M) having a smaller hydrophobic side chain.For reshaping the large pocket, substitution of L271 with glutamic acid (E), asparagine (N), or phenylalanine (F) and substitution of H254 with tyrosine (Y) were employed.Moreover, we explored replacing a non-catalytic residue D233 as it was in proximity to the predicted binding of (Rp)diastereomer.Although the substitution of D233 was not aimed at enhancing activity toward (Rp)-diastereomer, it was included to assess the efficacy of rational design.In total, we generated 22 PTE variants in addition to WT-PTE, G60A and In1W.These variants were then employed in screening for stereospecificity.

Screening for diastereoselective variants
Crude enzymes were used in initial screening.Selective hydrolysis was assessed through the relative hydrolysis of individual pure diastereomers of sofosbuvir precursor and remdesivir precursor (Figure 4).The results of In1W and G60A variants agreed with previous report.In1W exhibits a strong preference for (Sp)-rem and (Sp)-sof, while G60A was capable of hydrolyzing only (Rp)-sof and not (Rp)-rem (Xiang et al., 2019;Bigley et al., 2020).Nevertheless, crude WT enzyme showed higher activity for (Rp)-sof (79.6%) than the G60A variant.Among mutations in the small pocket, the substitution of I106A resulted in an enhancement of hydrolytic activity for (Sp)-rem, (Sp)-sof, and (Rp)-sof compared to either WT or G60A variant.I106V variant displayed superior hydrolysis of the (Rp)-sof compared to G60A however, it exhibited diminished activity against (Rp)-rem.No substitutions of amino acids in the large pocket resulted in hydrolysis improvement of (Rp)-isomer.Interestingly, W131M mutation in the leaving group pocket, demonstrated nearly absolute 100% relative hydrolysis of both (Rp)-sof and (Rp)-rem.The double-mutant G60A/W131M yielded a variant with a higher catalytic efficiency toward (Rp)-isomer over (Sp)-isomer for both remdesivir and sofosbuvir precursors.Conversely, the two-point mutation I106A/W131M led to non-selective hydrolysis.The variants exhibiting (Rp)-selective hydrolysis including G60A, I106A, W131M, G60A/I106A, G60A/W131M, and I106A/ W131M were chosen for purification and detailed analysis.

Enzymatic activity and kinetic characterization of the PTE variants
Kinetic parameters (K m , k cat and k cat /K m ) were used to characterize enzyme variants (Table 1).Michaelis constant (K m ) represents the substrate concentration needed for an enzyme to achieve half of its maximum velocity (Vmax).A high Km value indicates that a large amount of substrate is necessary to saturate the enzyme, implying that the enzyme has low affinity for the substrate.Catalytic activity is indicated by turnover number (k cat ).High k cat value suggests a rapid reaction rate.The ratio of (k cat /K m ) represents the catalytic efficiency.Furthermore, the ratio of the kinetic parameters of (Rp)-to (Sp)-(abbreviated Rp/Sp) was calculated for comparing substrate stereoselectivity.If the mutant enzyme shows a preference for the Rp substrate over the Sp substrate, the K m ratio (Rp/Sp) will be less than 1.Conversely, if Rp substrate is more efficient, the ratio of k cat or k cat /K m for Rp substrate compared to that for the Sp substrate will be greater than 1.
The results (Table 1) showed that both affinity and reaction rate were imperative to considered for optimizing the effectiveness of the enzyme in catalyzing enantioselective process.For example, the lower of K m value as well as the K m ratio of (Rp/Sp) in WT variants indicated the higher affinity toward the (Rp)diastereomers of both ProTide precursors compared to theG60A variant.Despite the better affinity, WT enzyme was unable to separate the racemic mixture during the reaction while G60A variant showed that capability (Xiang et al., 2019).This could be explained by the significantly change in the reaction rate.Mutation of G60A resulted in enhancement in the turnover rate specially toward (Rp)-diastereomers.As a results, the overall catalytic efficiency of G60A toward (Rp)-diastereomers was better than its effectiveness toward (Sp)-diastereomers.Most PTE variants except G60A exhibited improved affinity toward (Rp)-diastereomer as indicated by the lower K m ratio of (Rp/Sp).However, only the I106A/W131M showed a preference for both (Rp)-rem (K m ratio = 0.70) and (Rp)-sof (K m ratio = 0.11).Remarkably elevated turnover numbers were detected in variant W131M and I106A/W131M (Table 1).In W131M variant, the turnover number ratio toward (Rp/Sp)-rem (291.04) was highest compared to the relative ratio of other variants.Interestingly, while the introduction of a single point mutation at I106A showed marginal enhancement of turnover rate compared to W131M, the synergistic effect of combining I106A and W131M led to significant improvement.The catalytic efficiency ratio of (Rp/Sp)-rem, increased by 333.79-fold, while in (Rp/Sp)-sof, it surged to 995.58 (Table 1).The results suggested the W131M and I106A/W131M variants have potential for utilization in stereoselective hydrolysis.
As hydrolysis time is also an important factor for industrial use of a biocatalyst, we further investigated time course of hydrolysis for the preparation of pure diastereomer using W131M and I106A/W131M in comparison to WT (Figure 5).In WT, the kinetic properties (Table 1) showed a greater catalytic efficiency for sofosbuvir precursors compared to remdesivir precursors and a higher preference for (Rp)-diastereomer than (Sp)-diastereomer.Nevertheless, WT could not completely hydrolyze (Rp)-sof within 20 min of hydrolysis.At 20 min, the remaining (Rp)sof was 3%, while the remaining (Sp)-sof was 23% (Figure 5A).On the other hand, W131M showed complete hydrolysis of (Rp)-sof and (Rp)rem within 8 min and 16 min, respectively (Figure 5B).At the same time, it retained a high amount of the desired isomer, with 81% of (Sp)-sof at 8 min and 77% of (Rp)-rem at 16 min.
For I106A/W131M variant, the enzyme could completely hydrolyze the (Rp)-sof within 12 min while the remaining of the pure (Sp)-sof was over 90% (Figure 5C).However, I106A/W131M variant could not completely hydrolyze (Rp)-rem within 20 min.At 20 min, hydrolysis of I106A/W131M left 2% of (Rp)-rem remaining while 79% of (Sp)-rem was retained.The rapid hydrolysis achieved in a short period (<20 min) effectively depleted unwanted diastereomers, justifying the expediency of preparing pure diastereomers.These results imply that the enzymes developed in this study could facilitate the development of cost-effective bioprocesses for preparing ProTide drugs in a batch system.The mutant PTEs screening for determination of stereoselectivity.PTE exhibits proficiency in hydrolyzing diverse substrates, notable variations arise in their catalytic efficiencies, often attributed to differences in the size of the binding pocket (Chen-Goodspeed et al., 2001b;Bigley and Raushel, 2013;2019).Previous studies on the synthesis of sofosbuvir and remdesivir focused on selective degradation of Sp-isomers (Xiang et al., 2019;Bigley et al.,    Time course of enzymatic hydrolysis of the ProTide precursors with different PTE variants including (A) wild type PTE, (B) W131M, and (C) I106A/ W131M.The reactions were carried out as follows.The Enzymes were incubated at 30 °C with 60 µM of each diastereomeric substrate in a reaction mixture at optimized concentrations determined from prior experimentation.WT-PTE was used at concentrations of 7.99 µM for (Rp, Sp)-sof and 15.99 µM for (Rp, Sp)-rem, while mutant enzymes (W131M and I106A/W131M) were used at 35.57 nM and 0.72 µM for both precursors, respectively.Aliquots were taken every 2 min, quenched with methanol, and analyzed by HPLC using a CHIRALPAK IG-U column with gradient elution.Detection was performed at 270 nm, and the remaining substrates percentage were determined from the area under the curve.
The kinetic analysis showed an aspect that should be considered for future development.Although the affinity of the variants were improved as designing, the kinetic data suggested that flexibility had a major role in enzyme selectivity of PTE.High affinity for the ligand may resulted in decreasing of turnover speed (Kari et al., 2021) which could consequently resulted in reduction of catalytic efficiency and substrate selectivity.Moving forward, further studies leveraging advanced computational modeling and turnover number prediction hold the promise of unlocking new strategies to fine-tune biocatalysts for enhanced selectivity and efficiency in diverse applications.By navigating this intricate balance, biocatalyst designers can pave the way for transformative advancements in biotechnology and pharmaceutical synthesis.

FIGURE 1 FDA
FIGURE 1FDA-approved and clinical candidate phosphoramidate nucleotide prodrugs.

FIGURE 3
FIGURE 3 Three-dimensional structure of phosphotriesterase (PTE).(A) Interaction of PTE with paraoxon analog, diethyl 4-methylbenzylphosphonate (PDB ID: 1DPM).Receptor and ligand exposures as well as the physiochemical nature of the residues are explained in the color legend.(B) Distribution of molecular hydrophobicity at the active site.Ribbon diagrams and molecular surface colored according to hydrophobicity potential are shown.Color scheme: red, hydrophobic; blue, hydrophilic.The structure of ligand is presented by purple balls and sticks model.(C) Binding site of PTE.The amino acids at small pocket are presented in orange, large pocket in green and leaving pocket in pink.

TABLE 1
Kinetic parameters of purified wild type PTE (WT) and their variants.All values are means of three independent determinations.
(Continued on following page)Frontiers in Bioengineering and Biotechnology frontiersin.org

TABLE 1 (
Continued) Kinetic parameters of purified wild type PTE (WT) and their variants.All values are means of three independent determinations.