Abstract
The most commonly occurring intrinsically disordered proteins (IDPs) are polyampholytes, which are defined by the duality of low net charge per residue and high fractions of charged residues. Recent experiments have uncovered nuances regarding sequence-ensemble relationships of model polyampholytic IDPs. These include differences in conformational preferences for sequences with lysine vs. arginine, and the suggestion that well-mixed sequences form a range of conformations, including globules, conformations with ensemble averages that are reminiscent of ideal chains, or self-avoiding walks. Here, we explain these observations by analyzing results from atomistic simulations. We find that polyampholytic IDPs generally sample two distinct stable states, namely globules and self-avoiding walks. Globules are favored by electrostatic attractions between oppositely charged residues, whereas self-avoiding walks are favored by favorable free energies of hydration of charged residues. We find sequence-specific temperatures of bistability at which globules and self-avoiding walks can coexist. At these temperatures, ensemble averages over coexisting states give rise to statistics that resemble ideal chains without there being an actual counterbalancing of intra-chain and chain-solvent interactions. At equivalent temperatures, arginine-rich sequences tilt the preference toward globular conformations whereas lysine-rich sequences tilt the preference toward self-avoiding walks. We also identify differences between aspartate and glutamate containing sequences, whereby the shorter aspartate sidechain engenders preferences for metastable, necklace-like conformations. Finally, although segregation of oppositely charged residues within the linear sequence maintains the overall two-state behavior, compact states are highly favored by such systems. Significance Statement Intrinsically disordered regions (IDRs) of proteins, when tethered to folded domains, function either as flexible tails or as linkers between domains. Most IDRs are polyampholytes that comprise a mixture of oppositely charged residues. Recent measurements of tethered polyampholytes showed that tendency of arginine- and lysine-rich sequences to behave very differently from one another. Using computer simulations, we show that these differences are determined by differences in free energies of hydration, steric volumes, and other considerations. Further, the interplay between electrostatic attractions and favorable free energies of hydration creates distinct stable states for polyampholytic IDRs. These findings have implications for switch-like transitions and the regulation of effective concentrations of interaction motifs by IDRs.