research papers
Structure and stability of an apo thermophilic esterase that hydrolyzes polyhydroxybutyrate
aSchool of Chemistry and Biochemistry, Georgia Institute of Technology, 901 Atlantic Drive NW, Atlanta, GA 30332, USA, and bKimberly-Clark Corporation, 1400 Holcomb Bridge Road, Roswell, GA 30076, USA
*Correspondence e-mail: raquel.lieberman@chemistry.gatech.edu
Pollution from plastics is a global problem that threatens the biosphere for a host of reasons, including the time scale that it takes for most plastics to degrade. Biodegradation is an ideal solution for remediating bioplastic waste as it does not require the high temperatures necessary for thermal degradation and does not introduce additional pollutants into the environment. Numerous organisms can scavenge for bioplastics, such as polylactic acid (PLA) or poly-(R)-hydroxybutyrate (PHB), which they can use as an energy source. Recently, a promiscuous PHBase from the thermophilic soil bacterium Lihuaxuella thermophila (LtPHBase) was identified. LtPHBase can accommodate many substrates, including PHB granules and films and PHB block copolymers, as well as the unrelated polymers polylactic acid (PLA) and polycaprolactone (PCL). LtPHBase uses the expected Ser–His–Asp catalytic triad for hydrolysis at an optimal near 70°C. Here, the 1.75 Å resolution of apo LtPHBase is presented and its chemical stability is profiled. Knowledge of its substrate preferences was extended to different-sized PHB granules. It is shown that LtPHBase is highly resistant to unfolding, with barriers typical for thermophilic enzymes, and shows a preference for low-molecular-mass PHB granules. These insights have implications for the long-term potential of LtPHBase as an industrial PHB hydrolase and shed light on the evolutionary role that this enzyme plays in bacterial metabolism.
PDB reference: apo LtPHBase, 9byu
1. Introduction
Pollution from plastics is a global problem that jeopardizes the health of our biosphere, with massive amounts of plastic deposited annually into landfills, freshwater sources and the ocean. Only 21% of petroleum-derived plastic waste is successfully managed through recycling and other waste-reducing processes. Most plastics have a lengthy half-life and take decades or longer to degrade. In addition, petroleum-based plastics can decompose into microplastics (<5 µm), which can re-enter human metabolism via the consumption of livestock, marine life and drinking water, all of which pose health hazards (Idris et al., 2023; Andrady, 2011). Moreover, the SARS-CoV-2 pandemic has exacerbated the use of single-use disposable plastic materials in the form of personal protective equipment such as gloves (Idris et al., 2023, Kaushal et al., 2021).
Most synthetic plastics in use are petroleum-based (Mohanan et al., 2020). Although efforts are under way to more efficiently degrade such plastics, an alternative is the use of bioplastics to manufacture consumable goods. Bioplastics are made from biomass and other readily available materials that are non-petroleum-based and that are naturally sourced (Tokiwa et al., 2009). Common examples include polylactic acid (PLA) and poly(butylene adipate-co-terephthalate). These plastics are used in compostable cutlery and consumable goods that are commercially available today. Bioplastics degrade more quickly than fossil-fuel-based plastics (Tokiwa et al., 2009). However, bioplastics can still produce microplastic particles if not processed properly during recycling, and not all bioplastic materials are readily degraded (Pascoe Ortiz, 2023).
Efforts to discover bioplastic-degrading organisms and adapt them to industrial waste-management schemes have significantly accelerated in recent years. Many microorganisms, including fungi, algae, bacteria and insects, possess bioplastic-depolymerizing enzymes (Yang et al., 2020; Ping et al., 2017; Sánchez, 2020). These organisms can scavenge for bioplastics, such as PLA or PHB, as an energy source. Biodegradation is an ideal solution for remediating bioplastic waste as it does not require the high temperatures necessary for thermal degradation and does not introduce additional pollutants into the environment.
Recently, we identified a PHBase from the thermophilic soil bacterium Lihuaxuella thermophila (LtPHBase; UniProt A0A1H8IKU3). LtPHBase exhibits optimal near 70°C (Thomas et al., 2022). Like all PHBases, LtPHBase contains the canonical Ser–His–Asp catalytic triad (Ser121–His270–Asp197). The of LtPHBase with serendipitously bound ethylene glycol (EG) and 2-propanol (IPA) revealed a relatively shallow active site compared with other PHBases (Thomas et al., 2022), which we posited to explain its ability to accommodate many different substrates beyond PHB granules and films, such as PHB block copolymers and the unrelated polymers PLA and PCL.
Here, we present the LtPHBase lacking any bound ligands and profile its chemical stability. We also extend our knowledge of substrate preferences to different-sized PHB granules. We show that LtPHBase is highly resistant to unfolding, with barriers typical for thermophilic enzymes, and is an enzyme that prefers low-molecular-mass PHB granules. These insights have implications for the long-term potential of LtPHBase for use as an industrial PHB hydrolase and shed light on the evolutionary role that this enzyme plays in bacterial metabolism.
of2. Materials and methods
2.1. Chemicals and laboratory reagents
All chemicals and β-D-1-thiogalactopyranoside (IPTG), antibiotics and high-molecular-mass PHB polymer granules. Low-molecular-mass polyhydroxybutyrate polymer samples were obtained from Polysciences, Warrington, Pennsylvania, USA. The polymers were supplied as a powder/granulated powder and were used without further characterization. All laboratory supplies were purchased from Fisher Scientific. Chemically competent Escherichia coli cells were purchased from New England Biolabs, Beverly, Massachusetts, USA.
resins were purchased from Millipore–Sigma, St Louis, Missouri, USA, including all buffer components, media, isopropyl2.2. Expression, purification and activity assay
The enzyme was expressed and purified to et al., 2022). The final protein preparation was concentrated to 10 mg ml−1 in 5.0 mM Tris–HCl pH 9.0 and frozen in small aliquots at −20°C until use. Enzymatic activity was determined via direct measurement of the β-hydroxybutyrate liberated during the reaction using the MAK272 hydroxybutyrate assay kit from Millipore–Sigma, St Louis, Missouri, USA as described previously (Thomas et al., 2022). The substrate concentration was 0.1–50 µM, which spans the range of Km values for each substrate.
following the procedure described previously (Thomas2.3. Crystallization and structure determination
Crystallization conditions for apo LtPHBase were identified by sparse-matrix screening that avoided 2-propanol (IPA) and ethylene glycol (EG). Crystals were grown at 25°C via hanging-drop crystallization in VDX plates with a reservoir buffer consisting of 18% PEG 4000, 100 mM Tris–HCl pH 9.0, 50 mM MgCl2, 50 mM KCl. A 4 µl drop volume was utilized by mixing 2 µl protein stock and 2 µl reservoir buffer. Hexagonal rod-shaped crystals formed in 3–4 days with approximate dimensions of 1–2 mm in length and a cross-sectional diameter of 0.2–0.3 mm. Crystals were harvested, incubated briefly in reservoir buffer supplemented with 20% PEG 400 and flash-cooled in liquid nitrogen. Data were collected at 100 K on beamline 5.0.1 at the Advanced Light Source (ALS) using an oscillation angle of 0.25° and a total rotation of 180°. The diffraction data were processed using XDS and scaled using XSCALE (Kabsch, 2010) using CC1/2 > 0.9 as a cutoff for resolution. The apo structure of LtPHBase was solved by using PDB entry 8daj (Thomas et al., 2022) as a search model. Phenix.xtriage did not detect any evidence of (Liebschner et al., 2019). The apo LtPHBase structure was built in Coot (Emsley et al., 2010) and refined using phenix.refine (Liebschner et al., 2019; McCoy et al., 2007; Table 1). PDBePISA (Krissinel & Henrick, 2007) was used to identify crystal contacts and PyMOL (Schrödinger) was used to generate figures. The structure has been deposited in the PDB with accession code 9byu.
|
2.4. Substrate docking
Apo LtPHBase was submitted as the receptor and the hydroxybutyrate trimer from PDB entry 2d81 as the ligand to the online EDock server (Zhang et al., 2020). The image for the top docked pose was generated in PyMOL.
2.5. Thermodynamic unfolding
LtPHBase was chemically denatured in guanidinium hydrochloride (GdHCl; Millipore–Sigma, St Louis, Missouri, USA) at a concentration of 10 mM in 10 mM Tris–HCl pH 9.0. After mixing the protein with various amounts of the denaturant, the tubes were incubated at 25°C for 12 h to ensure that the samples were at equilibrium. Fluorescence emission intensity was measured using a Chirascan V100 spectropolarimeter (Applied Photophysics, Beverly, Massachusetts, USA) that was fitted with a CCD-array fluorescence detector, which simultaneously collects fluorescence emission intensity from 190 to 900 nm. The excitation wavelength of 280 nm was delivered via a fiber-optic cable connected to a light source. Emission spectra were analyzed in the region 300–450 nm. The fraction of unfolded protein was determined by fitting the fluorescence emission intensity to the relationship of Clarke & Fersht (1993),
where fu is the fraction of unfolded enzyme, m is the m value for cooperativity of the transition, x is the concentration of denaturant, d50 is the midpoint, R is the universal and K is the absolute temperature. The Ku and the free-energy change ΔGu for the unfolding reaction were calculated according to
and
The free energy in the absence of denaturant, , was calculated by a linear least-squares fit of the data (using QtiPlot; https://qtiplot.com/) in the linear portion of the unfolding curve to
2.6. Kinetics of unfolding
The kinetics of the unfolding/refolding reaction as a function of denaturant concentration were determined using an SF3 stopped-flow device coupled to the Chirascan V100 spectropolarimeter. For unfolding reactions, LtPHBase in 10 mM Tris–HCl pH 9.0 was rapidly mixed with GdHCl at concentrations ranging from 0.5 to 8 M. The for mixing in the SF3 is approximately 2 ms. Single- or double-exponential functions were fitted to the experimental data using both the supplied instrument software and the QtiPlot scientific graphics package (https://qtiplot.com/). A single exponential provided the best fit, as evaluated by the coefficient of determination (R2), the root-mean-squared error (RMSE) and χ2 per degree of freedom (χ2/DOF) for multiple tested fits. The chevron plot was analyzed by determining the slope of the line between 0.5 and 3.5 M GdHCl to obtain the value of mf and between 3.5 and 7.0 M GdHCl to obtain the value of mu. These slopes describe the dependence of the folding and unfolding reactions on GdHCl. Initially, the values of kf and ku in the absence of GdHCl, i.e. water only, were obtained by extrapolating the linear portions of the chevron curve to the y axis. Unfolding and refolding terms were also calculated by fitting the equation (Garg et al., 2022)
where kobs is the observed and are the folding and unfolding rate constants in water, mf and mu are the slopes of the chevron curve on either side of the minimum and x is the concentration of GdHCl.
The value of Cm, the GdHCl concentration at which 50% of the protein is unfolded, was determined as the GdHCl value at the minimum of the chevron curve and by calculating
The free energy of the folding/unfolding reaction is calculated by
where ΔG is the free energy and and are the folding and unfolding rate constants in water.
3. Results
3.1. of apo LtPHBase
In our previous LtPHBase structure (PDB entry 8daj; Thomas et al., 2022), two IPA molecules from the crystallization reservoir buffer and a molecule of EG from the cryoprotectant were bound serendipitously within the active site. These molecules, and a network of water molecules, formed a hydrogen-bonding network with surface loops and the catalytic triad. To test the hypothesis that LtPHBase might adopt different active-site loop conformations in the absence of bound ligands, we identified new crystallization conditions that do not include IPA and we used PEG 400 as the cryoprotectant instead of EG.
The 1.75 Å resolution apo LtPHBase structure was solved by using the IPA/EG-bound LtPHBase structure (Table 1). The overall root-mean-squared deviation (r.m.s.d.) is 0.2 Å, indicating a high level of agreement between the two structures. The two most notable differences are distal from the active site. Firstly, the loop comprising residues 221–225 is disordered in apo LtPHBase, whereas in the IPA/EG-bound structure the loop was readily modeled (Fig. 1a). Secondly, the five N-terminal residues adopt different conformations. These changes in the structure can be attributed to differences in crystal contacts (Figs. 1a and 1b) and suggest flexibility in these regions.
Subtle changes are also seen in the active site. Crystal contacts stabilize the active-site conformation observed in each structure (Fig. 1b). The original IPA/EG-bound structure is somewhat more occluded compared with the new apo structure in the region above Cys40 and Cys78, which appear to form a disulfide bond (Fig. 1b). Slight changes in the surface helix (residues 155–164) appear to be primarily responsible for this difference. Most notably, in the apo structure Met155 is oriented towards Phe81 and away from the Tyr159 side chain, which itself is in a new orientation compared with IPA/EG-bound LtPHBase (Fig. 1b). Arg91 is also found in an alternative conformation compared with that in IPA/EG-bound LtPHBase (Fig. 1c). In both structures the thermal B factors of Met155, Thr156 and Ser157 are considerably higher than the remaining helix and active-site environments (Fig. 1c), suggesting that this region is somewhat mobile and that additional loop conformations might be accessible in solution. Despite these differences, just like the IPA/EG-bound structure, a hydroxybutyrate trimer can still be docked computationally into the active site, indicating that the observed conformation is competent for catalysis (Fig. 1d). Taken together, the changes detected between the two structures imply some dynamics of surface features in solution, but the structures are otherwise nearly indistinguishable.
3.2. Thermodynamics and kinetics of LtPHBase unfolding
To further characterize apo LtPHBase, we next probed its room-temperature unfolding properties. Excitation at λex = 280 nm monitors the microenvironment of the five buried tryptophans in LtPHBase. The fluorescence emission intensity decreases, and the emission λ maximum is red-shifted from 328 to 341 nm, as a function of denaturant concentration. The unfolding transition begins in ∼2.8 M denaturant and is complete in 4.2 M GdHCl, with a midpoint of 3.4 M GdHCl. The unfolding reaction is best-fitted to a simple two-state model (N↔D; Fig. 2, Table 2). After the sample in 7.0 M GdHCl was dialyzed overnight versus 10 mM Tris–HCl pH 9.0 at 4°C, the enzyme regained activity, with a Km of 8.3 mM and a kcat of 3.2 s−1 when measured using the assay described in Section 2. These values are comparable to those for the native enzyme against the same substrate, indicating that chemical unfolding is fully reversible. The model fit shows that the enzyme stability in the absence of denaturant, , is 8.9 kcal mol−1.
|
Kinetic unfolding measurements further support a two-state unfolding transition for LtPHBase at room temperature. Single mixing experiments of LtPHBase and GdHCl indicated that LtPHBase could rapidly be denatured at GdHCl concentrations above the midpoint of unfolding determined by equilibrium unfolding measurements (not shown). Fluorescence emission intensity at 328 nm (excitation wavelength of 280 nm) was monitored over the course of the unfolding or refolding reaction. For the unfolding reactions, enzyme at a concentration of 10 µM in 10 mM Tris–HCl pH 9.0 was rapidly mixed with a 10× volume of GdHCl at concentrations between 0.5 and 8.0 M. When the folded enzyme is mixed with 8.0 M denaturant, rapid loss of fluorescence emission intensity at 328 nm is observed over 10 s. The curve is best-fitted to a single exponential, yielding a first-order for unfolding of 2.7 s−1 (Fig. 3a). Similarly, LtPHBase refolding is observed in a rapid mixing experiment by diluting unfolded LtPHBase in reaction buffer/7.0 M GdHCl with a tenfold excess of reaction buffer alone (final concentration of 0.7 M GdHCl), seen as a rapid increase in the fluorescence emission intensity over a 20 s time course. These kinetic data were also best-fitted to a single exponential, yielding a first-order for folding of 17.9 s−1 (Fig. 3b). The chevron plot (the logarithm of the observed unfolding/folding versus the denaturant concentration; Fig. 3c, Table 2) shows a linear dependence of kobs on both sides of the GdHCl midpoint. Comparison of the slopes of the folding reaction, mf (the left arm of the plot), and the unfolding reaction, mu (the right arm of the plot), indicates that the folding reaction is faster than the unfolding reaction. Compared with the average folding and unfolding rates of thermophilic proteins that unfold via a two-state mechanism (Glyakina & Galzitskaya, 2020), LtPHBase folds at a similar rate [ln(kf) = 4.22 compared with an average ln(kf) of 4.75 ± 1.2] but unfolds faster [ln(ku) = −8.55, average ln(ku) = −5.63 ± 1.12; range −12 to 4].
In summary, both kinetics and equilibrium experiments indicate that LtPHBase unfolding at room temperature is best described as a simple two-state process. The denaturant midpoint and free-energy values are within 1.1 kcal mol−1. This, and the symmetric shape of the chevron plot, indicates that there is not likely to be an intermediate form in the unfolding/refolding pathways and that the reaction is fully reversible.
3.3. Enzymatic activity versus PHB molecular mass
Given the shallow active site and broad substrate repertoire, we next probed the effect of PHB polymer size on the activity of LtPHBase. Enzyme kinetic data were acquired on PHB samples with molecular masses ranging from 500 to 350 000 Da (Table 3). LtPHBase hydrolyzes all substrates within this molecular-mass range; however, there is an inverse relationship between PHB mass and the catalytic parameters. Namely, the Km, which is indicative of substrate binding, ranges from an average of 3.5 to 19.2 mM and the catalytic efficiency, kcat, ranges from 11.2 to 0.82 s−1 over the PHB from 500 to 350 000 Da. In sum, LtPHBase prefers lower-molecular-mass PHB. Whether the LtPHBase reaction is endohydrolytic or exohydrolytic is unknown, but the enzyme rates could be higher for the lower-molecular-mass PHBs if the catalytic mechanism is solely endohydrolytic.
‡For kcat/Km, the values in parentheses are the calculated Gaussian error propagation: ΔZ = Z[(ΔX/X)2 + (ΔY/Y)2]1/2, where X is kcat and Y is Km. |
4. Discussion
A future that uses an environmentally circular bioplastics scheme would likely involve bioplastic degradation as part of industrial waste management. Such a process would further require a bioplastic-degrading enzyme that is able to function at and withstand the associated harsh conditions. To this end, we identified a promiscuous PHBase from the thermophilic soil bacterium L. thermophila (LtPHBase). LtPHBase can accommodate many different substrates (Thomas et al., 2022).
Despite the fact that other PHBase structures have active sites with a range of accessibilities (Hisano et al., 2006; Jendrossek et al., 2013), our study shows that the active site of LtPHBase has loops that adopt virtually the same position whether in the apo form or bound by IPA/EG. Thus, it appears that PHBases have diverged evolutionarily with regard to the degree of conformational change that is required to accommodate PHB substrate(s) into the binding cleft and ultimately into position by the catalytic triad. This mechanistic difference is most likely to drive the observed substrate preferences.
To our knowledge, LtPHBase is the first PHBase to be surveyed for PHB molecular-mass preference. Given that LtPHBase hydrolyzes multiple different substrates, we hypothesized that there would not be a substantial molecular-mass preference. Our experiments did not support this hypothesis. LtPHBase is likely to prefer lower-molecular-mass, non-film polymers due to the apparent lack of conformational access to the catalytic triad, as seen in our apo structure. The lower-molecular-weight polymer samples are more readily accommodated into the binding cleft of the enzyme, so that the labile ester bond is positioned within the catalytic triad. The 350 000 Da sample is a poorer substrate than the PHB film initially utilized to characterize the enzyme (Thomas et al., 2022). Yet, the enzyme is capable of hydrolyzing high-molecular-mass PHB polymers and films, so there is potentially a change in the active-site loop conformation that allows these larger substrates to be accommodated in the active site that is not captured in the current The importance of the active-site loops in LtPHBase are underscored by unsuccessful attempts to express a soluble enzyme with a shortened loop in this region (not shown).
Generally, thermophilic enzymes are kinetically slower than their mesophilic or psychrophilic counterparts at their optimum temperature due to the need of thermophiles to invest more structural et al., 2007). Direct kinetic comparisons across known PHBases are hampered by differences in the available PHB substrates, primarily molecular mass, average granule diameter or film dimensions, used in different laboratories for the enzymatic assays. Still, casual inspection of the rates for mesophilic forms of the enzyme [for example PHBases from Aeromonas caviae (Amir et al., 2024), Microbacterium paraoxydans (Sayyed et al., 2019) and Penicillium expansum (Gowda & Shivakumar, 2015)] indicate that LtPHBase is slower than these mesophilic PHBases. Thus, another curiosity regarding LtPHBase is the interplay between its slow kinetics and the aforementioned substrate preferences. One potential caveat of the current work is that crystallization was performed and the kinetics of folding and unfolding were determined at room temperature, which is ∼30°C below the optimal temperature of the enzyme kinetics. Although this may affect the absolute values of rate constants and loop positions in the the results are directionally correct because the enzyme is functional at room temperature. Future efforts will be directed towards solution dynamics of LtPHBase not accessible in the crystalline state or kinetics at room temperature, which may in turn underlie preferences for both different sizes and types of substrates.
during the electrostatic reorganization process (RocaThe high degree of thermostability of LtPHBase reported by us previously, as expected, is manifested in a high degree of resistance to chemical with an apparent transition midpoint of 3.4 M GdHCl. In general, resistance to is typical of thermophilic and hyperthermophilic enzymes (Feller, 2018). For example, the MTH1880 protein from Methanobacterium thermoautotrophicum has a chemical unfolding midpoint of 4.0 M GdHCl (Kim et al., 2016). Stabilization of the native structure is a hallmark of thermostable enzymes and a mechanism that is entropically driven (Pica & Graziano, 2016). Although there is not a singular mechanism that explains thermophilic enzymes, resistance to unfolding is often associated with more well packed hydrophobic cores and an increase in the number of electrostatic interactions relative to mesophilic counterparts (Reed et al., 2013). LtPHBase has an extensive hydrophobic core and a substantial number of electrostatic interactions. Kinetic stability is in part afforded by a slower unfolding rate and a faster folding rate (Mukaiyama & Takano, 2009), and in the case of LtPHBase the folding rate is faster than its chemical unfolding rate. Interestingly, this is a smaller difference in rates than has been seen in other thermophilic enzymes, but still results in an increase in the activation barrier to unfolding (Karshikoff et al., 2015; Wittung-Stafshede, 2004).
LtPHBase belongs to the secreted subfamily of PHBases, which most likely function to scavenge hydroxybutyrate (HB) from a few environmental biopolyester polymers. Presumably, the liberated HB is transported into the bacterium for use as an energy source. However, there are no other PHB-cycle enzymes that can be identified in the annotated genome or by BLAST analyses (Altschul et al., 1990) against the L. thermophila genome using homologous PHB-biosynthesis enzymes (not shown). L. thermophila is not unique in this regard. A BLAST search of LtPHBase against microbial genomes reveals PHBase homologs in organisms such as Bacillus aquiflavi and Shouchella shacheensis that similarly lack PHB-biosynthesis enzymes, although these are In the long term, beyond clarification of the relationship between the substrate scope and kinetics of LtPHBase for the potential use of PHB as a renewable biopolymer, is the solution of the mystery of why L. thermophila and other organisms have evolved to utilize an imported polymer as an energy source.
Acknowledgements
Data were collected on Advanced Light Source beamline 5.0.1 through a partnership with SER-CAT. The Berkeley Center for Structural Biology is supported by the Howard Hughes Medical Institute, Participating Research Team members and the National Institutes of Health, National Institute of General Medical Sciences, ALS-ENABLE grant P30 GM124169. The Advanced Light Source is a Department of Energy Office of Science User Facility under Contract No. DE-AC02-05CH11231. The PILATUS detector on beamline 2.0.1 was funded under NIH grant S10OD021832. The PILATUS detector on beamline 5.0.1 was funded under NIH grant S10OD026941. SER-CAT is supported by its member institutions, equipment grants (S10_RR25528, S10_RR028976 and S10_OD027000) from the National Institutes of Health, and funding from the Georgia Research Alliance. This research used resources of the Advanced Photon Source, a US Department of Energy (DOE) Office of Science user facility operated for the DOE Office of Science by Argonne National Laboratory under Contract No. DE-AC02-06CH11357.
Funding information
This work was partially funded by NSF grant 1817796 to RLL. GMT was supported by the GAANN P200A210014.
References
Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. (1990). J. Mol. Biol. 215, 403–410. CrossRef CAS PubMed Web of Science Google Scholar
Amir, M., Bano, N., Gupta, A., Zaheer, M. R. & Roohi (2024). Biodegradation, 35, 137–153. Google Scholar
Andrady, A. L. (2011). Mar. Pollut. Bull. 62, 1596–1605. CrossRef CAS PubMed Google Scholar
Clarke, J. & Fersht, A. R. (1993). Biochemistry, 32, 4322–4329. CrossRef CAS PubMed Web of Science Google Scholar
Emsley, P., Lohkamp, B., Scott, W. G. & Cowtan, K. (2010). Acta Cryst. D66, 486–501. Web of Science CrossRef CAS IUCr Journals Google Scholar
Feller, G. (2018). Semin. Cell Dev. Biol. 84, 129–137. CrossRef CAS PubMed Google Scholar
Garg, M., Sharma, D. & Kumar, R. (2022). Biophys. Chem. 290, 106892. CrossRef PubMed Google Scholar
Glyakina, A. V. & Galzitskaya, O. V. (2020). Biomolecules, 10, 197. CrossRef PubMed Google Scholar
Gowda, U. S. V. & Shivakumar, S. (2015). 3 Biotech, 5, 901–909. PubMed Google Scholar
Hisano, T., Kasuya, K.-I., Tezuka, Y., Ishii, N., Kobayashi, T., Shiraki, M., Oroudjev, E., Hansma, H., Iwata, T., Doi, Y., Saito, T. & Miki, K. (2006). J. Mol. Biol. 356, 993–1004. CrossRef PubMed CAS Google Scholar
Idris, S. N., Amelia, T. S. M., Bhubalan, K., Lazim, A. M. M., Zakwan, N., Jamaluddin, M. I., Santhanam, R., Amirul, A. A., Vigneswari, S. & Ramakrishna, S. (2023). Environ. Res. 231, 115988. CrossRef PubMed Google Scholar
Jendrossek, D., Hermawan, S., Subedi, B. & Papageorgiou, A. C. (2013). Mol. Microbiol. 90, 649–664. Web of Science CrossRef CAS PubMed Google Scholar
Kabsch, W. (2010). Acta Cryst. D66, 125–132. Web of Science CrossRef CAS IUCr Journals Google Scholar
Karshikoff, A., Nilsson, L. & Ladenstein, R. (2015). FEBS J. 282, 3899–3917. Web of Science CrossRef CAS PubMed Google Scholar
Kaushal, J., Khatri, M. & Arya, S. K. (2021). Clean. Eng. Technol. 2, 100083. CrossRef Google Scholar
Kim, H., Kim, S., Jung, Y., Han, J., Yun, J.-H., Chang, I. & Lee, W. (2016). PLoS One, 11, e0145853. CrossRef PubMed Google Scholar
Krissinel, E. & Henrick, K. (2007). J. Mol. Biol. 372, 774–797. Web of Science CrossRef PubMed CAS Google Scholar
Liebschner, D., Afonine, P. V., Baker, M. L., Bunkóczi, G., Chen, V. B., Croll, T. I., Hintze, B., Hung, L.-W., Jain, S., McCoy, A. J., Moriarty, N. W., Oeffner, R. D., Poon, B. K., Prisant, M. G., Read, R. J., Richardson, J. S., Richardson, D. C., Sammito, M. D., Sobolev, O. V., Stockwell, D. H., Terwilliger, T. C., Urzhumtsev, A. G., Videau, L. L., Williams, C. J. & Adams, P. D. (2019). Acta Cryst. D75, 861–877. Web of Science CrossRef IUCr Journals Google Scholar
McCoy, A. J., Grosse-Kunstleve, R. W., Adams, P. D., Winn, M. D., Storoni, L. C. & Read, R. J. (2007). J. Appl. Cryst. 40, 658–674. Web of Science CrossRef CAS IUCr Journals Google Scholar
Mohanan, N., Montazer, Z., Sharma, P. K. & Levin, D. B. (2020). Front. Microbiol. 11, 580709. CrossRef PubMed Google Scholar
Mukaiyama, A. & Takano, K. (2009). Int. J. Mol. Sci. 10, 1369–1385. CrossRef PubMed CAS Google Scholar
Pascoe Ortiz, S. (2023). PLoS Biol. 21, e3002045. CrossRef PubMed Google Scholar
Pica, A. & Graziano, G. (2016). Biopolymers, 105, 856–863. CrossRef CAS PubMed Google Scholar
Ping, L.-F., Chen, X.-Y., Yuan, X.-L., Zhang, M., Chai, Y.-J. & Shan, S.-D. (2017). Int. J. Biol. Macromol. 104, 1238–1245. CrossRef CAS PubMed Google Scholar
Reed, C. J., Lewis, H., Trejo, E., Winston, V. & Evilia, C. (2013). Archaea, 2013, 373275. Web of Science CrossRef PubMed Google Scholar
Roca, M., Liu, H., Messer, B. & Warshel, A. (2007). Biochemistry, 46, 15076–15088. CrossRef PubMed CAS Google Scholar
Sánchez, C. (2020). Biotechnol. Adv. 40, 107501. PubMed Google Scholar
Sayyed, R. Z., Wani, S. J., Alyousef, A. A., Alqasim, A., Syed, A. & El-Enshasy, H. A. (2019). PLoS One, 14, e0212324. CrossRef PubMed Google Scholar
Thomas, G. M., Quirk, S., Huard, D. J. E. & Lieberman, R. L. (2022). Protein Sci. 31, e4470. CrossRef PubMed Google Scholar
Tokiwa, Y., Calabia, B. P., Ugwu, C. U. & Aiba, S. (2009). Int. J. Mol. Sci. 10, 3722–3742. Web of Science CrossRef PubMed CAS Google Scholar
Wittung-Stafshede, P. (2004). Biochim. Biophys. Acta, 1700, 1–4. PubMed CAS Google Scholar
Yang, Y., Wang, J. & Xia, M. (2020). Sci. Total Environ. 708, 135233. CrossRef PubMed Google Scholar
Zhang, W., Bell, E. W., Yin, M. & Zhang, Y. (2020). J. Cheminform, 12, 37. CrossRef PubMed Google Scholar
This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.