As a versatile nanomaterial derived from renewable sources, nanocellulose has attracted considerable attention for its potential applications in various sectors, especially those focused on water treatment and remediation. Here, we have combined atomic force microscopy (AFM) and reactive molecular dynamics (RMD) simulations to characterize the interactions between cellulose nanofibers modified with carboxylate or phosphate groups and the protein foulant model bovine serum albumin (BSA) at pH 3.92, which is close to the isoelectric point of BSA. Colloidal probes were prepared by modification of the AFM probes with the nanofibers, and the nanofiber coating on the AFM tip was for the first time confirmed through fluorescence labeling and confocal optical sectioning. We have found that the wet-state normalized adhesion force is approximately 17.87 +/- 8.58 pN/nm for the carboxylated cellulose nanofibers (TOCNF) and about 11.70 +/- 2.97 pN/nm for the phosphorylated ones (PCNF) at the studied pH. Moreover, the adsorbed protein partially unfolded at the cellulose interface due to the secondary structure's loss of intramolecular hydrogen bonds. We demonstrate that nanocellulose colloidal probes can be used as a sensitive tool to reveal interactions with BSA at nano and molecular scales and under in situ conditions. RMD simulations helped to gain a molecular- and atomistic-level understanding of the differences between these findings. In the case of PCNF, partially solvated metal ions, preferentially bound to the phosphates, reduced the direct protein-cellulose connections. This understanding can lead to significant advancements in the development of cellulose-based antifouling surfaces and provide crucial insights for expanding the pH range of use and suggesting appropriate recalibrations.