2D-Cross Correlation Spectroscopy Coupled with
Molecular Fluorescence Spectroscopy for Analysis of
Molecular Structure Modification of Camel Milk and
Cow Milk Mixtures during Coagulation

Oumayma Boukria, El Mestafa El Hadrami, Shaxnoza Sultanova, Jasur Safarov, Françoise Leriche, Abderrahmane Ait Kaddour

► To cite this version:

Oumayma Boukria, El Mestafa El Hadrami, Shaxnoza Sultanova, Jasur Safarov, Françoise Leriche, et al.. 2D-Cross Correlation Spectroscopy Coupled with Molecular Fluorescence Spectroscopy for Analysis of Molecular Structure Modification of Camel Milk and Cow Milk Mixtures during Coagulation. Foods, 2020, 9 (6), pp.724. 10.3390/foods9060724 . hal-03081570

HAL Id: hal-03081570

https://vetagro-sup.hal.science/hal-03081570

Submitted on 18 Dec 2020

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

2D-Cross Correlation Spectroscopy Coupled with
Molecular Fluorescence Spectroscopy for Analysis of
Molecular Structure Modification of Camel Milk and
Cow Milk Mixtures during Coagulation

Oumayma Boukria <fnref n="1" />, El Mestafa El Hadrami <fnref n="1" /> ID, Shaxnoza Sultanova <fnref n="2" />, Jasur Safarov <fnref n="2" /> ID,
Françoise Leriche <fnref n="3" /> and Abderrahmane Aït-Kaddour <fnref n="3" />,* ID

Received: 2 May 2020; Accepted: 27 May 2020; Published: 2 June 2020

<figure><img src="image_2.png" /></figure>

check for
updates

**Abstract:** Synchronous fluorescence spectroscopy (SFS) coupled with two-dimensional correlation spectroscopy (2DCOS) was employed to monitor, at the molecular level, the coagulation of five mixture ratios of camel's milk (CaM) and cow's milk (CM) (100% CaM, 75% CaM:25% CM, 50% CaM:50% CM, 25% CaM:75% CM and 100% CM). The dissimilarities among the different formulations are highlighted on the synchronous 2DCOS-SFS. In addition, according to the cross-peak symbols in synchronous and asynchronous spectra, the rate of response modification in riboflavin, protein and vitamin A matched with common coagulation phenomena usually reported during chymosin coagulation (hydrolysis of κ-casein, destabilization of casein micelles and aggregation). This study demonstrated that 2DCOS-SFS is a successful strategy to discriminate milk mixtures and to monitor molecular structure modifications during coagulation process.

**Keywords:** milk; mixture; spectroscopy; 2D COS

**1. Introduction**

Nowadays, the production of cow's milk (CM) is dominant, but the original quality of non-bovine
milk, such as camel's milk (CaM), is now better described and is experiencing growing interest [1].
CaM, goat's milk, ewe's milk and buffalo's milk have received special attention because of their
recognition as potential functional foods from a nutritional point of view. It is admitted that CaM
presents a high nutritional quality; e.g., it has three times more vitamin C, minerals (e.g., K+, Cu2+ and
Mn2+), and essential and polyunsaturated Fatty Acids (FAs) than CM [2,3]. It is also considered
as exhibiting properties to manage chronic ailments [4], e.g., tuberculosis, jaundice and asthma.
Thus, mixing milks from different species appears as a good strategy to increase the consumption of
non-bovine milk and enable consumers and dairy companies to benefit from their nutritional and
technological advantages. This opportunity requires the availability of suitable methods to assess the
quality features of products but also to control the different stages of their manufacturing process.
Milk coagulation is an important step in the development of texture of different dairy products

<footnote n="1">¹ Applied Organic Chemistry Laboratory, Sciences and Techniques Faculty, Sidi Mohamed Ben Abedallah University, Route d'Immouzer, Fès BP 2202, Morocco; boukriaoumayma@hotmail.com (O.B.); elmestafa.elhadrami@usmba.ac.ma (E.M.E.H.)</footnote>

<footnote n="2">² Department of Food Engineering, Faculty of Mechanical Building, Tashkent State Technical University Named after Islam Karimov, University str. 2, Tashkent 100095, Uzbekistan; sh.sultanova@yahoo.com (S.S.); jasursafarov@yahoo.com (J.S.)</footnote>

<footnote n="3">³ Université Clermont-Auvergne, INRAE, VetAgro Sup, 63370 Lempdes, France;
francoise.leriche@vetagro-sup.fr</footnote>

<footnote>* Correspondence: abderrahmane.aitkaddour@vetagro-sup.fr; Tel.: +33-(0)4-73-98-13-78</footnote>
(e.g., cheeses, yogurts). Therefore, a good management of this operation is essential to obtain a cheese with a good texture. Among other factors, it was reported that texture features of cheese are affected by their structural properties, composition and fat distribution that are highly dependent on the making process [5]. The management of the coagulation is also important because it was reported that the microstructure of dairy products is closely related to other quality features of cheeses, like physicochemical properties, flavor, color, nutritional profile and bioavailability of nutrients [6].

Dairy products contain natural fluorophores like Maillard reaction products, amino acid groups (e.g., tryptophan, tyrosine), vitamins (riboflavin and vitamin A) and nicotinamide adenine dinucleotide (NADH), among others. Each fluorophore presents a specific excitation and emission spectrum that can be used to study structural changes at the molecular level and predict and identify the composition of food products. The presence of these fluorescence molecules is used as an advantage in fluorescence spectroscopy in order to perform non-destructive analysis of dairy products especially when coupled with chemometrics [7].

Different excitation or emission configurations can be used to study food and dairy products. Among these configurations, synchronous fluorescence spectroscopy (SFS) seems very interesting because it is a good compromise between classical fluorescence spectroscopy (fixing the excitation wavelength and recording emission spectrum and vice versa) and total excitation emission matrix (EEM). In the SFS configuration, spectra are recorded after scanning simultaneously the excitation and emission monochromators with a wavelength difference (λ) providing a good fingerprint identification of the product fluorophores [8]. In dairy product analysis, SFS was used to discriminate between milk samples [9], to determine the effects of mild heating, acidification and added minerals on structural changes of milk components [7,10,11] and to study effects of salts and heating on cheese molecular structure [12–14].

The understanding of the impact of the association of milk from different species on molecular structure, location and interactions of components (e.g., fat, proteins) during manufacture is essential for predicting quality attributes of dairy products and particularly cheese. However, it is generally difficult to describe the milk coagulation by using fluorescence spectroscopy from conventional one-dimension fluorescence spectra (i.e., wavelength vs. fluorescence intensity) or by using conventional chemometric analytical methods (e.g., principal component analysis and partial least squares discriminant analysis).

The concept of 2D correlation spectroscopy (2D COS) was first introduced by Noda [15,16] and later extended to the generalized 2D correlation spectroscopy [17]. The generalized 2D COS is a useful technique that has been developed to deal with signal fluctuations induced by different external perturbation factors, like time, or any other physical variables such as temperature or pressure [17–19]. Generalized 2D COS is particularly powerful to simplify spectral data analysis due to their general complexity generally assigned to the number of observed bands and their potential overlapping. This method relies on calculating synchronous and asynchronous changes in signal intensities at different wavelengths. This method can be used to unfold spectra on a second dimension in order to increase spectral resolution and to carefully analyze chemical and structural information concerning molecular events [20,21], spectral response sequences, relative intensity changes, and so on [22].

Recently, 2DCOS was widely investigated in various areas to deeply analyze fluorescence spectra of, e.g., tissue culture, recycled water systems, and humic substances [23–26]. Nonetheless, no study evaluated the potential of this method coupled with fluorescence spectroscopy to study the molecular structure modifications occurring during milk coagulation and to discriminate mixture of milk samples from different species during coagulation. Therefore, the present paper will focus on describing coagulation of milk mixtures from CaM and CM by a method combing synchronous fluorescence spectra and 2DCOS. This was performed in order to appraise the development of an accurate, rapid and feasible analytical method to monitor milk coagulation and to distinguish between different formulations of milk during coagulation because this can be useful for identification of milk adulteration.
## 2. Materials and Methods

### 2.1. Materials

Raw CM was purchased from a local dairy farm located in Marmilaht, France, while raw CaM was purchased from a local farm located at Fez, Morocco. Milk samples were frozen at -20 °C and forwarded to the laboratory in refrigerated condition. The mean biochemical composition was analyzed with the Gerber method for fat, Kjeldahl method for protein and colorimetric method for lactose according to the normative references (ISO 488:2008, ISO 8968-1:2014, and ISO 26462:2010). The average composition of CM was 3.10% fat, 3.21% protein and 4.16% lactose and 4.77% fat, 2.96% protein and 4.82% lactose for CaM. The samples were maintained at this temperature until analysis. Before the analysis, samples were defrosted at 4 °C for 24 h. Then, five different formulations were prepared by mixing both types of milk in different proportions. Before mixing CaM and CM, samples were warmed to 40 °C (±1 °C) in a water bath and gently mixed for 1 min by hand. The volume fractions (%) of CaM in the different formulations were 100% (i.e., CaM), 75% (i.e., CaM3:CM1, v/v), 50% (i.e., CaM1:CM1, v/v), 25% (i.e., CaM1:CM3, v/v) and 0% (i.e., CM). Based on the recommendation of the manufacturer, milk coagulations were realized at 40 °C (±1 °C) after adding 2.5 µL/mL of commercial rennet (CHY-MAX® M-Chr. Hansen) containing 50 mg/L of chymosin. After adding chymosin, the milk was gently mixed by hand for a few seconds and rapidly deposited in pre-warmed (40 °C) measurement cell of the spectral device.

### 2.2. Fluorescence Spectroscopy

The front face excitation synchronous fluorescence spectra were recorded by a FluoroMax-4 (Jobin-Yvon, Longjumeau, France) with a signal to noise ratio of 3000:1, equipped with a solid sample holder and a quartz cell (30 × 10 × 10 mm). During acquisition, the excitation wavelength (λex) and emission wavelength (λem) were scanned in a synchronous way with an 80 nm interval (Δλ = λex – λem) [7]. Spectra were recorded from 250 to 550 nm at 40 °C (±1 °C) every 5 min during coagulation. The temperature was controlled with a water bath equipped with a thermostat control [27]. Collections of spectral data were obtained by using the FluorEssence software (2010), version 3.5 (Horiba scientific, Jobin-Yvon, Longjumeau, France). For each formulation, three replicates were performed. Synchronous fluorescence spectra were preprocessed by Standard Normal Variate (SNV) method performed with the PLS_toolbox 8.5.2 (Eigenvector Research, Manson, WA, USA) in MATLAB 8.0.0.783 (R2012b) (The MathWorks, Natick, MA, USA).

### 2.3. Two-Dimensional Correlation Spectroscopy

The 2DCOS-synchronous fluorescence spectra were calculated after treatment of the collections of coagulation time-dependent dynamic spectra with 2DCOS analysis method by using a home software program in MATLAB. The 2DCOS spectra (synchronous and asynchronous) were obtained after applying to initial sample spectra an external perturbation, here coagulation time, and measuring a series of dynamic spectra and then processing these dynamic spectra by the Noda algorithm [17]. The reference spectra were chosen as the average spectrum, as previously reported by Noda and Ozaki [28]. In the present study, the 2DCOS spectral intensities are plotted as a function of wavelength variables. The two orthogonal axes of spectral variables define the 2D map or spectra, and the third axis, normal to the plane, which is the intensity.

## 3. Results and Discussion

### 3.1. Spectral Analysis and Comparison between Milk Formulations

A high part of the fluorescence chemical molecules of milk samples can be revealed and identified in the synchronous fluorescence spectra. The normalized mean emission synchronous fluorescence patterns of the different milk samples recorded during coagulation are presented in
Figure 1. Four bands can be observed centered on 297 nm (emission 377 nm), 321 nm (emission 401 nm), 360 nm (emission 440) and 449 nm (emission 529 nm). Those bands were previously identified, in milk and cheese [13,14], and were assigned to tryptophan, vitamin A and riboflavin compounds, respectively. The analysis of the band intensities and shapes suggested that physicochemical and molecular structure differences of the aforementioned molecules could be observed during coagulation of the different milk formulations. The analysis of the initial spectra is very useful to extract general differences between samples by observing the occurrence of fluorophores in the materials after observing, principally, the peak positions and, to a lesser extent, their shapes and intensities. However, when considering all the spectra (results not shown) recorded during coagulation of the five milk formulations (i.e., 5 formulations × 24 spectra recorded during coagulation time = a total of 120 spectra) the small differences among them make it difficult to discriminate samples only by using the conventional chemometric technique (e.g., principal component analysis) and especially when only one SFS offset (i.e., 80 nm) is used for spectral acquisition. In order to overcome this difficulty, the 2D COS method was used.

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Normalized mean synchronous fluorescence spectra recorded during coagulation of CM (cow's milk), CaM (camel's milk) and their mixtures.</figcaption></figure>

## 3.2. Analysis of 2D COS Synchronous and Asynchronous Spectra

The synchronous spectrum of 2DCOS presents in the diagonal line peaks, named auto-peaks. Those auto-peaks are the result of autocorrelation between fluorescence bands and the dynamic fluctuation. They highlight the sensitivity of fluorescence to the external perturbation, here, coagulation time. The cross-peaks located at the off diagonal revealed the relativity of intensity variations of a pair of group vibrations corresponding to their frequencies. A positive cross-peak (in red or green color area in Figure 2) indicates a simultaneous increase or decrease of different fluorescence bands under the coagulation time. The more the intensity changes between fluorescence peaks is coordinated, the stronger the cross-peak is. In contrast, a negative cross-peak (blue area color in Figure 2) represents the coordinated changes of band intensities in the opposite directions, increase versus decrease [15,17].

The Figure 2a, b exhibits the synchronous 2DCOS fluorescence spectra and their respective auto-peaks spectra. The regions from 250–280 nm and 500–550 nm were excluded from the analysis because they only describe baseline modifications. Figure 3 presented the intensity of auto-peaks identified in the synchronous spectra. Figure 2a,b indicates that three major auto-peaks can be clearly identified for the CaM during coagulation, while five auto-peaks are clearly visible for the other formulations containing CM (CM, CaM1:CM1, CaM1:CM3 and CaM3:CM1). For CaM, the auto-peaks were located at 298, 322 and 361 nm. For CM, peaks identified were centered at 298, 322, 361, 412 and 449 nm.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Synchronous spectra (a), auto-peak spectra (b) and asynchronous spectra (c).</figcaption></figure>

The peak cluster between 280 and 390 nm was present in the five formulations and the peaks were all positively auto correlated. The 449 nm peak, when present in the map, is positively auto correlated to the peak cluster located between 280 and 390 nm, suggesting synchronous modification of the peak intensities during coagulation time. Nonetheless, differences in intensities were observed between samples (Figure 3). The highest intensities were noted for CaM, then for CM, CaM1:CM1, and finally for CaM1:CM3. For the 449 nm band, the highest intensity was observed for CaM1:CM3, CaM3:CM1 and finally CM.
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Relative auto-peak intensity of CM (cow's milk), CaM (camel's milk) and their mixtures.</figcaption></figure>

Concerning the band maxima of proteins, it generally reflects the average exposure of their
tryptophan residues to the aqueous phase [29]. The comparison between the tryptophan position
band of CaM, CM and samples containing a mixture of the two revealed a shift to lower wavelength
(298 towards 296 nm). Spectral shifts have been observed because of several phenomena, such as
tertiary structure change, binding of ligands and protein-protein association. In a previous study,
fluorescence of tryptophan was used to monitor the structural modifications of proteins and their
physicochemical environment during the coagulation process [30]. The excitation peak of tryptophan
is highly sensitive to its local environment, and can be assigned to different protein conformational
changes [29] during coagulation and to the gel structure differences obtained depending on milk
formulation. Concerning the fluorescence of vitamin A (322 nm), it was reported that this band did
not exhibit any correlation to the milk coagulation process. As reported by Boubellouta et al. [11],
during the coagulation of CM, the modification of the band intensity can indicate a fluorescence transfer
between tryptophan residues of caseins and vitamin A located in the membrane of fat globules. This is
in agreement with an interaction of caseins with fat-globule membrane during the formation of the
protein network [31,32].

When compared to CM and CaM, the formulations containing both CaM and CM presented
peaks located at the same positions only for two bands (i.e., 322 and 449). Concerning, the band
at 361 nm, this one was identified at 362, 362 and 367 for CaM3:CM1, CaM1:CM1 and CaM1:CM3
samples, respectively. This highlighted a red shift of the 361 nm band when increasing the content of
CaM in the formulation. This fluorescence band could be assigned to the photochemical degradation
product derived from riboflavin, lumichrome; with excitation/emission maxima generally identified at
360/450 nm in a model system [33] and 370/430 in yoghurt [34].

The small peak located at 412 nm in CM shifted to lower wavelength (403 nm) for formulations
containing both CaM and CM. No literature was found that accounted with certainty for the fluorophores
at the origin of this band. However, Kikugawa and Beppu [35] and Wold et al. [36] associated this
peak to fluorescence from stable oxidation products, resulting from aldehydes and amino acids.
The peaks centered at 403 nm presented the highest intensities for CaM3:CM1 compared to CaM1:CM3,
CaM1:CM1 and finally CM samples. This peak was negatively auto correlated with the other peaks.
This suggested that when this peak increased the other ones decreased during coagulation time
(Figure 2b).
Concerning the 449 nm band, it was generally assigned to riboflavin and can be used to identify protein conformation and interaction changes since it interacts with proteins [37].

The number, position and intensity of the auto-peaks are different among the milk formulations and can be used to differentiate molecular structure difference during their coagulation. Moreover, the profiles of the components are in agreement with the first inspection of the fluorescence landscapes, and further indicate the fact that the obtained fluorescence signal from the samples arise from tryptophan, vitamin A and riboflavin compounds. Nonetheless, in addition of identifying the synchronous relationships between fluorescence bands during coagulation, a new peak was highlighted at 403/412 nm. The results emphasized that changes in the intensity of vitamin A and riboflavin at 322 and 482 nm, respectively, can provide elements concerning molecular structure modifications taking place in the milk micelle structure [11]. The results also bring to the forefront that tryptophan and riboflavin fluorophores are relevant tracer molecules to monitor milk coagulation.

The asynchronous 2DCOS map gives the sequence modifications of the spectral band intensity at different wavelengths or λ alongside the coagulation time (external perturbation). When the perturbation is applied, a cross-peak at variable λ1 and λ2 emerges in the asynchronous 2DCOS spectra if the spectral intensity at variable λ1 changes before or after variation at variable λ2. As the intensity of the spectral bands is varying simultaneously with itself, no peak is observed along the diagonal line in the asynchronous 2DCOS map. According to Noda's theory, the sign of the cross peaks located at λ1 and λ2 in both synchronous and asynchronous 2DCOS maps can be used to bring to the forefront the sequence of the spectral intensity variations at λ1 and λ2 [38]. If the same signs are observed in synchronous and asynchronous cross-peaks, this would highlight that the intensity change in λ1 occurs before λ2, while opposite signs in the synchronous and asynchronous cross-peaks highlight that changes in band intensity at λ2 occurs before λ1.

The asynchronous spectra further revealed the sequence of the molecular fluorescence modification during milk coagulation (structure, interaction and environment). Below the diagonal line of the asynchronous maps of the different milk formulations, five peaks were identified at 298, 322, 360, 403 and 450 nm. The different cross-peaks formed by those bands are presented in Table 1. According to Noda's theory, the asynchronous and synchronous maps features can be analyzed as the sequence of the fluorescence modifications following the order of 360 nm (riboflavin/lumichrome) > 450 nm (riboflavin) > 403 nm (aldehydes and amino acids) > 298 nm (tryptophan) > 322 nm (vitamin A). This is consistent with the milk coagulation procedure observed after adding chymosin, hydrolysis of κ-casein, destabilization of casein micelles and aggregation [39].

<table><caption>**Table 1.** The results of two-dimensional correlation spectroscopy (2DCOS) analysis of synchronous fluorescence spectra during milk coagulation: peak assignments and the signs of the cross-peaks in synchronous spectra and asynchronous spectra (shown in bracket).</caption><thead><tr><th rowspan="2">Peak Position (nm)</th><th colspan="5">Sign of Identified Peak</th></tr><tr><th>298 nm</th><th>321 nm</th><th>360 nm</th><th>403 nm</th><th>450 nm</th></tr></thead><tbody><tr><td>298 nm</td><td>+</td><td>+</td><td>+</td><td>-</td><td>+</td></tr><tr><td></td><td></td><td>(-)</td><td>(+)</td><td>(-)</td><td>(+)</td></tr><tr><td>322 nm</td><td></td><td>+</td><td>+</td><td>-</td><td>+</td></tr><tr><td></td><td></td><td></td><td>(+)</td><td>(-)</td><td>(+)</td></tr><tr><td>360 nm</td><td></td><td></td><td>+</td><td>-</td><td>+</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td>(-)</td></tr><tr><td>403 nm</td><td></td><td></td><td></td><td>+</td><td>+</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td>(+)</td></tr><tr><td>450 nm</td><td></td><td></td><td></td><td></td><td>+</td></tr></tbody></table>

## 4. Conclusions

Molecular fluorescence spectroscopy associated with 2DCOS spectroscopy was used to monitor modifications occurring during the coagulation of different milk formulations based on CM and CaM
mixtures. The method demonstrated that this strategy successfully discriminates milk mixtures and monitors molecular structure modifications occurring during coagulation process.

The dissimilarities among the different formulations are precisely noticed on the synchronous 2DCOS fluorescence spectra. In addition, according to the cross-peak symbols in synchronous and asynchronous spectra, the speed of modification in the fluorescence molecules (riboflavin, protein and vitamin A) during coagulation time corroborated with common coagulation phenomena usually reported during chymosin coagulation. This research gave a promising methodology to discriminate and monitor the coagulation of milk by fluorescence and 2DCOS in an efficient, effective and non-destructive way. This method could also be proposed in the near future to identify adulteration of CaM by CM before and during processing.

**Author Contributions:** O.B. co-conceived and designed the study, performed milk sampling, biochemical and spectral analysis, wrote the paper, analyzed the results and approved the final draft. A.A.-K. co-conceived and designed this study, wrote the paper, analyzed the results, and approved the final draft. F.L. co-conceived and designed this study and approved the final draft. J.S. provided camel milk samples for the project and approved the final draft. S.S. and J.S. provided camel milk samples for the project. E.M.E.H. co-conceived and designed the study and approved the final draft. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the Centre National pour la Recherche Scientifique et Technique (CNRST) of Morocco and the VetAgro Sup Institute located in France.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Faye, B.; Konuspayeva, G. The sustainability challenge to the dairy sector e The growing importance of non-cattle milk production worldwide. *Int. Dairy J.* **2012**, *24*, 50–56. [CrossRef]

2. Farah, Z.; Rettenmaier, R.; Atkins, D. Vitamin content of camel milk. *Int. J. Vitam. Nutr. Res.* **1992**, *62*, 30–33. [PubMed]

3. Sawaya, W.N.; Khalil, J.K. Chemical Composition and Nutritional Quality of Camel Milk. *J. Food Sci.* **1984**, *49*, 744–747. [CrossRef]

4. Khalesi, M.; Salami, M.; Moslehishad, M.; Winterburn, J.; Moosavi-Movahedi, A.A. Biomolecular content of camel milk: A traditional superfood towards future healthcare industry. *Trends Food Sci. Technol.* **2017**, *62*, 49–58. [CrossRef]

5. Lobato-Calleros, C.; Vernon-Carter, E.J.; Hornelas-Uribe, Y. Microstructure and texture of cheese analogs containing different types of fat. *J. Texture Stud.* **1998**, *29*, 569–586. [CrossRef]

6. Lamichhane, P.; Kelly, A.L.; Sheehan, J.J. Symposium review: Structure-function relationships in cheese. *J. Dairy Sci.* **2018**, *101*, 2692–2709. [CrossRef]

7. Boubellouta, T.; Dufour, É. Effects of Mild Heating and Acidification on the Molecular Structure of Milk Components as Investigated by Synchronous Front-Face Fluorescence Spectroscopy Coupled with Parallel Factor Analysis. *Appl. Spectrosc.* **2008**, *62*, 490–496. [CrossRef]

8. Dubayová, K.; Kušnír, J.; Makuchová, V. Synchronous fluorescence fingerprint as the identity card of food products. *Chem. Pap.* **1998**, *52*, 797.

9. Hammami, M.; Rouissi, H.; Salah, N.; Selmi, H.; Al-Otaibi, M.; Blecker, C.; Karoui, R. Fluorescence spectroscopy coupled with factorial discriminant analysis technique to identify sheep milk from different feeding systems. *Food Chem.* **2010**, *122*, 1344–1350. [CrossRef]

10. Boubellouta, T.; Galtier, V.; Dufour, É. Effects of added minerals (calcium, phosphate, and citrate) on the molecular structure of skim milk as investigated by mid-infrared and synchronous fluorescence spectroscopies coupled with chemometrics. *Appl. Spectrosc.* **2009**, *63*, 1134–1141. [CrossRef] [PubMed]

11. Boubellouta, T.; Galtier, V.; Dufour, É. Structural changes of milk components during acid-induced coagulation kinetics as studied by synchronous fluorescence and mid-infrared spectroscopy. *Appl. Spectrosc.* **2011**, *65*, 284–292. [CrossRef] [PubMed]

12. Loudiyi, M.; Aït-Kaddour, A. Evaluation of the effect of salts on chemical, structural, textural, sensory and heating properties of cheese: Contribution of conventional methods and spectral ones. *Crit. Rev. Food Sci. Nutr.* **2019**, *59*, 2442–2457. [CrossRef] [PubMed]
13. Loudiyi, M.; Karoui, R.; Rutledge, D.N.; Lavigne, R.; Montel, M.; Aït-kaddour, A. Contribution of fluorescence spectroscopy and independent components analysis to the evaluation of NaCl and KCl effects on molecular-structure and fat melting temperatures of Cantal-type cheese. *Int. Dairy J.* **2017**, *73*, 116–127. [CrossRef]

14. Loudiyi, M.; Karoui, R.; Rutledge, D.N. Fluorescence spectroscopy coupled with independent components analysis to monitor molecular changes during heating and cooling of Cantal-type cheeses with different NaCl and KCl contents. *J. Sci. Food Agric.* **2018**, *98*, 963–975. [CrossRef] [PubMed]

15. Noda, I. Two-dimensional infrared spectroscopy of synthetic and biopolymers. *Bull. Am. Phys. Soc.* **1986**, *31*, 520–524.

16. Noda, I. Two-dimensional infrared (2D IR) spectroscopy: Theory and applications. *Appl. Spectrosc.* **1990**, *44*, 550–561. [CrossRef]

17. Noda, I. Generalized two-dimensional correlation method applicable to infrared, Raman, and other types of spectroscopy. *Appl. Spectrosc.* **1993**, *47*, 1329–1336. [CrossRef]

18. Ozaki, Y.; Noda, I. Potential of generalised two-dimensional correlation spectroscopy in the near infrared region. *J. Near Infrared Spectrosc.* **1996**, *4*, 85–99. [CrossRef]

19. Wang, Y.; Murayama, K.; Myojo, Y.; Tsenkova, R.; Hayashi, N.; Ozaki, Y. Two-dimensional Fourier transform near-infrared spectroscopy study of heat denaturation of ovalbumin in aqueous solutions. *J. Phys. Chem. B* **1998**, *102*, 6655–6662. [CrossRef]

20. Axelrod, D.; Koppel, D.E.; Schlessinger, J.; Elson, E.; Webb, W.W. Mobility measurement by analysis of fluorescence photobleaching recovery kinetics. *Biophys. J.* **1976**, *16*, 1055. [CrossRef]

21. Koppel, D.E.; Axelrod, D.; Schlessinger, J.; Elson, E.L.; Webb, W.W. Dynamics of fluorescence marker concentration as a probe of mobility. *Biophys. J.* **1976**, *16*, 1315–1329. [CrossRef]

22. Ma, F.; Chen, J.; Wu, X.; Zhou, Q.; Sun, S. Rapid discrimination of Panax notoginseng of different grades by FT-IR and 2DCOS-IR. *J. Mol. Struct.* **2016**, *1124*, 131–137. [CrossRef]

23. Chen, N.; Chen, N.; Li, J.; Cao, C.; Wang, J. Rapid authentication of different ages of tissue-cultured and wild Dendrobium huoshanense as well as wild Dendrobium henanense using FTIR and 2D-COS IR. *J. Mol. Struct.* **2015**, *1101*, 101–108. [CrossRef]

24. Henderson, R.K.; Baker, A.; Murphy, K.R.; Hambly, A.; Stuetz, R.M.; Khan, S.J. Fluorescence as a potential monitoring tool for recycled water systems: A review. *Water Res.* **2009**, *43*, 863–881. [CrossRef]

25. Hur, J.; Jung, K.; Mee, Y. Characterization of spectral responses of humic substances upon UV irradiation using two-dimensional correlation spectroscopy. *Water Res.* **2011**, *45*, 2965–2974. [CrossRef]

26. Hur, J.; Lee, B. Comparing the Heterogeneity of Copper-Binding Characteristics for Two Different-Sized Soil Humic Acid Fractions Using Fluorescence Quenching Combined with 2D-COS. *Sci. World J.* **2011**, *11*, 1865–1876. [CrossRef]

27. Rysstad, G.; Abrahamsen, R.K. Formation of volatile aroma compounds and carbon dioxide in yogurt starter grown in cows and goats milk. *J. Dairy Res.* **1987**, *54*, 257–266. [CrossRef]

28. Noda, I.; Ozaki, Y. Principle of Two-Dimensional Correlation Spectroscopy. In *Two-Dimensional Correlation Spectroscopy: Applications in Vibrational and Optical Spectroscopy*; John Wiley & Sons, Ltd.: The Atrium, Southern Gate, Chichester, West Sussex, UK, 2004; ISBN 9780470012390.

29. Lakowicz, J.R. Quenching of fluorescence. In *Principles of Fluorescence Spectroscopy*; Springer Science+Business Media, LLC: New York, NY, USA, 1999; pp. 237–265. ISBN 9781475730630.

30. Herbert, S.; Riaublanc, A.; Bouchet, B.; Gallant, D.J.; Dufour, E. Fluorescence Spectroscopy Investigation of Acid- or Rennet-Induced Coagulation of Milk. *J. Dairy Sci.* **1999**, *82*, 2056–2062. [CrossRef]

31. Dufour, E.; Devaux, M.F.; Fortier, P.; Herbert, S. Delineation of the structure of soft cheeses at the molecular level by fluorescence spectroscopy—Relationship with texture. *Int. Dairy J.* **2001**, *11*, 465–473. [CrossRef]

32. Dufour, E.; Lopez, C.; Riaublanc, A.; Mouhous Riou, N. La spectroscopie de fluorescence frontale: Une approche non invasive de la structure et des interactions entre les constituants des aliments. *Agoral* **1998**, *10*, 209–215.

33. Fox, J.J.B.; Thayer, D.W. Radical oxidation of riboflavin. *Int. J. Vitam. Nutr. Res.* **1998**, *68*, 174–180. [PubMed]

34. Christensen, J.; Becker, E.M.; Frederiksen, C.S. Fluorescence spectroscopy and PARAFAC in the analysis of yogurt. *Chemom. Intell. Lab. Syst.* **2005**, *75*, 201–208. [CrossRef]

35. Kikugawa, K.; Beppu, M. Involvement of lipid oxidation products in the formation of fluorescent and cross-linked proteins. *Chem. Phys. Lipids* **1987**, *44*, 277–296. [CrossRef]
36. Wold, P.J.; Veberg, A.; Lundby, F.; Nikolai, A.; Moan, J. Influence of storage time and color of light on photooxidation in cheese : A study based on sensory analysis and fluorescence spectroscopy. *Int. Dairy J.* **2006**, *16*, 1218–1226. [CrossRef]

37. Abbas, K.; Karoui, R.; Aït-Kaddour, A. Application of synchronous fluorescence spectroscopy for the determination of some chemical parameters in PDO French blue cheeses. *Eur. Food Res. Technol.* **2012**, *234*, 457–465. [CrossRef]

38. Sun, S.Q.; Zhou, Q.; Qin, Z. *Atlas of Two-Dimensional Correlation Infrared Spectroscopy for Traditional Chinese Medicine Identification*; Chemical Industry Press: Beijing, China, 2003.

39. Dalgleish, D.G.; Brinkhuis, J.; Payens, T.A.J. The Coagulation of Differently Sized Casein Micelles by Rennet. *Eur. J. Biochem.* **1981**, *119*, 257–261. [CrossRef]

<figure><img src="image_6.png" /></figure>

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).