Vibrational circular dichroism (VCD) spectroscopy is utilized to reveal the detailed conformational distributions of the dominant serine species in aqueous solutions under three representative pH conditions of 1.0, 5.7, and 13.0, together with vibrational absorption (VA) spectroscopy, density functional theory (DFT), and molecular dynamics simulation. The experimental VA and VCD spectra of serine in H(2)O and D(2)O in the fingerprint region under three pH values are obtained. DFT calculations at the B3LYP/6-311++G(d,p) level are carried out for the protonated, zwitterionic, and deprotonated serine species. The lowest-energy conformers of all three species are identified and their corresponding VA and VCD spectra simulated. A comparison between the gas-phase simulations and the experimental VA and VCD spectra suggests that one or two of the most stable conformers of each species contribute predominantly to the observed data, although some discrepancies are noted. To account for the solvent effects, both the polarizable continuum model and the explicit solvation model are considered. Hydrogen-bonded protonated, zwitterionic, and deprotonated serine-(water)(6) clusters are constructed based on radial distribution function analyses and molecular dynamics snapshots. Geometry optimization and VA and VCD simulations are performed for these clusters at the B3LYP/6-311++G(d,p) level. Inclusion of the explicit water molecules is found to improve the agreement between theory and experiment noticeably in all three cases, thus enabling conclusive conformational distribution analyses of serine in aqueous solutions directly.