Arslan Yousaf‡
a,
Muhammad Hafeez Ullah‡a,
Haq Nawaz*a,
Muhammad Irfan Majeed
*a,
Nosheen Rashidb,
Abdulrahman Alshammaric,
Norah A. Albekairic,
Arslan Alia,
Munawar Hussaina,
Abu Bakar Salfia,
Muhammad Aamir Aslam*d,
Kinza Idreesd and
Allah Dittae
aDepartment of Chemistry, University of Agriculture Faisalabad, Faisalabad (38000), Pakistan. E-mail: haqchemist@yahoo.com; irfan.majeed@uaf.edu.pk
bDepartment of Chemistry, University of Education, Faisalabad Campus, Faisalabad (38000), Pakistan
cDepartment of Pharmacology and Toxicology, College of Pharmacy, King Saud University, Post Box 2455, Riyadh, 11451, Saudi Arabia
dInstitute of Microbiology, Faculty of Veterinary Sciences, University of Agriculture Faisalabad, Faisalabad (38000), Pakistan. E-mail: aamir.aslam@uaf.edu.pk
eInstitute for Experimental Molecular Imaging, RWTH Aachen University Hospital, Aachen 52074, Germany
First published on 10th February 2025
Acinetobacter baumannii (A. baumannii) is an emerging Gram-negative nosocomial pathogen responsible for infection on a global scale. It has the ability to develop biofilms on different surfaces, especially abiotic surfaces, which is considered a major contributor of its pathogenicity. Surface-enhanced Raman spectroscopy (SERS) holds great potential as an effective method for identifying and characterizing the biochemical composition of biofilm-forming species. In this study, cell mass samples from different strains of A. baumannii, categorized based on their biofilm-forming ability (strong, medium and non-biofilm forming) using a 96-well microtiter plate assay (MTP), were analyzed by SERS. The identified spectral features of the SERS spectra were used to characterize bacterial strains capable of producing biofilms. Silver nanoparticles (Ag-NPs) served as the SERS substrate to differentiate biofilm-forming strains of A. baumannii. Chemometric tools, such as principal component analysis (PCA) and partial least squares discriminant analysis (PLS-DA), were employed for the classification and differentiation of SERS spectra from bacterial strains with varying biofilm-producing capacities, achieving 100% sensitivity, 94.3% specificity, and an area under the curve (AUC) value of 0.81 through Monte Carlo cross-validation. Furthermore, K-fold (Leave-K-out cross-validation (LKOCV)) was applied to verify the robustness of the PLS-DA model, and the AUC value was found to be 0.90, with a sensitivity of 100% and specificity of 98%. These results demonstrate that the PLS-DA model is highly effective for the differentiation and classification of bacterial strains with varying capacities for biofilm production.
Biofilms are complex microbial communities encased in an extracellular matrix (ECM) composed of extracellular polymeric substances (EPS), which include proteins, DNA and polysaccharides.8,9 These communities can comprise various microorganisms, including protozoa, fungi, algae and both Gram-positive and Gram-negative bacteria. The EPS matrix plays a key role in biofilm formation by facilitating cell-to-cell communication (quorum sensing), protecting bacteria from environmental stresses, and preventing desiccation, as water can constitute up to 98.5% of the biofilm mass.10 Recent research highlights a communication process known as auto-induction, where cells exchange signals within and between bacterial species. The Gram-negative bacteria typically use acylated homoserine lactones as auto-inducers, while the Gram-positive bacteria employ processed oligopeptides for communication.11
Biofilms are typically found on aortic and pulmonary valves, intrauterine devices, central venous catheters, contact lenses, urinary catheters, and water lines in dental units. Among the major causes of dangerous infections in patients, one is biofilms that grow on prosthetic heart valves or surrounding tissue. Microorganisms have been shown to smoothly attach and grow on the stitching cuff fibers of valves.12 The phenomenon of biofilm development in A. baumannii is influenced by various factors, such as the existence of antibiotic resistance genes, the environment in which bacteria grow, and the density of the bacterial cells.6 Quorum sensing is employed by both Gram-positive and Gram-negative bacteria to regulate several physiological functions, including biofilm formation, pathogenicity, antibiotic synthesis, motility and conjugation sporulation.11 Quorum sensing controls common steps in the development of biofilms, which include initial attachment to the biotic or abiotic surface, formation of a microcolony, maturation and formation of the architecture, and lastly, detachment of the biofilm.13 Due to its ability to develop resistance to all currently available antimicrobial agents, A. baumannii has been classified as a “red alert” human pathogen.14A. baumannii has been added to the drug-resistant and antimicrobial resistance bacteria research priority list of the World Health Organization (WHO), due to its rising multidrug-resistant strains, mortality and high morbidity.15
Different techniques, including confocal laser scanning microscopy (CLSM),16 scanning electron microscopy (SEM),17 X-ray crystallography,18 transmission electron microscopy (TEM),19 Fourier transform infrared spectroscopy (FTIR)20 and nuclear magnetic resonance (NMR) spectroscopy,21 are widely used for the characterization of bacterial biofilms. Notably, CLSM, SEM, TEM and AFM are the most important microscopic techniques for the visualization of biofilm morphology from the nanometer scale to the micrometer scale on sample surfaces. However, these techniques are limited due to their low axial resolution, time-consuming sample preparation and high cost.22–24 Moreover, TEM, FTIR, X-ray crystallography and NMR provide information associated with biochemical composition; however, they have some drawbacks. TEM is a pre-treatment technique involving procedures, such as freezing and fixations, which may modify the biofilm structure.22,23 FTIR, NMR and IR are not suitable for aqueous environments, time-consuming experiments, and need extensive sample preparation.25–27 X-ray crystallography requires a large quantity of powdered sample, and its high energy can destroy the sample.28,29 In order to overcome these problems, the development of a rapid identification method is essential for the differentiation and classification of bacterial strains having different capacities of producing biofilms.30
Raman spectroscopy is a practical technique that can be used to overcome the limitations of the abovementioned techniques.31 This is a non-destructive, suitable analytical technique used for the rapid identification and characterization of samples in aqueous environments. However, it faces challenges due to its weak signal, which can cause sample degradation at high laser intensities and long exposure time.32 SERS has effectively resolved the challenges with Raman scattering by enhancing the Raman signal by up to 1014 through adherence of sample molecules in close proximity to a roughened surface of metallic nanostructures.33 Silver nanoparticles (Ag-NPs) with suitable size and morphology can serve as a SERS platform to increase the detection sensitivity.34 Due to its great sensitivity, low cast, and ease of sample preparation, SERS has great potential for quickly and accurately identifying bacterial species and the molecular structure of bacteria, including carbohydrates, proteins, nucleic acids and lipids, and can be responsible for their characteristic SERS spectral features.35
In the present study, SERS has been employed using Ag NPs as a substrate for the rapid, sensitive and cost-effective discrimination among three different bacterial strains having different capacities for producing biofilms, which are categorized as strong, medium and non-biofilm forming on the basis of a 96-well microtiter plate assay (MTP). Moreover, SERS spectral data of different biofilm forming and non-biofilm forming strains are compared to identify specific SERS features that are directly or indirectly related to biochemical changes that occur due to the development of biofilm matrix and the unique features that correspond to individual strains. Moreover, different chemometric techniques such as PCA are used to identify the individual strains and PLS-DA to classify the biofilm forming and non-biofilm forming strains.
PLS-DA is a supervised chemometric method that is applied in the current study on SERS spectra for calibration and validation of the preprocessed spectral dataset in order to categorize the respective bacterial strains into non, medium, and strong based on their biofilm forming ability.38 Monte Carlo and K-fold cross validation through LOOC (Leave one sample out cross validation) were used to build the PLS-DA model. Monte Carlo and LKOCV were generally implemented for the purpose to ensure maximum statistical relevancy, and also to select the minimum number of latent variables, but avoiding the over fitting under the given spectral data. It can also be used to demonstrate that the analysis is not biased, as all spectra collected from a sample are either included in the calibration or the test sets, but cannot be represented in both. Particularly, calibration was performed on 60% of data sets by selecting it randomly, and the remaining 40% of data sets were utilized for prediction.43 In order to determine the specificity and sensitivity of the classification of SERS data sets, the ROC curve was used to display the results of the SERS spectral group classification as false positive, true positive, false negative, and true negative.37
![]() | ||
Fig. 1 Mean plot of the SERS spectra of non, medium and strong biofilm-forming strains of A. baumannii. |
Raman shift (cm−1) | Assignments | Components | References |
---|---|---|---|
a ν = Stretching (vs = symmetric, vas = antisymmetric), δ = bending/scissoring, ρt = twisting. | |||
326 | β (C–C–O) ring vibration | Carbohydrates | 37 |
343 | C–C–O δ, C–OH ρt | Carbohydrates | 44 |
393 | S–S ν | Carbohydrates | 45 |
420 | (C–C–C deformation) | Carbohydrates | 38 |
467 | C–C–O δ, C–OH ρt, β(CCC), β(CCO), β(OCO) | Carbohydrates | 44 and 46 |
516 | S–S ν, disulfide | Proteins | 47 and 48 |
633 | C–C twisting of guanine ring region | DNA | 49 and 47 |
679 | Thymine | DNA | 50 |
699 | (Guanine (COO deformation) | DNA | 51 and 52 |
720 | Adenine, glycosidic ring mode of adenine | DNA | 53 and 54 |
729 | Adenine | DNA | 29 |
747 | (C–S) ν | Proteins | 55 |
762 | Ring breathing of tryptophan | Proteins | 56 |
861 | C–C ν, C–O–C 1,4-glycosidic linkage | Proteins | 57 and 53 |
CCH deformation in ring breathing of tyrosine | |||
953 | υ(CO) and C–N of proteins | Proteins | 58 |
990 | C–C ring breathing in phenylalanine | Proteins | 59 and 60 |
1030 | C–H δ | Proteins | 61 |
1070 | C–N ν. (protein) | Proteins | 62 |
1097 | PO2− vs | DNA/RNA | 29 |
1130 | (δ (CC, CO, –COH) | Carbohydrates | 38 |
1327 | δ(CH) | Carbohydrates | 7 |
1395 | COO−1vs | Proteins | 63 |
1445 | CH2 δ vibration (lipid band) | Lipids | 46 |
1550 | Tryptophan | Proteins | 64 |
1693 | Amide I α-helix | Proteins | 49 |
The major differentiating SERS bands in the mean spectral plot of bacterial strains include 420 (C–C–C deformation), 467 (C–OH ρt, β(CCC), β(CCO), β(OCO)), 516 (S–S ν, disulfide), 633 (the guanine ring region and the tyrosine amino acid chain include aromatic rings with a C–C configuration), 699 (guanine (COO deformation), 720 (adenine), 729 (the glycosidic ring conformation of adenine), 762 (ring breathing of tryptophan), 953 (υ(CO) and C–N of proteins), 990 (C–C ring breathing of phenylalanine), 1030 (C–H δ), 1070 (C–N ν of proteins), 1130 (CC, CO, –COH), 1327 (CH), 1395 (COO−1vs), 1445 (CH2δ vibration (lipid band)) and 1693 cm−1 (amide I α-helix).
The dashed lines are used to represent the intensity-based differences of the SERS bands in the non-biofilm and high biofilm-forming strains, indicating that all sample strains belong to the same species of A. baumannii; hence, they have some similar biomolecular features. These SERS features include 326 (β (C–C–O) ring vibration), 343 (C–C–O δ, C–OH ρt), 393 (S–S ν), 679 (thymine), 747 ((C–S) ν), 861 (C–C ν, C–O–C 1,4-glycosidic linkage, CCH deformation in the ring breathing of tyrosine, namely in the CCH group; Gram-positive bacterium cell wall (teichuronic acid, which is associated with protein), 1097 (PO2-symmetrical ν) and 1550 cm−1 (tryptophan).
The SERS features observed that are related to carbohydrates include those at 326, 343, 393, 420, 467, 1130 and 1327 cm−1. The SERS bands at 420 (C–C–C deformation), 467 (C–OH ρt, β(CCC), β(CCO), β(OCO), corresponding to a specific sugar profile), 1130 ((δ (CC, CO, –COH) (specific SERS feature of glucose) and 1327 cm−1 (belonging to δ(CH)) appeared in biofilm-forming bacterial strains, and were absent in non-biofilm-forming strains.65,66 The SERS biomarkers at 326 (β (C–C–O) ring vibration), 343 (C–C–O δ, C–OH ρt), and 393 cm−1 (S–S ν) are observed in both strong and non-biofilm forming strains, but the increasing intensities of these SERS biomarkers indicate the high contents of polysaccharides in the cell mass/pellets of biofilm-forming strains in comparison to those that do not form biofilms. The peak observed at 393 cm−1 related to S–S ν has lower intensity in the biofilm-forming A. baumannii strains than that in those A. baumannii strains which do not form biofilms. In the process of biofilm formation and accumulation, the concentration of polysaccharide intercellular adhesion (PIA) is high, which mediates the attachment of A. baumannii bacteria to the fibrous net that builds the biofilm mass.67–69
The SERS bands related to proteins are 516 (S–S ν, disulfide), 747 ((C–S) ν), 762 (ring breathing of tryptophan), 861 (C–C ν, C–O–C 1,4-glycosidic linkage; cell wall of Gram-positive bacteria associated with protein), 953 (υ(CO) and C–N of proteins), 990 (C–C ring breathing in phenylalanine), 1070 (C–N ν of proteins), 1395 (COO−1 symmetric ν), 1550 (tryptophan) and 1693 cm−1 (alpha-helix of Amide I), of which the peaks at 762 (ring breathing of tryptophan), 990 (C–C ring breathing in phenylalanine), 1395 (COO−1 symmetric ν) and 1693 cm−1 (alpha-helix of Amide I) are present only in bacterial strains involved in biofilm formation. Moreover, the SERS band at 953 cm−1 (υ(CO) is present only in the non-biofilm-forming A. baumannii strains and absent in the biofilm-forming A. baumannii strains (medium, strong) due to alterations in the protein structure.70 The SERS features, including 516 (S–S ν, disulfide), 747 ((C–S) ν), 861 (C–C ν, C–O–C 1,4-glycosidic linkage, CCH deformation in tyrosine ring breathing; cell wall of Gram-positive bacteria) and 1550 cm−1 (tryptophan), are observed in all bacterial cell mass samples with intensity-based differences. There are thirteen cell wall anchor proteins (CWA) related to virulence in A. baumannii, which may be helpful in biofilm development and are covalently bound to peptidoglycan in the cell wall. Moreover, multimeric cellular proteins, structural proteins and surface proteins (such as pili, flagella, fimbriae and appendages proteins) are also present in the biofilm forming process. The function of these proteins in the biofilm-forming process is to connect the cells with one another through various microorganisms.71,72
SERS features observed at 633 (C–C twisting of the guanine ring region), 679 (thymine), 699 (guanine (COO deformation), 720 (adenine), 729 (the glycosidic ring conformation of adenine) and 1097 cm−1 (PO2− symmetrical ν) are linked to DNA. The SERS band at 633 (the guanine ring region and the tyrosine amino acid chain include aromatic rings with a C–C configuration), 699 (guanine (COO deformation), 720 (adenine) and 729 cm−1 (the glycosidic ring conformation of adenine) are found only in bacterial strains capable of producing biofilms, and not in ones that are unable to do so. This indicates the high contents of DNA in biofilm-forming strains due to the double strand chromosomal plasmid DNA utilized in New Delhi metallo beta-lactamase (NDM) synthesis, which transforms and replicates its genes in new colonies that serve as hosts and play a crucial role in bacterial cell adhesion to surfaces and biofilm formation. This plasmid genetic transformation of information changes the non-biofilm-forming bacterial strains into biofilm-forming bacterial strains.73–75 Moreover, the SERS band at 679 (thymine) and 1097 cm−1 (PO2− vs) are observed in all samples with intensity-based differences.
![]() | ||
Fig. 2 PCA scatter plot of the score for the SERS spectral groups of non, medium and strong biofilm-forming A. baumannii strains. |
Fig. 3a illustrates a pairwise comparison of the PCA scatter score of the SERS spectra of non- and medium biofilm-forming strains of A. baumannii bacteria, in which the SERS spectral datasets are differentiated along PC-2 with variability of 23.59%, while 60.68% variance was obtained, as explained by PC-1. The spectra of the medium biofilm A. baumannii strains concentrated along the positive axis of PC-2 are shown with blue dots, while the red dots indicate the spectra of the non-biofilm forming A. baumannii strain, which are grouped along the negative axis of PC-2. This indicates that the changes in the biomolecular contents of the two types of A. baumannii bacterial strains through their SERS spectral features. The basis for distinction in the SERS spectral datasets of the bacterial cell mass (pellets) is explained by the pairwise PCA loadings in Fig. 3b, where the negative loadings represent the SERS spectral bands of bacterial strains which do not form a biofilm as their SERS spectra are clustered on the negative side as red dots in the PCA scatter plot, including characteristic SERS bands at 326, 343, 467, 729, 990, 1327, 1395 and 1693 cm−1. The positive loadings represent the spectra of samples of medium biofilm strains, which appear on the positive side as blue dots in the PCA scatter plot and are observed at 393, 633, 679, 762, 861 and 1097 cm−1.
![]() | ||
Fig. 3 Pairwise PCA: (a) scatter plot of scores and (b) loadings for SERS spectral groups of non-biofilm- versus medium biofilm-forming strains. |
Fig. 4a shows a PCA scatter score of the spectral data of cell mass samples of A. baumannii strains producing medium and strong biofilms, which clearly separates the spectral datasets of strains producing medium and strong biofilms, providing an explanation for 79.70% of the variation along PC-1. The spectra of the bacteria with a medium level of biofilm formation are clustered on the negative side of PC-1, which are represented by blue dots, whereas bacteria with a strong biofilm formation tendency are shown by green dots concentrated along the positive side of PC-1. The two clusters of spectral data are distant from each other, showing significant differences in the SERS spectral characteristics of these groups of samples at the biochemical level.
![]() | ||
Fig. 4 Pairwise PCA: (a) scatter plot of scores and (b) loadings of the SERS spectra of strong biofilm- versus medium biofilm-forming A. baumannii strains. |
Fig. 4b presents data regarding distinct biochemical characteristics in the manner of pairwise comparison of the SERS spectra of strong versus medium biofilm-forming A. baumannii strains as PCA loadings. The loadings on the positive side of PC-1 are observed at 393 (S–S ν of carbohydrates), 679 (thymine related to DNA), 720 (adenine), 953 (υ(CO) and C–N of proteins) and 1395 cm−1 (the protein has COO−1 symmetry in its structure), as the SERS spectral features of the strong biofilm (which are green dots) are plotted on the positive side of the PC-1 scatter plot. Meanwhile, the negative loading indicates the SERS bands of medium biofilm-forming strains that include 467 (C–C–O δ, C–OH ρt, β(CCC), β(CCO), β(OCO) in carbohydrates), 861 (C–C ν, C–O–C 1,4-glycosidic linkage, CCH deformation in ring breathing of tyrosine, namely in the CCH group; Gram-positive bacterium cell wall (teichuronic acid, which is associated with protein), 1097 (PO2− vs in DNA), and 1550 cm−1 (tryptophan linked to protein), which are represented with blue dots scattered on the negative side of the PC-1.
The pairwise PCA of the SERS spectra of the non-biofilm forming A. baumannii strains (red dots), on the negative axis, and the strong biofilm-forming A. baumannii strains (green dots), along the positive axis of PC-1, is shown in Fig. 5a. PC-1 exhibits a variance of 79.70%, while 16.53% is shown along PC-2. The SERS spectra of both types of bacterial strains are differentiated from each other on the basis of biomolecular changes, which can be associated with their capacity of producing the biofilm.
![]() | ||
Fig. 5 Pairwise PCA: (a) scatter plot of scores and (b) loadings of the SERS spectra of non-biofilm- versus strong biofilm-forming A. baumannii strains. |
Fig. 5b displays the PC-1 loadings for the PCA of the SERS spectra from the cell mass samples of the strong non-biofilm forming A. baumannii strains. The positive loadings correspond to spectra of strong biofilm-forming (green dots) strains grouped on the positive axis, and the negative loadings are associated with the spectra of the non-biofilm forming (red dots) strains clustered on the negative axis. The major SERS bands of the strong biofilm strains in the positive loadings include 393 (S–S ν in carbohydrates), 679 (related to guanosine, thymine in DNA), 720 (adenine), 953 (υ(CO) and C–N of proteins) and 1395 cm−1 (COO−1vs group that is commonly found in proteins), while the non-biofilm SERS features in negative loading are observed at 420, 467, 699, 861, 1097 and 1550 cm−1.
Fig. 6 displays the scatter score of the PLS-DA model showing the SERS spectra of the strong (green dots), medium (blue dots), and non-biofilm forming bacterial strains (red dots) on the positive axis, negative axis, and the negative side of the PLS-DA scatter score plot, respectively. This clearly differentiates between the SERS spectra of the positive biofilm strains (blue clustered dots, green clustered dots) and the negative biofilm strains (red clustered dots) due to their characteristic SERS spectral features, which can be associated with the increasing concentration levels of DNA, protein, lipids and carbohydrates in the high biofilm-forming strains.
![]() | ||
Fig. 6 PLS-DA scatter plots: (a) Monte Carlo and (b) K-fold cross-validation of the SERS spectra of non, medium and strong biofilm-forming strains of A. baumannii. |
The SERS features observed in the PLS-DA loading plot of the non-, medium and strong biofilm-forming strains of A. baumannii are shown in the ESI (Fig. S2),† corresponding to the SERS bands observed in the mean plot of Fig. 1. The SERS spectral features observed in the non-, medium and strong biofilm-forming strains include bands at 326, 343, 420, 516, 633, 679, 720, 861, 953, 990, 1097, 1130, 1327, 1445 and 1693 cm−1.
Fig. 7 displays the ROC curve of the PLS-DA model used to classify the SERS spectral data sets of the positive and negative biofilm strains. The ROC curve was drawn between 1-specificity (false positive rate) and sensitivity (true positive rate). The Monte Carlo validation method was employed, and the Area Under the Curve (AUC) value for the PLS-DA model is found to be 0.81 with a sensitivity of 100% and specificity of 94.3%. Moreover, LKOCV was applied to further verify the robustness of the classification model. The values of AUC (0.90), sensitivity (100%), and specificity (98%), respectively, are determined through LKOCV for the PLS-DA model. A value of AUC at close to 1 indicates that the model is a good fit with the high accuracy and validation, while a value at close to zero would suggest an inaccurate fit and bad performance of this model.
The biofilm-forming and non-biofilm-forming bacterial strains have common intensity-based different SERS spectral features that are labeled with dotted line, and some features that are present with the same intensity in the biofilm-forming and non-biofilm-forming strains are shown in Fig. 1. These features may complicate the chemometric model's ability to achieve 100% specificity. Although SERS provides high sensitivity, some biochemical features may be common across different strains, making it difficult for the chemometric models to differentiate between them with 100% accuracy. The chemometric model PLS-DA relies on spectral features to classify strains. This model is effective, but sometimes struggles with borderline cases where the spectra of the biofilm and non-biofilm strains have common features. This results in misclassification, which reduces the specificity of the model. Improvements can be achieved through optimized preprocessing, by minimizing the variation in the experimental conditions, and by addressing the class imbalance.76,77 These steps should help to enhance the model's performance and bring the specificity closer to 100%.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4ra06267a |
‡ The first two authors contributed equally to this work. |
This journal is © The Royal Society of Chemistry 2025 |