This study aims to present the deterministic response spectrum for various important district of Maharashtra state by employing deterministic seismic hazard analysis. Seismogenic sources and earthquake data are gathered from multiple resources within a radius of 400 km from the state boundary. The earthquake data are processed for homogenisation, declustering and completeness analysis. The homogenised earthquake data are used for the development of seismic source zone of the study region. The declustering of the homogenised earthquake catalogue is performed four times using four different declustering methods. The common earthquake events from the four declustered catalogues are used for the completeness analysis. The complete earthquake data is superimposed on the fault map, and these earthquake events are utilised for estimating the maximum magnitude potential of all seismogenic sources, which relies on three distinct scenarios. The three-scenario earthquake is used independently as an input of the attenuation model. In the current study, three region-specific attenuation models are used to estimate the hazard at a specific site. To address the epistemic uncertainty associated with seismogenic sources and earthquake data, a total of twelve branches are prepared in the logic tree. Finally, design ground acceleration values of ten significant district in Maharashtra are finalised at bedrock level, C-type, and D- type site and consecutively, a smooth response spectrum is developed for each district. The findings indicate that, except the response spectrum of Region 1 for the Jalgaon, Chandrapur, Latur, Pune, and Thane districts, the response spectra of all regions of all the districts are more conservative than the response spectrum developed using the zone factor mentioned in IS1893-Part 1 (2016).