Open Access Article
Liangyu Qian†
a,
Isaiah T. Dishner†
b,
Dana L. Carper
a,
Jerry M. Parks
a,
Vilmos Kertesz
a,
Nicholas T. Zolnierczuk
a,
Nduka D. Ogbonna
b,
Nikki A. Thiele
b,
John F. Cahill
a,
Jeffrey C. Foster
*b and
Joshua K. Michener
*a
aBiosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN, USA. E-mail: michenerjk@ornl.gov
bChemical Sciences Division, Oak Ridge National Laboratory, Oak Ridge, TN, USA. E-mail: fosterjc@ornl.gov
First published on 23rd June 2026
Nylons, a major class of synthetic polyamides, are widely used due to their excellent mechanical strength, thermal stability, and chemical resistance. Conventional nylon production relies on the polymerization of lactams or stoichiometric nylon salts. However, applying these approaches to unconventional precursors such as bio-derived glutaric acid produces polymers with low molecular weight and limited applications. To address these challenges, we demonstrated that chemically-synthesized nylon diads enable the production of higher-molecular-weight polyamides compared with traditional salts. We then identified a biosynthetic approach using amide synthetases to convert unprotected bifunctional substrates into nylon-relevant diads. Using a cofactor regeneration system, enzymatic diad synthesis was scaled to produce sufficient material for laboratory-scale characterization and solid-state polymerization. Amide synthetases demonstrated broad substrate scope, catalyzing the regioselective assembly of diverse nylon-relevant diacids, diamines, and ω-amino acids. This strategy offers a novel route to synthesize challenging nylon monomers and advances production of bioderived nylons.
Green foundation1. This work describes the biosynthesis of process-advantaged amide diads for improved polymerization of nylons from renewable feedstocks2. By increasing the molecular weights of bioderived nylons by more than 5-fold, amide diads overcome a key obstacle in developing sustainable polymers 3. In the future, amide synthetases can be applied in vivo using engineered bacteria for direct conversion of renewable feedstocks into amide diads |
Commercial nylons are typically produced by polymerizing either a single ω-amino acid (e.g., nylon-6, also known as PA6) or a stoichiometric nylon salt formed from a diamine and a diacid (e.g., nylon-66, or PA66). Efforts to improve nylon sustainability have led to interest in incorporating bioderived components such as glutarate and cadaverine.5–13 However, polymerization of nylon salts containing glutarate often yields polymers with low molecular weight, which has hindered their broader exploration and application.5,14 Despite the limited studies on glutarate-based nylons, PA55 composed of cadaverine and glutarate has been reported to possess desirable properties, including a high melting temperature and excellent thermal stability.14 PA65 remains largely unexplored but shows a distinct hydrogen bonding pattern and moderate crystallinity.5,15 More flexible synthetic routes that can readily incorporate non-conventional monomers could expand the scope of potential polyamides.
As an alternative synthesis route, PA66 diads generated during enzymatic PA66 hydrolysis have been polymerized to yield PA66 with higher molecular weight than that generated from polymerization of the corresponding nylon salt.16 However, it was unclear whether other nylon diads would similarly outperform their respective salts. Furthermore, while oligomeric products can be formed from polyamide hydrolysis, direct diad synthesis is challenging. Peptide synthesis via solid-phase methods is well-established, but it involves expensive protection and deprotection steps to control chain length, sequence and reactivity, as well as extensive use of organic solvents.17 Alternatively, enzymatic synthesis could potentially offer a more efficient and environmentally friendly route for generating these intermediates. A biosynthetic approach would be particularly useful if combined with a pathway for biosynthesis of nylon-relevant monomers, since production of zwitterionic diads would avoid the cost of separately synthesizing, neutralizing, and purifying diacids and diamines.18
Biology provides several possible routes for amide bond formation. Ribosomes can synthesize amide bonds between natural α-amino acids. However, while ribosomes can generate polyamides with high molecular weight and exquisite sequence control, expanding the monomer scope beyond canonical α-amino acids is challenging.19–25 Additionally, non-ribosomal peptide synthetases (NRPSs) are multimodular enzymes that offer greater flexibility in monomer scope and bond formation but are difficult to express and engineer.26–28 In addition to ribosomal and non-ribosomal peptide biosynthesis, recent work has demonstrated that enzyme promiscuity can also enable amide bond formation through distinct enzymatic mechanisms. Certain alcohol dehydrogenases have been shown to catalyze oxidative coupling of alcohols with ammonia or amines via hemiaminal intermediates, providing a redox-driven route for synthesis of primary and secondary amides.29 Additionally, amide synthetases can catalyze amide bond formation in small molecules by using ATP to adenylate a carboxylic acid, which is then condensed with an amine.30,31 These enzymes natively function in the biosynthesis of complex secondary metabolites such as antibiotics and siderophores and have been mainly repurposed for the production of pharmaceuticals.32–36 They are relatively small (400–500 amino acids), easy to express and purify, and active in heterologous hosts.
Enzymatic amide bond formation has emerged as a sustainable alternative to traditional synthetic approaches. For example, the amide synthetase McbA from the ANL superfamily of adenylating enzymes is known for its broad substrate specificity and has been applied to the chemoenzymatic synthesis of β-carboline derivatives and other pharmaceutical-type amides.33,37 However, McbA has been shown to have low activity with aliphatic substrates, suggesting that it may not accept nylon-relevant monomers.32 Similarly, the amide synthetase SfaB, a self-sufficient AMP ligase from the ANL superfamily that participates in an NRPS pathway, accepts a broad range of short-chain fatty acids and catalyzes amidation or thioesterification with various amine or thiol nucleophiles, but reportedly does not tolerate dicarboxylic acids such as those used in nylons.38 Therefore, the application of natural or engineered amide synthetases for the synthesis of nylon building blocks remains largely unexplored.
In this study, we showed that chemically-synthesized nylon diads can be used to generate a range of polyamides with higher molecular weights compared to polymerization of traditional salts. Motivated by these findings, we then demonstrated a simplified biosynthetic strategy using amide synthetases to produce nylon diad precursors from unprotected bifunctional substrates. By coupling amide bond formation with an ATP-regeneration system, we optimized and scaled up this biocatalytic process and established it as a viable method for synthesizing nylon diads. Furthermore, we showed that amide synthetases exhibit a broad substrate scope and can catalyze the regioselective assembly of a diverse range of nylon-relevant diacids, diamines, and ω-amino acids. This approach provides a facile route for the polymerization of challenging monomers and expands the potential for tailoring nylon properties through enzymatically derived precursors.
To test these hypotheses, we used isothermal thermogravimetric analysis (TGA) to monitor the change in sample mass under simulated SSP conditions. Diad and nylon salt samples were initially conditioned at 100 °C for 1 h to remove any residual solvent and/or moisture, and subsequently polymerized by heating to 220 °C for 8 h (Fig. S40–S45). As shown in Fig. 1D, most nylon salts lose a considerably larger portion of their initial mass than their diad counterparts over the course of the polymerization, even when accounting for differences in expected mass loss caused by the evolution of water. Moreover, the majority of this mass loss occurs at relatively low monomer conversion (i.e. within approximately the first 30 min). These results imply that evaporation of diamine and/or diacid components of the nylon salts creates a stoichiometric imbalance which limits maximum Xn in accordance with Carother's equation.40 Exceptionally high mass losses were observed during the polymerization of succinic acid-containing PA64 and PA54 samples regardless of whether the diad or nylon salt starting materials were used (Fig. S44 and S45). This mass loss was ultimately attributed to competing imidization pathways which can create volatile byproducts when imidization occurs adjacent to a polymer chain end. Imide formation was clearly visible from Fourier-transform infrared (FT-IR) spectroscopy of the PA64 and PA54 polymers, which showed an imide C
O stretch at ∼1770 cm−1 (Fig. S46). Thermal imidization also explains the relatively low Mw of PA64 and PA54 samples since imidization imbalances the ratio of carboxylic acid to amine chain ends and also acts as a form of ring-chain equilibrium which further limits Xn.41
To compare the physical properties of the polyamides prepared from diads with those prepared from the corresponding nylon salt, we used thermogravimetric analysis (TGA) and differential scanning calorimetry (DSC) to measure the degradation temperatures (Td), glass transition temperatures (Tg), melting temperatures (Tm), and enthalpies of melting (ΔHm) (Fig. S47–S52 and Table S2). No considerable difference was observed in Td for PA66, PA56, PA65, PA55, or PA54 prepared from either nylon salts or diads. However, a somewhat lower Td was observed for PA64 synthesized from MS compared to its corresponding salt (Fig. S51A). Additionally, PA64 prepared from MS showed an apparent two step decomposition profile, whereas the salt-derived polyamide decomposed in a single step. Despite this difference in thermal stability, both PA64 samples display nearly identical thermal properties, consistent with the hypothesis that the two samples are chemically identical. Likewise, similar agreement in Tg, Tm, and ΔHm for salt-derived and diad-derived polyamides were observed for PA66, PA65, and PA54. In combination, we found that diads are preferred feedstocks for nylon polymerization, since they yield polyamides that are chemically identical but have higher molecular weights compared to the corresponding salt-derived polyamides.
Encouragingly, DdaG and SfaB showed significant activities and distinct specificities for nylon diad synthesis. The molar yields after overnight incubation ranged from 3% to 80%, with DdaG demonstrating the highest activity in synthesizing MS (∼80%) and CS (∼60%) but low production of MG and MA. In contrast, SfaB exhibited lower activity than DdaG for production of MS and CS but produced higher amounts of MG (∼20%) and MA (∼3%). No detectable formation of triads, such as diacid triads (e.g., SMS and AMA) and diamine triads (e.g., MSM and MAM) or higher-order oligoamides were observed. This result indicated that those amide synthetases actively catalyze diad formation but do not extend to higher-order oligoamides. We next scaled up MA production with SfaB. As an additional safeguard, we sought to confirm the structure of enzymatically-produced MA. To this end, we developed a purification protocol wherein the crude reaction mixture was reacted with 2-acetyldimedone (dde-OH) to provide a dde-functionalized adduct that is easily isolated chromatographically, and subsequently deprotected via hydrazinolysis (Scheme S1 and Fig. S57 and S58). This allowed us to confirm the structure of enzymatically synthesized MA by comparing the 1H NMR spectra to the chemically synthesized standard (Fig. 2D). Scale-up of this process would require the development of a more efficient purification route.
We did not observe detectable product formation in reactions with SuCphA1 or McbA, which is consistent with findings of a previous report that McbA does not accept aliphatic substrates.32 Interestingly, previous studies showed that SfaB cannot adenylate diacids for subsequent condensation with hydroxylamine, based on colorimetric detection of pyrophosphate released during adenylation.38 However, our observations reveal that SfaB can act on a range of diacids (i.e., succinic acid, glutaric acid and adipic acid) in combination with cadaverine and hexamethylenediamine. These differing results suggest that SfaB also exhibits specificity towards the amine acceptor.
Since the reactions did not fully convert the substrates to products, we next investigated factors that might limit yield. To assess the potential for product inhibition, we conducted reactions in which we incubated DdaG with 5 mM succinate and 5 mM hexamethylenediamine, with or without the addition of 2.5 mM of the purified MS diad. The reactions produced similar increases in MS concentration (Fig. S95A), indicating that product inhibition was not significantly limiting the yield. Next, we tested for diad hydrolysis by incubating purified DdaG with the MS diad. Diad concentrations remained unchanged relative to the no-enzyme control (Fig. S95B), suggesting negligible hydrolysis. Finally, we incubated reactions for 24 h, measured the MS concentration, and then added an additional 5 μM of DdaG and let the reaction continue for an additional 24 h. The addition of fresh enzyme approximately doubled the final MS concentration (Fig. S95C), indicating the enzyme deactivation likely limits conversion. Together, these results indicate that the limited yield is unlikely to arise from strong product inhibition or hydrolytic equilibrium and, instead, is associated with additional factors such as enzyme deactivation that could be addressed through enzyme engineering.
DdaG and SfaB are ATP-dependent amide synthetases, and providing super-stoichiometric ATP to ensure effective conversion is cost-prohibitive. As an alternative, we implemented an ATP recycling system using a purified type 2-III polyphosphate kinase (PPK12).34,45,46 In addition, we included inorganic pyrophosphatase from E. coli (ecPPase) to prevent potential inhibition by pyrophosphate. We then performed reaction optimization using MS production by DdaG. We iteratively optimized the concentrations of the substrates, polyphosphate, magnesium, AMP, and pH (Fig. S61–65). Based on these results, the finalized optimal conditions for MS production were 50 mM succinic acid, 100 mM hexamethylenediamine, 100 mM MgCl2, 10 μM DdaG, 1 mg mL−1 PPK12, 5 mM AMP, 0.2 U ecPPase and 100 mM HEPES (pH 8). With these optimal conditions, we subsequently scaled up our reaction system with ATP regeneration from 50 μL to a 5 mL reaction volume. After 24 h incubation, a small volume (50 μL) of sample was analyzed by I.DOT/OPSI-MS, indicating a promising yield of ∼50% (∼5 g L−1) and demonstrating a more than 3-fold improvement from the initial conditions (Fig. 3). Although these experiments were performed with a single reaction, they demonstrate proof of principle that enzymatic catalysis using amide synthetases for nylon diad synthesis is both scalable and amenable to optimization. Future work should focus on further scale-up and process optimization to improve productivity, robustness, and applicability to specific reactions. In addition, we calculated the E-factor, a key metric in green chemistry, to evaluate and compare the environmental impact of enzymatic and chemical diad synthesis. The ATP-regenerated enzymatic approach exhibits a substantially lower E-factor (i.e., 4.3–8.7; excluding or including HEPES buffer, respectively) than the conventional chemical synthesis route (i.e., 63), highlighting the potential of biocatalysis as a greener alternative.
In vitro synthesis is well-suited for rapid lab-scale oligomer production to enable polymer synthesis and characterization. Longer term, in vivo strategies offer significant advantages for scale-up. Through metabolic engineering, intracellular ATP availability can be modulated to support ATP-dependent biosynthetic reactions without the need for external supplementation.47 Moreover, the microbial biosynthesis of nylon-relevant monomers has been achieved. For example, engineered Escherichia coli strains have been developed for the production of adipic acid at high titers.48 Similarly, de novo microbial synthesis of hexamethylenediamine and related diamines from glucose has been demonstrated.49 Collectively, these advances highlight the potential for future work of developing in vivo platforms for the synthesis of nylon-relevant precursors, in which heterologous expression of amide synthetases could enable intracellular diad formation from endogenously produced substrates.
We initially detected product formation using direct injection IDOT/OPSI-MS for untargeted characterization of enzymatic reactions.50 Using this approach, 41 out of 96 reaction combinations generated detectable products above the control using either DdaG or SfaB (Fig. 4 and S66). To more accurately identify these products, we confirmed 24 structures through analysis of fragmentation patterns by MS2 (Fig. S67–S90). Both enzymes displayed activities in ligating the carboxylic groups of the diacids with a wide spectrum of amine acceptors including linear diamines (cadaverine C, hexamethylenediamine M), an aromatic diamine (p-xylylenediamine X), a cyclic diamine (cis-1,4-cyclohexanediamine N), and the amine groups of ω-amino acids (6-aminohexanoic acid H).
![]() | ||
| Fig. 4 Amide synthetases have broad substrate range for potential nylon precursors. In vitro biochemical assays were conducted by incubating different diacids (i.e., succinic acid, glutaric acid, adipic acid and thiomalic acid) with diamines or ω-amino acids in the presence of enzymes (i.e., DdaG or SfaB) or no-enzyme control. Products were assayed using OPSI-MS. The log10MS signal intensity was cut off at 5.0 to filter out noise. C: cadaverine; M: hexamethylenediamine; X: p-xylylenediamine; N: cis-1,4-cyclohexanediamine; B: 4-aminobutyrate; V: 5-aminovalerate; H: 6-aminohexanoic acid; E: 4-amino-2-hydroxybutanoic acid; F: 4-amino-3-hydroxybutanoic acid. Sample size is N = 3. Data for CS, MS, CG, MG, CA, and MA shown in Fig. 2B are replotted here as unquantified MS intensities for comparison. | ||
Notable differences between DdaG and SfaB were observed in their preference for carboxylic acid donors, consistent with our quantification of their activity in producing more conventional nylon-relevant diads. DdaG exhibited excellent catalytic performance in ligating succinic acid (S) with a wide range of amine acceptors but showed lower activities with glutaric acid (G) and adipic acid (A) compared to SfaB. In contrast, SfaB was better using G and A as carboxylic donors with diverse amine acceptors. In addition, DdaG can also catalyze reactions involving terephthalic acid (T) and α-ketoglutaric acid (K) with multiple diamines, while SfaB was inactive in these reactions. Alternatively, SfaB was capable of activating 1,4-cyclohexanedicarboxylic acid (O) with specific diamines (i.e., M and X), whereas DdaG showed negligible activity (Fig. S66). It is also notable that both enzymes were unable to activate the carboxylic groups of ω-amino acids. Interestingly, we noticed both enzymes could activate diacids bearing additional functional groups such as thiomalic acid (P). The broad substrate scope of DdaG and SfaB suggested the potential of using amide synthetases to produce novel value-added oligoamides.
To obtain mechanistic insights into whether substrate binding interactions could contribute to the distinct diacid preferences of DdaG and SfaB, we used Boltz-251,52 to model representative enzyme–substrate complexes (i.e., DdaG-succinate/glutarate and SfaB-succinate/glutarate) and analyzed predicted affinity values across 100 models. For each enzyme–substrate pair, all 100 models had high Boltz confidence values (DdaG-glutarate > 0.89, DdaG-succinate > 0.89, SfaB-succinate > 0.91; SfaB-glutarate > 0.93). The predicted affinity trends were consistent with the experimental preferences. For SfaB, glutarate was favored over succinate in 82 of 100 paired models, with median predicted IC50 values of 626 and 851 μM, respectively. For DdaG, succinate was favored over glutarate in 99 of 100 paired models, with median predicted IC50 values of 265 and 519 μM, respectively (Table S3). We then used protein–ligand interaction profiler (PLIP)53 to analyze the highest-affinity model for each enzyme–substrate pair to identify favorable noncovalent interactions (Fig. S96). In both enzymes, the two substrates showed similar hydrophobic and salt-bridge interactions, whereas the preferred substrates showed additional hydrogen-bonding interactions (Tables S4 and S5). In DdaG, the succinate-bound model contained two more hydrogen bonds than the glutarate-bound model. In SfaB, the glutarate-bound model contained two hydrogen bonds that were absent in the succinate-bound model. Together, the predicted affinity values from Boltz-2 and subsequent PLIP analysis suggest that differences in substrate binding interactions may contribute to the observed diacid preferences of DdaG and SfaB. Future crystal structures of specific enzyme–substrate complexes would provide more definitive molecular insights for understanding the substrate specificities.
We also determined that DdaG regioselectively ligated S with the asymmetric polyamide spermidine. Of the two possible primary amines in spermidine, DdaG selectively coupled S with the primary amine nearer the middle amine, producing only a single diad product. This regioselectivity highlights the potential of amide synthetases to synthesize ordered polymers from asymmetric substrates, which is challenging for traditional chemical synthesis methods (Fig. S91).54
Given the distinct activities of these enzymes with diacids and their considerable tolerance for various diamines, we investigated the effect of diacid and diamine chain length on enzyme activity. We tested linear diacids and diamines with carbon lengths of 7 to 10, excluding C10 diamines due to solubility issues (Fig. S92). The results for DdaG were consistent with our previous findings, showing a preference for S as the diacid and comparable activity across the full range of diamine lengths. In contrast, SfaB, although less active in ligating S with longer diamines, showed activity with C7–C10 diacids and M.
To further characterize DdaG performance, we conducted kinetic studies using succinic acid and hexamethylenediamine as substrates. The kcat and Km values for DdaG with these substrates were determined to be 3.2 ± 0.2 min−1 and 8.8 ± 1.8 mM (Table 1 and Fig. S93), respectively. Compared to its native reaction with fumaric acid and 2,3-diaminopropionic acid,36 the kcat with non-native substrates is 4-fold lower and Km is 15-fold higher.
| Substrates | kcat (min−1) | Km (mM) | Note |
|---|---|---|---|
| Succinic acid and hexamethylenediamine | 3.2 ± 0.2 | 8.8 ± 1.8 | This study |
| Fumaric acid and 2,3-diaminopropionic acid | 11.6 ± 0.4 | 0.55 ± 0.07 | Hollenhorst et al.36 |
Finally, we tested amide synthetase activity using cell-free enzyme expression. When expressed in a cell-free system, DdaG showed high activity ligating its native substrates but low production of the PA64 diad (Fig. S94). Similarly, SfaB exhibited higher activity with a short chain fatty acid (i.e., 5-chlorovalerate) than glutarate. These results suggest that amide synthetases have substantial engineering potential to enhance their activity with non-native or non-preferred substrates, and a cell-free system could serve as an efficient platform for high-throughput enzyme engineering.32
His-PPK12 was produced as described previously.45 Expression and purification of His-DdaG, His-SfaB, His-McbA, His-PylC, His-SUMO-PylC and His-SuCphA1 were conducted following previously published methods. A colony or glycerol stock of BL21(DE3) E. coli containing plasmid DNA was used to inoculate 10 mL of LB medium supplemented with appropriate antibiotics. The culture was incubated at 37 °C and 250 rpm overnight. The overnight culture was then diluted 100-fold into 200 mL of Terrific Broth. The subculture was incubated at 37 °C at 250 rpm until an OD of 0.4–0.6 was reached. Protein expression was induced by adding IPTG to a final concentration of 0.1 mM. The temperature was decreased to 16 °C and the culture was grown with shaking overnight. The cells were harvested by centrifugation (4000 rpm, 40 min, 4 °C), resuspended in 5 mL of lysis buffer (50 mM HEPES, 300 mM NaCl, 10 mM imidazole, 10 mM MgCl2, 10% glycerol, pH 8.0), and sonicated on ice. Cell debris was removed by centrifugation (10
000 rpm, 30 min, 4 °C). The supernatant was then filtered with an MCE filter (0.22 μm) and loaded at 3 mL min−1 onto an ÄKTA Start FPLC (Cytiva Marlborough, MA) equipped with a 5 mL His-Trap™ column (Cytiva). The column was washed with wash buffer (50 mM HEPES, 300 mM NaCl, 30 mM imidazole, 10 mM MgCl2, 10% glycerol, pH 8.0) before eluting with elution buffer (50 mM HEPES, 300 mM NaCl, 250 mM imidazole, 10 mM MgCl2, 10% glycerol, pH 8.0) at 3 mL min−1. Purified proteins were concentrated with an Amicon® Ultra centrifugal filter at 10 kDa cutoff, and buffer exchanged into exchange buffer (50 mM HEPES, pH 8.0, 100 mM NaCl, 10 mM MgCl2, 10% glycerol). The final purified product was analyzed by SDS-PAGE and enzyme concentration was measured using a NanoDrop™ 1000 Spectrophotometer (Thermo Scientific) with exchange buffer as a blank, and the extinction coefficient of each protein was calculated using ProtParam from the ExPASy Proteomics Server. The molecular weight of each protein was calculated using ProtParam from the ExPASy Proteomics Server, and further confirmed by the SDS-PAGE method. Pure enzymes were stored at −80 °C for subsequent activity assays.
:
100 (v/v) in high performance liquid chromatography (HPLC) grade water with 500 nM propranolol acting as an internal standard for droplet capture. 40 µL of diluted reactions were transferred to a I.DOT S.100 96-well plate and analyzed without further separations. The I.DOT system was used to eject 20 nL of sample into the OPSI, into a flow of 75/25/0.1 (v/v/v) acetonitrile/water/formic acid and transported the sample to the electrospray ion source of the mass spectrometer. Dispensing throughput was 4 s per sample and each sample was measured in triplicate. Peak widths were ∼1.2 s wide.
For identification of products by high mass resolution, a Thermo Q-Exactive HF mass spectrometer (ThermoFisher Scientific) in positive ion mode was used for characterization. Scan settings were sheath gas = 80, auxiliary gas = 40, electrospray voltage = 4 kV, ion injection time = 50 ms, automatic gain control = 3e6, capillary temperature = 200 °C, mass/charge (m/z) range = 100–750 m/z, and OPSI solvent flow of 250 µL min−1. A 60
000 mass resolution was used for quantitative analyses while 240
000 resolution was used for compound identification. Tandem MS data were acquired using the Q-Exactive HF, but water adduct formation was commonly observed in the ion trap. To avoid water adducts in tandem mass spectra, tandem mass spectrometry data was collected using a Sciex 5600 TripleTOF time of flight mass spectrometer (Sciex) with the following settings: GS1 = 75, GS2 = 35, electrospray voltage = 5.5 kV, capillary temperature = 200 °C, and OPSI solvent flow of 150 µL min−1. Collision energies across 10–50 eV were evaluated for each ion.
Control of I.DOT settings (positioning, well selection, timing and droplet dispensing cycles), OPSI peak finding, data extraction from vendor file formats, and signal integration were enabled through custom software packages written in Delphi, Python, and C#, available upon request. For each droplet sampling event, the average mass spectrum from the resulting mass spectral peak was normalized to propranolol internal standard signal and background subtracted. Log–log matrix matched calibration curves were used for quantitation. Mass spectrometric signals typically spanned several orders of magnitude, thus, log–log calibrations resulted in reduced error across the full concentration range.
Isothermal TGA studies were performed on a TA Instruments Discovery TGA 55 using a platinum pan under a nitrogen atmosphere. Samples were initially conditioned at 100 °C for 1 h to remove any residual moisture or solvent, and subsequently heated to 220 °C for 8 h to simulate SSP conditions. Theoretical mass loss (Δmth) was calculated from the following equation:
O)[O–])C(
O)[O–] and O
C([O–])CCCC(
O)[O–], respectively. Multiple sequence alignments required by Boltz-2 were generated with the ColabFold server.56 A single contact restraint was applied during co-folding to ensure substrate binding in the substrate binding pocket. We applied a restraint with a maximum distance of 6 Å between a pocket residue (Thr328 for SfaB and Arg93 for DdaG) and one of the chemically equivalent carboxylate carbon atoms of each substrate.
For each enzyme–substrate pair, all 100 Boltz-2 predicted affinities were analyzed as an ensemble. Because affinity_pred_value is reported as log10(IC50), with IC50 in μM, affinity predictions were characterized on the log10-transformed scale. For each complex, we calculated the mean, median, standard deviation, interquartile range, minimum, and maximum of the log10(IC50 μM) values. Nonparametric bootstrap sampling was used to estimate 95% confidence intervals for the mean and median values. For each enzyme–substrate pair, 10
000 bootstrap datasets were generated by sampling 100 affinity values with replacement from the original set of 100 Boltz-2 affinity values. The mean and median were calculated for each bootstrap dataset, and the 2.5th and 97.5th percentiles of the resulting bootstrap distributions were used as the lower and upper bounds of the 95% confidence intervals, respectively. Back-transformed IC50 values are provided in the text for interpretability. Predicted IC50 values in μM were obtained by backtransforming the logarithmic values using 10x. Substrate preferences for each enzyme were assessed by a paired comparison of matched model replicates, with the substrate having the lower affinity value considered the more strongly bound substrate.
For simplicity, we selected the model of each enzyme–substrate complex with the most favorable ligand binding affinity for structural analysis with the protein ligand interaction profiler (PLIP).53 However, the affinity conclusions are based on the full ensembles.
This manuscript has been authored by UT-Battelle, LLC under Contract No. DE-AC05-00OR22725 with the U.S. Department of Energy. The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a non-exclusive, paid-up, irrevocable, world-wide license to publish or reproduce the published form of this manuscript, or allow others to do so, for United States Government purposes. DOE will provide public access to these results of federally sponsored research in accordance with the DOE Public Access Plan (https://energy.gov/downloads/doe-public-access-plan).
Footnote |
| † These authors contributed equally. |
| This journal is © The Royal Society of Chemistry 2026 |