Next Article in Journal
Isolation and Characterization of Galloylglucoses Effective against Multidrug-Resistant Strains of Escherichia coli and Klebsiella pneumoniae
Previous Article in Journal
The Potential of Lamiaceae Herbs for Mitigation of Overweight, Obesity, and Fatty Liver: Studies and Perspectives
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Terahertz Spectroscopic Analysis of Lactose in Infant Formula: Implications for Detection and Quantification

by
Sopanant Datta
1,
Kiattiwut Prasertsuk
2,
Nuttawat Khammata
3,
Patharakorn Rattanawan
2,
Jia Yi Chia
2,
Rungroj Jintamethasawat
2,
Thawatchart Chulapakorn
2,4 and
Taweetham Limpanuparb
1,*
1
Science Division, Mahidol University International College, Mahidol University, Salaya, Phutthamonthon, Nakhon Pathom 73170, Thailand
2
National Electronics and Computer Technology, National Science and Technology Development Agency, 112 Thailand Science Park, Khlong Luang 12120, Thailand
3
Department of Physics and Materials Science, Faculty of Science, Chiang Mai University, Chiang Mai 50200, Thailand
4
Department of Construction Sciences, Lund University, 22100 Lund, Sweden
*
Author to whom correspondence should be addressed.
Molecules 2022, 27(15), 5040; https://doi.org/10.3390/molecules27155040
Submission received: 7 July 2022 / Revised: 29 July 2022 / Accepted: 4 August 2022 / Published: 8 August 2022
(This article belongs to the Section Molecular Structure)

Abstract

:
Lactose plays a significant role in daily lives as a constituent of various food and pharmaceutical products. Yet, lactose intolerance conditions demand low-lactose and lactose-free products in the market. These increasing nutritional claims and labels on food products entail simple and reliable methods of analysis that can be used for meeting quality standards, nutritional claims and legal requirements. In this study, terahertz time–domain spectroscopy (THz-TDS) was employed to analyse α-lactose monohydrate qualitatively and quantitatively in food products. Both absorption spectra and absorption coefficient spectra were investigated for their prediction performance. Regression models for lactose quantification using peak area and height of the absorption peaks 0.53 and 1.37 THz were developed and assessed in infant formula samples. Satisfactory prediction results were achieved in ideal conditions with pure standards, but not in all predictions of infant formula samples. Reasons and further implications are discussed.

Graphical Abstract

1. Introduction

Amongst various populations, milk and milk products are considered as integral parts of the diet. Over 6 billion people worldwide consume dairy products [1]. Milk is a good source of protein, calcium, phosphorus, magnesium and other essential macro- and micronutrients. The disaccharide lactose is the principal sugar in milk, which can be broken down into its monomers glucose and galactose by the enzyme lactase in our gastrointestinal tract [2]. Lactose intolerance is a common condition in individuals with lactase deficiency, causing symptoms such as flatulence and diarrhoea [3]. This condition is generally treated with the avoidance of lactose-containing food products. To meet the dietary requirements of the lactose intolerant population, various low-lactose and lactose-free dairy products have been developed and sold in the market.
In addition to being found naturally in milk, lactose has important applications in the food industry. It is added to infant formulas to imitate the lactose content in human milk [2]. With its lack of flavour and low sweetness, it is used as a bulking agent that does not significantly add to the overall sweetness, a diluent or enhancer of flavour and a reducing sugar in the Maillard reaction to form flavour and colour compounds [4].
In the pharmaceutical industry, lactose is a widely used excipient of drugs. It is generally used as a filler and binder in oral solid dosage formulations, such as tablets and dry powder inhalers [5]. The most common form of lactose used as an excipient is crystalline α-lactose monohydrate [6].
The detection and quantification of lactose are vital in the food and pharmaceutical industries for product research and development, quality control, nutritional labelling and verification of nutritional claims. In these processes, it is necessary to handle appropriate equipment for the rapid, simple and real-time determination of lactose content in food and pharmaceutical products.
Commonly used and recognised analytical methods for lactose determination include spectrophotometric [AOAC 2006.06], spectroscopic [AOAC 972.16], polarimetric [AOAC 896.01], gravimetric [AOAC 930.28] and chromatographic methods. Conventional methods are modified and novel methods are developed to improve the sensitivity, selectivity and limit of detection [7]. These improvements can be achieved with chromatographic and spectroscopic methods (Table 1) such as the conventional high-performance liquid chromatography (HPLC) [8], hydrophilic interaction chromatography (HILIC) [9], high-performance thin layer chromatography (HPTLC) [10], high-performance anion exchange chromatography with pulsed amperometric detection (HPAEC-PAD) [7], near-infrared spectroscopy (NIR) [11] and nuclear magnetic resonance (NMR) [12]. In spite of the sensitivity, selectivity and accuracy achieved with these methods, they are limited by their high requirements for operation cost, time consumption and operator expertise [13]. These limitations can introduce difficulties in the routine analysis of samples to meet quality standards and labelling requirements.
Recently, terahertz (THz) spectroscopy has been demonstrated to have significant potential as a rapid, non-destructive and reliable method of analysis in the food and pharmaceutical industries. Terahertz waves, with a frequency range between the microwave and infrared regions (0.1–10 THz), are non-ionising, have high penetrability and can provide intermolecular information on vibrational activities of molecules. They can detect weak intermolecular interactions such as hydrogen bonds and van der Waals forces [14]. Terahertz time–domain spectroscopy (THz-TDS) has been used for the detection and identification of pesticides [15], microplastics [16], amino acids and sugars including lactose [11,17,18,19]. As this method measures crystal lattice vibrations, it detects crystalline forms of molecules [11,20,21]. Therefore, it can be used as a preliminary method of detection and quantification or used to measure the crystallinity of compounds in food and pharmaceutical samples.
In this study, THz-TDS is employed to detect and quantify α-lactose monohydrate in food samples. A quantitative model for lactose analysis is developed using two different well-known absorption peaks at 0.53 and 1.37 THz [22]. The resulting model is applied to dairy product samples to assess THz-TDS as a potential rapid and real-time method of the detection and quantification of compounds in food.

2. Results and Discussion

2.1. Spectra of Pure Samples

THz-TDS measurements were conducted with pure α-lactose monohydrate samples at concentrations of 0%, 1%, 3%, 5%, 10%, 15%, 20%, 49%, 80% and 100% (w/w). The spectra are in the range of 0.3–2.0 THz with a high signal-to-noise ratio. Measurements were conducted under ambient conditions (≈37% RH) and a nitrogen-purged atmosphere (≤7% RH) to explore the effects of moisture. Figure 1 shows the absorption spectra from measurements under ambient conditions and one of the replicate measurements in a nitrogen atmosphere. The absorption peaks observed are at the frequencies of 0.53, 1.20, 1.37 and 1.82 THz, in agreement with previous studies using THz-TDS [20,21,22,23,24,25] and density functional theory simulation [24]. The absorption peaks obtained from measurements in a nitrogen atmosphere and those obtained under ambient conditions are of the same frequency and magnitude, indicating that the lactose absorption peaks are independent of water absorption.
Amongst the absorption peaks, strong peaks are observed at 0.53 THz and 1.37 THz, with clear detection at concentrations of 10% (w/w) and higher. The absorbance value of an absorption peak varies with its corresponding lactose concentration. That is, absorption peaks in spectra of higher lactose concentrations are higher in absorbance. Nevertheless, the etalon feature or fringe pattern are present in the frequency domain spectra of all samples. This etalon feature arises as the transmitted THz radiation propagates through the air/sample and sample/air interfaces. This effect is more pronounced in pellets of a low thickness (≈1 mm), as used in this study [26]. The etalon pattern leads to difficulties in baseline removal, an increased limit of detection, and inaccuracies in the quantitative analysis.
Spectral artefacts induced by the etalon effect are not easily removed by ratioing to a reference spectrum. One of the techniques of etalon effect elimination is the use of a reference spectrum obtained from a combination of thicknesses and the refractive index of the sample [26]. Thus, absorption coefficient spectra are investigated in this study, as the spectra involve more parameters including the correction for sample thickness (see calculation methods in Section 3). The refractive index and absorption coefficient spectra of pure lactose samples are shown in Figure 2. In comparison to Figure 1, the etalon effect in the absorption coefficient spectra is considerably reduced, yet not totally eliminated. Both refractive indices and absorption coefficients increase with the increasing concentration of lactose. However, pronounced transitions can be observed for the lactose concentrations of 5% (w/w) and higher. In the absorption coefficient spectra, strong absorption peaks centred at 0.53, 1.20, 1.37 and 1.82 THz are observed. Dispersion regions at the same frequencies are observed in the refractive index spectra, whereby those at 1.82 THz are considerably less pronounced, which is in agreement with a previous study [25].

2.2. Quantification Model

The positive relationship between the α-lactose monohydrate concentration and peak area or peak height of the two major absorption peaks, 0.53 THz and 1.37 THz, were observed as shown in Figure 3. Peaks from the absorption spectra are considered in comparison to those from the absorption coefficient spectra. Gaussian fitting and baseline correction were performed to extract the peak area and peak height values from the spectra. The resulting models and predictions from different datasets (experimental replicates) and baseline correction methods can be found in the Supplementary Materials. A model of lactose quantification, using the absorption and absorption coefficient spectra as predictors, from a measurement performed in a nitrogen atmosphere is shown in Table 2.
For the absorption spectra, the plot of the lactose concentration with the peak area shows that the slope from the 1.37 THz absorption peak is higher than that from the 0.53 THz peak, denoting how lactose vibrational modes are more sensitive to 1.37 THz waves. Peak area of the 0.53 THz absorption peak results in a slightly more accurate prediction of lactose concentration, with r2 values of 0.9923 and 0.9697 for the 0.53 THz and 1.37 THz absorption peaks, respectively. This is also reflected in the RMSE values shown in Table 2. The better performance of the 0.53 THz peak is consistent with the analysis of a previous study, whereby the peak is found to have a more linear baseline and clear start and end points [27]. The plot of the lactose concentration with the peak height, on the other hand, shows very similar trendlines between the 0.53 THz and 1.37 THz absorption peaks. Similar to the peak area models, the model from the 0.53 THz absorption peak has a higher r2 value (and lower RMSE value) than the 1.37 THz absorption peak, 0.9334 and 0.8885, respectively. Both models from the peak area had a better prediction performance than those from the peak height, as the peak area takes into account the unequal peak broadenings and peak shape distortions. The result is consistent with a previous study on α-lactose monohydrate, whereby prediction using the peak area of the 0.53 THz absorption peak had a better performance than that using the peak height [19].
For the absorption coefficient spectra, plots of the 1.37 THz show higher slopes than those of 0.53 THz in a similar manner to the absorption spectral data. However, slightly more accurate predictions are achieved with the 1.37 THz peak, in contrast to the results from the absorption spectra. In addition, the peak height of both peaks resulted in higher r2 values than the peak area.
Comparing the results from the absorption and absorption coefficient spectral data, most models from the absorption coefficient spectra had better prediction accuracy than those from the absorption spectra. The more pronounced etalon effect in the absorption spectra could have led to greater difficulties in the baseline removal. This is responsible for the relatively high limit of detection (see Section 2.1). The use of the absorption coefficient spectra can reduce the etalon effect, thereby improving the limit of detection with high prediction accuracy.

2.3. Application on Diary Product Samples

The qualitative and quantitative determination of lactose using THz spectroscopy were performed on three commercial infant formula samples of the same brand. Sample 1 and Sample 2 are normal infant formula products and Sample 3 is a lactose-free product (see product labels in Figures S2–S4, Supplementary Materials). To assess the prediction performance of the quantification model, Sample 1 and Sample 2 were sent to two external accredited laboratories for lactose quantification (using AOAC official methods 980.13 and 982.14) and are used as reference values.
The absorption and absorption coefficient spectra obtained from measurement in a nitrogen atmosphere are shown together with pure lactose standard samples (0% and 15% (w/w) lactose) in Figure 4. Absorption peaks at 0.53 THz and 1.37 THz are clearly observed from the spectra of Sample 1 and Sample 2, signifying the presence of α-lactose monohydrate in these samples. These absorption peaks are not observed in the spectrum of the lactose-free Sample 3.
For the quantification performance, the mean predicted concentration of lactose using the absorption peak area are shown in Table 3. The lactose concentration obtained from the HPLC measurements in the external laboratories and carbohydrate content from the product nutritional labels are shown as a reference. In this table, the predicted concentrations of lactose anhydrous (molar mass = 342.3 g mol−1), which is ≈95%, by mass, that of α-lactose monohydrate (molar mass = 360.3 g mol−1), are used for comparison to align to the values from the product labels and HPLC.
According to product label and HPLC quantification, Sample 1 and Sample 2 have similar lactose concentrations (~60%). However, THz-TDS results differ significantly (~60% vs. ~20% for Samples 1 and 2, respectively). This discrepancy can be due to the distinctive formulations of the two samples (see Figures S2 and S3, Supplementary Materials), with different sources of lactose and preparation techniques. This can result in lactose of different forms and properties, to which THz-TDS is sensitive [25].
The models from the 0.53 THz peak in the absorption coefficient spectra resulted in a prediction closest to the result from the gold standard method of HPLC in Sample 1. However, the prediction from the other models were far from the HPLC result. The inaccuracy of prediction can be due to multiple reasons. First, there is a significant baseline shift in the spectra of the infant formula samples compared to the spectra of the pure standards (Figure 4), resulting from the infant formula matrix. It has been found that the magnitude of the baseline can be influenced by particle size [27]. The scattering effects in a nonlinear baseline also vary with the sample thickness [26]. Second, to achieve an accurate prediction of compound concentration, the selection of optimum conditions and parameters are necessary. For instance, the linearity of the data was improved with the use of absorption coefficient spectra, resulting in a better prediction performance compared to the other models.
While the prediction in Sample 1 was quite close to the HPLC result, the prediction in Sample 2 differed significantly. This can be attributed to the distinctive formulation and processing methods for the two samples, resulting in the presence of lactose in different forms. Lactose may exist in its amorphous form and crystalline forms of α-lactose monohydrate, α-lactose anhydrous and β-lactose anhydrous. The major form of lactose used in commercial infant formula is α-lactose monohydrate, due to its non-hygroscopicity [28]. However, it may not be at its maximum purity and may undergo mutarotation or conversion to other anomers during storage [28,29]. This can occur depending on storage conditions such as temperature and humidity. While HPLC resulted in similar concentration values between the two samples, THz-TDS results differed considerably due to its sensitivity to distinctive forms of a compound.
For Sample 3, the negative concentration value predicted from the 0.53 THz absorption peak could be due to an improper baseline fitting, as mentioned earlier. The other three predictors resulted in lactose concentrations close to 0% (0 is in the range of mean ± SD) or no detection of lactose, which is in line with the “lactose-free” label.
Overall, THz-TDS can be employed to verify the lactose-free claim in products. However, the accuracy of quantification in lactose-containing products can be influenced by its presence in multiple forms. Pioneer THz-TDS studies on mixtures of pure compounds with the utilisation of nonlinear quantification models were explored by the authors [30]. This may be applicable to lactose in the future.

3. Materials and Methods

Analytical grade α-lactose monohydrate and polyethylene (PE) powder of >98% purity were purchased from Sigma-Aldrich. Commonly available infant formula products were bought from a local supermarket. Chemicals and samples were stored in a dry cabinet at room temperature without any pre-treatment.
To obtain the absorption spectra of pure α-lactose monohydrate samples, twelve mixtures of α-lactose monohydrate standard and polyethylene were prepared at concentrations of 0%, 1%, 3%, 5%, 10%, 15%, 20%, 49%, 80% and 100% (w/w). The 0% and 100% (w/w) samples were prepared in duplicates. Mixtures were grounded to fine powders and mixed in a mortar until well-combined. They were formed into disc pellets of 13 mm diameter and approximately 100 mg weight. Each mixture was pressed into pellets with a hydraulic press under a pressure of 7 tons for 5 min. The thicknesses of the pellets were between 0.550 and 0.980 mm. The mass and thickness of each pellet are reported in Table S1, Supplementary Materials. Pellets removed from the mold were stored in a plastic (PE) bag inside a dry cabinet.
Samples of commercial infant formula were used as a test set to assess the quantification models derived from the spectra of pure samples. The samples were prepared to a final weight of 100 mg (without PE), in a similar manner to the pure α-lactose monohydrate standards.
Spectroscopic measurements were conducted with a commercial THz-TDS (TF4-1511, Toptica Photonics, Germany) with a resolution of 0.005 THz at room temperature. The system is housed in a custom-made cabinet (Figure S1, Supplementary Materials). Radiated THz waves were delivered by four parabolic mirrors, propagated through the samples and then to the THz waves receiver. The signal transmitted from the sample holder (without a sample) was used as the reference signal E ref t and the signal transmitted through the sample is subsequently obtained as the sample signal E sam t . The time domain pulse is transformed to the frequency domain with fast Fourier transform (FFT) based on E ref t and E sam t , Equation (1):
E ˜ ω = E ω e i φ ω = E t e i ω t φ ω d t
where E ω and φ ω are the amplitude and phase of the electric field, respectively, and E t is the time–domain waveform. The refractive index and absorption coefficient spectra were obtained using the optical constant extraction model by Dorney [31] and Duvillaret [32,33]. The refractive index n ω and absorption coefficient α ω of the samples were calculated by the following Equations (2) and (3):
n ω = 1 + φ ω c ω d
α ω = 2 d ln 4 n ω ρ ω n ω + 1 2
where d is the thickness of the sample, c = 3 × 108 ms−1 is the speed of electromagnetic wave in vacuum, φ ω is the phase difference between the sample signal and the reference signal, and ρ ω is the module ratio of the amplitude of sample signal E ω sam to the amplitude of reference signal E ω ref after Fourier transform.
All spectra in this study are averaged results of over 100 acquisitions to optimise signal-to-noise ratio. Sets of THz-TDS measurements were performed under ambient conditions (≈37% RH) and in nitrogen-purged atmosphere (≤7% RH). The first set and second set were performed at sampling rates of 250 and 500 times per measurement, respectively, with replicate measurements in nitrogen atmosphere. Results presented in this paper are from the second dataset. Results from all datasets can be found in the Supplementary Materials.
Fourier transform of raw time–domain data, curve-fitting (Gaussian and Lorentzian) and baseline correction for the computation of peak area and height were performed with Python 3.10 programming language and Jupyter IDE (version 2022). Visualisation of spectra and linear regression analysis for the quantification models were performed using Wolfram Mathematica 12.3. Raw data and source codes (Jupyter and Mathematica notebooks) for analysis are available in the Supplementary Materials.

4. Conclusions

In this pilot study, a lactose detection and quantification procedure by THz-TDS was presented and evaluated with infant formula samples. Quantification models from the peak area and peak height of the absorption and absorption coefficient spectra were compared. The peak areas of two absorption peaks at frequencies of 0.53 THz and 1.37 THz resulted in a high quantification performance in pure α-lactose monohydrate standard samples. Although satisfactory results were obtained in pure standards, the prediction in infant formula samples may not be reflective of the total lactose content in the product. These were due to baseline shifts arising from the infant formula matrix and potential mutarotation of lactose anomers. Optimum conditions and predictors are also necessary for accurate quantification. Since THz-TDS is sensitive to different forms of the analyte, it may not be suitable for the accurate in-line determination of sugars in processed food and pharmaceutical products. THz-TDS can be a promising technique for the rapid, non-destructive and in-line detection of compounds. Nevertheless, caution must be taken in the quantitative analysis when the analyte can be present in a myriad of forms.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules27155040/s1, Table S1: Pellet mass and thickness of pure α-lactose monohydrate standards and infant formula samples for THz-TDS measurement; Figure S1: Custom-made THz-TDS system in this study at National Electronics and Computer Technology Centre, Pathum Thani, Thailand (a zoom-in of the sample holder is framed in red); Figure S2: Product label of infant formula Sample 1; Figure S3: Product label of infant formula Sample 2; Figure S4: Product label of infant formula Sample 3; data and source codes.

Author Contributions

Conceptualisation, K.P. and T.L.; methodology, K.P., N.K., P.R., J.Y.C. and T.C.; software, S.D., K.P., N.K. and T.L.; validation, S.D.; formal analysis, S.D.; investigation, K.P., P.R., J.Y.C., R.J. and T.C.; resources, K.P.; data curation, K.P. and T.L.; writing—original draft preparation, S.D., K.P. and T.L.; writing—review and editing, K.P. and T.L.; visualisation, S.D., N.K. and T.L.; supervision, K.P.; project administration, K.P.; funding acquisition, K.P. and T.L. All authors have read and agreed to the published version of the manuscript.

Funding

S.D. and T.L were supported by a mid-career researcher development grant (Project No. NRCT5-RSA63015-22) jointly funded by the National Research Council of Thailand and Mahidol University, and MU-GPI grant (Contract No. 06/2565) by Mahidol University. J.Y.C. was funded by National Electronics and Computer Technology Center (Project No. P2150039). K.P. is grateful for equipment support by NSTDA’s Food and Feed Innovation Center (Project No. P1850906).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data are included in the Supplementary Materials.

Acknowledgments

The authors acknowledge the assistance of Kridtin Chinsukserm in the early stage of our investigation. We thank Kantapong Sucharitpongpan, Nichakarn Termsaithong and Alicha Wararuk for their assistance in the sample preparation. We thank Junming Ho (UNSW) for discussions on the possibility of THz-TDS spectral simulation by computational methods. The late Kittipong Kasamsook’s foundation of the Terahertz Technology Research Team is greatly appreciated.

Conflicts of Interest

The authors declare no conflict of interest. Funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Sample Availability Statement

Samples of the compounds are not available from the authors. Sample information and their sources are reported in the Supplementary Materials.

References

  1. Food and Agriculture Organizaion. Gateway to Dairy Production and Products. Available online: https://www.fao.org/dairy-production-products/products/en/ (accessed on 2 March 2022).
  2. Fox, P.F.; Uniacke-Lowe, T.; McSweeney, P.L.H.; O’Mahony, J.A. Lactose. In Dairy Chemistry and Biochemistry; Fox, P.F., Uniacke-Lowe, T., McSweeney, P.L.H., O’Mahony, J.A., Eds.; Springer: Cham, Switzerland, 2015; pp. 21–68. [Google Scholar]
  3. National Institute of Diabetes and Digestive and Kidney Diseases. Lactose Intolerance. Available online: https://www.niddk.nih.gov/health-information/digestive-diseases/lactose-intolerance (accessed on 2 March 2022).
  4. Fox, P.F. Lactose: Chemistry and Properties. In Advanced Dairy Chemistry; McSweeney, P., Fox, P.F., Eds.; Springer: New York, NY, USA, 2009; pp. 1–15. [Google Scholar]
  5. Hebbink, G.A.; Dickhoff, B.H.J. Chapter 5—Application of lactose in the pharmaceutical industry. In Lactose; Paques, M., Lindner, C., Eds.; Academic Press: Cambridge, MA, USA, 2019; pp. 175–229. [Google Scholar]
  6. Paterson, A. Production and uses of lactose. In Advanced Dairy Chemistry; McSweeney, P., Fox, P.F., Eds.; Springer: New York, NY, USA, 2009; pp. 105–120. [Google Scholar]
  7. Churakova, E.; Peri, K.; Vis, J.S.; Smith, D.W.; Beam, J.M.; Vijverberg, M.P.; Stor, M.C.; Winter, R.T. Accurate analysis of residual lactose in low-lactose milk: Comparing a variety of analytical techniques. Int. Dairy J. 2019, 96, 126–131. [Google Scholar] [CrossRef]
  8. Eadala, P.; Waud, J.; Matthews, S.B.; Green, J.T.; Campbell, A.K. Quantifying the ‘hidden’ lactose in drugs used for the treatment of gastrointestinal conditions. Aliment. Pharmacol. Ther. 2009, 29, 677–687. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Garballo-Rubio, A.; Soto-Chinchilla, J.; Moreno, A.; Zafra-Gómez, A. Determination of residual lactose in lactose-free cow milk by hydrophilic interaction liquid chromatography (HILIC) coupled to tandem mass spectrometry. J. Food Compos. Anal. 2018, 66, 39–45. [Google Scholar] [CrossRef]
  10. Martínez-Aviñó, A.; Molins-Legua, C.; Campíns-Falcó, P. Combining high performance thin layer chromatography with minispectrometer-fiber optic probe-coupled to smartphone for in place analysis: Lactose quantification in several matrices. J. Chromatogr. A 2022, 1661, 462694. [Google Scholar] [CrossRef]
  11. Warnecke, S.; Wu, J.X.; Rinnan, Å.; Allesø, M.; van den Berg, F.; Jepsen, P.U.; Engelsen, S.B. Quantifying crystalline α-lactose monohydrate in amorphous lactose using terahertz time domain spectroscopy and near infrared spectroscopy. Vib. Spectrosc. 2019, 102, 39–46. [Google Scholar] [CrossRef]
  12. Monakhova, Y.B.; Kuballa, T.; Leitz, J.; Andlauer, C.; Lachenmeier, D.W. NMR spectroscopy as a screening tool to validate nutrition labeling of milk, lactose-free milk, and milk substitutes based on soy and grains. Dairy Sci. Technol. 2012, 92, 109–120. [Google Scholar] [CrossRef] [Green Version]
  13. Pokrzywnicka, M.; Koncki, R. Disaccharides determination: A review of analytical methods. Crit. Rev. Anal. Chem. 2018, 48, 186–213. [Google Scholar] [CrossRef]
  14. Bo, Y.; Fang, J.; Zhang, Z.; Xue, J.; Liu, J.; Hong, Z.; Du, Y. Terahertz and Raman Spectroscopic Investigation of Monohydrate Cocrystal of Antitubercular Isoniazid with Protocatechuic Acid. Pharmaceutics 2021, 13, 1303. [Google Scholar] [CrossRef]
  15. Hua, Y.; Zhang, H. Qualitative and quantitative detection of pesticides with terahertz time-domain spectroscopy. IEEE Trans. Microw. Theory Tech. 2010, 58, 2064–2070. [Google Scholar] [CrossRef]
  16. Im, J.; Goo, T.; Kim, J.; Choi, S.; Hong, S.J.; Bahk, Y.-M. Detection of Microplastic in Salts Using Terahertz Time-Domain Spectroscopy. Sensors 2021, 21, 3161. [Google Scholar] [CrossRef]
  17. Taday, P.F. Applications of terahertz spectroscopy to pharmaceutical sciences. Philos. Trans. R. Soc. London. Ser. A Math. Phys. Eng. Sci. 2004, 362, 351–364. [Google Scholar] [CrossRef] [PubMed]
  18. Yamauchi, S. Terahertz time-domain spectroscopy to identify and evaluate anomer in lactose. Am. J. Anal. Chem. 2013, 4, 756. [Google Scholar] [CrossRef] [Green Version]
  19. Wang, Y.; Zhao, Z.; Qin, J.; Liu, H.; Liu, A.; Xu, M. Rapid in situ analysis of l-histidine and α-lactose in dietary supplements by fingerprint peaks using terahertz frequency-domain spectroscopy. Talanta 2020, 208, 120469. [Google Scholar] [CrossRef] [PubMed]
  20. McIntosh, A.I.; Yang, B.; Goldup, S.M.; Watkinson, M.; Donnan, R.S. Crystallization of amorphous lactose at high humidity studied by terahertz time domain spectroscopy. Chem. Phys. Lett. 2013, 558, 104–108. [Google Scholar] [CrossRef]
  21. Yamauchi, S.; Hatakeyama, S.; Imai, Y.; Tonouchi, M. Nondestructive evaluation of crystallized-particle size in lactose-powder by terahertz time-domain spectroscopy. Opt. Eng. 2013, 53, 031203. [Google Scholar] [CrossRef] [Green Version]
  22. Peretti, R.; Mitryukovskiy, S.; Froberger, K.; Mebarki, M.A.; Eliet, S.; Vanwolleghem, M.; Lampin, J.-F. THz-TDS time-trace analysis for the extraction of material and metamaterial parameters. IEEE Trans. Terahertz Sci. Technol. 2018, 9, 136–149. [Google Scholar] [CrossRef] [Green Version]
  23. Hou, L.; Shi, W.; Dong, C.; Yang, L.; Wang, Y.; Wang, H.; Hang, Y.; Xue, F. Probing trace lactose from aqueous solutions by terahertz time-domain spectroscopy. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2021, 246, 119044. [Google Scholar] [CrossRef]
  24. Huang, H.; Shao, S.; Wang, G.; Ye, P.; Su, B.; Zhang, C. Terahertz spectral properties of glucose and two disaccharides in solid and liquid states. iScience 2022, 25, 104102. [Google Scholar] [CrossRef]
  25. Baranova, A.; Lykina, A.; Antonova, D.; Smolyanskaya, O. Optical Properties of Crystalline Lactose Fluidized with Dilutions of Various Substances in the Terahertz Frequency Range. Pharmaceutics 2022, 14, 32. [Google Scholar] [CrossRef]
  26. Namkung, H.; Kim, J.; Chung, H.; Arnold, M.A. Impact of pellet thickness on quantitative terahertz spectroscopy of solid samples in a polyethylene matrix. Anal. Chem. 2013, 85, 3674–3681. [Google Scholar] [CrossRef]
  27. Smith, G.; Hussain, A.; Bukhari, N.I.; Ermolina, I. Quantification of residual crystallinity in ball milled commercially sourced lactose monohydrate by thermo-analytical techniques and terahertz spectroscopy. Eur. J. Pharm. Biopharm. 2015, 92, 180–191. [Google Scholar] [CrossRef] [PubMed]
  28. Portnoy, M.; Barbano, D.M. Lactose: Use, measurement, and expression of results. J. Dairy Sci. 2021, 104, 8314–8325. [Google Scholar] [CrossRef]
  29. Altamimi, M.J.; Wolff, K.; Nokhodchi, A.; Martin, G.P.; Royall, P.G. Variability in the α and β anomer content of commercially available lactose. Int. J. Pharm. 2019, 555, 237–249. [Google Scholar] [CrossRef] [Green Version]
  30. Loahavilai, P.; Datta, S.; Kingkan, C.; Prasertsuk, K.; Jintamethasawat, R.; Rattanawan, P.; Chia, J.Y.; Thanapirom, C.; Limpanuparb, T. Chemometric Analysis of Ternary Mixture of Caffeine, Quinic Acid and Nicotinic Acid by Terahertz Spectroscopy. ACS Omega 2022. [Google Scholar] [CrossRef]
  31. Dorney, T.D.; Baraniuk, R.G.; Mittleman, D.M. Material parameter estimation with terahertz time-domain spectroscopy. JOSA A 2001, 18, 1562–1571. [Google Scholar] [CrossRef] [Green Version]
  32. Duvillaret, L.; Garet, F.; Coutaz, J.-L. A reliable method for extraction of material parameters in terahertz time-domain spectroscopy. IEEE J. Sel. Top. Quantum Electron. 1996, 2, 739–746. [Google Scholar] [CrossRef] [Green Version]
  33. Duvillaret, L.; Garet, F.; Coutaz, J.-L. Influence of noise on the characterization of materials by terahertz time-domain spectroscopy. JOSA B 2000, 17, 452–461. [Google Scholar] [CrossRef] [Green Version]
Figure 1. THz absorption spectra of α-lactose monohydrate at varying concentrations under (a) ambient condition (≈37% RH) and (b) nitrogen-rich (≤7% RH) atmosphere. Frequencies of absorption peaks (0.53 THz, 1.20 THz, 1.37 THz, 1.82 THz) are signified by vertical dashed lines. Absorbance values (arbitrary unit) are the ratio between the intensity of the substance and the reference, as described in Equation (1) under Section 3.
Figure 1. THz absorption spectra of α-lactose monohydrate at varying concentrations under (a) ambient condition (≈37% RH) and (b) nitrogen-rich (≤7% RH) atmosphere. Frequencies of absorption peaks (0.53 THz, 1.20 THz, 1.37 THz, 1.82 THz) are signified by vertical dashed lines. Absorbance values (arbitrary unit) are the ratio between the intensity of the substance and the reference, as described in Equation (1) under Section 3.
Molecules 27 05040 g001
Figure 2. (a) THz refractive index (n) and (b) absorption coefficient (α in cm−1) spectra of α-lactose monohydrate at varying concentrations under nitrogen-rich atmosphere. Frequencies of absorption peaks (0.53 THz, 1.20 THz, 1.37 THz, 1.82 THz) are signified by vertical dashed lines. The calculation of refractive index and absorption coefficient are described in Equations (2) and (3) under Section 3.
Figure 2. (a) THz refractive index (n) and (b) absorption coefficient (α in cm−1) spectra of α-lactose monohydrate at varying concentrations under nitrogen-rich atmosphere. Frequencies of absorption peaks (0.53 THz, 1.20 THz, 1.37 THz, 1.82 THz) are signified by vertical dashed lines. The calculation of refractive index and absorption coefficient are described in Equations (2) and (3) under Section 3.
Molecules 27 05040 g002
Figure 3. The relationship between α-lactose monohydrate concentration and (a) absorbance peak area, (b) absorbance peak height, (c) absorption coefficient peak area and (d) absorption coefficient peak height of the absorption peaks at 0.53 THz and 1.37 THz, from measurements in nitrogen atmosphere. Simple linear regression analysis [19] was applied on the datasets and results are reported in Table 2.
Figure 3. The relationship between α-lactose monohydrate concentration and (a) absorbance peak area, (b) absorbance peak height, (c) absorption coefficient peak area and (d) absorption coefficient peak height of the absorption peaks at 0.53 THz and 1.37 THz, from measurements in nitrogen atmosphere. Simple linear regression analysis [19] was applied on the datasets and results are reported in Table 2.
Molecules 27 05040 g003
Figure 4. THz (a) absorption and (b) absorption coefficient spectra of infant formula samples, 15% (w/w) lactose and 0% (w/w) lactose standard samples.
Figure 4. THz (a) absorption and (b) absorption coefficient spectra of infant formula samples, 15% (w/w) lactose and 0% (w/w) lactose standard samples.
Molecules 27 05040 g004
Table 1. Comparison of common analytical methods for lactose quantification.
Table 1. Comparison of common analytical methods for lactose quantification.
MethodPrincipleSample PreparationReagents
chromatographic: HPLC, HILIC, HPTLC, HPAECseparation of compounds based on polarity, electrical charge or molecular sizesample purification, solid removalsolvents
spectrophotometricconversion of lactose into its monomers for further reactions, followed by spectrophotometric quantificationhydrolysis of lactose, addition of reagentssolvent,
enzyme,
oxidising agents
spectroscopic:
NMR
determination of molecular structure through chemical shifts when magnetic field is appliedsample purification, solid removalsolvents
spectroscopic:
NIR, THz
determination of molecular structure through molecular vibrationsgrounding, water removal-
Table 2. Equation of trendline and RMSE of different prediction models for lactose quantification.
Table 2. Equation of trendline and RMSE of different prediction models for lactose quantification.
PredictorAbsorption PeakEquation of TrendlineRMSE
AbsorbancePeak area0.53 THz3.43 × 10−4 x + 0.00210.0011
1.37 THz7.06 × 10−4 x + 0.00410.0047
Peak height0.53 THz6.54 × 10−3 x + 0.08400.0664
1.37 THz6.46 × 10−3 x + 0.09900.0870
Absorption coefficient Peak area0.53 THz1.98 × 10−2 x − 0.02880.0840
1.37 THz6.71 × 10−2 x − 0.24680.2127
Peak height0.53 THz6.08 × 10−1 x − 0.67941.8840
1.37 THz1.06 x − 2.60183.1344
Table 3. Mean predicted concentration of lactose anhydrous in infant formula samples using 0.53 THz and 1.37 THz absorption peak area from THz-TDS measurements conducted in nitrogen atmosphere, mean concentration obtained from HPLC measurements in external laboratories and product nutritional labelling.
Table 3. Mean predicted concentration of lactose anhydrous in infant formula samples using 0.53 THz and 1.37 THz absorption peak area from THz-TDS measurements conducted in nitrogen atmosphere, mean concentration obtained from HPLC measurements in external laboratories and product nutritional labelling.
Mean Concentration ± SD in % (w/w)
Sample 1Sample 2Sample 3
Absorbance0.53 THz39.0 ± 1.61.8 ± 1.0–19.4 ± 1.5
1.37 THz45.1 ± 1.511.4 ± 0.70.0 ± 0.6
Absorption
coefficient
0.53 THz62.5 ± 4.423.6 ± 9.21.8 ± 10.1
1.37 THz68.4 ± 1.620.8 ± 0.5not detected *
HPLC
(from external laboratories)
61.97 (n = 2)59.97 (n = 3)not done
Product label:
carbohydrate content
6058assumed to be 0
(lactose-free)
*Absorption coefficient <–0.1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Datta, S.; Prasertsuk, K.; Khammata, N.; Rattanawan, P.; Chia, J.Y.; Jintamethasawat, R.; Chulapakorn, T.; Limpanuparb, T. Terahertz Spectroscopic Analysis of Lactose in Infant Formula: Implications for Detection and Quantification. Molecules 2022, 27, 5040. https://doi.org/10.3390/molecules27155040

AMA Style

Datta S, Prasertsuk K, Khammata N, Rattanawan P, Chia JY, Jintamethasawat R, Chulapakorn T, Limpanuparb T. Terahertz Spectroscopic Analysis of Lactose in Infant Formula: Implications for Detection and Quantification. Molecules. 2022; 27(15):5040. https://doi.org/10.3390/molecules27155040

Chicago/Turabian Style

Datta, Sopanant, Kiattiwut Prasertsuk, Nuttawat Khammata, Patharakorn Rattanawan, Jia Yi Chia, Rungroj Jintamethasawat, Thawatchart Chulapakorn, and Taweetham Limpanuparb. 2022. "Terahertz Spectroscopic Analysis of Lactose in Infant Formula: Implications for Detection and Quantification" Molecules 27, no. 15: 5040. https://doi.org/10.3390/molecules27155040

APA Style

Datta, S., Prasertsuk, K., Khammata, N., Rattanawan, P., Chia, J. Y., Jintamethasawat, R., Chulapakorn, T., & Limpanuparb, T. (2022). Terahertz Spectroscopic Analysis of Lactose in Infant Formula: Implications for Detection and Quantification. Molecules, 27(15), 5040. https://doi.org/10.3390/molecules27155040

Article Metrics

Back to TopTop