Protein oligomerization regulates many critical physiological processes, and its dysregulation can contribute to dysfunction and diseases. Elucidating the assembly pathways and quantifying their underlying thermodynamic and kinetic parameters are crucial for a comprehensive understanding of biological processes and for advancing therapeutics targeting abnormal protein oligomerization. Established binding assays, with limited mass precision, often rely on simplified models for data interpretation. In contrast, high-resolution native mass spectrometry (nMS) can directly determine the stoichiometry of biomolecular complexes in vitro. However, quantification is hindered by the fact that the relative abundances of gas-phase ions generally do not reflect solution concentrations due to nonuniform response factors. Recently, slow mixing mode (SLOMO)-nMS, which can quantify the relative response factors of interacting species, has been demonstrated to reliably measure the affinity (Kd) of binary biomolecular complexes. Here, we introduce an extended form of SLOMO-nMS that enables simultaneous quantification of the thermodynamics in multistep association reactions. Application of this method to homo-oligomerization of concanavalin A and insulin confirmed the reliability of the assay and uncovered details about the assembly processes that had previously resisted elucidation. Results acquired using SLOMO-nMS implemented with charge detection shed new light on the binding of recombinant human angiotensin-converting enzyme 2 and the SARS-CoV-2 spike protein. Importantly, new assembly pathways were uncovered, and the affinities of these interactions, which regulate host cell infection, were quantified. Together, these findings highlight the tremendous potential of SLOMO-nMS to accelerate the characterization of protein assembly pathways and thermodynamics and, in so doing, enhance fundamental biological understanding and facilitate therapeutic development. https://orcid.org/0000-0002-3389-7112.