Sequence charge decoration dictates coil-globule transition in intrinsically disordered proteins

Coarse-grained residue-based models of disordered protein condensates: utility and limitations of simple charge pattern parameters

Biomolecular condensates undergirded by phase separations of proteins and nucleic acids serve crucial biological functions. To gain physical insights into their genetic basis, we study how liquid-liquid phase separation (LLPS) of intrinsically disordered proteins (IDPs) depends on their sequence charge patterns using a continuum Langevin chain model wherein each amino acid residue is represented by a single bead. Charge patterns are characterized by the "blockiness" measure κ and the "sequence charge decoration" (SCD) parameter. Consistent with random phase approximation (RPA) theory and lattice simulations, LLPS propensity as characterized by critical temperature Tcr* increases with increasingly negative SCD for a set of sequences showing a positive correlation between κ and -SCD. Relative to RPA, the simulated sequence-dependent variation in Tcr* is often-though not always-smaller, whereas the simulated critical volume fractions are higher. However, for a set of sequences exhibiting an anti-correlation between κ and -SCD, the simulated Tcr*'s are quite insensitive to either parameter. Additionally, we find that blocky sequences that allow for strong electrostatic repulsion can lead to coexistence curves with upward concavity as stipulated by RPA, but the LLPS propensity of a strictly alternating charge sequence was likely overestimated by RPA and lattice models because interchain stabilization of this sequence requires spatial alignments that are difficult to achieve in real space. These results help delineate the utility and limitations of the charge pattern parameters and of RPA, pointing to further efforts necessary for rationalizing the newly observed subtleties.

Preface: Special Topic on Single-Molecule Biophysics

Single-molecule measurements are now almost routinely used to study biological systems and processes. The scope of this special topic emphasizes the physics side of single-molecule observations, with the goal of highlighting new developments in physical techniques as well as conceptual insights that single-molecule measurements bring to biophysics. This issue also comprises recent advances in theoretical physical models of single-molecule phenomena, interpretation of single-molecule signals, and fundamental areas of statistical mechanics that are related to single-molecule observations. A particular goal is to illustrate the increasing synergy between theory, simulation, and experiment in single-molecule biophysics.

Charge fluctuation effects on the shape of flexible polyampholytes with applications to intrinsically disordered proteins

Random polyampholytes (PAs) contain positively and negatively charged monomers that are distributed randomly along the polymer chain. The interaction between charges is assumed to be given by the Debye-Huckel potential. We show that the size of the PA is determined by an interplay between electrostatic interactions, giving rise to the polyelectrolyte effect due to net charge per monomer (σ) and an effective attractive PA interaction due to charge fluctuations, δσ. The interplay between these terms gives rise to non-monotonic dependence of the radius of gyration, R _g , on the inverse Debye length, κ, when PA effects are important ( δ σ σ > 1 ). In the opposite limit, R _g decreases monotonically with increasing κ. Simulations of PA chains, using a charged bead-spring model, further corroborate our theoretical predictions. The simulations unambiguously show that conformational heterogeneity manifests itself among sequences that have identical PA parameters. A clear implication is that the phases of PA sequences, and by inference intrinsically disordered proteins (IDPs), cannot be determined using only the bare PA parameters (σ and δσ). The theory is used to calculate the changes in R _g on N, the number of residues for a set of IDPs. For a certain class of IDPs, with N between 24 and 441, the size grows as R _g ∼ N ^0.6, which agrees with data from small angle X-ray scattering experiments.

Conformational preferences and phase behavior of intrinsically disordered low complexity sequences: insights from multiscale simulations

While many proteins and protein regions utilize a complex repertoire of amino acids to achieve their biological function, a subset of protein sequences are enriched in a reduced set of amino acids. These so-called low complexity (LC) sequences, specifically intrinsically disordered variants of LC sequences, have been the focus of recent investigations owing to their roles in a range of biological functions, specifically phase separation. Computational studies of LC sequences have provided rich insights into their behavior both as individual proteins in dilute solutions and as the drivers and modulators of higher-order assemblies. Here, we review how simulations performed across distinct resolutions have provided different types of insights into the biological role of LC sequences.

Modulating charge patterning and ionic strength as a strategy to induce conformational changes in intrinsically disordered proteins

We present an analytical theory to describe conformational changes as a function of salt for polymers with a given sequence of charges. We apply this model to describe Intrinsically Disordered Proteins (IDPs) by explicitly accounting for charged residues and their exact placement in the primary sequence while approximating the effect of non-electrostatic interactions at a mean-field level by effective short-range (two body and three-body) interaction parameters. The effect of ions is introduced by treating electrostatic interactions within Debye-Huckle approximation. Using typical values of the short-range mean-field parameters derived from all-atom Monte Carlo simulations (at zero salt), we predict the conformational changes as a function of salt concentration. We notice that conformational transitions in response to changes in ionic strength strongly depend on sequence specific charge patterning. For example, globule to coil transition can be observed upon increasing salt concentration, in stark contrast to uniformly charged polyelectrolyte theories based on net charge only. In addition, it is possible to observe non-monotonic behavior with salt as well. Drastic differences in salt-induced conformational transitions is also evident between two doubly phosphorylated sequences-derived from the same wild type sequence-that only differ in the site of phosphorylation. Similar effects are also predicted between two sequences derived from the same parent sequence differing by a single site mutation where a negative charge is replaced by a positive charge. These effects are purely a result of charge decoration and can only be understood in terms of metrics based on specific placement of charges, and cannot be explained by models based on charge composition alone. Identifying sequences and hot spots within sequences-for post translational modification or charge mutation-using our high-throughput theory will yield fundamental insights into design and biological regulation mediated by phosphorylation and/or local changes in salt concentration.

The Unfolded State of the C-Terminal Domain of L9 Expands at Low but Not at Elevated Temperatures

The temperature dependence of the overall dimensions of the denatured state ensemble (DSE) of proteins remains unclear. Some studies indicate compaction of the DSE at high temperatures, whereas others argue that dimensions do not decrease. The degree of compaction or expansion in the cold-denatured state has been less studied. To investigate the temperature dependence of unfolded state dimensions, small angle x-ray scattering measurements were performed in native buffer in the absence of denaturant for a designed point mutant of the C-terminal domain of L9, a small cooperatively folded α-β protein, at 14 different temperatures over the range of 5-60°C. The I98A mutation destabilizes the domain such that both the DSE and the folded state are populated at 25°C in the absence of denaturant or extreme pH. Thermal unfolding as well as cold unfolding can thus be observed in the absence of denaturant, allowing a direct comparison of these regimes for the same protein using the same technique. The temperature of maximal stability, T_s, is 30°C. There is no detectable change in R_g of the unfolded state as the temperature is increased above T_s, but a clear expansion is detected as the temperature is decreased below T_s. The R_g of the DSE populated in buffer was found to be 27.8 ± 1.7 Å at 5°C, 21.8 ± 1.9 Å at 30°C, and 21.7 ± 2.0 Å at 60°C. In contrast, no significant temperature dependence was observed for the value of R_g measured in 6 M guanidine hydrochloride. The small angle x-ray scattering data reported here indicate clear differences between the cold- and thermal-unfolded states and show that there is no significant compaction at elevated temperatures.