A model gel of whey protein isolate (WPI) was prepared by cold gelation with calcium. This system was modified by the addition of free cysteine residues (Cys) at different steps of the process. The WPI cold-set gels obtained were then subjected to heat treatment at 90 degrees C. First, the effect of Cys addition on the heat-induced aggregation of WPI was studied through Atomic Force Microscopy (AFM) and infrared spectroscopy (ATR-FTIR), while Cys' effect on cold gelation was observed by AFM, Confocal Laser Scanning Microscopy (CLSM) and oscillatory rheology (amplitude sweeps). The impact of heating on the microstructure and the viscoelastic properties of the WPI cold-set gels were finally investigated through several techniques, including DSC, ATR-FTIR, CLSM, cryo-SEM, and rheological measurements (temperature sweeps). When added during the first step of cold gelation, Cys modified heat-induced aggregation of WPI, resulting in the formation of a denser gel network with a fractal dimension (Df) of 2.8. However, the addition of Cys during the second step of cold gelation led to the formation of highly branched clusters of WPI and a looser gel network was observed (Df = 2.4). In this regard, the use and limitations of oscillatory rheology and the "Kraus model" to determine the Df of WPI cold-set gels was discussed. The viscoelastic properties and the microstructure of the WPI cold-set gels were irreversibly modified by heating. Gels were stiffer, more brittle, and coarser after heat treatment. New disulfide bonds and calcium bridges formed, as well as H-bonded beta-sheets, all contributing to the formation of the final gel network structure.