On the basis of molecular dynamics simulation, we model the ester carbonyl stretch FTIR signals of methyl acetate in D2O and DMSO. An ab initio map is constructed at the B3LYP/6-311++G** level to relate the carbonyl stretch frequency to the external electric field. Using this map, fluctuating Hamiltonian of the carbonyl stretch is constructed from the MD simulation trajectory. The IR spectra calculated based on this Hamiltonian are found to be in good agreement with the experiment. For methyl acetate in D2O, hydrogen bonding on alkoxy oxygen causes a blue shift of frequency, while that on carbonyl oxygen causes a red shift. Two peaks observed in FTIR signals originate from the balance of these two effects. Furthermore, in both D2O and DMSO solutions, correlations are found between the instantaneous electric field on C=O and the frequencies. Broader line width of the signal in D2O suggests a more inhomogeneous electric field distribution due to the complicated hydrogen-bonding environment.