Quantifying charge state heterogeneity for proteins with multiple ionizable residues

A pH-Centric Model of Nucleolar Activity and Regulation

The nucleolus is essential for the efficient and accurate production of ribosomal subunits, which are crucial for assembling ribosomes-the cellular machinery responsible for protein synthesis. Emerging insights into its liquid-like nature have shed new light on the role of its unique biophysical properties in the activity and regulation of this organelle. In this perspective, we examine recent insights into nucleolar biophysical homeostasis, with a focus on its regulation as an acidic biomolecular condensate. We review current evidence on how nucleolar composition and biochemical activities could generate and maintain a proton gradient. Additionally, we propose an integrative model explaining how nucleolar acidity contributes to homeostasis at a molecular level, providing a unified framework for its role in health and disease.

MEDOC: A Fast, Scalable, and Mathematically Exact Algorithm for the Site-Specific Prediction of the Protonation Degree in Large Disordered Proteins

Intrinsically disordered regions are found in most eukaryotic proteins and are enriched with positively and negatively charged residues. While it is often convenient to assume that these residues follow their model-compound pKa values, recent work has shown that local charge effects (charge regulation) can upshift or downshift side chain pKa values with major consequences for molecular function. Despite this, charge regulation is rarely considered when investigating disordered regions. The number of potential charge microstates that can be populated through acid/base regulation of a given number of ionizable residues in a sequence, N, scales as ∼2N. This exponential scaling makes the assessment of the full charge landscape of most proteins computationally intractable. To address this problem, we developed "multisite extent of deprotonation originating from context" (MEDOC) to determine the degree of protonation of a protein based on the local sequence context of each ionizable residue. We show that we can drastically reduce the number of parameters necessary to determine the full, analytical Boltzmann partition function of the charge landscape at both global and site-specific levels. Our algorithm applies the structure of the q-canonical ensemble, combined with novel strategies to rapidly obtain the minimal set of parameters, thereby circumventing the combinatorial explosion of the number of charge microstates even for proteins containing a large number of ionizable amino acids. We apply MEDOC to several sequences, including a global analysis of the distribution of pKa values across the entire DisProt database. Our results show differences in the distribution of predicted pKa values for different amino acids and good agreement with NMR-measured pKa values in proteins.

Liquid-liquid phase separation driven by charge heterogeneity

Despite the intrinsic charge heterogeneity of proteins plays a crucial role in the liquid-liquid phase separation (LLPS) of a broad variety of protein systems, our understanding of the effects of their electrostatic anisotropy is still in its early stages. We approach this issue by means of a coarse-grained model based on a robust mean-field description that extends the DLVO theory to non-uniformly charged particles. We numerically investigate the effect of surface charge patchiness and net particle charge on varying these features independently and with the use of a few parameters only. The effect of charge anisotropy on the LLPS critical point is rationalized via a thermodynamic-independent parameter based on orientationally averaged pair properties, that estimates the particle connectivity and controls the propensity of the liquid phase to condensate. We show that, even though directional attraction alone is able to lower the particle bonding valence-thus shifting the critical point towards lower temperatures and densities-directional repulsion significantly and systematically diminishes the particle functionality, thus further reducing the critical parameters. This electrostatically-driven shift can be understood in terms of the additional morphological constraints introduced by the directional repulsion, that hinder the condensation of dense aggregates.

Phosphorylation of disordered proteins tunes local and global intramolecular interactions

Protein post-translational modifications, such as phosphorylation, are important regulatory signals for diverse cellular functions. In particular, intrinsically disordered protein regions (IDRs) are subject to phosphorylation as a means to modulate their interactions and functions. Toward understanding the relationship between phosphorylation in IDRs and specific functional outcomes, we must consider how phosphorylation affects the IDR conformational ensemble. Various experimental techniques are suited to interrogate the features of IDR ensembles; molecular simulations can provide complementary insights and even illuminate ensemble features that may be experimentally inaccessible. Therefore, we sought to expand the tools available to study phosphorylated IDRs by all-atom Monte Carlo simulations. To this end, we implemented parameters for phosphoserine (pSer) and phosphothreonine (pThr) into the OPLS version of the continuum solvent model, ABSINTH, and assessed their performance in all-atom simulations compared with published findings. We simulated short (<20 residues) and long (>80 residues) phospho-IDRs that, collectively, survey both local and global phosphorylation-induced changes to the ensemble. Our simulations of four well-studied phospho-IDRs show near-quantitative agreement with published findings for these systems via metrics including changes to radius of gyration, transient helicity, and persistence length. We also leveraged the inherent advantage of sequence control in molecular simulations to explore the conformational effects of diverse combinations of phospho-sites in two multiphosphorylated IDRs. Our results support and expand on previous observations that connect phosphorylation to changes in the IDR conformational ensemble. Herein, we describe phosphorylation as a means to alter sequence chemistry, net charge and charge patterning, and intramolecular interactions, which can collectively modulate the local and global IDR ensemble features.

Emergent microenvironments of nucleoli

In higher eukaryotes, the nucleolus harbors at least three sub-phases that facilitate multiple functionalities including ribosome biogenesis. The three prominent coexisting sub-phases are the fibrillar center (FC), the dense fibrillar component (DFC), and the granular component (GC). Here, we review recent efforts in profiling sub-phase compositions that shed light on the types of physicochemical properties that emerge from compositional biases and territorial organization of specific types of macromolecules. We highlight roles played by molecular grammars which refers to protein sequence features including the substrate binding domains, the sequence features of intrinsically disordered regions, and the multivalence of these distinct types of domains / regions. We introduce the concept of a barcode of emergent physicochemical properties of nucleoli. Although our knowledge of the full barcode remains incomplete, we hope that the concept prompts investigations into undiscovered emergent properties and engenders an appreciation for how and why unique microenvironments control biochemical reactions.

A Few Charged Residues in Galectin-3's Folded and Disordered Regions Regulate Phase Separation

Proteins with intrinsically disordered regions (IDRs) often undergo phase separation to control their functions spatiotemporally. Changing the pH alters the protonation levels of charged sidechains, which in turn affects the attractive or repulsive force for phase separation. In a cell, the rupture of membrane-bound compartments, such as lysosomes, creates an abrupt change in pH. However, how proteins' phase separation reacts to different pH environments remains largely unexplored. Here, using extensive mutagenesis, NMR spectroscopy, and biophysical techniques, it is shown that the assembly of galectin-3, a widely studied lysosomal damage marker, is driven by cation-π interactions between positively charged residues in its folded domain with aromatic residues in the IDR in addition to π-π interaction between IDRs. It is also found that the sole two negatively charged residues in its IDR sense pH changes for tuning the condensation tendency. Also, these two residues may prevent this prion-like IDR domain from forming rapid and extensive aggregates. These results demonstrate how cation-π, π-π, and electrostatic interactions can regulate protein condensation between disordered and structured domains and highlight the importance of sparse negatively charged residues in prion-like IDRs.

Biomolecular Condensates are Characterized by Interphase Electric Potentials

Biomolecular condensates form via processes that combine phase separation and reversible associations of multivalent macromolecules. Condensates can be two- or multiphase systems defined by coexisting dense and dilute phases. Here, we show that solution ions partition asymmetrically across coexisting phases defined by condensates formed by intrinsically disordered proteins or homopolymeric RNA molecules. Our findings were enabled by direct measurements of the activities of cations and anions within coexisting phases of protein and RNA condensates. Asymmetries in ion partitioning between coexisting phases vary with protein sequence, macromolecular composition, salt concentration, and ion type. The Donnan equilibrium set up by the asymmetrical partitioning of solution ions generates interphase electric potentials known as Donnan and Nernst potentials. Our measurements show that the interphase potentials of condensates are of the same order of magnitude as membrane potentials of membrane-bound organelles. Interphase potentials quantify the degree to which microenvironments of coexisting phases are different from one another. Importantly, and based on condensate-specific interphase electric potentials, we reason that condensates are akin to capacitors that store charge. Interphase potentials should lead to electric double layers at condensate interfaces, thereby explaining recent observations of condensate interfaces being electrochemically active.

MEDOC: A fast, scalable and mathematically exact algorithm for the site-specific prediction of the protonation degree in large disordered proteins

Intrinsically disordered regions are found in most eukaryotic proteins and are enriched in positively and negatively charged residues. While it is often convenient to assume these residues follow their model-compound pKa values, recent work has shown that local charge effects (charge regulation) can upshift or downshift sidechain pKa values with major consequences for molecular function. Despite this, charge regulation is rarely considered when investigating disordered regions. The number of potential charge microstates that can be populated through acid/base regulation of a given number of ionizable residues in a sequence, N , scales as~ 2 N . This exponential scaling makes the assessment of the full charge landscape of most proteins computationally intractable. To address this problem, we developed MEDOC (Multisite Extent of Deprotonation Originating from Context) to determine the degree of protonation of a protein based on the local sequence context of each ionizable residue. We show that we can drastically reduce the number of parameters necessary to determine the full, analytic, Boltzmann partition function of the charge landscape at both global and site-specific levels. Our algorithm applies the structure of the q-canonical ensemble, combined with novel strategies to rapidly obtain the minimal set of parameters, thereby circumventing the combinatorial explosion of the number of charge microstates even for proteins containing a large number of ionizable amino acids. We apply MEDOC to several sequences, including a global analysis of the distribution of pKa values across the entire DisProt database. Our results show differences in the distribution of predicted pKa values for different amino acids, in agreement with NMR-measured distributions in proteins.

Beyond monopole electrostatics in regulating conformations of intrinsically disordered proteins

Conformations and dynamics of an intrinsically disordered protein (IDP) depend on its composition of charged and uncharged amino acids, and their specific placement in the protein sequence. In general, the charge (positive or negative) on an amino acid residue in the protein is not a fixed quantity. Each of the ionizable groups can exist in an equilibrated distribution of fully ionized state (monopole) and an ion-pair (dipole) state formed between the ionizing group and its counterion from the background electrolyte solution. The dipole formation (counterion condensation) depends on the protein conformation, which in turn depends on the distribution of charges and dipoles on the molecule. Consequently, effective charges of ionizable groups in the IDP backbone may differ from their chemical charges in isolation-a phenomenon termed charge-regulation. Accounting for the inevitable dipolar interactions, that have so far been ignored, and using a self-consistent procedure, we present a theory of charge-regulation as a function of sequence, temperature, and ionic strength. The theory quantitatively agrees with both charge reduction and salt-dependent conformation data of Prothymosin-alpha and makes several testable predictions. We predict charged groups are less ionized in sequences where opposite charges are well mixed compared to sequences where they are strongly segregated. Emergence of dipolar interactions from charge-regulation allows spontaneous coexistence of two phases having different conformations and charge states, sensitively depending on the charge patterning. These findings highlight sequence dependent charge-regulation and its potential exploitation by biological regulators such as phosphorylation and mutations in controlling protein conformation and function.

Biomolecular condensates are characterized by interphase electric potentials

Biomolecular condensates form via processes that combine phase separation and reversible associations of multivalent macromolecules. Condensates can be two- or multi-phase systems defined by coexisting dense and dilute phases. Here, we show that solution ions can partition asymmetrically across coexisting phases defined by condensates formed by intrinsically disordered proteins or homopolymeric RNA molecules. Our findings were enabled by direct measurements of the activities of cations and anions within coexisting phases of protein and RNA condensates. Asymmetries in ion partitioning between coexisting phases vary with protein sequence, condensate type, salt concentration, and ion type. The Donnan equilibrium set up by asymmetrical partitioning of solution ions generates interphase electric potentials known as Donnan and Nernst potentials. Our measurements show that the interphase potentials of condensates are of the same order of magnitude as membrane potentials of membrane-bound organelles. Interphase potentials quantify the degree to which microenvironments of coexisting phases are different from one another. Importantly, and based on condensate-specific interphase electric potentials, which are membrane-like potentials of membraneless bodies, we reason that condensates are mesoscale capacitors that store charge. Interphase potentials lead to electric double layers at condensate interfaces. This helps explain recent observations of condensate interfaces being electrochemically active.

Phosphorylation of disordered proteins tunes local and global intramolecular interactions

Protein post-translational modifications, such as phosphorylation, are important regulatory signals for diverse cellular functions. In particular, intrinsically disordered protein regions (IDRs) are subject to phosphorylation as a means to modulate their interactions and functions. Toward understanding the relationship between phosphorylation in IDRs and specific functional outcomes, we must consider how phosphorylation affects the IDR conformational ensemble. Various experimental techniques are suited to interrogate the features of IDR ensembles; molecular simulations can provide complementary insights and even illuminate ensemble features that may be experimentally inaccessible. Therefore, we sought to expand the tools available to study phosphorylated IDRs by all-atom Monte Carlo simulations. To this end, we implemented parameters for phosphoserine (pSer) and phosphothreonine (pThr) into the OPLS version of the continuum solvent model, ABSINTH, and assessed their performance in all-atom simulations compared to published findings. We simulated short (< 20 residues) and long (> 80 residues) phospho-IDRs that, collectively, survey both local and global phosphorylation-induced changes to the ensemble. Our simulations of four well-studied phospho-IDRs show near-quantitative agreement with published findings for these systems via metrics including changes to radius of gyration, transient helicity, and persistence length. We also leveraged the inherent advantage of sequence control in molecular simulations to explore the conformational effects of diverse combinations of phospho-sites in two multi-phosphorylated IDRs. Our results support and expand on prior observations that connect phosphorylation to changes in the IDR conformational ensemble. Herein, we describe phosphorylation as a means to alter sequence chemistry, net charge and charge patterning, and intramolecular interactions, which can collectively modulate the local and global IDR ensemble features.

Direct prediction of intermolecular interactions driven by disordered regions

Intrinsically disordered regions (IDRs) are critical for a wide variety of cellular functions, many of which involve interactions with partner proteins. Molecular recognition is typically considered through the lens of sequence-specific binding events. However, a growing body of work has shown that IDRs often interact with partners in a manner that does not depend on the precise order of the amino acid order, instead driven by complementary chemical interactions leading to disordered bound-state complexes. Despite this emerging paradigm, we lack tools to describe, quantify, predict, and interpret these types of structurally heterogeneous interactions from the underlying amino acid sequences. Here, we repurpose the chemical physics developed originally for molecular simulations to develop an approach for predicting intermolecular interactions between IDRs and partner proteins. Our approach enables the direct prediction of phase diagrams, the identification of chemically-specific interaction hotspots on IDRs, and a route to develop and test mechanistic hypotheses regarding IDR function in the context of molecular recognition. We use our approach to examine a range of systems and questions to highlight its versatility and applicability.

Macromolecular condensation organizes nucleolar sub-phases to set up a pH gradient

Nucleoli are multicomponent condensates defined by coexisting sub-phases. We identified distinct intrinsically disordered regions (IDRs), including acidic (D/E) tracts and K-blocks interspersed by E-rich regions, as defining features of nucleolar proteins. We show that the localization preferences of nucleolar proteins are determined by their IDRs and the types of RNA or DNA binding domains they encompass. In vitro reconstitutions and studies in cells showed how condensation, which combines binding and complex coacervation of nucleolar components, contributes to nucleolar organization. D/E tracts of nucleolar proteins contribute to lowering the pH of co-condensates formed with nucleolar RNAs in vitro. In cells, this sets up a pH gradient between nucleoli and the nucleoplasm. By contrast, juxta-nucleolar bodies, which have different macromolecular compositions, featuring protein IDRs with very different charge profiles, have pH values that are equivalent to or higher than the nucleoplasm. Our findings show that distinct compositional specificities generate distinct physicochemical properties for condensates.

Electric Potential at the Interface of Membraneless Organelles Gauged by Graphene

Eukaryotic cells contain membrane-bound and membrane-less organelles that are often in contact with each other. How the interface properties of membrane-less organelles regulate their interactions with membranes remains challenging to assess. Here, we employ graphene-based sensors to investigate the electrostatic properties of synapsin 1, a major synaptic phosphoprotein, either in a single phase or after undergoing phase separation to form synapsin condensates. Using these graphene-based sensors, we discover that synapsin condensates generate strong electrical responses that are otherwise absent when synapsin is present as a single phase. By introducing atomically thin dielectric barriers, we show that the electrical response originates in an electric double layer whose formation governs the interaction between synapsin condensates and graphene. Our data indicate that the interface properties of the same protein are substantially different when the protein is in a single phase versus within a biomolecular condensate, unraveling that condensates can harbor ion potential differences at their interface.

Intrinsically disordered regions are poised to act as sensors of cellular chemistry

Intrinsically disordered proteins and protein regions (IDRs) are abundant in eukaryotic proteomes and play a wide variety of essential roles. Instead of folding into a stable structure, IDRs exist in an ensemble of interconverting conformations whose structure is biased by sequence-dependent interactions. The absence of a stable 3D structure, combined with high solvent accessibility, means that IDR conformational biases are inherently sensitive to changes in their environment. Here, we argue that IDRs are ideally poised to act as sensors and actuators of cellular physicochemistry. We review the physical principles that underlie IDR sensitivity, the molecular mechanisms that translate this sensitivity to function, and recent studies where environmental sensing by IDRs may play a key role in their downstream function.

Counterion-controlled phase equilibria in a charge-regulated polymer solution

We study phase equilibria in a minimal model of charge-regulated polymer solutions. Our model consists of a single polymer species whose charge state arises from protonation-deprotonation processes in the presence of a dissolved acid, whose anions serve as screening counterions. We explicitly account for variability in the polymers' charge states. Homogeneous equilibria in this model system are characterised by the total concentration of polymers, the concentration of counter-ions and the charge distributions of polymers which can be computed with the help of analytical approximations. We use these analytical results to characterise how parameter values and solution acidity influence equilibrium charge distributions and identify for which regimes uni-modal and multi-modal charge distributions arise. We then study the interplay between charge regulation, solution acidity and phase separation. We find that charge regulation has a significant impact on polymer solubility and allows for non-linear responses to the solution acidity: Re-entrant phase behaviour is possible in response to increasing solution acidity. Moreover, we show that phase separation can yield to the coexistence of local environments characterised by different charge distributions.

Driving forces of the complex formation between highly charged disordered proteins

Highly disordered complexes between oppositely charged intrinsically disordered proteins present a new paradigm of biomolecular interactions. Here, we investigate the driving forces of such interactions for the example of the highly positively charged linker histone H1 and its highly negatively charged chaperone, prothymosin α (ProTα). Temperature-dependent single-molecule Förster resonance energy transfer (FRET) experiments and isothermal titration calorimetry reveal ProTα-H1 binding to be enthalpically unfavorable, and salt-dependent affinity measurements suggest counterion release entropy to be an important thermodynamic driving force. Using single-molecule FRET, we also identify ternary complexes between ProTα and H1 in addition to the heterodimer at equilibrium and show how they contribute to the thermodynamics observed in ensemble experiments. Finally, we explain the observed thermodynamics quantitatively with a mean-field polyelectrolyte theory that treats counterion release explicitly. ProTα-H1 complex formation resembles the interactions between synthetic polyelectrolytes, and the underlying principles are likely to be of broad relevance for interactions between charged biomolecules in general.

Pervasive, conserved secondary structure in highly charged protein regions

Understanding how protein sequences confer function remains a defining challenge in molecular biology. Two approaches have yielded enormous insight yet are often pursued separately: structure-based, where sequence-encoded structures mediate function, and disorder-based, where sequences dictate physicochemical and dynamical properties which determine function in the absence of stable structure. Here we study highly charged protein regions (>40% charged residues), which are routinely presumed to be disordered. Using recent advances in structure prediction and experimental structures, we show that roughly 40% of these regions form well-structured helices. Features often used to predict disorder-high charge density, low hydrophobicity, low sequence complexity, and evolutionarily varying length-are also compatible with solvated, variable-length helices. We show that a simple composition classifier predicts the existence of structure far better than well-established heuristics based on charge and hydropathy. We show that helical structure is more prevalent than previously appreciated in highly charged regions of diverse proteomes and characterize the conservation of highly charged regions. Our results underscore the importance of integrating, rather than choosing between, structure- and disorder-based approaches.

Dissecting the phase separation and oligomerization activities of the carboxysome positioning protein McdB

Across bacteria, protein-based organelles called bacterial microcompartments (BMCs) encapsulate key enzymes to regulate their activities. The model BMC is the carboxysome that encapsulates enzymes for CO2 fixation to increase efficiency and is found in many autotrophic bacteria, such as cyanobacteria. Despite their importance in the global carbon cycle, little is known about how carboxysomes are spatially regulated. We recently identified the two-factor system required for the maintenance of carboxysome distribution (McdAB). McdA drives the equal spacing of carboxysomes via interactions with McdB, which associates with carboxysomes. McdA is a ParA/MinD ATPase, a protein family well studied in positioning diverse cellular structures in bacteria. However, the adaptor proteins like McdB that connect these ATPases to their cargos are extremely diverse. In fact, McdB represents a completely unstudied class of proteins. Despite the diversity, many adaptor proteins undergo phase separation, but functional roles remain unclear. Here, we define the domain architecture of McdB from the model cyanobacterium Synechococcus elongatus PCC 7942, and dissect its mode of biomolecular condensate formation. We identify an N-terminal intrinsically disordered region (IDR) that modulates condensate solubility, a central coiled-coil dimerizing domain that drives condensate formation, and a C-terminal domain that trimerizes McdB dimers and provides increased valency for condensate formation. We then identify critical basic residues in the IDR, which we mutate to glutamines to solubilize condensates. Finally, we find that a condensate-defective mutant of McdB has altered association with carboxysomes and influences carboxysome enzyme content. The results have broad implications for understanding spatial organization of BMCs and the molecular grammar of protein condensates.

Phase Transitions of Associative Biomacromolecules

Multivalent proteins and nucleic acids, collectively referred to as multivalent associative biomacromolecules, provide the driving forces for the formation and compositional regulation of biomolecular condensates. Here, we review the key concepts of phase transitions of aqueous solutions of associative biomacromolecules, specifically proteins that include folded domains and intrinsically disordered regions. The phase transitions of these systems come under the rubric of coupled associative and segregative transitions. The concepts underlying these processes are presented, and their relevance to biomolecular condensates is discussed.

Combining Experiments and Simulations to Examine the Temperature-Dependent Behavior of a Disordered Protein

Intrinsically disordered proteins are a class of proteins that lack stable folded conformations and instead adopt a range of conformations that determine their biochemical functions. The temperature-dependent behavior of such disordered proteins is complex and can vary depending on the specific protein and environment. Here, we have used molecular dynamics simulations and previously published experimental data to investigate the temperature-dependent behavior of histatin 5, a 24-residue-long polypeptide. We examined the hypothesis that histatin 5 undergoes a loss of polyproline II (PPII) structure with increasing temperature, leading to more compact conformations. We found that the conformational ensembles generated by the simulations generally agree with small-angle X-ray scattering data for histatin 5, but show some discrepancies with the hydrodynamic radius as probed by pulsed-field gradient NMR spectroscopy, and with the secondary structure information derived from circular dichroism. We attempted to reconcile these differences by reweighting the conformational ensembles against the scattering and NMR data. By doing so, we were in part able to capture the temperature-dependent behavior of histatin 5 and to link the observed decrease in hydrodynamic radius with increasing temperature to a loss of PPII structure. We were, however, unable to achieve agreement with both the scattering and NMR data within experimental errors. We discuss different possible reasons for this including inaccuracies in the force field, differences in conditions of the NMR and scattering experiments, and issues related to the calculation of the hydrodynamic radius from conformational ensembles. Our study highlights the importance of integrating multiple types of experimental data when modeling conformational ensembles of disordered proteins and how environmental factors such as the temperature influence them.

Modulation of Functional Phosphorylation Sites by Basic Residues in the Unique Domain of c-Src

In contrast to the well-studied canonical regulatory mechanisms, the way by which the recently discovered Src N-terminal regulatory element (SNRE) modulates Src activity is not yet well understood. Phosphorylation of serine and threonine residues modulates the charge distribution along the disordered region of the SNRE and may affect a fuzzy complex with the SH3 domain that is believed to act as an information transduction element. The pre-existing positively charged sites can interact with the newly introduced phosphate groups by modulating their acidity, introducing local conformational restrictions, or by coupling various phosphosites into a functional unit. In this paper, we use pH-dependent NMR measurements combined with single point mutations to identify the interactions of basic residues with physiologically important phosphorylated residues and to characterize the effect of these interactions in neighbor residues, thus providing insight into the electrostatic network in the isolated disordered regions and in the entire SNRE. From a methodological point of view, the linear relationships observed between the mutation-induced pKa changes of the phosphate groups of phosphoserine and phosphothreonine and the pH-induced chemical shifts of the NH groups of these residues provide a very convenient alternative to identify interacting phosphate groups without the need to introduce point mutations on specific basic residues.

Uncovering the Contributions of Charge Regulation to the Stability of Single Alpha Helices

The single alpha helix (SAH) is a recurring motif in biology. The consensus sequence has a di-block architecture that includes repeats of four consecutive glutamate residues followed by four consecutive lysine residues. Measurements show that the overall helicity of sequences with consensus E4 K4 repeats is insensitive to a wide range of pH values. Here, we use the recently introduced q-canonical ensemble, which allows us to decouple measurements of charge state and conformation, to explain the observed insensitivity of SAH helicity to pH. We couple the outputs from separate measurements of charge and conformation with atomistic simulations to derive residue-specific quantifications of preferences for being in an alpha helix and for the ionizable residues to be charged vs. uncharged. We find a clear preference for accommodating uncharged Glu residues within internal positions of SAH-forming sequences. The stabilities of alpha helical conformations increase with the number of E4 K4 repeats and so do the numbers of accessible charge states that are compatible with forming conformations of high helical content. There is conformational buffering whereby charge state heterogeneity buffers against large-scale conformational changes thus making the overall helicity insensitive to large changes in pH. Further, the results clearly argue against a single, rod-like alpha helical conformation being the only or even dominant conformation in the ensembles of so-called SAH sequences.

Pervasive, conserved secondary structure in highly charged protein regions

Understanding how protein sequences confer function remains a defining challenge in molecular biology. Two approaches have yielded enormous insight yet are often pursued separately: structure-based, where sequence-encoded structures mediate function, and disorder-based, where sequences dictate physicochemical and dynamical properties which determine function in the absence of stable structure. Here we study highly charged protein regions (>40% charged residues), which are routinely presumed to be disordered. Using recent advances in structure prediction and experimental structures, we show that roughly 40% of these regions form well-structured helices. Features often used to predict disorder-high charge density, low hydrophobicity, low sequence complexity, and evolutionarily varying length-are also compatible with solvated, variable-length helices. We show that a simple composition classifier predicts the existence of structure far better than well-established heuristics based on charge and hydropathy. We show that helical structure is more prevalent than previously appreciated in highly charged regions of diverse proteomes and characterize the conservation of highly charged regions. Our results underscore the importance of integrating, rather than choosing between, structure- and disorder-based approaches.

The biophysics of disordered proteins from the point of view of single-molecule fluorescence spectroscopy

Intrinsically disordered proteins (IDPs) and regions (IDRs) have emerged as key players across many biological functions and diseases. Differently from structured proteins, disordered proteins lack stable structure and are particularly sensitive to changes in the surrounding environment. Investigation of disordered ensembles requires new approaches and concepts for quantifying conformations, dynamics, and interactions. Here, we provide a short description of the fundamental biophysical properties of disordered proteins as understood through the lens of single-molecule fluorescence observations. Single-molecule Förster resonance energy transfer (FRET) and fluorescence correlation spectroscopy (FCS) provides an extensive and versatile toolbox for quantifying the characteristics of conformational distributions and the dynamics of disordered proteins across many different solution conditions, both in vitro and in living cells.

A conceptual framework for understanding phase separation and addressing open questions and challenges

Macromolecular phase separation is being recognized for its potential importance and relevance as a driver of spatial organization within cells. Here, we describe a framework based on synergies between networking (percolation or gelation) and density (phase separation) transitions. Accordingly, the phase transitions in question are referred to as phase separation coupled to percolation (PSCP). The condensates that result from PSCP are viscoelastic network fluids. Such systems have sequence-, composition-, and topology-specific internal network structures that give rise to time-dependent interplays between viscous and elastic properties. Unlike pure phase separation, the process of PSCP gives rise to sequence-, chemistry-, and structure-specific distributions of clusters that can form at concentrations that lie well below the threshold concentration for phase separation. PSCP, influenced by specific versus solubility-determining interactions, also provides a bridge between different observations and helps answer questions and address challenges that have arisen regarding the role of macromolecular phase separation in biology.

Competing interactions give rise to two-state behavior and switch-like transitions in charge-rich intrinsically disordered proteins

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 intrachain 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 side chain 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.