Coarse-grained simulations highlight residues and sequence features that are essential drivers of biomolecular phase separation . Residue-level CG simulations have been successful in describing protein phase behavior , capturing trends in experimental phase separation propensities . Since residue-level models provide efficient access to long time-and large lengthscales, it becomes possible to connect sequence to protein and material properties as well as to condensate morphologies . It has been shown with the CALVADOS and Mpipi residue-level CG models that realistic simulations of disordered protein sequences at the IDRome scale are possible , providing insights into the rules that govern the global conformational ensembles of IDPs. Generally, an accurate description of phase separation requires multi-chain simulations. In exciting recent works, von Bülow et al. and Ginell et al. have extended the IDRome-wide simulations to investigate trends in phase separation and co-condensation of proteins . Pesce et al. showed how sequence patterning determines phase behavior at a fixed protein sequence composition . The trade-off between between the tendency of a protein sequence to interact and the viscosity of the dense protein phase was demonstrated by An et al., who then used active learning to find sequences that form dynamic condensates . Simulations can also reveal how different protein sequences mix or demix and how differences in interactions can give rise to a diversity of condensate morphologies . Depending on whether the interactions are compatible or not, mixtures of proteins form either distinct condensates , or partially or fully engulfed condensates in which one phase surrounds another. For example, Driver and Onck demonstrated how orthogonal hydrophobic interactions and cation-controlled interactions determine the morphology of condensates of IDRs of FET family transcription factors . Phopshorylation of the disordered C-terminal domain (CTD) of RNA polymerase II can lead to demixing between condensate dominated by hydrophobic interactions, which is surrounded by a condensate of negatively-charged phosphorylated CTD and a positively charged IDR from an elongation factor . Computing interfacial tensions from simulations enables one to rationalize the formation of different morphologies and to decipher how molecular scale interactions shape these morphologies .