Open Access Article
Harshal V
Barkale
a,
Bappa
Maiti
b and
Nilanjan
Dey
*a
aDepartment of Chemistry, Birla Institute of Technology and Science Pilani, Hyderabad campus, Hyderabad, Telangana 500078, India. E-mail: nilanjandey.iisc@gmail.com; nilanjan@hyderabad,bits-pilani.ac.in
bMichael Smith Laboratories, University of British Columbia, Vancouver V6T 1Z4, Canada
First published on 22nd July 2025
Herein, we report the design and synthesis of two pyrimidine derivatives with different donor substituents, pyrene (1) and anthracene (2), that show aggregation in an aqueous medium. The aggregates of compound 1 were found to be pH-sensitive as well as thermoresponsive. The energy-optimized structure of compound 1 revealed a nearly planar conformation with a dihedral angle of ∼0.2° between the donor and the acceptor moieties, which appeared to be ∼50.8° for compound 2. A 2.5-fold and 5-fold turn-on fluorescence response was observed after the addition of melamine for compounds 1 and 2, respectively. Importantly, compound 1 was successfully applied for the quantification of melamine in real milk products, with high recovery rates, and the results were validated using liquid chromatography-mass spectrometry (LC-MS). Accurate estimation of melamine levels was achieved with a limit of detection (LOD) of 0.8 ppm, well below regulatory thresholds. Furthermore, we developed chemically modified paper strips capable of detecting melamine directly in adulterated milk samples without any instrumentation, offering a practical and low-cost on-site screening method. To gain mechanistic insight, molecular dynamics (MD) simulations were employed to analyze the stability, interaction energies, and aggregation tendencies of the compounds. The self-assembled structures of 1 after 120 ns showed an anti-parallel arrangement with the pyrene moiety of one molecule π-stacked with the pyrimidine unit of another molecule. In contrast, compound 2 exhibited face-to-face π-stacking among anthracene moieties and multiple hydrogen bonds (1.8–2.2 Å) among pyrimidine residues. Finally, various machine learning models were used to predict melamine intensity from concentration vs. fluorescence data. The linear regression model (R2 = 0.9959) delivered the best performance, reinforcing the efficiency of linear approaches in this dataset.
In such cases, integration of machine learning-based algorithms to analyse the data set obtained from traditional solution-based techniques would be particularly interesting, as they can improve the analytical efficacy of the system by enhancing the sensitivity of probe molecules, accuracy, and real-time processing etc. ML algorithms manage large datasets, filter noise, and improve analyte detection in complex environments, making them crucial for environmental monitoring and diagnostics applications. ML-driven sensors offer predictive analysis and personalized solutions, while also reducing costs and optimizing design, establishing machine learning as a key force in developing scalable, efficient sensing technologies.5 Furthermore, the incorporation of polyaromatic scaffolds into the molecular backbone will result in building up self-assembled architectures. These moieties, through non-covalent interactions such as hydrophobic interactions or π–π stacking, guide the self-assembly of molecules into ordered structures. The extent of self-assembly can affect the morphology of the materials, surface area, and arrangement of functional groups, directly enhancing the sensing efficiency by improving signal transduction and analyte recognition.6 Also, it has often been observed that probe molecules in close-packed self-assembled structures can provide a suitable binding cleft for target analytes through multipoint interactions.7
In this paper, we have designed two easily-synthesizable pyrimidine derivatives with polyaromatic scaffolds, namely anthracene and pyrene, as signalling moieties that can form self-assembled nanostructures in an aqueous medium. It was found that both compounds showed “turn-on” fluorescence response exclusively in the presence of melamine, presumably through the formation of a complementary hydrogen bonding network.8 Despite such similarities, the degree of response (sensitivity of the system) was found to be dependent on the electronic nature, conformation and self-assembly behavior of the probe molecules.
Upon spectroscopic investigation, we would like to utilize this system for the determination (detection and quantification) of melamine in various real-life samples (milk and other milk products). Additionally, we aim to develop chemically modified coated paper strips for on-site detection of melamine.
The presence of polycyclic aromatic scaffolds such as pyrene and anthracene improves the self-assembly properties of the molecules via π–π stacking and hydrophobic interactions. Moreover, these moieties are known for their susceptibility to the microenvironment. Changes in local pH, temperature, viscosity, etc. can affect the degree of aggregation, as noted by the distinct monomer-to-excimer ratio. The incorporation of the pyrimidine unit in the molecule is crucial as it can form a long-range supramolecular architecture via extensive donor–acceptor type complementary hydrogen bonding interaction.
Moreover, connecting such electron-deficient pyrimidine units with electron-rich pyrene/anthracene moieties will enhance charge transport within the molecular framework, leading to improved sensing performance.
Furthermore, FMO (Frontier molecular orbital) analysis of compounds 1 and 2 indicated that the HOMO orbitals were mostly concentrated on the polyaromatic scaffolds (pyrene for 1 and anthracene for 2), while the LUMO orbitals were located on the pyrimidine site along with donor moiety (Fig. 1c).12 Such a kind of localization of FMOs at different parts of the probe molecules indicated the charge-transfer (CT) nature of the system, where the electron density can be transferred from the donor (pyrene or anthracene) to the acceptor (pyrimidine) moiety upon excitation.13 The dipole moment of compound 1 (4.71 D) is found to be greater than that of compound 2 (2.78 D). The greater dipole moment in the pyrene-pyrimidine moiety implies a stronger interaction between the donor (pyrene) and acceptor (pyrimidine moiety).14
This can enhance the intramolecular charge transfer (ICT) properties, making the molecule more efficient in electronic and optical applications.15 The electrostatic potential map (ESP) provides insights into the electron density distribution and identifies the most nucleophilic and electrophilic regions of the molecule. The blue-colored surface indicated a high electron density region, while the low electron density area was defined by a red-colored area. On the contrary, areas on the molecular surface where the enclosed negative electron charge and the positive core charge balanced each other (overall neutral charge), are mostly designated by green. From this, we can conclude that the carbonyl groups (C
O) present in the pyrimidine units are electron-deficient and act as hydrogen bond acceptors, while the –NH groups are electron-rich, serving as hydrogen bond donors.16
The UV-visible spectra of compound 1 exhibited two distinct absorption bands in the 365 (2 × 104 M−1 cm−1) and 440 nm (4 × 104 M−1 cm−1) regions in a THF medium, while in the aqueous medium, compound 1 experienced a hypochromic shift with substantial broadening (Fig. 2d). Notably, residual absorbance observed in the longer wavelength region was presumably due to the Mie scattering effect, caused by the presence of colloidal nanoscopic aggregates in the water medium. Additionally, the red-shifted absorption maxima (470 nm) in the aqueous medium indicated the charge-transfer nature of the aggregates, induced by hydrogen bonding interactions.18,19
Excitation spectra were recorded at both the 385 nm and 472 nm emission maxima in a water medium, revealing a distinct signature, indicating the presence of at least two different photoactive ground state species20 (Fig. 2c). The spectrum corresponding to the 385 nm emission band appeared broad and blue-shifted compared to the excitation spectrum at the 472 nm emission band. The excitation spectra correspond to the 385 nm emission band, resembling the excitation spectrum of the compound observed in the THF medium. This indicated that the emission band centered at 385 nm could have originated from monomer species. Thus, it can be argued that the broad and red-shifted peak observed at ∼472 nm could be assigned to aggregated species.21,22
The dynamic light scattering experiment suggested that compound 1 formed aggregates with an average hydrodynamic diameter of 308 ± 5.8 nm in the aqueous medium (Fig. 3c), while the analysis of the FE-SEM images indicated oval-shaped morphology (Fig. 3c inset). Additionally, time-correlated single-photon counting (TCSPC) experiments were conducted for compound 1 in both water and THF medium to understand the stability and decay kinetics of excited-state photoactive species (Fig. 3d). In the THF medium, compound 1 exhibited multiexponential decay with an average lifetime of 5.66 ns, while in the water medium, it showed a relatively shorter lifetime (3.88 ns). This shorter decay time could be attributed to the formation of π-stacked self-assembled structures in the aqueous medium.23
The pH and thermoresponsive properties highlight the sensor's environmental sensitivity, making it suitable for use in varied sample conditions and offering tunable performance in real-world applications. A temperature-variation study of compound 1 in the aqueous medium showed an enhancement of fluorescence intensity along with a 5 nm blue shift in emission maxima (Fig. 3a). This could be attributed to the association of aggregates, induced by thermal agitation at high temperatures.24 At high temperature, there could be a change in the nature of the aggregates due to the dissociation of the hydrogen-bonding network. The π-stacked self-assemblies, formed via hydrophobic interactions at high temperatures, showed enhancement of excimer emission. Additionally, the blue-shifted maxima observed at elevated temperature indicated a hydrophobic microenvironment around the probe molecules. Fluorescence spectra of compound 1 at different pH conditions exhibited distinct features (Fig. 3b). At acidic pH (pH 4), the pyrimidine unit is protonated, which promotes intermolecular hydrogen bonding and π–π stacking of the pyrene units. This leads to increased aggregation and consequently, a higher intensity of the aggregation band relative to the monomer band. As the pH increases to neutral (pH 7.4), partial deprotonation occurs, reducing hydrogen bonding and aggregation. This results in a relative increase in the monomer emission band, reducing the prominence of the excimer band. At basic pH (9.4), complete deprotonation introduces significant electrostatic repulsion, disrupting aggregation and stabilizing the monomeric form. This leads to a further increase in the monomer emission band relative to the aggregation band. Thus, it could be concluded that the changing microenvironment (pH and temperature) showed a significant influence on the aggregation of compound 1.25
Interestingly when the same amount of melamine was added to compound 2 the band at 470 nm was increased approximately 5-fold (Fig. 5d). Compound 2 shows a better response (approximately 4 times) as compared to compound 1 (Fig. 5c). This may be attributed to a different conformation as compared to pyrene, which allows it to have stronger interactions with melamine through π–π stacking. This enhances the efficiency of energy transfer, leading to a more noticeable fluorometric change.27,28 Given that selectivity is crucial in real-life sample analysis, the response of compound 1 to other competitive analytes was also examined (Fig. 4d). The extents of interaction with other analytes, as defined by changes in FL intensity at the 470 nm band, were found to be significantly less compared to melamine. The enhanced selectivity arises primarily from the complementary hydrogen bonding interactions between the pyrimidine unit of the probe and the triamine structure of melamine. Specifically, melamine contains three exocyclic –NH2 groups capable of acting as both hydrogen bond donors and acceptors, which align precisely with the C
O and –NH groups of the pyrimidine scaffold to form stable, directional, and multivalent hydrogen-bonded assemblies. In contrast, common milk components such as sugars (e.g., lactose, glucose, etc.) primarily form weak hydrogen bonds with –OH groups, lacking the defined geometry and complementarity required for strong probe-analyte complexation. Proteins are too large and structurally complex, with interactions that are often sterically hindered and non-specific. This indicated the robustness of the present system in estimating melamine in real-life samples. In the time-correlated single-photon counting (TCSPC) experiment, the average fluorescence lifetime of compound 1 monitored at 470 nm increased from 3.88 to 3.93 ns in the presence of melamine (Fig. 4c). This nearly unchanged lifetime suggests that the intrinsic deactivation pathways, both radiative and non-radiative, of the excited state remain largely unaffected after the addition of melamine.29
To evaluate the effects of the local surroundings, we have examined the sensing efficacy of the probe molecules at different temperatures and pH conditions. pH variation studies revealed that at a neutral pH (7.4), the maximum change was observed after the addition of melamine (Fig. 4a), though melamine-driven changes in fluorescence signal were observed in a wide range of pH (2–10), the maximum response was noticed at neutral pH conditions (6–8), because at neutral pH, both the melamine and pyrimidine moiety are in their optimal protonation states, facilitating efficient hydrogen bonding.30
This results in stabilized excited states and an enhanced fluorescence response. Probably at extreme pH conditions (basic as well as acidic), the shift in tautomeric equilibrium diminishes the hydrogen bonding efficiency of the pyrimidine unit, leading to a mild response towards melamine. To determine the stability of the hydrogen-bonded adduct with melamine, we recorded the fluorescence spectra at different temperatures (20 °C to 90 °C) in the aqueous medium (Fig. S4, ESI†). With an increase in temperature the intensity at the excimer band increases; this may be attributed to the partial weakening of hydrogen bonds between the melamine and pyrimidine residues, which might reduce the rigidity of the system.31 This increased flexibility could promote aggregation of the pyrene units, enhancing excimer formation.
000 synthetic data points by interpolating between the minimum and maximum values of the original dataset.
Out of many regression algorithms, we have selected support vector regression (SVR) with both linear and quadratic kernels, linear regression, and AdaBoost regression. The rationale for choosing these models stems from their diverse capabilities to capture complex relationships in the data. SVR, in particular, is adept at handling non-linear interactions commonly observed in spectroscopic data.32 Each model was trained on the augmented dataset, with the training process involving fitting the models to the synthetic intensity data derived from the concentration values. We evaluated the performance of these models using several statistical metrics, including mean squared error (MSE), coefficient of determination (R2), and root mean squared error (RMSE).33
These metrics provide quantitative assessments of the model's predictive accuracy, with lower values indicating superior performance as shown in Table 1.
| Method | Mean squared error | Root mean squared error | R-squared |
|---|---|---|---|
| SVR linear | 0.0042 | 0.0649 | 0.9724 |
| Linear regression | 0.0006 | 0.0249 | 0.9959 |
| SVR quadratic | 0.0103 | 0.1015 | 0.9326 |
| AdaBoost regression | 0.0008 | 0.0295 | 0.9942 |
To visualize the predictive capabilities of each model, we generated plots of actual versus predicted intensity values. These plots depicted the relationship between the observed and predicted intensities, complete with regression lines that illustrate the model fit to the data. The analysis revealed varying degrees of accuracy, with SVR quadratic (Fig. 6c) and AdaBoost regression (Fig. 6d) consistently demonstrating the best performance in predicting normalized intensity.
![]() | ||
| Fig. 6 (a) Plots of actual versus (Fig. 1) predicted intensity values of (a) SVR linear, (b) linear regression, (c) SVR quadratic, and (d) AdaBoost regression. | ||
The linear regression model (R2 = 0.9959) stands out with the best performance metrics (Fig. 6b), reinforcing that a linear approach is highly effective for predicting normalized intensity based on concentration in this dataset.34 The simplicity and interpretability of linear regression make it the most suitable model given the strong linear relationship. The other models (Fig. 6a), while also performing well, do not surpass the linear model's effectiveness, and more complex models like SVR quadratic (Fig. 6c) may introduce unnecessary complexity without providing significant benefits.
By employing machine learning techniques in this research, we underscore the significant potential for these methodologies to enhance predictive analysis in chemical sensing applications. Our findings contribute to the growing body of literature highlighting the effectiveness of advanced computational methods in analytical chemistry, particularly in accurately predicting intensity linked to melamine concentration.
O (of pyrimidine unit) stretching frequency to the lower energy region, from 1728 to 1660 cm−1 (Fig. 7e), which indicated the formation of hydrogen bonds with melamine –NH units. Additionally, hydrogen bonding interaction involving the pyrimidine–NH group was also evident from FT-IR analysis, as the corresponding stretching band (3390 cm−1) became broad.35,36 Though all the spectroscopic studies were performed in an aqueous medium, we have conducted 1H-NMR studies in DMSO-d6 due to the limited solubility of the probe molecules in D2O. The 1H NMR spectra of compound 1 with melamine in DMSO-d6 were recorded to further investigate the interaction mode between the probe and melamine (Fig. 7b). The peak labelled as ‘a’ shifted from 9.1 ppm to 10.6 ppm, and the protons labelled as ‘b’ and ‘c’ nearly disappeared upon interaction with melamine. This is indicative of hydrogen bonding interactions between melamine and the pyrimidine moiety.37
As expected, the aromatic proton adjacent to the pyrimidine moiety exhibited a downfield shift from approximately 8.4 to 8.6 ppm upon interaction with melamine, indicating a change in the chemical shift due to the formation of a complex. In contrast, the protons located farther from the pyrimidine unit did not display significant shifts, suggesting minimal change in their environment. Additionally, upon the addition of melamine, the –NH resonances at around 11.2 and 11.5 ppm were no longer observed, likely due to hydrogen bonding interactions between the pyrimidine residue and melamine, which resulted in the disappearance of these signals.38,39
To evaluate the hydrogen-bonding interaction between the probe and melamine, we examined the interaction of the compound with melamine in the presence of a known chaotropic agent, urea. The fluorescence response of compound 1 towards melamine was reduced significantly in the presence of 5 M urea (Fig. 7d). This observation could indirectly prove the essential role of hydrogen bonding interactions in the formation of supramolecular co-assembly with melamine (Fig. 1b).40
Additionally, dynamic light scattering (DLS) experiments revealed the formation of larger colloidal particles with an average diameter of 638.2 ± 5.7 nm in the presence of melamine (Fig. 7c). Consequently, FESEM images of the 1 melamine adduct displayed larger spherical aggregated structures (Fig. 7c inset). Changes in aggregation were also evident from FL-lifetime analysis. These spectroscopic results suggest that the formation of hydrogen bonds between compound 1 and melamine leads to notable changes in both the IR and NMR spectra, providing strong evidence for the mode of interaction and the resultant supramolecular assembly formation, which is responsible for the observed turn-on fluorescence response (Fig. 7a).
Furthermore, we attempted to study the interaction of 1 and 2 with melamine. In both cases, the addition of melamine resulted in more aggregated structures with the presence of a greater number of hydrogen bonds (Fig. 8c and d).
The polyaromatic residues were found to be in close proximity, which could lead to an increase in the excimer emission. This was also supported by experimental observation, where we witnessed a turn-on fluorescence response for both 1 and 2 upon the addition of melamine. A closer view indicated that the pyrene units in 1.melamine assembly mostly preferred the slipped conformation. At the same time, the anthracene moieties in the case of 2 showed a well-π-stacked conformation, though some conformations are T/L-shaped. The larger fluorescence response (F/F0) for compound 2 with melamine might be attributed to such differences in the arrangements of the polyaromatic residues in the aggregated state. The average number of hydrogen bonds in the case of 1 was found to be ∼22–25 after 50 ns, which appeared to be ∼12–15 for compound 2 (Fig. 8f). These hydrogen bonds were mostly observed between the pyrimidine residues and melamine units. The differences in the conformations of compounds 1 and 2 (planar vs. twisted) might be one of the contributing factors to dictate the hydrogen bonding interaction with melamine.43 A significant reduction in the electrostatic interaction energy for both 1 (−411.42 kJ mol−1) and 2 (−158.48 kJ mol−1) in the presence of melamine suggested that the self-assembly formation was mostly dominated by hydrogen bonding and π-stacking interaction. Also, smaller VdW energy values (−163.73 kJ mol−1 for 1 and −102.18 kJ mol−1 for 2) observed in both cases suggested the formation of large-sized aggregates with the addition of melamine, which was also confirmed by the DLS experiment.
The UV-visible spectra of compounds 1 and 2 in an aqueous medium exhibit distinct absorption features corresponding to their electronic transitions (Fig. S2, ESI†). A broad absorption band centered around 470 nm for compound 1 is indicative of an n–π* transition, while for compound 2, this band appears at approximately 440 nm. The observed 30 nm blue shift in compound 2 can be attributed to its less extended π-conjugation system compared to the planar structure of pyrene in compound 1.44
As already discussed, compound 2 adopts a twisted molecular conformation, whereas compound 1 assumes a planar arrangement. This structural difference significantly influences their self-assembly behavior, as demonstrated by dynamic light scattering (DLS) and field-emission scanning electron microscopy (FE-SEM) analyses (Fig. 5b). The DLS studies showed the formation of particles with a larger hydrodynamic diameter (485.7 ± 18.2 nm) for compound 2 than that of compound 1. Compound 2 exhibits a higher degree of aggregation than compound 1, which can be attributed to its twisted molecular arrangement, promoting enhanced intermolecular interactions and aggregation. The difference in the charge transfer for both compounds also results from the difference in the dihedral angle.45 The fluorescence spectra of compound 2 in THF and aqueous medium display similar features to those observed for compound 1 (Fig. 5a). In THF, compound 2 exhibits characteristic monomer emission, indicative of its molecularly dispersed state. In contrast, in an aqueous medium, compound 2 demonstrates fluorescence associated with aggregation, consistent with the behavior observed for compound 1. In the time-correlated single-photon counting (TCSPC) experiment, the average fluorescence lifetime of compound 1 monitored at 470 nm is 3.88 ns, and for compound 2 the average lifetime is around 3.45 ns (Fig. S5, ESI†). This extended lifetime may be attributed to the enhanced aggregation of compound 2 in the aqueous medium.
Overall, the distinct conformation, extended conjugation, and charge transfer in the probe molecules result in different analytical performances for the melamine.
Furthermore, fluorescence spectral analysis revealed no significant changes in the emission profile of the probe molecule after the pre-treatment process. This observation confirms the absence of background interference attributable to the chloroacetic acid treatment, ensuring the integrity of the analytical procedure and the specificity of melamine detection in the treated milk samples.
The samples were prepared and analyzed according to the procedure as follows.48 A 4.0 mL raw milk sample was placed into a 10 mL centrifuge tube, and 1.2 mL of 300 g L−1 trichloroacetic acid was added and mixed with a vortex for 5 min to deposit protein in the sample matrix. The mixture was centrifugated at 3500 rpm for 10 min. The supernatants were transferred into another centrifuge tube and adjusted to pH 7.4 with a small amount of 6 M NaOH, then filtered with a 0.22 μm filter. This filtrate was directly titrated with probe 1 for the detection of melamine (Fig. 9b). The fluorescence intensity of compound 1 at the 470 nm band increased with increasing volume of spiked milk samples (Fig. 9c). A nearly 2-fold increase in the emission intensity was observed at 470 nm after the addition of the 450 μM melamine spiked milk sample, which is almost in the same range as observed during titration studies. The minimum detectable concentration (LOD) of melamine was found to be less than the permissible limit. Calculations of percentage recovery values clearly stated that the present method could achieve a quantitative estimation of melamine.
This linear fluorescence response suggests that the present method could be utilized for quantitatively detecting melamine even in milk and other dairy products. Across all cases, the recovery values range between 96.2 and 105.7% with a relative standard deviation (RSD) of less than 5% (Fig. 9e). Furthermore, the same protocol was employed for the determination of melamine in other milk products such as milk powder and infant formula. In all cases, the estimated values were found to be in the same range as other reported values (Fig. 9d). Furthermore, the melamine content of each sample was also independently estimated by the LC-MS method.49 An excellent agreement was observed between the values obtained from this method and the method described in this work (Table S1, ESI†).
A comparative overview of conventional melamine detection techniques, including their detection limits, analytical methods, and limitations, is provided in Table S2 (ESI†) to highlight the advantages of our ML-assisted fluorescent probe.
O and –NH groups of the pyrimidine units and melamine's –NH groups altered the charge-transfer and aggregation behavior of the probes. The method was successfully applied for quantitative melamine detection in milk samples, with validation using alternative analytical methods.
Molecular dynamics (MD) simulations revealed that compound 1 forms anti-parallel π-stacked self-assemblies via pyrene–pyrimidine interactions, whereas compound 2 forms face-to-face anthracene stacking with multiple hydrogen bonds (1.8–2.2 Å) among pyrimidine residues, contributing to its enhanced stability. Additionally, we developed chemically modified paper strips for rapid, on-site melamine detection in milk. This portable method eliminates the need for complex instruments, making it accessible to non-specialists. Future efforts will focus on extending detection to other food adulterants and integrating the sensor into a user-friendly, portable device for field use.
Footnote |
| † Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d5ma00470e |
| This journal is © The Royal Society of Chemistry 2025 |