Cooperative nature of residue interactions in biomolecular condensates

Cooperativity of residue interactions is prevalent in proteins and plays an essential role in protein folding. However, whether cooperativity exists in intrinsically disordered proteins (IDPs) and their assembled biomolecular condensates remains elusive. Here, we reveal the presence of residue interaction cooperativity in various IDP condensates using atomistic molecular dynamics simulations. In all these condensates, cooperative residue interactions manifest as packing between finite residue motifs (3 to 5 residues) with lifetimes three orders of magnitude longer than residue contacts. Notably, these higher-order interactions (HOIs) exhibit a significant nonadditive free energy, which originates from a local low-dielectric environment. The underlying sequence patterns like “Y-X-Q”, “D-X-R”, and “S-Y-X” are identified. Furthermore, HOIs dominate the cohesive energy of IDPs and condensates can be treated as HOI-crosslinked networks. This model can successfully reproduce experimental results of condensate physical properties, including saturation concentration, viscoelasticity, and interfacial tension.