A new multistate empirical valence bond model (MS-EVB3) is developed for proton solvation and transport in aqueous solutions. The new model and its quantum version (qMS-EVB3) are based on the MS-EVB2 model [Day et al., J. Chem. Phys. 2002, 117, 5839] and recently developed flexible water models-the SPC/ Fw model [Wu et al. J. Chem. Phys. 2006, 124, 24503] and the qSPC/Fw model [Paesani et al. J. Chem. Phys. 2006, 125, 184507]-for classical and quantum simulations, respectively. Using ab initio data as benchmarks, the binding energies and optimized geometries calculated with the new model for protonated water clusters, as well as the potential energy surface for proton shuttling between water molecules in a cluster environment, are improved in comparison to the MS-EVB2 model. For aqueous solutions, classical and quantum molecular dynamics simulations with the MS-EVB3 model yield a more accurate description of the solvation structure and diffusive dynamics of the excess proton. New insight is also provided into the proton solvation and hopping dynamics in water, as well as the "amphiphilic" nature of the hydrated proton that has been predicted to give rise to its enhanced concentration at aqueous interfaces and an effectively lower pH of the air-water interface.