A search for neutral Higgs bosons of the Minimal Supersymmetric Standard Model (MSSM) is performed using a data sample corresponding to an integrated luminosity of up to 13.3 fb from proton--proton collisions at TeV recorded by the ATLAS detector at the LHC. The heavy resonance is assumed to decay to a pair with at least one lepton decaying hadronically. The search is performed in the mass range of 0.2-1.2 TeV. The data are in good agreement with the background predicted by the Standard Model, hence results are given as upper limits on the production cross section times branching fraction as a function of the Higgs boson mass. The results are interpreted in MSSM benchmark scenarios. The most stringent MSSM parameter space constraints from the search exclude at 95% confidence level for ~\GeV\ and for ~\GeV\ in the MSSM scenario.