As the water molecules found at the interface of lipid bilayers exhibit distinct structural and reorientation dynamics compared to water molecules found in bulk, the fluctuations in their hydrogen bond (HB) network are expected to be different from those generated by the bulk water molecules. The research presented here aims to gain an insight into temperature-dependent fluctuations of a HB network of water molecules found in an interfacial layer of multilamellar liposomes (MLVs) composed of anionic 1,2-dimyristoyl-sn-glycero-3-phospho-L-serine (DMPS) lipids. Besides suspending DMPS lipids in phosphate buffer saline (PBS) of different pH values (6.0, 7.4, and 8.0), the changes in HB network fluctuations were altered by the incorporation of a non-polar flavonoid molecule myricetin (MCE) within the hydrocarbon chain region. By performing a multivariate analysis on the water combination band observed in temperature-dependent FTIR spectra, the results of which were further mathematically analyzed, the temperature-dependent fluctuations of interfacial water molecules were captured; the latter were the greatest for DMPS in PBS with a pH value of 7.4 and in general were greater for DMPS multibilayers in the absence of MCE. The presence of MCE made DMPS lipids more separated, allowing deeper penetration of water molecules towards the non-polar region and their restricted motion that resulted in decreased fluctuations. The experimentally observed results were supported by MD simulations of DMPS (+MCE) lipid bilayers.