Hydrophobic interaction chromatography (HIC) is used as a critical polishing step in the downstream processing of biopharmaceuticals. Normally the process development of HIC is a cumbersome and time-consuming task, and the mechanical models can provide a powerful tool to characterize the process, assist process design and accelerate process development. However, the current estimation of model parameters relies on the inverse method, which lacks an efficient and logical parameter estimation strategy. In this study, a parameter-by-parameter (PbP) method based on the theoretical derivation and simplifying assumptions was proposed to estimate the Mollerup isotherm parameters for HIC. The method involves three key steps: (1) linear regression (LR) to estimate the salt-protein interaction parameter and the equilibrium constant; (2) linear approximation (LA) to estimate the stoichiometric parameter and the maximum binding capacity; and (3) inverse method to estimate the protein-protein interaction parameter and the kinetic coefficient. The results indicated that the LR step should be used for dilution condition (loading factor below 5%), while the LA step should be conducted when the isotherm is in the transition or nonlinear regions. Six numerical experiments were conducted to implement the PbP method. The results demonstrated that the PbP method developed allows for the systematic estimation of HIC parameters one-by-one, effectively reducing the number of parameters required for inverse method estimation from six to two. This helps prevent non-identifiability of structural parameters. The feasibility of the PbP-HIC method was further validated by real-world experiments. Moreover, the PbP method enhances the mechanistic understanding of adsorption behavior of HIC and shows a promising application to other stoichiometric displacement model-derived isotherms.