Abstract
Many polar organisms produce antifreeze proteins (AFPs) and ice-binding proteins (IBPs) to protect themselves from ice formation. As IBPs protect cells and organisms, the potential of IBPs as natural or biological cryoprotective agents (CPAs) for the cryopreservation of animal cells, such as oocytes and sperm, has been explored to increase the recovery rate after freezing–thawing. However, only a few IBPs have shown success in cryopreservation, possibly because of the presence of protein denaturants, such as dimethyl sulfoxide, alcohols, or ethylene glycol, in freezing buffer conditions, rendering the IBPs inactive. Therefore, we investigated the thermal and chemical stability of FfIBP isolated from Antarctic bacteria to assess its suitability as a protein-based impermeable cryoprotectant. A molecular dynamics (MD) simulation identified and generated stability-enhanced mutants (FfIBP_CC1). The results indicated that FfIBP_CC1 displayed enhanced resistance to denaturation at elevated temperatures and chemical concentrations, compared to wildtype FfIBP, and was functional in known CPAs while retaining ice-binding properties. Given that FfIBP shares an overall structure similar to DUF3494 IBPs, which are recognized as the most widespread IBP family, these findings provide important structural information on thermal and chemical stability, which could potentially be applied to other DUF3494 IBPs for future protein engineering.
Similar content being viewed by others
Introduction
Cryopreservation is an important technique in modern biotechnology and has many applications in areas such as in vitro fertilization (IVF), stem cell research, and tissue engineering. Cryopreservation allows biological samples to be stored for long periods of time without compromising their viability, which is critical for research and clinical applications. Many organisms living in cold areas produce ice-binding proteins (IBPs) or antifreeze proteins (AFPs) as a defense mechanism against the damage caused by ice formation at subzero temperatures. (Although the term AFPs is more specific in explaining proteins with thermal hysteresis (TH) activity and/or ice recrystallization inhibition (IRI) activity1,2, the term IBPs has a broader conceptual scope, accurately representing the function of the proteins and encompassing the previously identified IBPs and AFPs. Therefore, in the present study, we used both AFP and IBP interchangeably to encompass the previously identified IBPs and AFPs3.) Glycosylated AFPs (AFGPs) and AFPs were first discovered in the blood serum of polar fish, such as the Antarctic notothenioid4, winter flounder5,6, and ocean pout7. AFPs have also been identified in insects, plants, and bacteria and categorized into different types based on their structural features8.
AFPs and IBPs can inhibit ice growth by binding to specific planes of ice and lowering the freezing temperature through a mechanism called absorption-inhibition9,10. Their binding properties are measured via the gap between the melting and freezing temperatures, referred to as thermal hysteresis (TH) activity. Previous research has primarily focused on the interactions between ice and IBPs (or AFPs) and the effects of their binding, such as lowering the freezing point and inhibiting the growth of ice crystals, which are unique abilities of biological molecules that are considered crucial characteristics for protecting host cells. Given that AFPs and IBPs can protect host cells from freezing, research has explored their potential applicability in preserving important mammalian cells such as gametes and embryos11. For instance, AFPIII, one of the most well-studied AFPs, has been reported to improve the preservation of sperm, oocytes, and embryos from various organisms and can be used with various freezing methods, including vitrification, chilling, and cryopreservation. Glycosylated AFGPs have also been investigated for the cryopreservation of sperm, oocytes, and embryos.
Dimethyl sulfoxide (DMSO), glycerol, ethylene glycol, methanol, and 1,2-propanediol are common membrane-permeable cryoprotective agents (CPAs) utilized in the cryopreservation or vitrification of mammalian cells and tissues, known as membrane-permeable cryoprotectants11,12,13. Among these, DMSO and glycerol are well-recognized CPAs primarily owing to their ability to reduce the formation of intracellular ice during cell freezing. Therefore, 10% DMSO or 20% glycerol is commonly utilized in cryopreservation procedures to maximize cell recovery after the freeze–thaw process, even in studies evaluating the effectiveness of AFPs or IBPs as CPAs14. Despite the discovery and study of various AFPs and IBPs that protect organisms from freezing, little is known about their chemical tolerance to CPAs. Therefore, understanding how these denaturants affect the stability and functionality of proteins under cryoprotective conditions is crucial for the development of improved cryopreservation methods and the design of proteins with enhanced stability in such environments.
The FfIBP from Flavobacterium frigoris PS1 exhibits high TH activity and relatively strong recrystallization inhibition activity15. The FfIBP structure is composed of parallel β-strands in a helical pattern, which is mainly stabilized by inner hydrophobic interactions and intramolecular hydrogen bonds with long α415. Additionally, the FfIBP has a capping head region that consists of two small turns and β-strands, which contain a disulfide bond. This disulfide bond increases the overall rigidity of the protein and reduces residual movement at the ice-binding site, thereby augmenting TH activity and protein stability16. These findings suggest that the introduction of additional disulfide bonds may further enhance TH activity and protein stability. Consequently, the present study aimed to investigate the stability and chemical resistance of FfIBP in the context of CPAs and to establish functionally improved IBPs using molecular dynamics simulations.
Results
FfIBP mutagenesis guided by circular dichroism and molecular dynamics
Most hyperactive IBPs and AFPs contain more than one disulfide bond, which is one of the factors contributing to protein stability and TH activity17. Our previous study demonstrated that the structural rigidity of FfIBP is positively correlated with its thermal stability and TH activity15,16. Based on these findings, we aimed to enhance the TH activity and thermal stability of IBP by introducing an additional disulfide bond in the structurally unstable regions in addition to the pre-existing disulfide bond.
To determine the ideal locations for introducing disulfide bonds into FfIBP, circular dichroism (CD) analysis was conducted to monitor changes in the secondary structure of the protein at various temperatures. Evaluating the propensity of a secondary structure during the unfolding process provides insight into the early steps of thermal unfolding18. The temperature range used was 20–90 °C, with a 2 °C increment per step. The results indicated that the structure of FfIBP remained stable between 20 and 44 °C but began to unfold at 46 °C (Fig. 1a). The prediction of secondary structures at these temperatures indicated that the amino acid propensities of FfIBP for the β-sheet structure slightly increased. However, compared with that at 20 °C, the ends of α-helices (Helix2) decreased from 5.7% to 1.9% (50 °C), and the turn structure also decreased in concert with increasing temperatures (Fig. 1b). This result indicates that backbone disorder is initiated in the helical region of FfIBP by heat and that the helical and capping head regions of FfIBP may be promising for introducing disulfide bonds and maximizing protein stability (Fig. 1c)19.
Molecular dynamics (MD) simulations were performed at 200 nm to gain insights into the local conformational fluctuations of FfIBP. The results demonstrated that the capping head region (Ser116–Ala136) and the bottom region (Gly160–Lys165) were the most flexible areas, whereas the residues in the body region were relatively stable, with an average root-mean-square-fluctuation (RMSF) value of 0.42 Å (Fig. 2a). The top loop of the capping head region (Asp118–Ala123) exhibited an average RMSF value of 0.85 Å. Another flexible region close to the capping head region is the edge of long helix4 (Asp128–Leu132). Specifically, the root mean square deviation (RMSD) of Ser130 increased by 1.0 Å. The bottom of the body region (Gly160–Lys165) also displayed a high degree of fluctuation starting from the short-turn region (Fig. 2b, c).
Based on the CD and MD simulation data, four mutants annotated as FfIBP_CC1 (A101C-A120C), FfIBP_CC2 (S96C-L132C), FfIBP_CC3 (M249C-N270C), and FfIBP_CC4 (G160C-S176C), possibly forming disulfide bonds, were generated (Fig. 2d). The criteria for choosing residues to mutate was i) the high degree of local fluctuation of the residues, ii) the Cβ distance between two residues, which should be close to 4 Å for disulfide bond formation, and iii) not close to the ice-binding site. The mutation site for FfIBP_CC1 was located in the top loop of the capping head region, which contained a largely flexible residue Gly121. The FfIBP_CC2 is located at the interface between the body domain and the long α-helix (α4). The thermal stability analysis suggests that the α-helix is likely vulnerable to high temperatures and is important for regular spacing of the β-strands that are connected to the β-sheet on the IBS of FfIBP15. The formation of the disulfide bond at FfIBP_CC3 is anticipated to increase the stability of the global structure as the hyperactive IBPs contain disulfide bond(s) for the intra-interaction inside of the body domain. Therefore, we generated the FfIBP_CC3, which was located in the internal body domain. Finally, the most flexible region identified by MD simulation was selected as the 4th disulfide mutant (FfIBP_CC4). To confirm the disulfide formation, we solved the crystal structures of the FfIBP_CC1, CC2, CC3, and CC4 at 2.0, 2.1, 2.0, and 2.0 Å, respectively (Supplementary Table S1). Among these mutants, only FfIBP_CC1 and FfIBP_CC4 formed disulfide bonds within the intended area (Fig. 2e).
Structural superimposition indicated that wildtype and mutants are almost identical, as indicated by an RMSD value average of ~ 0.14 Å, including the orientation of the side chains throughout the structure except for the mutation regions (Supplementary Fig. S1). When the RMSD values of alpha carbon mutants were compared with the wildtype, relatively high RMSD values exceeding 0.4 Å were observed around the mutation sites of FfIBP_CC2, FfIBP_CC3, and FfIBP_CC4. However, FfIBP_CC1 did not exhibit any significant difference (Supplementary Fig. S2). Moreover, MD simulations using the FfIBP_CC1 structure indicated that the disulfide bond formation between C101 and C120 stabilized the fluctuation of residues in the capping head region (Fig. 1A and Supplementary Fig. S3).
Thermal stability of FfIBP and its derivatives
Given that thermostability is associated with flexible parts of the protein structure20 and the addition of disulfide bonds generally influences protein stability in general21,22, a label-free thermal stability assay on FfIBP and its mutants was conducted to study protein unfolding characteristics (Supplementary Fig. S4). Consistent with our hypothesis, FfIBP_CC1, which has an additional disulfide bond on the capping head region, exhibited increased thermal stability, with a melting temperature (Tm) of 65 °C compared with wildtype FfIBP (60 °C) (Fig. 3a). In contrast, FfIBP_CC2 exhibited a drastic decrease in the melting temperature by approximately 8.3 °C. The mutation at the Ser96 and Lue132 to cysteines likely distorted the conformation of the residual interactions between the β helical body region and α4, which was vulnerable during thermal melting. FfIBP_CC3 exhibited a ΔTm of − 0.8 °C compared to the wildtype, which is consistent with the hypothesis that the body region is heat-insensitive (Fig. 3a). Unexpectedly, in the case of FfIBP_CC4, the Tm decreased by − 5 °C, possibly due to conformational distortion caused by disulfide bond formation and a decrease in protein stability (Supplementary Fig. S1c). Thus, it is evident that introducing an additional disulfide bond in the capping read region directly or indirectly stabilizes the IBP. Therefore, only FfIBP_CC1 was selected for further experiments.
Activity of FfIBP and FfIBP_CC1
To investigate the impact of the additional disulfide bond on ice-binding properties, TH activity and ice re-recrystallization inhibition (IRI) activity assays were measured using wild-type FfIBP and FfIBP_CC1. TH activity was measured as described previously. Briefly, the sample was placed on a cold stage and the temperature was decreased and raised slowly to obtain a single ice crystal. The temperature gap between the melting and freezing points of the solution was considered the TH value.
The results indicated that FfIBP_CC1 exhibited a TH value of 2.3 ± 0.2 °C at a concentration of 40 µM15,16. The TH activity of FfIBP_CC1 was akin to that of wild-type FfIBP at equivalent concentrations (Fig. 3b). Moreover, when observing the ice morphology at the freezing point, the ice busted perpendicular to the C-axis with a dendritic morphology, which was similar in the presence of both proteins in solution (Fig. 3c). Additionally, the mean grains size (MGS) values for the IRI activity of both FfIBP wildtype and FfIBP_CC1 were measured after a 30-min incubation at − 6 °C, called splat cooling23. The data indicated that FfIBP_CC1 significantly inhibited the growth of ice crystals, even at a low concentration of 0.4 µM, which is comparable to that of wildtype FfIBP (Fig. 4). IRI activity did not increase, even at a higher concentration of 0.4 µM IBPs. These results suggest that the additional disulfide did not interfere with the ice-binding properties of FfIBP.
Stability of IBPs
DMSO, glycerol, ethylene glycol, methanol, and 1,2-propanediol are well-recognized CPAs13,24. However, these compounds also negatively affect proteins, potentially causing denaturation, which can result in the loss of protein activity and function. Chemical stability analysis of wild-type FfIBP and FfIBP_CC1 with the mentioned chemicals demonstrated that FfIBP_CC1 still maintained a high melting temperature (Tm), indicating that disulfide bond formation remained stable and was effective in preserving protein stability in the presence of these chemicals (Fig. 5).
These chemicals can be categorized into three groups based on their melting temperatures. In group A, which included DMSO and ethylene glycol, the Tm of FfIBP_CC1 remained relatively unchanged, even with an increase in the concentration of the chemicals, implying that the IBPs remained stable even with 30% of the chemicals in group A. In group B, which contained methanol, urea, and propylene glycol, the melting temperature patterns with increased chemical concentrations were similar, and the gaps in the melting temperatures of the two proteins were maintained. The last group, consisting of ethanol and glycerol, exhibited drastic changes in the melting temperatures in concert with an increased chemical concentration. Notably, glycerol was the only chemical that increased the melting temperature of both FfIBP and FfIBP_CC1. These findings indicate that the disulfide bond formation remains intact in IBPs in various protein denaturants and enhances the chemical stability of FfIBP and FfIBP_CC1. Furthermore, to assess the stability of IBPs over a broad pH range, the CD was used to monitor changes in the secondary structure under different pH conditions and incubation was extended at low temperatures. The results revealed that the secondary structures of both IBPs remained unchanged from pH 4.0 to pH 9.0, and 7 days of incubation at 0 °C, indicating that the FfIBP and FfIBP_CC1 are resistant to a wide pH range and low temperatures (Supplementary Fig. S5).
IRI activity in the presence of CPAs
To assess whether the ice-binding properties of FfIBP and FfIBP_CC1 were retained in the presence of CPAs, the IRI activity was measured in the presence of DMSO 10%, Gibco, CellBanker, and commercially available CPAs (Fig. 6). The minimum protein concentration showing IRI activity (0.4 µM) and a 4 µM concentration was mixed with the chemicals and incubated for 30 min at − 8 °C. The results demonstrated that the IRI activity of FfIBP and FfIBP_CC1 remained effective at a concentration of 0.4 µM and 4 µM in 10% DMSO. Moreover, in the presence of the Gibco and CellBanker, FfIBP and FfIBP_CC1 showed similar IRI activity at a concentration of 0.4 µM, indicating that the ice-binding property of FfIBP, as well as that of FfIBP_CC1, is functional in the CPAs.
Discussion
The introduction of a new disulfide bond could lead to a reduction or loss of enzymatic activity because of interference with the precise stereochemistry or dynamic movement of residues in the active site25. However, this approach has been used in protein engineering to enhance protein stability and prevent thermal damage21,22,26. AFPs and IBPs interact with solid ice, and the rigidity or stability of IBPs has been suggested as a critical factor in extending the binding time and increasing the strength of the bond between the protein and ice16. In this context, engineering disulfide bonds in IBPs can be a useful strategy to generate stable and functional analogs and understand the structure–activity relationships of IBPs.
To develop a mutant IBP with improved TH activity and thermal stability, we first identified the structurally unstable regions of FfIBP by denaturation of the secondary structure and MD simulations. We designed a mutant IBP by introducing additional disulfide bonds to stabilize these regions. Four sets (FfIBP_CC1, FfIBP_CC2, FfIBP_CC3, and FfIBP_CC4) of mutation sites for disulfide bond formation in the flexible region of the FfIBP were introduced. Successful formation of the disulfide bond at the intended location in FfIBP_CC1 and FfIBP_CC4 was confirmed by X-ray crystallography, whereas the rest did not form a disulfide bond. The failure to form a disulfide bond may be due to dihedral angle constraints and the inappropriate distance between the two cysteines. Contrary to our expectations, FfIBP_CC4 experienced a decrease of 5 °C compared to the WT form. The slight side chain rotation of the C176, observed in the structural comparison, does not explain this destabilization. However, we postulate that the forced angle change due to the disulfide bond negatively affects neighboring residues involved in the hydrophobic core, such as Trp174, Leu178, and Val198 (Supplementary Figure S1). It is worth noting that the initial 350 nm/330 nm ratio for FfIBP_CC4 is lower than that of other mutants, including the wild-type (Supplementary Figure S4). This difference may be attributed to the effect of disulfide formation on Trp174 in FfIBP_CC4. A similar destabilization upon disulfide bond formation was observed in the titin polyprotein27. These findings highlight the importance of selecting the appropriate cysteine mutant residues to facilitate disulfide bond formation and enhance thermal stability in IBPs. In addition, compared to the wildtype FfIBP, only mutant FfIBP_CC1 exhibited an improvement in thermal stability by 5 °C.
The stability increase was also observed in the thermal stability measured by CD. FfIBP_CC1 exhibited greater thermal stability by about 3 °C in comparison to wildtype CD spectra (Supplementary Figure S6). However, the melting temperature of FfIBP measured by CD spectra was only 47.8 °C, whereas Tycho showed it to be 60 °C. Initially, we speculated that the salt increases the protein stability because we removed only NaCl in the optimized buffer condition for CD, and many proteins showed salt concentration-dependent stability. However, Tm data at different salt concentrations showed no significant difference (Supplementary Figure S7). Next, we hypothesized that this discrepancy could be due to the difference in measurement methods. CD spectra used 210 nm wavelength to calculate the Tm value, which reflects the helix structure of the protein28, whereas Tycho measured the ratio of tryptophan fluorescence at 350 and 330 nm. The three tryptophans are located inside the β-helical body region of FfIBP. Moreover, this indicates that the conformational/unfolding transition of FfIBP is initiated in the helix region, and the β-helical body region is the last step in the unfolding process at higher temperature.
Interestingly, unlike wildtype, FfIBP_CC1 underwent a two-step denaturation. The initial denaturation was similar to that of FfIBP until 66 °C. When the temperature exceeded 68 °C, the spectra of delta epsilon tended to be flattened. The secondary structure prediction indicated that the α-helices (Helix1) were completely denatured, and the β-sheets were reduced beyond this temperature (Supplementary Figure S6b). Although further research is required to determine the detailed denaturation process, it is believed that the disulfide bond could affect the protein’s conformational changes during temperature increases.
In addition, FfIBP_CC1 exhibited increased thermal stability, even when exposed to high concentrations of chemical denaturants such as DMSO, MeOH, Urea, propylene glycol, EtOH, and glycerol. Usually, chemical denaturants exhibit a high affinity for water and disrupt the hydrogen bonds or hydrophobic interactions that stabilize the folded conformation of the protein by attacking the flexible local areas of the proteins29. Thus, disulfide bond formation in FfIBP likely contributes to minimizing the structural perturbations induced by these chemicals and heat inputs. Moreover, FfIBP and FfIBP_CC1 exhibited effective IRI activity in the presence of DMSO 10%, Gibco, and CellBanker, which are commercial cryoprotectants, indicating that the proteins were intact and functional in these chemicals.
The structural stability of proteins is greatly influenced by non-covalent bonding, including hydrogen bonds, electrostatic interactions, and hydrophobicity. Among these, the hydrophobic effect within the protein core appears to be a crucial factor for protein stability30, as evidenced by studies such as the methyl introduction study31 that aimed to increase protein thermal stability and the comparison between mesophilic and thermophilic proteins32. The addition of disulfide bonds to proteins has been observed to decrease thermal unfolding and denaturation by increasing the activation energy barrier to the unfolded states, thus protecting the hydrophobic core30. Moreover, the hydrophobic effect is suggested to be also mainly associated with the cold denaturation of proteins33. Therefore, disulfide bond formation in FfIBP (FfIBP_CC1) may also contribute to resistance to cold denaturation, although we did not observe the loss of both FfIBP and FfIBP_CC1 during freezing–thawing repetition (Supplementary Figure S4). Nevertheless, we anticipate that FfIBP_CC1 would exhibit beneficial effects on cell protection as a cryoprotective agent during the cryopreservation process because of the improved chemical stability34.
Taken together, our findings suggest that FfIBP is a resilient protein that maintains its ice-binding properties, even in the presence of denaturants. Furthermore, the incorporation of disulfide bonds into the capping head region of FfIBP enhances its stability and resistance to thermal and chemical denaturation. Given that FfIBP shares a similar structure with DUF3494 IBPs, the most widespread IBP family, our strategy can be applied to engineer other IBPs for use as natural cryopreservation agents35. Overall, this study provides new insights into the physicochemical stability of IBPs and suggests novel approaches for enhancing their stability.
Methods
Site-directed mutagenesis
Desired mutants were generated using site-directed mutagenesis. The primers used for mutagenesis are listed in Supplementary Table S2. The mutated sequences were confirmed via DNA sequencing.
Overexpression and purification of FfIBP and mutants
The constructs of pColdI-FfIBP mutants were transformed into Escherichia coli Rosetta-gamiTM 2 (DE3) (Novagen, Darmstadt, Germany) and grown in Luria–Bertani (LB) medium containing ampicillin 100 µg/mL. The culture was initially grown at 37 °C until the optical density at 600 nm reached 0.6. Subsequently, the temperature was lowered to 15 °C, and the culture was induced with 0.5 mM IPTG for protein expression. After a 72-h incubation, the cells were harvested, resuspended, and lysed using a sonication buffer (30 mM Tris–HCl pH 8.5, 150 mM NaCl, and 5 mM imidazole). The supernatant was separated via centrifugation and loaded onto a pre-equilibrated Ni–NTA resin agarose column (Qiagen, Hilden, Germany). The bound protein was eluted and treated with factor Xa protease (New England Biolabs, Ipswich, MA, USA) to remove the 6×His tag. Subsequently, proteins were purified using a Superdex 200 prep-grade column (GE Healthcare, Chicago, IL, USA). The purified protein was collected and concentrated in a buffer consisting of 30 mM Tris–HCl (pH 8.5 and 300 mM NaCl using an Amicon Ultra Centrifugal Filter (Merck, Darmstadt, Germany).
Crystallization and data collection
The crystallization conditions for FfIBP mutants were obtained from a previous report36, which was originally used for wild-type FfIBP. The hanging-drop vapor diffusion method was employed by mixing 1.5 μL of the protein solution (approximately 20 mg/mL for each protein) with 1.5 μL of the crystallization solution (0.1 M sodium citrate pH 4.4, 3 M sodium chloride). X-ray diffraction datasets were collected at the BL-5C beamline of the Pohang Light Source operated by the Pohang Accelerator Laboratory in Pohang, Republic of Korea. XDS software was used for indexing, integrating, and scaling the dataset37.
Structure determination and refinement
The MOLREP program from the CCP4i suite was used to determine the crystal structure of the mutants through a molecular replacement method using the wildtype FfIBP structure as a search model. The structure was refined using a combination of the REFMAC5 and Phenix. Refine software. MolProbity was used to check model quality38, and all structural representations were generated using PyMOL39. The mutant coordinates and structural factors were deposited in the Protein Data Bank RCSB under the accession codes 8X0Z, 8X1L, 8X1O, and 8X1P. Detailed data collection and refinement statistics are provided in Supplementary Table S1.
MD simulation
Protein structures were subjected to MD simulations using GROMACS 2022 with the CHARMM Force Field40,41. The process and syntax are referred to in the GROMACS documentation. After the removal of water molecules from the crystallography, the FfIBP structure was used as the initial coordinate. The structure was located in the center of the cubic system, solvated with a TIP4P water model, and neutralized with 0.15 M sodium chloride. The sides of the cubic box had dimensions of 15 nm (3,375,000 Å3). One thousand steps of steepest descent (SD) followed by 50,000 steps of the Adopted Basis Newton–Raphson (ABNR) algorithms were applied for minimization. After 100 ps of isochoric (NVT) and isobaric isothermal equilibration (NPT), the MD simulation was performed at 300 K with a 2 fs time step for 100 ns in an NPT ensemble under periodic boundary conditions (PBC) in all three dimensions. The long-range electrostatic interactions were calculated using the particle-mesh Ewald (PME) method42. The RMSD, RMSF, and distances between residues from the trajectory files were analyzed using GROMACS. Pymol39 and VMD43 programs were used for analysis and figure generation.
TH activity
The TH activity of the IBPs was measured as previously described15,36. The sample was placed in the holes of an aluminum disk on the stage and quickly frozen at approximately − 40 °C using a nanoliter osmometer (Otago Osmometers, Dunedin, New Zealand). The temperature was slowly increased until only one ice crystal remained (melting temperature). Subsequently, the temperature was slowly lowered again (at a rate of ~ 0.5 °C per min) until the ice crystal started to grow (freezing temperature). The shape of the ice crystals at the freezing point was observed under a microscope, and photographs were taken using a Canon Power Shot 620. The measurements were performed in a buffer containing 30 mM Tris–HCl (pH 8.5) and 150 mM NaCl at various sample concentrations.
IRI assay
IRI was measured using a modified “splat cooling” method23. Protein samples were diluted to the desired concentration in a buffer containing 30 mM Tris–HCl (pH 8.5) and 150 mM NaCl. A 10 µL of the droplet was released from a height of 1.20 m onto a 16 mm glass coverslip resting on a polished aluminum surface that was situated on a pool of liquid nitrogen. Upon impact with the cold coverslip, the droplet instantaneously froze into a wafer-like structure with a diameter of approximately 12 mm and a thickness of 10 µm. The coverslip with the frozen sample was then placed on a Linkam TMHS6000 cold stage (Linkam Scientific Instruments, Surrey, UK) and held at − 6 °C or − 8 °C for 30 min. The IRI assay was performed using an Olympus BX51 microscope equipped with crossed polarizers, and images were captured using an Olympus DP71 CCD camera. The IRI activity was evaluated by measuring the largest dimension of individual ice grains in various fields of view. The MGS was calculated by averaging the values of the 30 largest ice grains. A smaller MGS indicates a higher level of IRI activity, as it indicates the effective prevention of ice crystals from growing during the constant continuous process at − 6 °C or − 8 °C.
Thermal shift assay
Label-free thermal shift assays were performed using a Tycho NT.6 instrument (NanoTemper Technologies GmbH, Munich, Germany). Protein samples were diluted to a concentration of 1 mg/mL in a buffer containing 20 mM Tris–HCl at pH 8.0. Following the incubation, the samples were heated in a glass capillary from 35 to 95 °C at a rate of 30 °C per min. During this process, the intrinsic fluorescence emanating from tryptophan and tyrosine residues was recorded at wavelengths of 330 and 350 nm, respectively. Subsequently, the ratio of fluorescence at 350/330 nm and the melting temperature were calculated using the internal evaluation capabilities of a Tycho instrument. The data presented in this study represent the results of three independent experiments.
Circular dichroism spectroscopy to analyze the unfolding of IBP
CD spectroscopy was performed using a Chirascan CD spectrometer (Applied Photophysics Co., Leatherhead, UK) and a Peltier temperature controller. Protein samples (0.1 mg/mL) in 10 mM Tris–HCl (pH 8.5) were loaded into quartz cuvettes with a path length of 1 mm (Hellma, Plainview, NY, USA), and the absorbance was measured from 190 to 250 nm with 0.5 nm steps and 0.5 s/step at a 1 mm path length. The temperature was set from 20 to 90 °C with a 2 °C increment per step. The measurements were repeated three times, and the buffer spectra were subtracted from the average spectra to obtain the final trace. The analyzed data were obtained using the interactive circular dichroism analysis software CDNN28,44 and expressed as the Delta epsilon, given in M−1 cm−1. The proportion of α-helix, β-strand, turn, and other (mostly coil) structures at different temperatures was calculated using BeStSel software45.
Stability of IBPs at different pHs and low temperature
To evaluate the impact of various pH levels on the stability of FfIBP, its structural denaturation profile was determined using CD within a pH range of 4–9. At pH 4 and 5, FfIBP and FfIBP_CC1 were diluted in sodium acetate buffer (10 mM sodium acetate buffer); for pH 6 and 7, sodium phosphate buffer (10 mM sodium phosphate buffer) was used; for pH 8 and 9, Tris–HCl buffer (10 mM Tris–HCl) was used.
To determine the reversibility of cold denaturation, FfIBP and FfIBP_CC1 were incubated at 0 °C for 1, 2, 3, and 7 days. Following each incubation period, the proteins were centrifuged at 16,000×g for 10 min, and the supernatant was analyzed using CD to evaluate their stability46. In both experiments conducted under various pH conditions and at low temperatures, proteins were used at a final concentration of 0.1 mg/mL and assessed by scanning at 190–250 nm as described above.
CPA information
Two commercial CPAs were examined to determine their impact on protein stability and functionality: Gibco™ Synth-a-Freeze™ Cryopreservation Medium (Gibco, Thermo Fisher Scientific, Waltham, MA, USA) and CELLBANKER 1 (Amsbio Inc., Cambridge, M). Gibco™ Synth-a-Freeze™ Cryopreservation Medium is employed for the cryopreservation of a wide variety of mammalian cells, including keratinocyte cell, embryonic stem cell, neural stem cell, and mesenchymal stem cell. CELLBANKER 1 is designed for a broad spectrum of mammalian cells, such as CHO, Jurkat, CEM, and K562. Both media contain 10% DMSO as cryoprotectant.
Data availability
The structure datasets generated during the current study are available in the RCSB PDB Protein Data Bank repository. The PDB accession numbers are as follows: 8X0Z for FfIBP_CC1, 8X1L for FfIBP_CC2, 8X1O for FfIBP_CC3, and 8X1P for FfIBP_CC4.
References
Davies, P. L. Ice-binding proteins: A remarkable diversity of structures for stopping and starting ice growth. Trends Biochem. Sci. 39, 548–555 (2014).
Białkowska, A., Majewska, E., Olczak, A. & Twarda-clapa, A. Ice binding proteins: Diverse biological roles and applications in different types of industry. Biomolecules 10, 274 (2020).
Kim, H. J. et al. Marine antifreeze proteins: Structure, function, and application to cryopreservation as a potential cryoprotectant. Mar. Drugs 15, 27 (2017).
DeVries, A. L. & Wohlschlag, D. E. Freezing resistance in some antarctic fishes. Science 163, 1073–1075 (1969).
Cheng, A. & Merz, K. M. Ice-binding mechanism of winter flounder antifreeze proteins. Biophys. J. 73, 2851–2873 (1997).
Marshall, C. B., Chakrabartty, A. & Davies, P. L. Hyperactive antifreeze protein from winter flounder is a very long rod-like dimer of α-helices. J. Biol. Chem. 280, 17920–17929 (2005).
Hew, C. L. et al. Multiple genes provide the basis for antifreeze protein diversity and dosage in the ocean pout Macrozoarces americanus. J. Biol. Chem. 263, 12049–12055 (1988).
DeVries, A. L. & Cheng, C. H. C. Antifreeze proteins and organismal freezing avoidance in polar fishes. Fish Physiol. 22, 155–201 (2005).
Pertaya, N., Marshall, C. B., Celik, Y., Davies, P. L. & Braslavsky, I. Direct visualization of spruce budworm antifreeze protein interacting with ice crystals: Basal plane affinity confers hyperactivity. Biophys. J. 95, 333–341 (2008).
Knight, C. A. & Wierzbicki, A. Adsorption of biomolecules to ice and their effects upon ice growth. 2. A discussion of the basic mechanism of ‘antifreeze’ phenomena. Cryst. Growth Des. 1, 439–446 (2001).
Robles, V., Valcarce, D. G. & Riesco, M. F. The use of antifreeze proteins in the cryopreservation of gametes and embryos. Biomolecules 9, 1–12 (2019).
Zhang, T. & Rawson, D. M. Feasibility studies on vitrification of intact zebrafish (Brachydanio rerio) embryos. Cryobiology 33, 1–13 (1996).
Murray, K. A. & Gibson, M. I. Chemical approaches to cryopreservation. Nat. Rev. Chem. 6, 579–593 (2022).
Carpenter, J. F. & Hansen, T. N. Antifreeze protein modulates cell survival during cryopreservation: Mediation through influence on ice crystal growth. Proc. Natl. Acad. Sci. U. S. A. 89, 8953–8957 (1992).
Do, H., Kim, S. J., Kim, H. J. & Lee, J. H. Structure-based characterization and antifreeze properties of a hyperactive ice-binding protein from the Antarctic bacterium Flavobacterium frigoris PS1. Acta Crystallogr. Sect. D Biol. Crystallogr. 70, 1061–1073 (2014).
Hwang, J. et al. Importance of rigidity of ice-binding protein (FfIBP) for hyperthermal hysteresis activity and microbial survival. Int. J. Biol. Macromol. 204, 485–499 (2022).
Graham, L. A., Qin, W., Lougheed, S. C., Davies, P. L. & Walker, V. K. Evolution of hyperactive, repetitive antifreeze proteins in beetles. J. Mol. Evol. 64, 387–398 (2007).
Wijma, H. J., Floor, R. J. & Janssen, D. B. Structure- and sequence-analysis inspired engineering of proteins for enhanced thermostability. Curr. Opin. Struct. Biol. 23, 588–594 (2013).
Singh, N. & Li, W. Recent advances in coarse-grained models for biomolecules and their applications. Int. J. Mol. Sci. 20, 3774 (2019).
Vihinen, M. Relationship of protein flexibility to thermostability. Protein Eng. 1, 477–480 (1987).
Kazlauskas, R. Engineering more stable proteins. Chem. Soc. Rev. 47, 9026–9045 (2018).
Liu, T. et al. Enhancing protein stability with extended disulfide bonds. Proc. Natl. Acad. Sci. U. S. A. 113, 5910–5915 (2016).
Knight, C. A., Wen, D. & Laursen, R. A. Nonequilibrium antifreeze peptides and the recrystallization of ice. Cryobiology 32, 23–34 (1995).
Yang, H. et al. Temporal and concentration effects of methanol on cryopreservation of zebrafish (Danio rerio) sperm. Zebrafish 17, 233–242 (2020).
Matsumura, M. & Matthews, B. W. Control of enzyme activity by an engineered disulfide bond. Science 243, 792–794 (1989).
Wetzel, R., Perry, L. J., Baase, W. A. & Becktel, W. J. Disulfide bonds and thermal stability in T4 lysozyme. Proc. Natl. Acad. Sci. U. S. A. 85, 401–405 (1988).
Manteca, A. et al. The influence of disulfide bonds on the mechanical stability of proteins is context dependent. J. Biol. Chem. 292, 13374–13380 (2017).
Greenfield, N. J. Using circular dichroism spectra to estimate protein secondary structure. Nat. Protoc. 1, 2876–2890 (2007).
Wallqvist, A., Covell, D. G. & Thirumalai, D. Hydrophobic interactions in aqueous urea solutions with implications for the mechanism of protein denaturation. J. Am. Chem. Soc. 120, 427–428 (1998).
Rahban, M. et al. Thermal stability enhancement: Fundamental concepts of protein engineering strategies to manipulate the flexible structure. Int. J. Biol. Macromol. 214, 642–654 (2022).
Pace, C. N. Contribution of the hydrophobic effect to globular protein stability. J. Mol. Biol. 226, 29–35 (1992).
Vijayabaskar, M. S. & Vishveshwara, S. Comparative analysis of thermophilic and mesophilic proteins using Protein energy networks. BMC Bioinform. 11, 1–7 (2010).
Dias, C. L. et al. The hydrophobic effect and its role in cold denaturation. Cryobiology 60, 91–99 (2010).
Lee, J. et al. Effects of three different types of antifreeze proteins on mouse ovarian tissue cryopreservation and transplantation. PLoS ONE 10, 1–16 (2015).
Vance, T. D. R., Bayer-Giraldi, M., Davies, P. L. & Mangiagalli, M. Ice-binding proteins and the ‘domain of unknown function’ 3494 family. FEBS J. 286, 855–873 (2019).
Do, H., Lee, J. H., Lee, S. G. & Kim, H. J. Crystallization and preliminary X-ray crystallographic analysis of an ice-binding protein (FfIBP) from Flavobacterium frigoris PS1. Acta Crystallogr. Sect. F Struct. Biol. Cryst. Commun. 68, 806–809 (2012).
Kabsch, W. X. D. S. xds. Acta Crystallogr. Sect. D Biol. Crystallogr. 66, 125–132 (2010).
Williams, C. J. et al. MolProbity: More and better reference data for improved all-atom structure validation. Protein Sci. 27, 293–315 (2018).
DeLano, L. & Pymol, W. An open-source molecular graphics tool. {CCP4} Newsl Protein Crystallogr. 40, 1–8 (2002).
Abraham, M. J. et al. Gromacs: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 1–2, 19–25 (2015).
Afgan, E. et al. The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2018 update. Nucleic Acids Res. 46, W537–W544 (2018).
Petersen, H. G. Accuracy and efficiency of the particle mesh Ewald method. J. Chem. Phys. 103, 3668–3679 (1995).
Humphrey, W., Dalke, A. & Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 14, 33–38 (1996).
Louis-Jeune, C., Andrade-Navarro, M. A. & Perez-Iratxeta, C. Prediction of protein secondary structure from circular dichroism using theoretically derived spectra. Proteins Struct. Funct. Bioinform. 80, 374–381 (2012).
Micsonai, A. et al. BeStSel: Webserver for secondary structure and fold prediction for protein CD spectroscopy. Nucleic Acids Res. 50, W90–W98 (2022).
Sanfelice, D. & Temussi, P. A. Cold denaturation as a tool to measure protein stability. Biophys. Chem. 208, 4–8 (2016).
Acknowledgements
This research was part of the project titled “Development of potential antibiotic compounds using polar organism resources” (KIMST 20200610, Korea Polar Research Institute (KOPRI) Grant PM24030), funded by the Ministry of Oceans and Fisheries, Korea. (J.H.L.). This study was supported by the KOPRI; Grant Number PE24120). We would like to thank the staff at the X-ray core facility of the Korea Basic Science Institute (KBSI; Ochang, Korea) and BL-5C of the Pohang Accelerator Laboratory (PAL, Pohang, Korea) for their kind help with X-ray diffraction data collection.
Author information
Authors and Affiliations
Contributions
Y.N. assisted with data curation, formal analysis, methodology, validation, visualization, and writing – original draft; D.N. assisted with methodology, validation, and visualization; T.H. assisted with data curation, formal analysis, and methodology; B.K. assisted with methodology, validation, and visualization; J.L. assisted with conceptualization, funding acquisition, investigation, project administration, supervision, and writing – review, and editing; H.D. assisted with conceptualization, investigation, supervision, writing – original draft, writing – review, and editing.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Nam, Y., Nguyen, D.L., Hoang, T. et al. Engineered ice-binding protein (FfIBP) shows increased stability and resistance to thermal and chemical denaturation compared to the wildtype. Sci Rep 14, 3234 (2024). https://doi.org/10.1038/s41598-024-53864-w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-024-53864-w
Keywords
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.