Protein−polysaccharide complexes constructed via self-assembly methods are often used to develop novel biomaterials for a wide range of applications in biomedicine, food, and biotechnology. The objective of this work was to investigate theoretically and to demonstrate via constant-pH Monte Carlo simulations that the complexation phenomenon between insulin (INS) and the cationic polyelectrolyte chitosan (CS) is mainly driven by an electrostatic mechanism.Experimental results obtained from FTIR spectra and ζ-potential determinations allowed us to complement the conclusions.The characteristic absorption bands for the complexes could be assigned to a combination of signals from CS amide I and INS amide II. The second peak corresponds to the interaction between the polymer and the protein at the level of amide II. INS− CS complexation processes not expected when INS is in its monomeric form, but for both tetrameric and hexameric forms, incipient complexation due to charge regulation mechanism took place at pH 5. The complexation range was observed to be 5.5 < pH < 6.5. In general, when the number of INS units increases in the simulation process, the solution pH at which the complexation can occur shifts toward acidic conditions. CS's chain interacts more efficiently, i.e. in a wider pH range, with INS aggregates formed by the highest monomer number. The charge regulation mechanism can be considered as a previous phase toward complexation (incipient complexation) caused by weak interactions of a Coulombic nature. 36 biotechnology. 2,3 Self-assembly is a ubiquitous and natural 37 process that lead to complex nanostructures with outstanding 38 functionalities. 4 In particular, the implementation of pH-