Intracellular bodies such as nucleoli, Cajal bodies, and various signaling assemblies, represent membraneless organelles, or condensates, that form via liquidliquid phase separation (LLPS) 1,2 . Biomolecular interactions, particularly homotypic interactions mediated by self-associating intrinsically disordered protein regions (IDRs), are thought to underlie the thermodynamic driving forces for LLPS, forming condensates that can facilitate the assembly and processing of biochemically active complexes, such as ribosomal subunits within the nucleolus. Simplified model systems 3-6 have led to the concept that a single fixed saturation concentration (C sat ) is a defining feature of endogenous LLPS 7-9 , and has been suggested as a mechanism for intracellular concentration buffering 2,7,8,10 . However, the assumption of a fixed C sat remains largely untested within living cells, where the richly multicomponent nature of condensates could complicate this simple picture. Here we show that heterotypic multicomponent interactions dominate endogenous LLPS, and give rise to nucleoli and other condensates that do not exhibit a fixed C sat . As the concentration of individual components is varied, their partition coefficients change, in a manner that can be used to extract thermodynamic interaction energies, that we interpret within a framework we term polyphasic interaction thermodynamic analysis (PITA). We find that heterotypic interactions between protein and RNA components stabilize a variety of archetypal intracellular condensates, including the nucleolus, Cajal bodies, stress granules, and P bodies. These findings imply that the composition of condensates is finely tuned by the thermodynamics of the underlying biomolecular interaction network. In the context of RNA processing condensates such as the nucleolus, this stoichiometric self-tuning manifests in selective exclusion of fully-assembled RNP complexes, providing a thermodynamic basis for vectorial ribosomal RNA (rRNA) flux out of the nucleolus. The PITA methodology is conceptually straightforward and readily implemented, and it can be broadly utilized to extract thermodynamic parameters from microscopy images. These approaches pave the way for a deep understanding of the thermodynamics of multicomponent intracellular phase behavior and its interplay with nonequilibrium activity characteristic of endogenous condensates.To determine the thermodynamics of LLPS for intracellular condensates we first focused on the liquid granular component (GC) of nucleoli within HeLa cells, in particular on the protein Nucleophosmin (NPM1), which is known to be a key driver of nucleolar phase separation 11,12 . Under typical endogenous expression levels, we estimate NPM1 concentration in the nucleoplasm to be roughly C dil~4 μM; from simple binary phase separation models (i.e. Regular solution theory) 13 , this apparent saturation concentration, C sat , is expected to be fixed even under varied NPM1 expression levels (Fig. 1C, Supplemental Text). Consistent with previous studies 11 , ...