Abstract

Mixed-mode chromatography (MMC), which combines features of ion exchange chromatography (IEC) and hydrophobic interaction chromatography (HIC), is an interesting method for protein separation and purification. The design of MMC processes is challenging as adsorption equilibria are influenced by many parameters, including ionic strength and the presence of different salts in solution. Systematic studies on the influence of those parameters in MMC are rare. Therefore, in the present work, the influence of four salts, namely, sodium chloride, sodium sulfate, ammonium chloride, and ammonium sulfate, on the adsorption of lysozyme on the mixed-mode resin Toyopearl MX-Trp-650M at pH 7.0 and 25°C was studied systematically in equilibrium adsorption experiments for ionic strengths between 0 mM and 3000 mM. For all salts, a noticeable adsorption strength was observed over the entire range of studied ionic strengths. An exponential decay of the loading of the resin with increasing ionic strength was found until approx. 1000 mM. For higher ionic strengths, the loading was found to be practically independent of the ionic strength. At constant ionic strength, the highest lysozyme loadings were observed for ammonium sulfate, the lowest for sodium chloride. A mathematical model was developed that correctly describes the influence of the ionic strength as well as the influence of the studied salts. The model is the first that enables the prediction of adsorption isotherms of proteins on mixed-mode resins in a wide range of technically interesting conditions, accounting for the influence of the ionic strength and four salts of practical relevance.

1. Introduction

Chromatography is an important separation and purification technique in biotechnological processes. Besides long-established techniques, such as ion exchange chromatography (IEC), hydrophobic interaction chromatography (HIC), and affinity chromatography (AC), hybrid variants known as mixed-mode chromatography (MMC) exist and have attracted much interest in recent years [1, 2]. MMC combines at least two interaction modes, which can act individually or mutually. The most commonly used MMC resins are combinations of IEC and HIC resins [3]. When operated at low ionic strength in solution, the interactions of the target molecules with the ionic ligands of such a MMC resin are usually dominant, i.e., it acts basically as an IEC resin. Higher ionic strengths weaken the ionic interactions and lead to an increased attraction between the target molecules and the hydrophobic part of the resin [4], i.e., it becomes more similar to a HIC resin. Such mixed-mode resins provide wider operating ranges than IEC resins in terms of the salt concentration [5]. This higher salt tolerance compared to IEC may enable direct processing of high conductivity feedstocks [6]. Moreover, MMC can achieve capturing and purification of target molecules from complex solutions in a single process step by simultaneous operation of the different interaction modes [7, 8]. Mixed-mode resins can offer high binding capacities for proteins [9] as well as improved selectivities compared to single-mode resins [10, 11]. There are several experimental studies on the application of mixed-mode resins in the literature, e.g., [1215]. Typically, MMC is thereby applied for challenging separation processes, such as the purification of recombinant proteins, including monoclonal antibodies [15].

The conceptual design of chromatographic downstream processes for proteins is complex as the adsorption equilibria are influenced by many parameters, including ionic strength, type of salt present in the solution, temperature, and pH value, and is usually done by experimental high-throughput screenings (HTS) [16]. The higher complexity of mixed-mode resins compared to conventional resins additionally hampers the identification of suitable or even optimized process conditions. The consideration of all relevant parameters by HTS is a time-consuming and cost-intensive procedure. While there are several studies on MMC process development based on HTS, e.g., [1719], there is still a lack of systematic experimental studies and suitable models to describe the adsorption of proteins in MMC for a broad range of process parameters [20].

Experimental studies on the influence of salts on the adsorption of proteins on mixed-mode resins are mostly limited to sodium chloride, e.g., [2124]. Only few studies, e.g., [25, 26], additionally considered the influence of ammonium sulfate on the protein adsorption for ionic strengths up to 1500 mM. Furthermore, attempts to develop mathematical models for the adsorption of proteins on mixed-mode resins are scarce: Ghose et al. [27] and Gao et al. [28] both developed an adsorption isotherm model for MMC based on a modified Langmuir equation. Their models consider the effect of the pH value and the influence of sodium chloride up to concentrations of 500 mM. Similarly, the adsorption isotherm model proposed by Nfor et al. [29], which is based on Mollerup’s generalized thermodynamic framework for protein adsorption in IEC and HIC [3033], describes the influence of sodium chloride concentrations up to 2000 mM and the influence of the pH value on the adsorption of different proteins. To our knowledge, no model that considers the influence of other salts than sodium chloride on the adsorption of proteins on mixed-mode resins has been reported up to now.

In recent studies, our group developed a mathematical model for describing the influence of different salts and salt mixtures on the protein loading of a hydrophobic resin [3437]. The model is based on a Taylor series expansion in single ion molarities and enables the prediction of equilibrium adsorption isotherms for any mixture of the considered salts up to high ionic strengths. Moreover, this adsorption model was applied for predicting elution profiles of proteins for conceptual process design in HIC [38].

The first objective of the present work is the comprehensive and systematic experimental analysis of the influence of different salts and the ionic strength in industrially relevant ranges on the protein adsorption on a mixed-mode resin. Therefore, the influence of four important salts, namely, sodium chloride, sodium sulfate, ammonium chloride, and ammonium sulfate, on the adsorption of the model protein lysozyme on the mixed-mode resin Toyopearl MX-Trp-650M at 25°C and pH 7.0 was studied experimentally in a systematic way for ionic strengths up to 3000 mM. A temperature of 25°C and a pH value of 7.0 were chosen as these values lie well within the ranges in which lysozyme is stable [3941]. The conditions studied here are therefore typical for adsorption processes with lysozyme in practice. We furthermore compare the results of this work with the adsorption of lysozyme on a single-mode HIC resin at the same conditions [35]. However, also studying the influence of different temperatures and pH values on the adsorption of lysozyme on mixed-mode resins will be of interest in future work. The second objective of the present work is the development of a mathematical model that describes the influence of the four studied salts on the adsorption of lysozyme. Such a model is presented here; it is based on the same principles as our previous model used in [3437]. However, not the single ions are considered but the salts. The model that is presented here is easy to use and based on only few parameters. It not only enables the correlation of experimental adsorption data but also allows the prediction of unmeasured adsorption isotherms of lysozyme on Toyopearl MX-Trp-650M for the entire range of technically relevant ionic strengths of the four studied salts at 25°C and pH 7.0. The proposed model approach can also be transferred to other salts, temperatures, pH values, proteins, or mixed-mode resins in a straightforward manner. It is therefore a valuable tool for the conceptual process design in MMC as it can, e.g., be applied for predicting elution profiles of proteins in MMC.

2. Materials and Methods

2.1. Materials

Hen egg white lysozyme (Lys, ) was obtained from Sigma-Aldrich with a purity over 90%. The mixed-mode resin Toyopearl MX-Trp-650M was obtained from Tosoh Bioscience. Toyopearl MX-Trp-650M combines ionic and hydrophobic interactions as the tryptophan ligand of the resin possesses weak carboxyl cation exchange groups as well as indole hydrophobic functional groups. The salts used for the buffer preparation, namely, monosodium phosphate (NaH2PO4 · 2 H2O), disodium phosphate (Na2HPO4 · 2 H2O), sodium chloride (NaCl), sodium sulfate (Na2SO4), ammonium chloride (NH4Cl), and ammonium sulfate ((NH4)2SO4), were of analytical grade. Ammonium chloride was obtained from Bernd Kraft; all other salts were obtained from Carl-Roth. 1 N sodium hydroxide (NaOH) and 1 N hydrochloric acid (HCl), which were used for the adjustment of the pH value of the buffer solutions, were obtained from Carl-Roth. For all solutions, high-purity water, produced with a Milli-Q system from Merck Millipore, was used as solvent.

2.2. Batch Adsorption Experiments

Batch adsorption experiments were carried out on a fully automated liquid handling station Freedom EVO 200 from Tecan. The equilibrium measurements were carried out with 50.9 μl resin, prepared with a Media Scout ResiQuot device from Atoll, and 500 μl lysozyme stock solution in Riplate deepwell plates from Ritter. Preliminary lysozyme solutions were prepared gravimetrically with an AX205 scale from Mettler Toledo by dissolving lysozyme in 25 mM sodium phosphate buffer. Prior to the addition of lysozyme, the pH value of the buffer was adjusted to 7.0 with 1 N sodium hydroxide or 1 N hydrochloric acid. The pH value was measured with a 780 pH meter from Metrohm. Besides the lysozyme solutions, additional solutions of the studied salt (NaCl, Na2SO4, NH4Cl, or (NH4)2SO4) were prepared gravimetrically using the same scale by dissolving the salts in 25 mM sodium phosphate buffer. The pH values of the salt solutions were adjusted as described above. The lysozyme stock solutions for the adsorption experiments were prepared by mixing the preliminary lysozyme solution and salt solution in a 1:1 () ratio yielding ionic strengths of the studied salts varying between 0 mM and 3000 mM in the final lysozyme stock solution. The temperature during the adsorption experiments was always 25°C. The equilibration time on a tempered shaker from Inheco was 6 h. The lysozyme concentration in solution after equilibration was analyzed by UV absorption at 280 nm using an Infinite M200 plate reader from Tecan. As blank value, the same procedure was performed three times with lysozyme-free solutions for each studied salt and the mean value of these measurements was used for calibrating the concentration determination of the lysozyme-containing solutions. The equilibrium experiments and the analytics were performed as described in detail by Werner et al. [42]. The relative error of the equilibrium loading of the resin determined in these experiments is about 10%, except for very small loadings [42].

2.3. Correlation of Individual Adsorption Isotherms

The experimental results of the adsorption isotherms were correlated with the semi-empirical approach of Oberholzer and Lenhoff [43, 44]: where is the lysozyme concentration in solution, is the lysozyme loading of the resin, and is the adsorption equilibrium constant. and are lumped fitting parameters, which were made dimensionless by introducing mg . Equation (1) was only used here for correlating individual isotherms. This could have been done with any mathematical correlation that yields an adequate description of the data without substantially influencing the results of this work.

2.4. Data Processing

In the present work, MATLAB was used for all data processing and modeling steps. For better comparability of the influence of the four studied salts (NaCl, Na2SO4, NH4Cl, and (NH4)2SO4), the results are discussed in terms of ionic strength , which is defined using the molarities of the ions. The small concentration of buffer salts ( mM) that was present in all solutions was not taken into account for calculating the ionic strength throughout this work.

3. Results and Discussion

3.1. Experimental Adsorption Isotherms

Figure 1 shows the experimental equilibrium adsorption data obtained in this work for solutions containing sodium chloride, sodium sulfate, ammonium chloride, or ammonium sulfate together with the corresponding individual correlations, cf. Equation (1). Lysozyme carries a positive net charge at the studied pH value of 7.0 [41]. Therefore, as expected, the highest lysozyme loadings were observed when no salt was added to the solution ( mM), i.e., if pure 25 mM sodium phosphate buffer was used as solvent. The positively charged lysozyme adsorbs on the negatively charged carboxyl groups of the mixed-mode resin. Shielding effects of salts are nearly negligible as only very small ion concentrations of the buffer salts are present in the solution here. The addition of all studied salts leads to a decrease of the lysozyme loading with increasing ionic strength up to approx. 1000 mM. This behavior can be explained by the shielding of the negatively charged groups of the resin, i.e., the cation exchange part of the resin, by the cations in the solution as well as by the shielding of the positively charged sites on the lysozyme by the anions. The cations compete with lysozyme for the negatively charged positions of the resin. With increasing ionic strength, the cation exchange sites of the mixed-mode resin are increasingly occupied by the present cations, such that these sites are not accessible for the adsorption of lysozyme anymore. Moreover, the shielding of lysozyme by anions decreases the attraction between lysozyme and the negatively charged resin. Overall, at low ionic strengths, the cation exchange mechanism is the dominant contribution to the capacity of the mixed-mode resin, but this contribution decreases with increasing ionic strength. Very high ionic strengths lead to a complete occupation of the cation exchange sites of the resin with the cations in solution, such that ionic interactions between lysozyme and the resin are less pronounced in this regime. At the same time, an increase in ionic strength results in the exposure of hydrophobic regions on the lysozyme surface. The hydrate shells surrounding the lysozyme molecules are degraded by the ions in solution interacting with the surrounding water molecules, also known as salting-out effect [45]. Hence, the exposure of the hydrophobic surface regions of lysozyme can lead to an enhanced interaction with the hydrophobic indole groups of the resin. Therefore, at high ionic strengths, namely, above about  mM, one can expect that the lysozyme loading increases with increasing ionic strength due to favored hydrophobic interactions. However, the lysozyme loadings show no significant dependence on : all differences among the experimental data taken in that region are within the experimental uncertainties. This was observed for all investigated salts, cf. Figure 1. This behavior can be attributed to the superposition of the two discussed opposing effects: the increasing ionic strength leads to a weakening of the cation exchange mechanism of the adsorption whereas it favors the hydrophobic interactions by a salting-out effect, which leads to a better adsorption. For all salts, a significantly higher capacity was found in the cation exchange region (at very low ) than in the hydrophobic interaction region (at high ), which is in agreement with other reports on mixed-mode resins in the literature [25, 28, 29]. Furthermore, the adsorption isotherms tend to be almost linear in the hydrophobic interaction region whereas a considerable curvature was observed in the cation exchange region for all studied salts.

In Figure 2, the experimental results for all studied salts at a constant ionic strength of 250 mM are depicted.

Solutions containing sulfate ions show higher lysozyme loadings than solutions containing chloride ions. The rather large size of sulfate ions might hamper the accumulation around the positively charged lysozyme and thereby lead to a weaker shielding effect compared to the rather small chloride ions. Hence, in solutions containing sulfate ions, lysozyme adsorbs significantly stronger at the ion exchange ligands than it adsorbs in solutions of the same ionic strength but containing chloride ions. The highest loadings were observed for ammonium sulfate, the lowest for sodium chloride. This was found for all studied ionic strengths for which significant differences between the lysozyme loadings were obtained, i.e., for rather low ionic strengths up to approx. 1000 mM, at which ionic adsorption is assumed predominant. This agrees well with the previous presumption: ammonium ions are significantly larger than sodium ions and, hence, can be expected to lead to a weaker shielding effect than sodium ions. The differences in the influence of sodium chloride and ammonium sulfate on the lysozyme loading of Toyopearl MX-Trp-650M correspond well with the findings of Wu et al. [25], who studied the influence of these two salts on the loading of the same mixed-mode resin with human serum albumin (HSA).

The results of the correlation of the individual experimental isotherms with Equation (1) are summarized in Table S.1 in the Supplementary Materials. For describing the influence of the studied salts and the ionic strength on the adsorption of lysozyme on the studied mixed-mode resin, we have developed a model along the lines of an approach previously introduced by our group [3437].

3.2. Model of Loading as a Function of Ionic Strength

For a better comparability of the adsorption of lysozyme for the different ionic strengths and studied salts, the adsorption isotherm data were reduced in the same way as in our previous work [34]: the lysozyme loading of the mixed-mode resin, as given by the individual correlations, cf. Equation (1), at a constant lysozyme concentration was considered and compared for the studied conditions. As an example, the results for the loading for  mM (cf. the vertical dashed lines in Figures 1 and 2) as a function of the ionic strength are shown in Figure 3 for all studied salts.

As described before, we could expect that the lysozyme loading decreases in the cation exchange region due to the cations in solution competing with the lysozyme molecules for the negatively charged positions of the resin, whereas the loading increases in the hydrophobic interaction region with increasing ionic strength due to, among others, the salting-out effect. However, the experimental results indicate an exponential decay of ( mM) with increasing ionic strength up to approx.  mM and a rather constant behavior of ( mM) at a low level for  mM for all four salts.

The exponential decay in the cation exchange region, where high lysozyme loadings were observed, indicates that the effect of adding a certain amount of salt on the adsorption is directly proportional to the amount of salt already present in the solution in this region. Here, the shielding effect of additional ions is strong if only few ions are already present in the solution, whereas it is weaker if there are already many ions in the solution. It can be speculated that the finding also holds for other proteins and MMC resins, and, more importantly, that it also holds for IEC resins. Results from other studies, e.g., [21, 26, 4648], support this assumption.

While basically no protein adsorption can be expected on a single-mode IEC resin at high ionic strengths [48], the results of the present work demonstrate that there is at least some adsorption of lysozyme in the entire studied range of ionic strengths (up to 3000 mM). This is shown in Figure 4, which constitutes enlarged representations of the results shown in Figure 3. Furthermore, in Figure 4, results for the adsorption of lysozyme on a single-mode HIC resin from a previous work of our group [35] are shown. For all salts except for sodium chloride, at high ionic strengths ( mM), the adsorption of lysozyme on the studied mixed-mode resin is comparable to the adsorption on the pure HIC resin. This demonstrates the salt-tolerant nature of the mixed-mode resin.

Based on the experimental results, an exponential correlation was used in this work for describing the dependence of the lysozyme loading at constant lysozyme concentration on the ionic strength in solution for each salt : where , , and are model parameters. If no salt is present in the solution, i.e., for  mM, equals . For very high salt concentrations, the influence of salt type and concentration vanishes, cf. Figure 3, and equals the difference . Hence, the parameters and do not depend on the salt and its concentration in solution. The parameter , in contrast, is salt-specific and describes how strong a specific salt influences the adsorption of lysozyme on the mixed-mode resin. The resulting parameters obtained by fitting Equation (2) to the data of this work are reported in Table S.2 in the Supplementary Materials. Sodium chloride shows the greatest effect on the loading in the cation exchange region, which might be caused by a strong shielding effect of the comparably small sodium and chloride ions, leading to a large value of and therefore to the smallest adsorption of lysozyme with sodium chloride compared to the other salts, cf. Equation (2). By contrast, ammonium sulfate has the smallest effect, which might be caused by a rather weak shielding effect of the comparably large ammonium and sulfate ions, leading to a small value of and therefore to the strongest adsorption of lysozyme with ammonium sulfate compared to the other salts, cf. Equation (2). Hence, the larger the value of , the stronger is the negative effect of the respective salt on the adsorption of lysozyme on the MMC resin and therefore the lower is the adsorption of lysozyme at a specific ionic strength.

As depicted in Figures 3 and 4, Equation (2) represents the dependence of the lysozyme loading on the ionic strength very well for  mM and all four studied salts: the correlations with Equation (2) (lines) reproduce the data from the individual adsorption isotherms (symbols) within their uncertainty. The procedure described above was carried out for lysozyme concentrations in the interval  mM with a step size of 0.01 mM and additionally for  mM to obtain information on highly diluted solutions, for which adsorption isotherms are usually linear and can be described by the adsorption equilibrium constant alone. At high lysozyme concentrations like  mM, the equilibrium loading of the mixed-mode resin is nearly reached in some cases. Figure S.1 in the Supplementary Materials illustrates exemplarily the results for  mM and  mM. Good correlations with Equation (2) were found for all lysozyme concentrations. The results for the parameters , , and obtained from fitting Equation (2) for  mM are shown in Table S.2 in the Supplementary Materials.

3.3. Dependence of Model Parameters on Lysozyme Concentration

The results presented in Table S.2 in the Supplementary Materials indicate a dependence of the parameters , , and on the lysozyme concentration . This is also illustrated in Figure 5.

If no salt is present in the solution, i.e., for  mM, the lysozyme loading is given by , cf. Equation (2). Hence, equals the equilibrium adsorption isotherm of lysozyme if pure 25 mM sodium phosphate buffer without additional salts is used as solvent. Therefore, the dependence of on was modeled with the semi-empirical correlation function of Oberholzer and Lenhoff, cf. Equation (1), which was already used for correlating the individual isotherms in this work: where , , and correspond to the parameters , , and from Equation (1), respectively. The obtained parameters , , and are summarized in Table 1.

At very high ionic strengths, the lysozyme loading is described by the difference in Equation (2). At these conditions, the equilibrium adsorption isotherms show a rather linear behavior as illustrated in Figure 1. Therefore, a linear correlation was introduced for the dependence of on here: where denotes the adsorption equilibrium constant of the considered linear isotherms at high ionic strengths and was fitted to the data for and shown in Figure 5, cf. also Table S.2 in the Supplementary Materials. The obtained number for is given in Table 1. The results for depicted in Figure 5 almost coincide with those for . This is due to the rather small numbers of the difference compared to the numbers of , i.e., the small lysozyme loadings at high ionic strengths compared to those for  mM.

The results in Figure 5 indicate an exponential dependence of the salt-dependent parameters on for all four salts. Therefore, a correlation function similar to Equation (2) was deployed for describing for each salt as a function of : where , , and are the salt-dependent model parameters, which were fitted to the data for , cf. Table S.2 in the Supplementary Materials and symbols in Figure 5. The obtained parameters are summarized in Table 1. As represented by the lines in Figure 5, an excellent correlation of as a function of was obtained with Equation (5). If no lysozyme is present in the solution, i.e.,  mM, equals the parameter . For highly concentrated lysozyme solutions, is given by the difference . The parameter describes the exponential decay of over the lysozyme concentration.

The lines in Figure 5 denote the model results for the dependence of , , and on obtained with Equation (3), (4), or (5), respectively. An excellent description for all parameters was obtained.

3.4. Prediction of Equilibrium Adsorption Isotherms

Inserting Equations (3)–(5) in Equation (2) yields:

Equation (6) can be considered as a new adsorption isotherm model. Using only the parameters from Table 1, Equation (6) enables the prediction of adsorption isotherms for lysozyme on the mixed-mode resin Toyopearl MX-Trp-650M at 25°C and pH 7.0 under the influence of sodium chloride, sodium sulfate, ammonium chloride, or ammonium sulfate for ionic strengths up to  mM.

In Figure 6, isotherms obtained with this model (lines) are compared to the experimental data obtained in this work (symbols).

For all four salts and all considered ionic strengths, a very good agreement between model and experimental data was observed. The model correctly describes the significantly higher lysozyme loadings at lower ionic strengths due to strong ionic interactions and the rather low loadings at higher ionic strengths that can be attributed to occupied cation exchange sites of the resin and enhanced hydrophobic interactions, which are in general weaker than ionic interactions, for all salts. For ionic strengths between 1000 mM and 3000 mM, the model shows no significant dependence of the adsorption isotherms on the ionic strength or salt. Additionally, the model correctly describes the increasing linearity of the adsorption isotherms with increasing ionic strength. Hence, the proposed model can be used for the prediction of adsorption isotherms within the trained parameter range and therefore constitutes the basis of simulating column experiments and elution profiles with mixed-mode resins.

4. Conclusions

In the present work, we have systematically studied the influence of different salts and the ionic strength on the adsorption of lysozyme on the mixed-mode resin Toyopearl MX-Trp-650M at pH 7.0 and 25°C. Substantial adsorption of lysozyme was found for the entire range of studied ionic strengths, which demonstrates the salt tolerance of the mixed-mode resin. For all studied salts (sodium chloride, sodium sulfate, ammonium chloride, and ammonium sulfate), the lysozyme loading of the resin decreases exponentially with increasing ionic strength up to approx. 1000 mM. This observation can be explained by physical considerations: the effect of adding a certain amount of salt is directly proportional to the amount of salt already present in the solution. Hence, at low ionic strengths, where the cation exchange mechanism of the mixed-mode resin is dominant, each ion in solution has a stronger influence on the adsorption by shielding the ionic ligands of the resin than at higher ionic strengths. For a constant ionic strength in this region, the strongest influence was found for sodium chloride, which might be explained by a strong shielding effect of the rather small sodium and chloride ions. We speculate that the observed exponential dependence of the protein loading on the ionic strength is generic for chromatography in which ion exchange mechanisms are dominant. At higher ionic strengths, i.e., from approx. 1000 mM up to 3000 mM, where hydrophobic interactions are predominant, no significant influence of different salts or the ionic strength was observed. For all salts except sodium chloride, the found adsorption in this region is similar to the adsorption of lysozyme on a single-mode HIC resin that was studied in a previous work. Rather linear adsorption isotherms were observed in the hydrophobic interaction region.

We furthermore present a mathematical model for the adsorption of lysozyme on Toyopearl MX-Trp-650M. The model successfully describes the influence of the different studied salts, the ionic strength, and the lysozyme concentration in solution on the lysozyme loading. The model allows the prediction of equilibrium adsorption isotherms for a wide range of process parameters, e.g., for in principle any practically relevant ionic strength, including the transition region from the cation exchange region to the hydrophobic interaction region. Hence, also adsorption isotherms for not studied conditions can be predicted with the presented model approach based on only few parameters. The proposed model approach can be transferred to other salts, temperatures, pH values, proteins, or mixed-mode resins in a straightforward manner and ultimately be applied for simulation, conceptual design, and optimization of separation processes of proteins with mixed-mode chromatography.

Data Availability

The adsorption data used to support the findings of this study are included within the Supplementary Materials.

Conflicts of Interest

The authors declare no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Supplementary Materials

Additional model results and parameters are included within the Supplementary Materials. (Supplementary Materials)