Perplexing cooperative folding and stability of a low-sequence complexity, polyproline 2 protein lacking a hydrophobic core

Collapse and Protein Folding: Should We Be Surprised That Biothermodynamics Works So Well?

A complete understanding of protein function and dynamics requires the characterization of the multiple thermodynamic states, including the denatured state ensemble (DSE). Whereas residual structure in the DSE (as well as in partially folded states) is pertinent in many biological contexts, here we are interested in how such structure affects protein thermodynamics. We examine issues related to chain collapse in light of new developments, focusing on potential complications arising from differences in the DSE's properties under various conditions. Despite some variability in the degree of collapse and structure in the DSE, stability measurements are remarkably consistent between two standard methods, calorimetry and chemical denaturation, as well as with hydrogen-deuterium exchange. This robustness is due in part to the DSEs obtained with different perturbations being thermodynamically equivalent and hence able to serve as a common reference state. An examination of the properties of the DSE points to it as being a highly expanded ensemble with minimal amounts of stable hydrogen bonded structure. These two features are likely to be critical in the broad and successful application of thermodynamics to protein folding. Our review concludes with a discussion of the impact of these findings on folding mechanisms and pathways.

Amino Acid Residue-Specific Ramachandran Distributions Derived from a Simple Mean Field Potential

Protein dynamics in the unfolded state, in the context of early stage protein folding or intrinsically disordered proteins (IDPs), is not well understood. The discovery of IDPs, and their sequence-dependent dynamics, has led to many computational and experimental investigations regarding the conformational preferences of short oligopeptides and individual amino acid residues in the unfolded state. As proteins consist of sequences of amino acid residues, characterizing the intrinsic conformational preferences of the individual residues in the unfolded state is crucial for understanding the emergent conformations of peptides and proteins. While advances have been made in understanding conformational preferences, the atomistic mechanisms driving these preferences remain unresolved. In this work, we show that the distributions of atomic overlaps between backbone and side chain atoms in Ramachandran space are unique for amino acid residue mimetic structures alanine, valine, leucine, and isoleucine in Ramachandran space indicating unique intrapeptide energy landscapes for each residue. We then construct a mean field potential consisting of only an empirical peptide backbone-water and average intrapeptide Lennard-Jones contributions to explore their influence on the conformational preferences. With this fairly simple model, we were able to produce Ramachandran distributions that qualitatively agree with previously reported experimental and computational predictions about the conformational preferences of these amino acid residues in the unfolded state in water. Our results indicate these conformational preferences are the result of the balance between pPII-stabilizing backbone-water interactions and repulsive side chain-backbone interactions where the latter will depend uniquely on the atomic makeup and geometry of the side chain.

Hydrogen bonding patterns and cooperativity in polyproline II helical bundles

Hydrogen bond cooperativity (HBC) plays an important role in stabilizing protein assemblies built by α-helices and β-sheets, the most common secondary structures. However, whether HBC exists in other types of protein secondary structures such as polyproline II (PPII) helices remains unexplored. This is intriguing, since PPII systems as assembling blocks are continuously emerging across multiple fields. Here, using a combination of computational chemistry tools and molecular modeling corroborated by experimental observables, we characterize the distinct H-bonding patterns present in PPII helical bundles and establish that HBC stabilizes intermolecular PPII helices as seen in other protein assemblies such as amyloid fibrils. In addition to cooperative interactions in canonical CO···HN H-bonds, we show that analogous interactions in non-canonical CO···HαCα H-bonds are relevant in Gly-rich PPII bundles, thus compensating for the inability of glycine residues to create hydrophobic cores. Our results provide a mechanistic explanation for the assembly of these bundles.

Architectonic principles of polyproline II helix bundle protein domains

Glycine rich polyproline II helix assemblies are an emerging class of natural domains found in several proteins with different functions and diverse origins. The distinct properties of these domains relative to those composed of α-helices and β-sheets could make glycine-rich polyproline II helix assemblies a useful building block for protein design. Whereas the high population of polyproline II conformers in disordered state ensembles could facilitate glycine-rich polyproline II helix folding, the architectonic bases of these structures are not well known. Here, we compare and analyze their structures to uncover common features. These protein domains are found to be highly tolerant of distinct flanking sequences. This speaks to the robustness of this fold and strongly suggests that glycine rich polyproline II assemblies could be grafted with other protein domains to engineer new structures and functions. These domains are also well packed with few or no cavities. Moreover, a significant trend towards antiparallel helix configuration is observed in all these domains and could provide stabilizing interactions among macrodipoles. Finally, extensive networks of Cα-H···OC hydrogen bonds are detected in these domains. Despite their diverse evolutionary origins and activities, glycine-rich polyproline II helix assemblies share architectonic features which could help design novel proteins.

How hydrophobicity, side chains, and salt affect the dimensions of disordered proteins

Despite the generally accepted role of the hydrophobic effect as the driving force for folding, many intrinsically disordered proteins (IDPs), including those with hydrophobic content typical of foldable proteins, behave nearly as self-avoiding random walks (SARWs) under physiological conditions. Here, we tested how temperature and ionic conditions influence the dimensions of the N-terminal domain of pertactin (PNt), an IDP with an amino acid composition typical of folded proteins. While PNt contracts somewhat with temperature, it nevertheless remains expanded over 10-58°C, with a Flory exponent, ν, >0.50. Both low and high ionic strength also produce contraction in PNt, but this contraction is mitigated by reducing charge segregation. With 46% glycine and low hydrophobicity, the reduced form of snow flea anti-freeze protein (red-sfAFP) is unaffected by temperature and ionic strength and persists as a near-SARW, ν ~ 0.54, arguing that the thermal contraction of PNt is due to stronger interactions between hydrophobic side chains. Additionally, red-sfAFP is a proxy for the polypeptide backbone, which has been thought to collapse in water. Increasing the glycine segregation in red-sfAFP had minimal effect on ν. Water remained a good solvent even with 21 consecutive glycine residues (ν > 0.5), and red-sfAFP variants lacked stable backbone hydrogen bonds according to hydrogen exchange. Similarly, changing glycine segregation has little impact on ν in other glycine-rich proteins. These findings underscore the generality that many disordered states can be expanded and unstructured, and that the hydrophobic effect alone is insufficient to drive significant chain collapse for typical protein sequences.

The relevance of short peptides for an understanding of unfolded and intrinsically disordered proteins

Over the last thirty years the unfolded state of proteins has attracted considerable interest owing to the discovery of intrinsically disordered proteins which perform a plethora of functions despite resembling unfolded proteins to a significant extent. Research on both, unfolded and disordered proteins has revealed that their conformational properties can deviate locally from random coil behavior. In this context results from work on short oligopeptides suggest that individual amino acid residues sample the sterically allowed fraction of the Ramachandran plot to a different extent. Alanine has been found to exhibit a peculiarity in that it has a very high propensity for adopting polyproline II like conformations. This Perspectives article reviews work on short peptides aimed at exploring the Ramachandran distributions of amino acid residues in different contexts with experimental and computational means. Based on the thus provided overview the article discussed to what extent short peptides can serve as tools for exploring unfolded and disordered proteins and as benchmarks for the development of a molecular dynamics force field.

New insight into the mechanism by which antifreeze peptides regulate the physiological function of Streptococcus thermophilus subjected to freezing stress

INTRODUCTION

Antifreeze peptides regulate the physiological functions of frozen cells and even their apoptosis; however, the mechanisms by which antifreeze peptides regulate these processes remain unclear, although the interactions between cell membranes and ice are well known to be important in this process.

OBJECTIVES

Our study aims to investigate how antifreeze peptides regulate cell physiological functions during the freezing process.

METHODS

We investigated the cryoprotective effect of rsfAFP on the physiological functions of S. thermophilus under freezing stress by measuring cellular metabolism activity, intracellular enzyme activity, cell membrane characterization, and cell apoptosis. The mechanism by which rsfAFP impacts S. thermophilus physiological functions under freezing stress was investigated using multispectral techniques and cryo-TEM.

RESULTS

We show that a recombinant antifreeze peptide (rsfAFP) interacts with the extracellular capsular polysaccharides and peptidoglycan of Streptococcus thermophilus and ice to cover the outer layer of the membrane, forming a dense protective layer that regulates the molecular structure of extracellular ice crystals, which results in reduced extracellular membrane damage, depressed apoptosis and increased intracellular metabolic activity. This interaction mechanism was indicated by the fact that S. thermophilus better maintained its permeability barrier, membrane fluidity, membrane structural integrity, and cytoplasmic membrane potential during freezing stress with rsfAFP treatment.

CONCLUSION

These results provide new insights into the mechanism by which rsfAFP regulates frozen cellphysiological functionsand apoptosis under freezing stress.

Proline C-H Bonds as Loci for Proline Assembly via C-H/O Interactions

Proline residues within proteins lack a traditional hydrogen bond donor. However, the hydrogens of the proline ring are all sterically accessible, with polarized C-H bonds at Hα and Hδ that exhibit greater partial positive character and can be utilized as alternative sites for molecular recognition. C-H/O interactions, between proline C-H bonds and oxygen lone pairs, have been previously identified as modes of recognition within protein structures and for higher-order assembly of protein structures. In order to better understand intermolecular recognition of proline residues, a series of proline derivatives was synthesized, including 4R-hydroxyproline nitrobenzoate methyl ester, acylated on the proline nitrogen with bromoacetyl and glycolyl groups, and Boc-4S-(4-iodophenyl)hydroxyproline methyl amide. All three derivatives exhibited multiple close intermolecular C-H/O interactions in the crystallographic state, with H⋅⋅⋅O distances as close as 2.3 Å. These observed distances are well below the 2.72 Å sum of the van der Waals radii of H and O, and suggest that these interactions are particularly favorable. In order to generalize these results, we further analyzed the role of C-H/O interactions in all previously crystallized derivatives of these amino acids, and found that all 26 structures exhibited close intermolecular C-H/O interactions. Finally, we analyzed all proline residues in the Cambridge Structural Database of small-molecule crystal structures. We found that the majority of these structures exhibited intermolecular C-H/O interactions at proline C-H bonds, suggesting that C-H/O interactions are an inherent and important mode for recognition of and higher-order assembly at proline residues. Due to steric accessibility and multiple polarized C-H bonds, proline residues are uniquely positioned as sites for binding and recognition via C-H/O interactions.

Insight into polyproline II helical bundle stability in an antifreeze protein denatured state

The use of polyproline II (PPII) helices in protein design is currently hindered by limitations in our understanding of their conformational stability and folding. Recent studies of the snow flea antifreeze protein (sfAFP), a useful model system composed of six PPII helices, suggested that a low denatured state entropy contributes to folding thermodynamics. Here, circular dichroism spectroscopy revealed minor populations of PPII like conformers at low temperature. To get atomic level information on the conformational ensemble and entropy of the reduced, denatured state of sfAFP, we have analyzed its chemical shifts and {1H}-15N relaxation parameters by NMR spectroscopy at four experimental conditions. No significant populations of stable secondary structure were detected. The stiffening of certain N-terminal residues at neutral versus acidic pH and shifted pKa values leads us to suggest that favorable charge-charge interactions could bias the conformational ensemble to favor the formation the C1-C28 disulfide bond during nascent folding, although no evidence for preferred contacts between these positions was detected by paramagnetic relaxation enhancement under denaturing conditions. Despite a high content of flexible glycine residues, the mobility of the sfAFP denatured ensemble is similar for denatured α/β proteins both on fast ps/ns as well as slower μs/ms timescales. These results are in line with a conformational entropy in the denatured ensemble resembling that of typical proteins and suggest that new structures based on PPII helical bundles should be amenable to protein design.

Do polyproline II helix associations modulate biomolecular condensates?

Biomolecular condensates are microdroplets that form inside cells and serve to selectively concentrate proteins, RNAs and other molecules for a variety of physiological functions, but can contribute to cancer, neurodegenerative diseases and viral infections. The formation of these condensates is driven by weak, transient interactions between molecules. These weak associations can operate at the level of whole protein domains, elements of secondary structure or even moieties composed of just a few atoms. Different types of condensates do not generally combine to form larger microdroplets, suggesting that each uses a distinct class of attractive interactions. Here, we address whether polyproline II (PPII) helices mediate condensate formation. By combining with PPII-binding elements such as GYF, WW, profilin, SH3 or OCRE domains, PPII helices help form lipid rafts, nuclear speckles, P-body-like neuronal granules, enhancer complexes and other condensates. The number of PPII helical tracts or tandem PPII-binding domains can strongly influence condensate stability. Many PPII helices have a low content of proline residues, which hinders their identification. Recently, we characterized the NMR spectral properties of a Gly-rich, Pro-poor protein composed of six PPII helices. Based on those results, we predicted that many Gly-rich segments may form PPII helices and interact with PPII-binding domains. This prediction is being tested and could join the palette of verified interactions contributing to biomolecular condensate formation.

Glycine rich segments adopt polyproline II helices: Implications for biomolecular condensate formation

Many intrinsically disordered proteins contain Gly-rich regions which are generally assumed to be disordered. Such regions often form biomolecular condensates which play essential roles in organizing cellular processes. However, the bases of their formation and stability are still not completely understood. Based on NMR studies of the Gly-rich H. harveyi "snow flea" antifreeze protein, we recently proposed that Gly-rich sequences, such as the third "RGG" region of Fused in Sarcoma (FUS) protein, may adopt polyproline II helices whose association might stabilize condensates. Here, this hypothesis is tested with a polypeptide corresponding to the third RGG region of FUS. NMR spectroscopy and molecular dynamics simulations suggest that significant populations of polyproline II helix are present. These findings are corroborated in a model peptide Ac-RGGYGGRGGWGGRGGY-NH₂, where a peak characteristic of polyproline II helix is observed using CD spectroscopy. Its intensity suggests a polyproline II population of 40%. This result is supported by data from FTIR and NMR spectroscopies. In the latter, NOE correlations are observed between the Tyr and Arg, and Arg and Trp side chain hydrogens, confirming that side chains spaced three residues apart are close in space. Taken together, the data are consistent with a polyproline II helix, which is bent to optimize interactions between guanidinium and aromatic moieties, in equilibrium with a statistical coil ensemble. These results lend credence to the hypothesis that Gly-rich segments of disordered proteins may form polyproline II helices which help stabilize biomolecular condensates.

Therapeutic peptides for the treatment of cystic fibrosis: Challenges and perspectives

Cystic fibrosis (CF) is the most common amongst rare genetic diseases, affecting more than 70.000 people worldwide. CF is characterized by a dysfunctional chloride channel, termed cystic fibrosis conductance regulator (CFTR), which leads to the production of a thick and viscous mucus layer that clogs the lungs of CF patients and traps pathogens, leading to chronic infections and inflammation and, ultimately, lung damage. In recent years, the use of peptides for the treatment of respiratory diseases, including CF, has gained growing interest. Therapeutic peptides for CF include antimicrobial peptides, inhibitors of proteases, and modulators of ion channels, among others. Peptides display unique features that make them appealing candidates for clinical translation, like specificity of action, high efficacy, and low toxicity. Nevertheless, the intrinsic properties of peptides, together with the need of delivering these compounds locally, e.g. by inhalation, raise a number of concerns in the development of peptide therapeutics for CF lung disease. In this review, we discuss the challenges related to the use of peptides for the treatment of CF lung disease through inhalation, which include retention within mucus, proteolysis, immunogenicity and aggregation. Strategies for overcoming major shortcomings of peptide therapeutics will be presented, together with recent developments in peptide design and optimization, including computational analysis and high-throughput screening.

Piecewise All-Atom SMD Simulations Reveal Key Secondary Structures in Luciferase Unfolding Pathway

Although the folding of single-domain proteins is well characterized theoretically and experimentally, the folding of large multidomain proteins is less well known. Firefly luciferase, a 550 residue three-domain protein, has been commonly used as a substrate to study chaperone reactions and as a model system for the study of folding of long polypeptide chains, including related phenomena such as cotranslational folding. Despite being characterized by various experimental techniques, the atomic-level contributions of various secondary structures of luciferase to its fold's mechanical stability remain unknown. Here, we developed a piecewise approach for all-atom steered molecular dynamics simulations to examine specific secondary structures that resist mechanical unfolding while minimizing the amount of computational resources required by the large water box of standard all-atom steered molecular dynamics simulations. We validated the robustness of this approach with a small NI3C protein and used our approach to elucidate the specific secondary structures that provide the largest contributions to luciferase mechanostability. In doing so, we show that piecewise all-atom steered molecular dynamics simulations can provide novel atomic resolution details regarding mechanostability and can serve as a platform for novel mutagenesis studies as well as a point for comparison with high-resolution force spectroscopy experiments.

Glycine in Water Favors the Polyproline II State

Conformational preferences of amino acid residues in water are determined by the backbone and side-chain properties. Alanine is known for its high polyproline II (pPII) propensity. The question of relative contributions of the backbone and side chain to the conformational preferences of alanine and other amino acid residues in water is not fully resolved. Because glycine lacks a heavy-atom side chain, glycine-based peptides can be used to examine to which extent the backbone properties affect the conformational space. Here, we use published spectroscopic data for the central glycine residue of cationic triglycine in water to demonstrate that its conformational space is dominated by the pPII state. We assess three commonly used molecular dynamics (MD) force fields with respect to their ability to capture the conformational preferences of the central glycine residue in triglycine. We show that pPII is the mesostate that enables the functional backbone groups of the central residue to form the most hydrogen bonds with water. Our results indicate that the pPII propensity of the central glycine in GGG is comparable to that of alanine in GAG, implying that the water-backbone hydrogen bonding is responsible for the high pPII content of these residues.

Water as a Good Solvent for Unfolded Proteins: Folding and Collapse are Fundamentally Different

The argument that the hydrophobic effect is the primary effect driving the folding of globular proteins is nearly universally accepted (including by the authors). But does this view also imply that water is a "poor" solvent for the unfolded states of these same proteins? Here we argue that the answer is "no," that is, folding to a well-packed, extensively hydrogen-bonded native structure differs fundamentally from the nonspecific chain collapse that defines a poor solvent. Thus, the observation that a protein folds in water does not necessitate that water is a poor solvent for its unfolded state. Indeed, chain-solvent interactions that are marginally more favorable than nonspecific intrachain interactions are beneficial to protein function because they destabilize deleterious misfolded conformations and inter-chain interactions.

Predicting Conformational Properties of Intrinsically Disordered Proteins from Sequence

Intrinsically disordered proteins (IDPs) can adopt a range of conformations from globules to swollen coils. This large range of conformational preferences for different IDPs raises the question of how conformational preferences are encoded by sequence. Global compositional features of a sequence such as the fraction of charged residues and the net charge per residue engender certain conformational biases. However, more specific sequence features such as the patterning of oppositely charged residues, expansion driving residues, or residues that can undergo posttranslational modifications can also influence the conformational ensembles of an IDP. Here, we outline how to calculate important global compositional features and patterning metrics that can be used to classify IDPs into different conformational classes and predict relative changes in conformation for sequences with the same amino acid composition. Although increased effort has been devoted to determining conformational properties of IDPs in recent years, quantitative predictions of conformation directly from sequence remain difficult and often inaccurate. Thus, if quantitative predictions of conformational properties are desired, then sequence-specific simulations must be performed.

Sequence Effects on Size, Shape, and Structural Heterogeneity in Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs) lack well-defined three-dimensional structures, thus challenging the archetypal notion of structure-function relationships. Determining the ensemble of conformations that IDPs explore under physiological conditions is the first step toward understanding their diverse cellular functions. Here, we quantitatively characterize the structural features of IDPs as a function of sequence and length using coarse-grained simulations. For diverse IDP sequences, with the number of residues ( NT) ranging from 20 to 441, our simulations not only reproduce the radii of gyration ( Rg) obtained from experiments, but also predict the full scattering intensity profiles in excellent agreement with small-angle X-ray scattering experiments. The Rg values are well-described by the standard Flory scaling law, Rg = Rg0 NTν, with ν ≈ 0.588, making it tempting to assert that IDPs behave as polymers in a good solvent. However, clustering analysis reveals that the menagerie of structures explored by IDPs is diverse, with the extent of heterogeneity being highly sequence-dependent, even though ensemble-averaged properties, such as the dependence of Rg on chain length, may suggest synthetic polymer-like behavior in a good solvent. For example, we show that for the highly charged Prothymosin-α, a substantial fraction of conformations is highly compact. Even if the sequence compositions are similar, as is the case for α-Synuclein and a truncated construct from the Tau protein, there are substantial differences in the conformational heterogeneity. Taken together, these observations imply that metrics based on net charge or related quantities alone cannot be used to anticipate the phases of IDPs, either in isolation or in complex with partner IDPs or RNA. Our work sets the stage for probing the interactions of IDPs with each other, with folded protein domains, or with partner RNAs, which are critical for describing the structures of stress granules and biomolecular condensates with important cellular functions.

Commonly used FRET fluorophores promote collapse of an otherwise disordered protein

The dimensions that unfolded proteins, including intrinsically disordered proteins (IDPs), adopt in the absence of denaturant remain controversial. We developed an analysis procedure for small-angle X-ray scattering (SAXS) profiles and used it to demonstrate that even relatively hydrophobic IDPs remain nearly as expanded in water as they are in high denaturant concentrations. In contrast, as demonstrated here, most fluorescence resonance energy transfer (FRET) measurements have indicated that relatively hydrophobic IDPs contract significantly in the absence of denaturant. We use two independent approaches to further explore this controversy. First, using SAXS we show that fluorophores employed in FRET can contribute to the observed discrepancy. Specifically, we find that addition of Alexa-488 to a normally expanded IDP causes contraction by an additional 15%, a value in reasonable accord with the contraction reported in FRET-based studies. Second, using our simulations and analysis procedure to accurately extract both the radius of gyration (R_g) and end-to-end distance (R_ee) from SAXS profiles, we tested the recent suggestion that FRET and SAXS results can be reconciled if the R_g and R_ee are "uncoupled" (i.e., no longer simply proportional), in contrast to the case for random walk homopolymers. We find, however, that even for unfolded proteins, these two measures of unfolded state dimensions remain proportional. Together, these results suggest that improved analysis procedures and a correction for significant, fluorophore-driven interactions are sufficient to reconcile prior SAXS and FRET studies, thus providing a unified picture of the nature of unfolded polypeptide chains in the absence of denaturant.

The Singular NMR Fingerprint of a Polyproline II Helical Bundle

Polyproline II (PPII) helices play vital roles in biochemical recognition events and structures like collagen and form part of the conformational landscapes of intrinsically disordered proteins (IDPs). Nevertheless, this structure is generally hard to detect and quantify. Here, we report the first thorough NMR characterization of a PPII helical bundle protein, the Hypogastrura harveyi "snow flea" antifreeze protein (sfAFP). J-couplings and nuclear Overhauser enhancement spectroscopy confirm a natively folded structure consisting of six PPII helices. NMR spectral analyses reveal quite distinct Hα2 versus Hα3 chemical shifts for 28 Gly residues as well as ¹³Cα, ¹⁵N, and ¹HN conformational chemical shifts (Δδ) unique to PPII helical bundles. The ¹⁵N Δδ and ¹HN Δδ values and small negative ¹HN temperature coefficients evince hydrogen-bond formation. ¹H-¹⁵N relaxation measurements reveal that the backbone structure is generally highly rigid on ps-ns time scales. NMR relaxation parameters and biophysical characterization reveal that sfAFP is chiefly a dimer. For it, a structural model featuring the packing of long, flat hydrophobic faces at the dimer interface is advanced. The conformational stability, measured by amide H/D exchange to be 6.24 ± 0.2 kcal·mol^-1, is elevated. These are extraordinary findings considering the great entropic cost of fixing Gly residues and, together with the remarkable upfield chemical shifts of 28 Gly ¹Hα, evidence significant stabilizing contributions from CαHα ||| O═C hydrogen bonds. These stabilizing interactions are corroborated by density functional theory calculations and natural bonding orbital analysis. The singular conformational chemical shifts, J-couplings, high hNOE ratios, small negative temperature coefficients, and slowed H/D exchange constitute a unique set of fingerprints to identify PPII helical bundles, which may be formed by hundreds of Gly-rich motifs detected in sequence databases. These results should aid the quantification of PPII helices in IDPs, the development of improved antifreeze proteins, and the incorporation of PPII helices into novel designed proteins.

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.

Response to Comment on “Innovative scattering analysis shows that hydrophobic disordered proteins are expanded in water”

Best et al. claim that we provide no convincing basis to assert that a discrepancy remains between FRET and SAXS results on the dimensions of disordered proteins under physiological conditions. We maintain that a clear discrepancy is apparent in our and other recent publications, including results shown in the Best et al. comment. A plausible origin is fluorophore interactions in FRET experiments.

Calcium-Mediated Allostery of the EGF Fold

The epidermal growth factor (EGF)-like domain is one of the most abundant disulfide-containing domains in nature and is involved in many cellular processes critical to life. Although many EGF-like domains participate in calcium-dependent functions by responding to the local calcium concentration, little is known about how this responsiveness is programmed at the molecular level. Here, we reveal the structural and environmental determinants underpinning the folding of a synthetic analogue of the EGF-A domain (from the low-density lipoprotein receptor). We show that calcium sensitivity is enabled by an allosteric folding pathway, in which calcium binding is connected to the peptide core through local inter-residue interactions. In the absence of calcium, the fold favors disorder because the inherently weak core is insufficient to stabilize the active form, resulting in substantial loss in activity of 2 orders of magnitude. The EGF-A fold, which can freely transition between active and disordered states, is volatile, and we found it to be intolerant of mutations, unlike other disulfide-rich peptides that have been used as stabilizing frameworks. This volatility is beneficial for modularity/plasticity and appears to have evolved for such a purpose, allowing cellular pathways to sense and respond to environmental cues.

Collapse Transitions of Proteins and the Interplay Among Backbone, Sidechain, and Solvent Interactions

Proteins can collapse into compact globules or form expanded, solvent-accessible, coil-like conformations. Additionally, they can fold into well-defined three-dimensional structures or remain partially or entirely disordered. Recent discoveries have shown that the tendency for proteins to collapse or remain expanded is not intrinsically coupled to their ability to fold. These observations suggest that proteins do not have to form compact globules in aqueous solutions. They can be intrinsically disordered, collapsed, or expanded, and even form well-folded, elongated structures. This ability to decouple collapse from folding is determined by the sequence details of proteins. In this review, we highlight insights gleaned from studies over the past decade. Using a polymer physics framework, we explain how the interplay among sidechains, backbone units, and solvent determines the driving forces for collapsed versus expanded states in aqueous solvents.

Peptide Solubility Limits: Backbone and Side-Chain Interactions

We calculate the solubility limit of pentapeptides in water by simulating the phase separation in an oversaturated aqueous solution. The solubility limit order followed by our model peptides (GGRGG > GGDGG > GGGGG > GGVGG > GGQGG > GGNGG > GGFGG) is found to be the same as that reported for amino acid monomers from experiment (R > D > G > V > Q > N > F). Investigation of dynamical properties of peptides shows that the higher the solubility of a peptide is, the lower the time spent by the peptide in the aggregated cluster is. We also demonstrate that fluctuations in conformation and hydration number of peptide in monomeric form are correlated with the solubility of the peptide. We considered energetic mechanisms and dynamical properties of interbackbone CO-CO and CO···HN interactions. Our results confirm that CO-CO interactions more than the interbackbone H-bonds are important in peptide self-assembly and association. Further, we find that the stability of H-bonded peptide pairs arises mainly from coexisting CO-CO and CO···HN interactions.