We develop an accurate simulation-based inference framework for high-mass (≳107M⊙) black-hole binaries observable by LISA. The method is implemented within the DINGO gravitational-wave parameter-estimation code, extending its application from ground-based detectors to the LISA band. We train a normalizing-flow model using aligned-spin higher-mode waveform models and a low-frequency approximation of the detector response. After sampling, we importance-sample to the true posterior. We validate performance on simulated signals spanning the signal-to-noise regimes relevant for LISA observations and benchmark our new DINGO implementation against standard methods. We report robust agreement in the inferred posterior distributions up to signal-to-noise ratios of ∼500. At higher signal-to-noise ratios of ∼1000, we observe a reduction in sampling efficiency, while still yielding unbiased and tightly localized posteriors that can be used as a starting point for follow-up with traditional this http URL trained flow can generate 20 thousand posterior samples in less than a minute, establishing DINGO as a promising neural inference framework for rapid full-parameter estimation of massive black-hole binaries in the LISA band. The likelihood-free nature of this approach allows for straightforward generalizations, including a time-dependent detector response, non-stationary noise artifacts such as gaps and glitches, and low-latency parameter estimations.