We study the linear and nonlinear optical properties of a well-known acid base indicator, bromophenol blue (BPB), in aqueous solution by employing static and integrated approaches. In the static approach, optical properties have been calculated using time-dependent density functional theory (TD-DFT) on the fully relaxed geometries of the neutral and different unprotonated forms of BPB. Moreover, both closed and open forms of BPB were considered. In the integrated approach, the optical properties have been computed over many snapshots extracted from molecular dynamics simulation using a hybrid time-dependent density functional theory/molecular mechanics approach. The static approach suggests closed neutral double right arrow anionic interconversion as the dominant mechanism for the red shift in the absorption spectra of BPB due to a change from acidic to basic pH. It is found by employing an integrated approach that the two interconversions, namely open neutral double right arrow anionic and open neutral double right arrow dianionic, can contribute to the pH- dependent shift in the absorption spectra of BPB. Even though both static and integrated approaches reproduce the pH-dependent red shift in the absorption spectra of BPB, the latter one is suitable to determine both the spectra and spectral broadening. Finally, the computed static first hyperpolarizability for various protonated and deprotonated forms of BPB reveals that this molecule can be used as a nonlinear optical probe for pH sensing in addition to its highly exploited use as an optical probe.