The initial recovery steps for recombinant protein were carried out using aqueous two-phase systems in the presence of guanidine hydrochloride. To investigate the phase behavior of such systems, equilibrium data was obtained for molecular weight fractions of PEG (3000 and 6000) + sodium sulfate (pH 7) + guanidine hydrochloride (5% and 10% mass) + water at 25 degrees C. Displacement of the binodal toward higher concentrations through an increase in guanidine hydrochloride concentration was clearly observed. A slightly larger two-phase was obtained as the PEG molecular weight increased from PEG 3000 to 6000. The experimental binodal data of the system was successfully correlated with Bleasdale equations (high R-2 and low SE). The effects of molecular weight, guanidine hydrochloride concentration, tie-line length, and slope of tie-line on the partition behaviors of guanidine hydrochloride were examined at constant pH. The compositions of each aqueous phase were calculated using two thermodynamic models for the activity coefficient. The calculated results showed that the extended Pitzer and extended UNIQUAC models can predict the data in the quaternary system well. The fitted binary interaction parameters of the models were reported. (C) 2016 Elsevier B.V. All rights reserved.