  Biologically functional liquid-liquid phase separation of intrinsically
disordered proteins (IDPs) is driven by interactions encoded by their amino
acid sequences. Little is currently known about the molecular recognition
mechanisms for distributing different IDP sequences into various cellular
membraneless compartments. Pertinent physics was addressed recently by applying
random-phase-approximation (RPA) polymer theory to electrostatics, which is a
major energetic component governing IDP phase properties. RPA accounts for
charge patterns and thus has advantages over Flory-Huggins and Overbeek-Voorn
mean-field theories. To make progress toward deciphering the phase behaviors of
multiple IDP sequences, the RPA formulation for one IDP species plus solvent is
hereby extended to treat polyampholyte solutions containing two IDP species.
The new formulation generally allows for binary coexistence of two phases, each
containing a different set of volume fractions $(\phi_1,\phi_2)$ for the two
different IDP sequences. The asymmetry between the two predicted coexisting
phases with regard to their $\phi_1/\phi_2$ ratios for the two sequences
increases with increasing mismatch between their charge patterns. This finding
points to a multivalent, stochastic, "fuzzy" mode of molecular recognition that
helps populate various IDP sequences differentially into separate phase
compartments. An intuitive illustration of this trend is provided by
Flory-Huggins models, whereby a hypothetical case of ternary coexistence is
also explored. Augmentations of the present RPA theory with a relative
permittivity $\epsilon_{\rm r}(\phi)$ that depends on IDP volume fraction
$\phi=\phi_1+\phi_2$ lead to higher propensities to phase separate, in line
with the case with one IDP species we studied previously. ...
