- Open Access
Interactive association between processing induced molecular structure changes and nutrient delivery on a molecular basis, revealed by cutting-edge vibrational biomolecular spectroscopy
Journal of Animal Science and Biotechnology volume 10, Article number: 85 (2019)
This study was conducted to determine protein molecular structure profiles and quantify the relationship between protein structural features and protein metabolism and bioavailability of blend pelleted products (BPP) based on co-products (canola or carinata) from processing with different proportions of pulse pea screenings and lignosulfonate chemical compound.
The protein molecular structures were determined using the non-invasive advanced vibrational molecular spectroscopy (ATR-FT/IR) in terms of chemical structure and biofunctional groups of amides (I and II), α-helix and β-sheet.
The results showed that increasing the level of the co-products in BPP significantly increased the spectral intensity of the amide area and amide height. The products exhibited similar protein secondary α-helix to β-sheet ratio. The protein molecular structure profiles (amides I and II, α-helix to β-sheet) were highly associated with protein degradation kinetics and intestinal digestion. In conclusion, the non-invasive vibrational molecular spectroscopy (ATR-FT/IR) could be used to detect inherent structural make-up characteristics in BPP.
The molecular structural features related to protein biopolymer were highly associated with protein utilization and metabolism.
Due to the high worldwide demands of oils and fuel in industries, bio-energy processing (bio-fuel, bio-oil, and bio-ethanol) resulted in huge amounts of co-products such as canola meal, new carinata meal, and distiller’s dried grains with solubles [1,2,3]. There are some studies that have investigated the application of the bio-coproduct of canola meal in ruminant or monogastric food/feed research [4,5,6]. Nevertheless, to our knowledge, there is limited information that could be found in the literature on the new bio-coproduct of carinata meal from bio-fuel processing when it is blended with other foodstuff such as peas screenings to optimize physicochemically or the biopolymer functions as well as nutritive value.
“Wet” chemistry analysis is usually used for feed evaluation, however, this technique usually damages the inherent structure of food/feed samples . The biofunctions and nutritive value, as well as fermentation features, have been reported to be influenced by the different inherent molecular structure make-up and conformation . The mid-infrared spectral region (ca. 4000–800 cm−1) has strong characteristic vibrational transitions compared with near-infrared region, especially in the wavelength range between ca. 1800 and 800 cm−1, which is called the “fingerprint region” [8, 9].
Vibrational spectroscopies such as Fourier transform infrared spectroscopy with attenuated total reflectance (ATR-FT/IR) is capable to detect the molecular structure of biomaterials. The ATR-FT/IR spectroscopy is a direct, rapid, non-destructive, and non-invasive bioanalytical technique used to reveal the infrared spectrum of absorptions or emissions of liquid, gas, or solids . The ATR-FT/IR spectroscopy has advantages such as revealing the molecules structural changes of different types of food/feed and determining the nutrient utilization and bioavailability . Moreover, this technique could recognize the molecular structure of different food crop varieties, food/feed ingredients, and studying the effects of food/feed processing on protein and carbohydrate-related molecular structures [4, 12,13,14]. However, there is no systematic study that has been conducted to determine how the blending process induces changes in protein intrinsic molecular structures and how these inherent structure changes influence protein metabolism and utilization. Therefore, the current study was performed to 1) investigate the differences among eight blend pellet products (BPP) from the bio-energy processing (new carinata meal vs. canola meal) with different proportions of pea screenings and lignosulfonate compound in terms of protein molecular structure; and 2) to quantify the protein inherent molecular structure changes in relation to protein profile, Cornell Net Carbohydrate, and Protein System (CNCPS) protein sub-fractions, energy values, protein digestion (rumen and intestine), and the metabolizable protein supply. The hypothesis of this study was that vibrational Fourier transform infrared spectroscopy could be used to detect an interactive association of processing induced molecule structural changes in biofunctional groups of protein amides I and II, alpha-helix and beta-sheet in the BPPs with protein metabolism and utilization.
All experimental procedures used in this study were approved by the University of Saskatchewan Animal Care Committee (UCACS Protocol No. 19910012) and were conducted in accordance with the Canadian Council of Animal Care guidelines .
Sample preparation for chemical and molecular structure studies
The experiment was performed at the SRP biomolecular spectroscopy lab, University of Saskatchewan (Saskatoon, Canada). The co-products of canola meal and new carinata meal from bio-fuel and bio-oil processing were used as BPPs by adding different levels of pulse co-products (pea screenings) and lignosulfonate compound. Eight blends were formulated; the BPP1 to BPP4 were based on new co-products of carinata meal (Agrisoma; Saskatoon, Canada) with different levels of pulse pea screenings and lignosulfonate; and the BPP5 to BPP8 were based on processing co-product of canola meal (Cargill Animal Nutrition, Clavet, Canada) with different levels of lignosulfonate and pulse pea screenings. The composition of the BPPs (on DM basis) is as follow: BPP1: lignosulfonate 0% + carinata meal 50% + pea screenings 50.0% DM. BPP2: lignosulfonate 4.8% + carinata meal 47.6% + pea screenings 47.6% DM. BPP3: lignosulfonate 0% + carinata meal 75% + pea screenings 25% DM. BPP4: lignosulfonate 4.8% + carinata meal 71.4% + pea screenings 23.8% DM. BPP5: lignosulfonate 0% + canola meal 50% + pea screenings 50.0% DM. BPP6: lignosulfonate 4.8% + canola meal 47.6% + pea screenings 47.6% DM. BPP7: lignosulfonate 0% + canola meal 75% + pea screenings 25% DM. BPP8: lignosulfonate 4.8% + canola meal 71.4% + pea screenings 23.8% DM. The use of different of the proportion of canola meal and pea screenings were made to obtain the best of amino acid profile and nutritional value. Adding lignosulfonate at 5%DM was based on a previous study  that showed a beneficial effect on canola meal by increasing the ruminal undegradable and the lactational performance of high producing dairy cows.
Pulse pea screenings were sourced from ILTA Grain Company (Surrey, Canada), while the lignosulfonate was obtained from Ameri-bond (Canada). The BPP was processed at two different times to make two different batches (n = 2) for each BPP. The blending and pelleting were conducted at Canadian Feed and Research Centre (CFRC, North Battleford, the University of Saskatchewan, Canada). For the pellet processing, the following procedure was followed to obtain the BPP: 1) mixing the combinations in the Scott Equipment model TSM 363 (New Prague, MN, USA) for 2 min.; 2) heating the different combinations by using Colorado Mill Equipment ECO-R30 (Cañon City, USA) at 65 °C and pelleting through a 3.6-mm diameter die such that the residence time of the blends in the die did not overtake 15 s; and 3) cooling at room temperature.
Detection of blending and pelleting impact on biofunctions and nutritive value using wet chemistry and biological techniques
The detailed chemical compositions, nutrient profiles, degradation kinetic profile, intestinal digestion, and the true nutrients supply, as well as metabolizable protein and feed milk value, were reported previously by Guevara-Ouendo et al. . The chemical profile, CNCPS fractions, and energy values of the BPPs (n = 16; carinata meal or canola meal with different combinations of pulse peas and lignosulfonate) are summarized in Table 1. These chemical and nutritive data were used for the correlation and regression studies on the association between the molecular structures and nutrition.
The CP was analyzed according to AOAC official method 984.13 . The protein subfractions i.e. neutral detergent insoluble crude protein (NDICP), non–protein nitrogen (NPN), and acid detergent insoluble crude protein (ADICP) were determined using the methods described by . The soluble crude protein (SCP) was analyzed by incubating samples with bicarbonate-phosphate buffer then filtrating through Whatman filter paper . For energy profiles, total digestible nutrient (TDN), metabolizable energy (ME), digestible energy (DE), and net energy (NE) were used for estimating the available energy in BPP. The truly digestible crude protein (tdCP), total digestible nutrient at a maintenance level (TDN1×), digestible energy at a production level of intake (DE3×), metabolizable energy at a production level of intake (ME3×), net energy for lactation at a production level (NEL3×) were determined by using a summative approach of the NRC .
The in situ degradation kinetics and the intestinal digestion of CP were performed according to . Degradation characteristics of CP were determined by applying the first-order kinetic model described by . The results were calculated using the NLIN procedure of SAS 9.4 and iterative least-squares regression (Gausse Newton method):
where U is the undegradable fraction (%); D is the potentially degradable fraction (%); Kd is the degradation rate (%/h), and T0 is the lag time (h).
The bypass crude protein (BCP) was determined according to the NRC model:
The truly nutrient supply was estimated using the NRC model. In this model, the MP (g/kg DM) is calculated based on the following equation :
where ME is the metabolizable protein; AMCP is the absorbable microbial protein; ARUP is the truly absorbable rumen undegraded feed protein, and AECP is the truly absorbable endogenous protein in the small intestine.
The milk value (FMV) was predicted based on the metabolizable protein content of BPP, where the efficiency of utilizing of metabolizable protein is assumed to be 0.67 and protein composition in milk is assumed to be 33 g protein/kg of milk .
Detection of blending and pelleting impact on protein molecular structure changes and collecting molecular spectra related to the protein primary and secondary structural components
For the molecular spectral analysis, the samples were grounded to pass a 0.12-mm sieve (Retsch ZM200, Rose, Scientific Ltd., Canada) for ATR-FT/IR spectroscopic analysis. Five subsamples from every sample were spectroscopically scanned. The molecular spectral data of samples were collected and corrected for the background spectrum using ATR-FT/IR molecular spectroscopy (JASCO 4200, JASCO International Co. Ltd., Tokyo, Japan). The spectra were generated in the mid-IR (ca. 4000–800 cm−1; Additional file 1: Figure S1) and the fingerprint region (ca. 1800–800 cm−1) with a spectral resolution of 4 cm−1. The ATR- FT/IR spectral was processed by using OMNIC 7.3 (Spectra-Tech, Madison, WI). The regions of specific interest in this study included the primary molecular protein structural (amide I and amide II) and the secondary molecular protein structural (α-helix and β-sheet) in the mid-IR. The structural spectral features on the protein were determined by analyzing the absorption peak parameter such as spectral region, baseline, peak, height, and area according to the published methods .
The univariate spectral analyses of protein structure comprised the primary and the secondary protein structures. The primary protein structures included amide I and II. The baseline of protein spectral was centered at ca. 1480–1730 cm−1 (Fig. 1). The baseline of the amide I area was centered at ca. 1569–1730 cm−1. The baseline of the amide II area was centered at ca. 1480–1569 cm−1. The peak height of the amide I was centered at ca. 1638–1649 cm−1, while the peak height of amide II was centered at ca. 1533–1540 cm−1. The secondary protein structures of the amide I region were determined by using the 2nd Derivative Function and Fourier Self–Deconvolution function in OMNIC 7.4 Software (Spectra Tech, Madison, WI) according to published methods [24, 25]. The secondary protein structures mainly comprised α-helix and β-sheet. The peak height of α-helix was centered at ca. 1647–1653 cm−1, while the peak height of β-sheet was centered at ca. 1625–1631 cm−1.
Principal Component Analysis (PCA) was performed using the Unscrambler 10.3 (CAMO Software AS, Oslo, Norway) for clustering any variation among BPP products. The raw data were prepossessed using baseline offset. The Savitzky-Golay algorithm was used to calculate the second derivative. Principal components with high eigenvalues were selected for further analysis. Two-dimensional score plots were used to display the PCA scores among the treatments. Loading plots for the important principal components were used to display the relations among PC components with IR variable data.
The data of functional groups in the protein region (ca. 1480–1730 cm−1) were analyzed by SAS 9.4 (SAS Institute, Inc., Cary, NC, USA). The experiment was designed using the randomized complete block design (RCBD) with blending/pellet processing time run as a random block effect.
The correlation between the biofunctional groups related to protein region (amide I, II peak highest and areas, α-helix, β-sheet and their ratio) and the chemical profiles of protein, energy values, degradation kinetics parameters, intestinal digestive characteristics of protein, and the truly absorbed protein supply was analyzed by using the PROC CORR procedure in SAS 9.4 (SAS Institute, Inc., Cary, NC, USA). Rank correlation with the SPEARMAN option and normality test with the UNIVARIATE option was used in the correlation study.
Multiple regression analysis (with model variable selection method) was performed to select the best biofunctional groups that would explain the nutritive values of BPP using the PROC REG procedure of SAS with a reversed stepwise option. The following model where be used for the multiple regression with model variable selection: model Y = spectral parameter 1 + spectral parameter 2 + spectral parameter 3 + spectral parameter 4 + ... + error. The model used a “STEPWISE” option with variable selection criteria: “SLENTRY = 0.05, SLSTAY = 0.05”. All variables left in the final prediction models were significant at the 0.05 level. Residual analysis was performed using the Univariate procedure of SAS with Normal and Plot options. Collinearity detection was conducted using the VIF option of SAS to eliminate the influence of correlated dependent variables to the prediction of independent variables.
Effect of blending/pelleting on protein molecular spectral intensities in blend pelleted products
The protein molecular spectral profiles include the primary and secondary structures of protein of different BPPs by using ATR-FT/IR vibrational spectroscopy is shown in Table 2. The results showed that BPP7 and BPP3 had the highest (P < 0.01) amide I peak height (averaged 0.307 IU), while the BPP2 and BPP6 had the lowest values (averaged 0.259 IU). The amide I area was significantly different (P < 0.01) among BPPs, where the BPP7 and BPP3 had the highest amide I area (averaged 21.2 IU), whereas the BPP2 and BPP6 had the lowest values.
Our results showed that the ratio of amide I to amide II was higher (P < 0.05) in BPP based on the new co-product of carinata meal compared with co-product of canola meal (2.12 vs. 1.98). Furthermore, adding the pulse pea screenings decreased (P < 0.05) the ratio of amide I to amide II in BPP based on the new co-product of carinata meal from 2.19 to 2.05, and in BPP based on the co-product of canola from 1.95 to 2.02.
The secondary structures such as α-helix and β-sheet and their ratio of BBPs are presented in Table 2. It has been found that the α-helix, β-sheet height ratio was the same for all treatments. This study showed that the ratio of α-helix to β-sheet decreased (P < 0.05) from 1.17 to 1.13 by decreasing the level of co-products in the BPPs.
Correlation analysis between protein molecular spectral features and nutrition profiles in the blend pelleted products
The correlation analysis between the vibrational spectral features and protein profiles, protein subfractions and the energy values of BPPs are shown in Table 3. The CP had positive correlations with amide I area (r = 0.70, P < 0.05), total amide area (r = 0.58, P = 0.02), and amide I height (r = 0.55, P = 0.03) in BPPs. However, CP exhibited a negative correlation with the ratio of α-helix to β-sheet (r = − 0.55, P = 0.03).
The results in the current study showed NDICP had a negative correlation with amide II height (r = − 0.50, P = 0.05) and a positive correlation with the ratio of amide I to amide II height (r = 0.52; P = 0.04). The content of ADICP was found to be positively correlated with the amide II height (r = 0.56, P = 0.02) or the amide II area (r = 0.51, P = 0.04). The concentration of NPN was negatively correlated with the amide I area (r = − 0.80, P = 0.01), amide I height (r = − 0.71, P < 0.05), amide I to amide II area ratio (r = − 0.68, P < 0.05), amide I to amide II height ratio (r = − 0.56, P < 0.05), β-sheet height (r = − 0.67, P < 0.05), and α-helix height (r = − 0.61, P < 0.05). However, there was a positive correlation between NPN and the ratio of α-helix to β-sheet (r = 0.64, P = 0.01).
For the truly digestible crude protein, the results showed correlations between tdCP and amide I area (r = 0.68, P = 0.01), amide I height (r = 0.53, P < 0.05), the amide area (r = 0.56, P = 0.02). There was no correlation (P > 0.05) between the energy values and the molecular structure features related to protein region.
For protein sub fractions partitioned by the CNCPS model, the results showed that the slowly degradable protein (PB2 fraction) was positively associated with the amide I to amide II height ratio (r = 0.50, P = 0.05) and negatively related to the amide II height (r = − 0.52, P = 0.04). The PC fraction was observed to be positively associated with the amide II height (r = 0.56, P = 0.02) and amide II area (r = 0.51, P = 0.04).
In the current study, it has been found that the protein degradation rate was negatively correlated with the amide I area (r = − 0.36, P < 0.01) and the amide I to amide II height ratio (r = − 0.36, P = 0.04), while the slowly degradable fraction of protein was positively correlated with the amide I height, amide I area and amide area. The amide II area and the ratio of amide I to amide II height were found to be negatively related to the EDCP of BPPs (P < 0.05; Table 4). Additionally, the undegradable CP fraction of BPPs was found to be negatively correlated with the heights or areas of amide I and amide II. The in vitro digestion of BCP (% dIDP) and the total intestinal digestibility of CP (IADP %) have been found to be positively correlated with the amide I to amide II height ratio (r = 0.55, P < 0.05). The results showed that the α-helix to β-sheet ratio was correlated with the slowly degradable fraction of CP (r = − 0.44, P = 0.01) and the undegradable fraction of CP (r = 0.45, P = 0.01).
The correlation between protein molecular structure and the truly absorbable protein supply of BPP is shown in Table 5. The data showed significant correlation between AMCP and the amide II height (r = − 0.53, P < 0.01), the amide II area (r = − 0.50, P < 0.01), the amide area (r = − 0.43, P < 0.01), and the helix to β-sheet ratio (r = 0.55, P < 0.01). For the truly absorbable rumen undegraded protein in the small intestine, ARUP has exhibited a positive correlation with the amide I to amide II height ratio (r = 0.62, P < 0.01). The MP had a positive correlation with amide I to amide II height ratio (r = 0.61, P < 0.01) and a negative correlation with the amide II height (r = − 0.43, P = 0.02). The DPB had positive correlations with amide II height (r = 0.70, P < 0.01), the total amide area (r = 0.65, P < 0.01), amide II area (r = 0.64, P < 0.01), amide I height (r = 0.63, P < 0.01), amide I area (r = 0.57, P < 0.01), β-sheet height (r = 0.79, P < 0.01), and α-helix height (r = 0.77, P < 0.01) but a negative correlation with α-helix to β-sheet ratio (r = − 0.62, P < 0.01). The FMV had a positive correlation with amide I to amide II height ratio (r = 0.45, P = 0.01), while negatively correlated with the amide II height (r = − 0.39, P = 0.03).
Model variable selection analysis to choose the most important protein spectral parameters to predict the protein nutrient profile and protein utilization and metabolism
The multiple regressions analysis is shown in Table 6. The equations of protein profiles showed that CP could be predicted from amide I area and α-helix height, taking 79% of the total variance. The SCP could be predicted from amide I to amide II ratio area and α-helix height and 68% of the total variance was taken by it. Amide I to amide II ratio could also predict NDICP and NPN. The amide II height could be used to predict the ADICP, PB2, and PC. For truly digestible nutrients, tdCP could be predicted from amide I area and α-helix height, while the Amide I height and β-sheet height could predict TDN3× with a total variance of 82%.
Table 7 shows that the amide I area and the α-helix height could be used to estimate the Kd and the undegradable fraction of CP with 46% and 62% of the total variance, respectively. The amide I height and β-sheet height were the best spectral variables to predict the EDCP with 67% of the variance. The amide I, II height ratio, Amide area, β-sheet height and α-helix to β-sheet ratio could be used to predict the intestinal digestion of CP with 90% of the variance. The results in Table 8 shows that the amide I to amide II height ratio was the best spectral feature in estimating the ARUP, MP, and FMV of BPPs, while the amide I height and β-sheet height would be used to estimate the DBP and AMCP.
Recently, the advanced vibrational spectroscopic techniques have been established to quantitatively estimate the primary and secondary molecular make-up of protein [26, 27]. Generally, the amide I and amide II bands are used to detect the information about protein concentration [28, 29]. However, the amide I is used more frequently than amide II to reveal the molecular structure of the protein, since the amide II originates from complex of vibrations that includes numerous functional groups such as ligneous compounds . Our results showed that amide I area and peak height are the highest in the BPP7and BPP3 and lowest in the BPP2 and BPP6. These results are in agreement with the results obtained from the wet chemistry analysis . The high CP or amide I area values are attributed to the high inclusion level of co-product (canola or carinata) in those BPPs. The total amide area was highly sensitive to the changes in blend BPPs composition, where the amide area increased with increasing the co-products levels or with decreasing pulse pea screenings levels in BPPs. For example, the BPP3 or BPP7 (carinata or canola meal 75% DM + pea screenings 25% DM) has a higher (P < 0.05) total amide area than BPP1 or BPP5 (carinata or canola meal 50% DM + pea screenings 50.0% DM).
The ratio of amide I to amide II has been found to be influenced by the heat-related processing of food/feed or by gene transformation of plant forage [32,33,34]. Previous studies noted that the amide I and II ratio had a positive correlation with the metabolizable protein . Based on the current study results, the high amide I to amide II ratio of BPPs could be a consequence of the high inclusion level of co-products or adding co-product of carinata meal to BPPs. In agreement with these findings, Guevara-Oquendo et al.  reported a higher indigestible protein content in the BPP based on the co-product of canola meal (1.5% CP) than that of BPP based on the new co-product of carinata meal (3.2% CP). The low PC in BPP based on the new co-product of carinata meal is attributed to the greater content of NDICP and a lower content ADICP in the co-product of carinata meal compared with the co-product of canola meal. The NDICP is slowly degraded in the rumen and largely contributes to escaping protein from ruminal degradation . Thus, a large amount of NDICP could reach the small intestine and hence, increasing the MP supply to dairy cows . On the other hand, the ADICP reflects the amount of protein that is completely indigestible in the gastrointestinal tract . Therefore, increasing the concentration of ADICP could limit the total tract digestibility of protein.
The secondary structures such as α-helix and β-sheet and their ratio are commonly used to detect the information about the protein’s molecular makeup [28, 36]. In the current study, all BPPs underwent the same processing, thus it is not surprising that the ratio was not changed among all BPPs. Yu and Samadi  found that the ratio of α-helix to β-sheet was altered by the moist heating of soybean and canola seeds. The alteration in the ratio of α-helix to β-sheet by thermal treatment was also reported by . The changes in the secondary structure of the protein are possibly related to the denaturation of α-helix and β-sheet during thermal treatment. The current results showed that the ratio of α-helix to β-sheet was decreased with decreasing the level of co-products in the BPPs, which would reflect a reduction in the MP supply in BPPs. In agreement with these results, Guevara et al.  found that MP supply has been reduced by decreasing the inclusion level of the co-products of carinata or canola meal in BPPs.
The PCA analysis in the current study was used to reduce the number of variables. The PCA was performed on the molecular structure related to protein region (ca. 1480–1730 cm−1). The first two PCs derived from the PCA classification of these spectra described 94% of the variance in the BPPs (Fig. 2a, b). Most of the BPPs based on the co-product of canola meal such as BPP6, BPP7 and BPP8 were clearly separated from the BPPs based on the new co-product of carinata meal by the PC2 which accounted for 5% of the variance. The BPP1 had exhibited the least negative values in PC2, while the BPP7 and BPP8 had the highest positive values. The PC1 which account for 89% of the variations among BPPs in terms of the molecular structure features did not cluster most of the BPPs. The overlapping between BPPs in the PC1 would indicate that these pellets had similar molecular structure features in the amide region. The loading point plots were used to determine the most important regions responsible for the clustering (Fig. 4). The amide I peak at ca. 1650 cm−1 was heavily loaded in PC1 and PC2, which separated the negative scores of spectra that belong to BPPs based on the new co-product of carinata meal from the positive score of the spectra that related to the BPP6, BPP7, and BPP8 (Fig. 3). These findings indicate that the amide region at ca. 1650 cm−1 of PC2 was the most important parameters for discriminating the BPPs. These data demonstrated that the amide I peak at ca. 1650 cm−1 for BPPs based on the new co-product of carinata meal was lower than that of the BPP6, BPP7, and BPP8. These data are in agreement with the univariate analysis (Table 2) that showed BPPs based on the co-product of canola meal were significantly higher in the amide I peak height compared with BPPs based on the new co-product of carinata meal. Based on these findings, the amide I band which is sensitive to small differentiation in molecular structure and hydrogen bonding motifs is important in the determination of protein structural and conformational changes.
In order to obtain clear and precise peak positions of protein bands by FT/IR-ATR, the raw spectra were processed by taking the second derivative (Fig. 4 a,b), which gives a negative peak for each band and shoulder in the absorption spectrum, and hence allows us to identify the individual peaks among complex spectra. The PCA score plot demonstrated that the clusters of all BPPs were overlapped along PC1 (51%) and PC2 (17%). The PCA loading plots of PC1 and PC2 are shown in Fig. 5. The loading plot showed that the variations along PC2 could be explained by the positive loading in the Amide II region (centered at ca. 1548 cm−1; N-H (60%) bending and C-N (40%) stretching vibrations: proteins α-helix), which separated the negative score of some BPPs based on the co-product of canola meal from the positive score of BPPs based on the new co-product of carinata meal.
The correlation analysis between the vibrational spectral features and protein profiles, protein subfractions and the predicted energy values of BPPs. Our results for the correlation between CP and primary structure and secondary structure are in agreement with previous studies that reported positive correlations of CP with amide I area and amide I height of the co-product of carinata meal or canola meal samples . Furthermore, the correlation between NDICP and primary structure are similar to the previous studies that showed NDICP had a negative correlation with the amide II height and a positive correlation with the ratio of amide I to amide II height . Our results showed there was no correlation (P > 0.05) between the metabolic and net energy values by the NRC-model and the molecular structure features related to protein region. In line with findings,  did not detect any association between the energy values and the molecular structure characteristics of the protein.
For the CNCPS fractions, the current study results are in agreement with  who reported a positive correlation between PB2 subfraction with the amide I to amide II height ratio. However, there was no association between the CNCPS fractions and α-helix to β-sheet ratio. In agreement with observations, Huang et al.  reported no correlation between α-helix to β-sheet ratio and the protein subfractions estimated by the CNCPS model.
For the in situ degradation kinetic parameters, a previous study noted that the ratio of amide I to amide II was highly correlated with the in situ protein degradation kinetic parameters . These associations have been found to be affected by the enzymatic digestion of protein . It has been observed the changes in the ratio of α-helix to β-sheet ratio would induce alterations in the protein molecular makeup . The thermal treatment has been found to decrease the solubility of protein and increase the ADICP and NDICP as a consequence of protein denaturation during the heating process. Furthermore, the heat treatment could increase the cross-linkages among the amino acids in the polypeptide chain and reduce the sugars and finally decrease the solubility of CP . Previous studies found that applying the heat treatment and increasing the heating time, caused an increase in the α-helix: β-sheet ratio in flaxseed and bio-ethanol co-products, respectively [34, 36]. In the current study, the α-helix to β-sheet ratio was correlated with the slowly degradable fraction of CP (r = − 0.44, P = 0.01) and the undegradable fraction of CP (r = 0.45, P = 0.01). These findings are in agreement with a previous study that found strong correlations between the α-helix to the β-sheet ratio of camelina seeds and the in situ protein degradation parameters . The correlation coefficient values in this study were lower than that reported by Khan et al.  due to the diversity in protein origin in BPPs which applied adverse effects on the accuracy of predictions.
In conclusion, the results in the current study indicated that the molecular structure features related to the protein of the blend pelleted products based on the co-product of canola or carinata meal could be revealed by the ATR-FT/IR spectroscopy. The univariate analysis showed differences in the absorption of the functional groups in the intestines related to the primary structure of the protein. The secondary structure of protein i.e. α-helix, β-sheet height ratio did not affect by BPP because all ingredients underwent to the same processing condition. The amide I to amide II height ratio was the best spectral parameter to estimate the changes in the protein degradation and the metabolizable protein of BPP.
Availability of data and materials
The datasets during and/or analyzed during the current study available from the corresponding author on reasonable request.
Two-dimensional score plots
Acid detergent insoluble crude protein
Truly absorbable endogenous protein in the small intestine
Absorbable microbial protein
Truly absorbable rumen undegraded feed protein
Fourier transform infrared spectroscopy with attenuated total reflectance
Bypass crude protein
Blend pelleted products
Cornell Net Carbohydrate and Protein System
Potentially degradable fraction (%)
- DE3× :
Digestible energy at the production level of intake
Feed milk value
Degradation rate (%/h)
Passage rate (%/h)
- ME3× :
Metabolizable energy at the production level of intake
Neutral detergent insoluble crude protein
- NEL3× :
Net energy for lactation at the production level
Principal Component Analysis
Residue present at t h incubation (%)
Soluble crude protein
Lag time (h)
Truly digestible crude protein
Total digestible nutrient
- TDN1× :
Total digestible nutrient at the maintenance level
Undegradable fraction (%)
Canola Council of Canada. Canola Meal Feed Industry Guide, 5th edition. Publication of Canola Council of Canada, Winnipeg, . 2015. https://www.canolacouncil.org/media/516716/2015_canola_meal_feed_industry_guide.pdf. (Accessed 30 June 2019).
Ban Y, Yu P. Structural and nutritional characterization of macromolecular complexes in new bioenergy feedstock by infrared radiation with advanced molecular spectroscopy and spectral chemometrics. Appl Spectrosc Rev. 2016;51:822–38. https://doi.org/10.1080/05704928.2016.1191017.
Xin H, Yu P. Chemical profile, energy values, and protein molecular structure characteristics of biofuel/bio-oil co-products (Carinata meal) in comparison with canola meal. J Agric Food Chem. 2013;61:3926–33. https://doi.org/10.1021/jf400028n.
Huang X, Zhang H, Yu P, Miller S, Sweat J, Yu P. Structural changes on a molecular basis of canola meal by conditioning temperature and time during pelleting process in relation to physiochemical (energy and protein) properties relevant to ruminants. PLoS One. 2017. https://doi.org/10.1371/journal.pone.0170173.
Theodoridou K, Yu P. Effect of processing conditions on the nutritive value of canola meal and presscake. Comparison of the yellow and brown-seeded canola meal with the brown-seeded canola presscake. J Sci Food Agric. 2013;93:1986–95. https://doi.org/10.1002/jsfa.6004.
Heendeniya RG, Christensen DA, Maenz DD, McKinnon JJ, Yu P. Protein fractionation byproduct from canola meal for dairy cattle. J Dairy Sci. 2012;95:4488–44500. https://doi.org/10.3168/jds./jds.2011-5029.
Yu P, Gamage IH, Zhang X. New approaches and recent advances on characterization of chemical functional groups and structures, physiochemical property, and nutritional values in feedstocks and by-products: advanced Spectro analytical and modeling investigations. Appl Spectrosc Rev. 2014;49:585–602. https://doi.org/10.1080/05704928.2013.879064.
Yu P. Application of advanced synchrotron radiation-based Fourier transform infrared (SR-FTIR) microspectroscopy to animal nutrition and feed science: a novel approach. Br J Nutr. 2004;92:869–85. https://doi.org/10.1079/BJN20041298.
Liu N. Ruminal nutrient availability and inherent structural features of six barley varieties using in situ technique and mid-IR spectroscopy. MSc thesis: University of Saskatchewan; 2009.
Smith BC. Fundamentals of Fourier transform infrared spectroscopy. Baseline. Boca Raton, FL, USA.: CRC Press; 2009.
Reynolds CK, Harmon DL, Cecava MJ. Absorption and delivery of nutrients for milk protein synthesis by portal-drained viscera. J Dairy Sci. 1994;77:2787–808. https://doi.org/10.3168/jds./jds.S0022-0302(94)77220-9.
Abeysekara S, Christensen DA, Yu P. Characterizations of structural, biochemical, and nutritive profiles in silage among cool-season corn cultivars in relation to heat units (aCHU, dCHU) with curvilinear response and multivariate analyses. J Agric Food Chem. 2013;61:12315–26. https://doi.org/10.1021/jf403245c.
Khan NA, Peng Q, Xin H, Yu P. Vibrational spectroscopic investigation of heat-induced changes in functional groups related to protein structural conformation in camelina seeds and their relationship to digestion in dairy cows. Anim Prod Sci. 2015;55:201–6. https://doi.org/10.1071/AN14400.
Xin H, Yu P. Using ATR-FT/iR to detect carbohydrate-related molecular structure features of carinata meal and their in situ residues of ruminal fermentation in comparison with canola meal. Spectrochim Acta - Part A Mol Biomol Spectrosc. 2013;114:599–606. https://doi.org/10.1016/j.saa.2013.05.056.
Canadian Council on Animal Care. Guidelines on: the care and use of farm animals in research, teaching and testing. Ottawa: ON; 2009.
Wright CF, von Keyserlingk MAG, Swift ML, Fisher LJ, Shelford JA, Dinn NE. Heat- and lignosulfonate-treated canola meal as a source of ruminal undegradable protein for lactating dairy cows. J Dairy Sci. 2010; 88:238–243. doihttps://doi.org/10.3168/jds./jds.S0022-0302(05)72681-3.
Guevara-Oquendo VH, Christensen DA, McKinnon JJ, Rahman MM, Yu P. Potential nitrogen to energy synchronization, rumen degradation kinetics, and intestinal digestibility of blend pelleted products of new co-products from bio-fuel processing, pulse screenings and lignosulfonate compound in dairy cows. Anim Feed Sci Technol. 2018;236:196–207. https://doi.org/10.1016/j.anifeedsci.2017.11.011.
AOAC. Official Methods of Analysis of AOAC International. Assoc Off Anal Chem Int. 1990. p. 1058–1059.
Licitra G, Hernandez TM, Van Soest PJ. Standardization of procedures for nitrogen fractionation of ruminant feeds. Anim Feed Sci Technol. 1996;57:347–58. https://doi.org/10.1016/0377-8401(95)00837-3.
Roe MB, Sniffen CJ, Chase LE, Cornell U. Techniques for measuring protein fractions in feedstuffs. Proc. 1990 Cornell Nutr. Conf. Feed Manuf. Proc. Cornell Nutr. Conf. Ithaca, NY; 1990.
NRC. Nutrient requirement of dairy cattle. 7th rev. National Research Council, Washington, DC: National Academy Press; 2001.
Yu P, Meier JA, Christensen DA, Rossnagel BG, McKinnon JJ. Using the NRC-2001 model and the DVE/OEB system to evaluate nutritive values of Harrington (malting-type) and Valier (feed-type) barley for ruminants. Anim Feed Sci Technol. 2003;107:45–60. https://doi.org/10.1016/S0377-8401(03)00062-2.
Ørskov ER, Mcdonald I. The estimation of protein degradability in the rumen from incubation measurements weighted according to rate of passage. J Agric Sci. 1979;92:499–503. https://doi.org/10.1017/S0021859600063048.
Yu P. Protein secondary structures (alpha-helix and beta-sheet) at a cellular level and protein fractions in relation to rumen degradation behaviours of protein: a new approach. Br J Nutr. 2005;94:655–65. https://doi.org/10.1079/BJN20051532.
Theodoridou K, Yu P. Application potential of ATR-FT/IR molecular spectroscopy in animal nutrition: revelation of protein molecular structures of canola meal and presscake, as affected by heat-processing methods, in relationship with their protein digestive behavior and utilization for dairy cattle. J Agric Food Chem. 2013;61:5449–58. https://doi.org/10.1021/jf400301y.
Khan NA, Booker H, Yu P. Molecular structures and metabolic characteristics of protein in brown and yellow flaxseed with altered nutrient traits. J Agric Food Chem. 2014;62:6556–64. https://doi.org/10.1021/jf501284a.
Yu P. Molecular chemical structure of barley proteins revealed by ultra-spatially resolved synchrotron light sourced FTIR microspectroscopy: comparison of barley varieties. Biopolymers. 2007;85:308–17. https://doi.org/10.1002/bip.20661.
Damiran D, Yu P. Molecular basis of structural makeup of hulless barley in relation to rumen degradation kinetics and intestinal availability in dairy cattle: a novel approach. J Dairy Sci. 2011;94:5151–9. https://doi.org/10.3168/jds./jds.2011-4564.
Peng Q, Khan NA, Wang Z, Yu P. Moist and dry heating-induced changes in protein molecular structure, protein subfractions, and nutrient profiles in camelina seeds. J Dairy Sci. 2014;97:446–57. https://doi.org/10.3168/jds./jds.2013-7298.
Jackson M, Mantsch HH. The use and misuse of FTIR spectroscopy in the determination of protein structure. Crit Rev Biochem Mol Biol. 1995;30:95–120. https://doi.org/10.3109/10409239509085140.
Liu B, McKinnon JJ, Thacker P, Yu P. Molecular structure and metabolic characteristics of the proteins and energy in triticale grains and dried distillers grains with solubles for dairy cattle. J Agric food Chem. Am Chem Soc. 2012;60:10064–1074. https://doi.org/10.1021/jf302382b.
Yu P. Synchrotron IR microspectroscopy for protein structure analysis: Potential and questions. Spectroscopy. 2006;20:229–51. https://doi.org/10.1155/2006/263634.
Doiron KJ, Yu P, Christensen CR, Christensen DA, McKinnon JJ. Detecting molecular changes in Vimy flaxseed protein structure using synchrotron FTIRM and DRIFT spectroscopic techniques: structural and biochemical characterization. Spectroscopy. 2009;23:307–22. https://doi.org/10.3233/SPE-2009-0393.
Russell JB, O’Connor JD, Fox DG, Van Soest PJ, Sniffen CJ. A net carbohydrate and protein system for evaluating cattle diets: I. Ruminal fermentation J Anim Sci. 1992;70:3551–61. https://doi.org/10.2527/1992.70113562x.
Yu P, Nuez-Ortín WG. Relationship of protein molecular structure to metabolizable proteins in different types of dried distillers grains with solubles: a novel approach. Br J Nutr. 2010;104:1429–37. https://doi.org/10.1017/S0007114510002539.
Samadi, Yu P. Dry and moist heating-induced changes in protein molecular structure, protein subfraction, and nutrient profiles in soybeans. J Dairy Sci. 2011;94:6092–102. https://doi.org/10.3168/jds.2011-4619.
Doiron KJ, Yu P. Recent research in flaxseed (oil seed) on molecular structure and metabolic characteristics of protein, heat processing-induced effect and nutrition with advanced synchrotron-based molecular techniques. Crit Rev Food Sci Nutr. 2017;57:8–17. https://doi.org/10.1080/10408398.2013.764513.
Ban Y, Prates LL, Yu P. Investigating molecular structures of bio-fuel and bio-oil seeds as predictors to estimate protein bioavailability for ruminants by advanced nondestructive vibrational molecular spectroscopy. J Agric Food Chem. 2017:9147–57. https://doi.org/10.1021/acs.jafc.7b02239.
Xin H, Yu P. Chemical profile, energy values, and protein molecular structure characteristics of biofuel/bio-oil co-products (Carinata meal) in comparison with canola meal. J Agric Food Chem. 2013;61:3926–33. https://doi.org/10.1021/jf400028n.
Yu P. Protein molecular structures, protein subfractions, and protein availability affected by heat processing: a review. Am J Biochem Biotechnol. 2007;3:70–90. https://doi.org/10.3844/ajbbsp.2007.66.86.
Licitra G, Hernandez TM, Van Soest PJ. Standardization of procedures for nitrogen fractionation of ruminant feeds. Anim Feed Sci Technol. 1996;57:347–58. https://doi.org/10.1016/0377-8401(95)00837-3.
Williams PC, Starkey PM. Influence of feed ingredients upon the prediction of protein in animal feed-mixes by near-infrared reflectance spectroscopy. J Sci Food Agric. Wiley Ltd 1980; 31:1201–1213. https://doi.org/10.1002/jsfa.2740311112.
The authors would like to thank Zhiyuan Niu (University of Saskatchewan), John Smillie and Rex Newkirk (Canadian Feed Research Centre- Feed Processing), Morgan Hobin (Rayner Research and Teaching Facility) and Lisette Mascarenhas and Constance Chiremba (SK Pulse Growers) for their assistance and/or technical support.
This research was financially supported by the grants from SaskPulse Growers, Natural Sciences and Engineering Research Council of Canada (NSERC), the SaskCanola, the Ministry of Agriculture Strategic Research Chair Program, and SaskMilk.
Ethics approval and consent to participate
All animal experimental procedures used in this study were approved by the University of Saskatchewan Animal Care Committee (UCACS Protocol No. 19910012) and were conducted in accordance with the Canadian Council of Animal Care guidelines (CCAC, 2009).
Consent for publication
The authors declare that they have no competing interest.
About this article
Cite this article
Ismael, A., Guevara-Oquendo, V.H., Refat, B. et al. Interactive association between processing induced molecular structure changes and nutrient delivery on a molecular basis, revealed by cutting-edge vibrational biomolecular spectroscopy. J Animal Sci Biotechnol 10, 85 (2019) doi:10.1186/s40104-019-0384-z
- Alpha-helix and beta-sheet
- Amides (I and II)
- Biofunctional groups
- Chemical structure
- Protein metabolism and bioavailability
- Protein molecular structure