We develop a "hybrid" model of the proton inelastic structure functions applicable for a wide region of invariant masses of produced states 𝑊 and invariant momentum transfer 𝑄 including deep inelastic scattering (DIS), nucleon resonance production as well as the region close to inelastic threshold. In the DIS region, we compute the structure functions in terms of the parton distributions together with higher-twist corrections which are determined from a global QCD fit. The resonant part is addressed in terms of the Breit-Wiegner contributions from five states: Δ(1232) resonance, 𝑁 (1440) Roper resonance, and three effective resonances describing the second and third resonance regions. The couplings of the nucleon resonances to photon are described in terms of helicity amplitudes. The nonresonant background is treated in terms of a smooth extrapolation of the DIS structure functions to low-𝑊 and low-𝑄 values with a proper inelastic threshold behavior and the real photon limit 𝑄 2 = 0 constrained by the photoproduction data. We independently treat the transverse, 𝐹 𝑇 , and the longitudinal, 𝐹 𝐿 , structure function and fix the model parameters from a global analysis of the electron-proton differential cross section data together with data on the total photoproduction cross section. We demonstrate a very good performance of the model by comparing our predictions with data on cross sections and the structure functions 𝐹 2 and 𝑅 = 𝐹 𝐿 /𝐹 𝑇 in the resonance and the DIS region.