We describe a reconstruction method that allows direct recovery of tissue chromophore concentrations and acoustic velocity using tomographic multi-spectral photoacoustic measurements. The method is tested in detail using extensive simulated data. Criteria for choosing optimum wavelength sets are analyzed for simultaneous reconstruction of tissue functional and acoustic properties. © 2010 IOS Press and the authors. All rights reserved.