Effects of denaturants and osmolytes on proteins are accurately predicted by the molecular transfer model

Mechanical Unfolding of Network Nodes Drives the Stress Response of Protein-Based Materials

Biomaterials synthesized from cross-linked folded proteins have untapped potential for biocompatible, resilient, and responsive implementations, but face challenges due to costly molecular refinement and limited understanding of their mechanical response. Under a stress vector, these materials combine the gel-like response of cross-linked networks with the mechanical unfolding and extension of proteins from well-defined 3D structures to unstructured polypeptides. Yet the nanoscale dynamics governing their viscoelastic response remains poorly understood. This lack of understanding is further exacerbated by the fact that the mechanical stability of protein domains depends not only on their structure, but also on the direction of the force vector. To this end, here we propose a coarse-grained network model based on the physical characteristics of polyproteins and combine it with the mechanical unfolding response of protein domains, obtained from single molecule measurements and steered molecular dynamics simulations, to explain the macroscopic response of protein-based materials to a stress vector. We find that domains are about 10-fold more stable when force is applied along their end-to-end coordinate than along the other tethering geometries that are possible inside the biomaterial. As such, the macroscopic response of protein-based materials is mainly driven by the unfolding of the node-domains and rearrangement of these nodes inside the material. The predictions from our models are then confirmed experimentally using force-clamp rheometry. This model is a critical step toward developing protein-based materials with predictable response and that can enable applications for shape memory and energy storage and dissipation.

The Action of Chemical Denaturants: From Globular to Intrinsically Disordered Proteins

Proteins perform their many functions by adopting either a minimal number of strictly similar conformations, the native state, or a vast ensemble of highly flexible conformations. In both cases, their structural features are highly influenced by the chemical environment. Even though a plethora of experimental studies have demonstrated the impact of chemical denaturants on protein structure, the molecular mechanism underlying their action is still debated. In the present review, after a brief recapitulation of the main experimental data on protein denaturants, we survey both classical and more recent interpretations of the molecular basis of their action. In particular, we highlight the differences and similarities of the impact that denaturants have on different structural classes of proteins, i.e., globular, intrinsically disordered (IDP), and amyloid-like assemblies. Particular attention has been given to the IDPs, as recent studies are unraveling their fundamental importance in many physiological processes. The role that computation techniques are expected to play in the near future is illustrated.

Residue-Dependent Transition Temperatures and Denaturant Midpoints in the Folding of a Multidomain Protein

As a consequence of the finite size of globular proteins, it is expected that there should be dispersions in the global melting temperature (T_m) and the denaturation midpoint (C_m). Thermodynamic considerations dictate that the dispersions, ΔT_m in T_m, and ΔC_m in C_m, should decrease with N, the number of residues in the protein. We performed coarse-grained simulations of the self-organized polymer (SOP) model of the multidomain protein adenylate kinase (ADK) with N = 214 in order to calculate thermal and denaturation unfolding titration curves. The results show that ΔT_m/T_m and ΔC_m/C_m are nonzero and follow the previously established ( Phys. Rev. Lett. 2004, 93, 268107) thermodynamic 1/N scaling for proteins accurately. For ADK, the dispersions are small (≈0.004), which implies that the melting temperature is more or less unique, which is unlike in BBL (N = 40) where ΔT_m/T_m ≈ 0.03.

Salt-Induced Transitions in the Conformational Ensembles of Intrinsically Disordered Proteins

Salts modulate the behavior of intrinsically disordered proteins (IDPs) and influence the formation of membraneless organelles through liquid-liquid phase separation (LLPS). In low ionic strength solutions, IDP conformations are perturbed by the screening of electrostatic interactions, independent of the salt identity. In this regime, insight into the IDP behavior can be obtained using the theory for salt-induced transitions in charged polymers. However, salt-specific interactions with the charged and uncharged residues, known as the Hofmeister effect, influence IDP behavior in high ionic strength solutions. There is a lack of reliable theoretical models in high salt concentration regimes to predict the salt effect on IDPs. We propose a simulation methodology using a coarse-grained IDP model and experimentally measured water to salt solution transfer free energies of various chemical groups that allowed us to study the salt-specific transitions induced in the IDPs conformational ensemble. We probed the effect of three different monovalent salts on five IDPs belonging to various polymer classes based on charged residue content. We demonstrate that all of the IDPs of different polymer classes behave as self-avoiding walks (SAWs) at physiological salt concentration. In high salt concentrations, the transitions observed in the IDP conformational ensembles are dependent on the salt used and the IDP sequence and composition. Changing the anion with the cation fixed can result in the IDP transition from a SAW-like behavior to a collapsed globule. An important implication of these results is that a suitable salt can be identified to induce condensation of an IDP through LLPS.

A Transfer Free Energy Based Implicit Solvent Model for Protein Simulations in Solvent Mixtures: Urea-Induced Denaturation as a Case Study

We developed a method for implicit solvent molecular dynamics simulations of proteins in solvent mixtures (model with implicit solvation thermodynamics, MIST). The MIST method introduces experimental group transfer free energies to the generalized Born formulation for generating molecular trajectories without the need for developing rigorous explicit-solvent force fields for multicomponent solutions. As a test case, we studied the urea-induced denaturation of the Trp-cage miniprotein in water. We demonstrate that our method allows efficient exploration of the conformational space of the protein in only a few hundreds of nanoseconds of all-atom unbiased simulations. Furthermore, selective implementation of the transfer free energies of specific peptide groups, backbone, and side chains enables us to decouple their specific energetic contributions to the conformational changes of the protein. The approach herein developed can readily be extended to the investigation of complex matrices as well as to the characterization of protein aggregation. The MIST method is implemented in Plumed (ver. 2.8) as a separate module called SASA.

Cooperativity and Folding Kinetics in a Multidomain Protein with Interwoven Chain Topology

Although a large percentage of eukaryotic proteomes consist of proteins with multiple domains, not much is known about their assembly mechanism, especially those with intricate native state architectures. Some have a complex topology in which the structural elements along the sequence are interwoven in such a manner that the domains cannot be separated by cutting at any location along the sequence. Such proteins are multiply connected multidomain proteins (MMPs) with the three-domain (NMP, LID, and CORE) phosphotransferase enzyme adenylate kinase (ADK) being an example. We devised a coarse-grained model to simulate ADK folding initiated by changing either the temperature or guanidinium chloride (GdmCl) concentration. The simulations reproduce the experimentally measured melting temperatures (associated with two equilibrium transitions), FRET efficiency as a function of GdmCl concentration, and the folding times quantitatively. Although the NMP domain orders independently, cooperative interactions between the LID and the CORE domains are required for complete assembly of the enzyme. Kinetic simulations show that, on the collapse time scale, multiple interconnected metastable states are populated, attesting to the folding heterogeneity. The network of kinetically connected states reveals that the CORE domain folds only after the NMP and LID domains, reflecting the interwoven nature of the chain topology.

Forming a Double-Helix Phase of Single Polymer Chains by the Cooperation between Local Structure and Nonlocal Attraction

Double-helix structures, such as DNA, are formed in nature to realize many unique functions. Inspired by this, researchers are pursuing strategies to design such structures from polymers. A key question is whether the double helix can be formed from the self-folding of a single polymer chain without specific interactions. Here, using Langevin dynamics simulation and theoretical analysis, we find that a stable double-helix phase can be achieved by the self-folding of single semiflexible polymers as a result of the cooperation between local structure and nonlocal attraction. The critical temperature of double-helix formation approximately follows T^{cri}∼ln(k_{θ}) and T^{cri}∼ln(k_{τ}), where k_{θ} and k_{τ} are the polymer bending and torsion stiffness, respectively. Furthermore, the double helix can exhibit major and minor grooves due to symmetric break for better packing. Our results provide a novel guide to the experimental design of the double helix.

Macromolecular Crowding Is More than Hard-Core Repulsions

Cells are crowded, but proteins are almost always studied in dilute aqueous buffer. We review the experimental evidence that crowding affects the equilibrium thermodynamics of protein stability and protein association and discuss the theories employed to explain these observations. In doing so, we highlight differences between synthetic polymers and biologically relevant crowders. Theories based on hard-core interactions predict only crowding-induced entropic stabilization. However, experiment-based efforts conducted under physiologically relevant conditions show that crowding can destabilize proteins and their complexes. Furthermore, quantification of the temperature dependence of crowding effects produced by both large and small cosolutes, including osmolytes, sugars, synthetic polymers, and proteins, reveals enthalpic effects that stabilize or destabilize proteins. Crowding-induced destabilization and the enthalpic component point to the role of chemical interactions between and among the macromolecules, cosolutes, and water. We conclude with suggestions for future studies. Expected final online publication date for the Annual Review of Biophysics, Volume 51 is May 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Environmental factors modulating protein conformations and their role in protein aggregation diseases

The adverse physiological conditions have been long known to impact protein synthesis, folding and functionality. Major physiological factors such as the effect of pH, temperature, salt and pressure are extensively studied for their impact on protein structure and homeostasis. However, in the current scenario, the environmental risk factors (pollutants) have gained impetus in research because of their increasing concentrations in the environment and strong epidemiologic link with protein aggregation disorders. Here, we review the physiological and environmental risk factors for their impact on protein conformational changes, misfolding, aggregation, and associated pathological conditions, especially environmental risk factors associated pathologies.

Temperature and Guanidine Hydrochloride Effects on the Folding Thermodynamics of WW Domain and Variants

We used simulations based on an all-atom Go model to calculate the folding temperatures (T_fs) and free energies (ΔGs) of two variants of the WW domain, which is a small all-β-sheet protein. The results, without adjusting any parameter, are in good agreement with experiments, thus validating the simulations. We then used the molecular transfer model to predict the changes in the ΔGs and T_fs as the guanidine hydrochloride concentration is varied. The predictions can be readily tested in experiments.

Integrating single-molecule spectroscopy and simulations for the study of intrinsically disordered proteins

Over the last two decades, intrinsically disordered proteins and protein regions (IDRs) have emerged from a niche corner of biophysics to be recognized as essential drivers of cellular function. Various techniques have provided fundamental insight into the function and dysfunction of IDRs. Among these techniques, single-molecule fluorescence spectroscopy and molecular simulations have played a major role in shaping our modern understanding of the sequence-encoded conformational behavior of disordered proteins. While both techniques are frequently used in isolation, when combined they offer synergistic and complementary information that can help uncover complex molecular details. Here we offer an overview of single-molecule fluorescence spectroscopy and molecular simulations in the context of studying disordered proteins. We discuss the various means in which simulations and single-molecule spectroscopy can be integrated, and consider a number of studies in which this integration has uncovered biological and biophysical mechanisms.

Atomic view of cosolute-induced protein denaturation probed by NMR solvent paramagnetic relaxation enhancement

The cosolvent effect arises from the interaction of cosolute molecules with a protein and alters the equilibrium between native and unfolded states. Denaturants shift the equilibrium toward the latter, while osmolytes stabilize the former. The molecular mechanism whereby cosolutes perturb protein stability is still the subject of considerable debate. Probing the molecular details of the cosolvent effect is experimentally challenging as the interactions are very weak and transient, rendering them invisible to most conventional biophysical techniques. Here, we probe cosolute-protein interactions by means of NMR solvent paramagnetic relaxation enhancement together with a formalism we recently developed to quantitatively describe, at atomic resolution, the energetics and dynamics of cosolute-protein interactions in terms of a concentration normalized equilibrium average of the interspin distance, [Formula: see text], and an effective correlation time, τ_c The system studied is the metastable drkN SH3 domain, which exists in dynamic equilibrium between native and unfolded states, thereby permitting us to probe the interactions of cosolutes with both states simultaneously under the same conditions. Two paramagnetic cosolute denaturants were investigated, one neutral and the other negatively charged, differing in the presence of a carboxyamide group versus a carboxylate. Our results demonstrate that attractive cosolute-protein backbone interactions occur largely in the unfolded state and some loop regions in the native state, electrostatic interactions reduce the [Formula: see text] values, and temperature predominantly impacts interactions with the unfolded state. Thus, destabilization of the native state in this instance arises predominantly as a consequence of interactions of the cosolutes with the unfolded state.

Electrostatic Interactions Explain the Higher Binding Affinity of the CR3022 Antibody for SARS-CoV-2 than the 4A8 Antibody

A structural understanding of the mechanism by which antibodies bind SARS-CoV-2 at the atomic level is highly desirable as it can tell the development of more effective antibodies to treat Covid-19. Here, we use steered molecular dynamics (SMD) and coarse-grained simulations to estimate the binding affinity of the monoclonal antibodies CR3022 and 4A8 to the SARS-CoV-2 receptor-binding domain (RBD) and SARS-CoV-2 N-terminal domain (NTD). Consistent with experiments, our SMD and coarse-grained simulations both indicate that CR3022 has a higher affinity for SARS-CoV-2 RBD than 4A8 for the NTD, and the coarse-grained simulations indicate the former binds three times stronger to its respective epitope. This finding shows that CR3022 is a candidate for Covid-19 therapy and is likely a better choice than 4A8. Energetic decomposition of the interaction energies between these two complexes reveals that electrostatic interactions explain the difference in the observed binding affinity between the two complexes. This result could lead to a new approach for developing anti-Covid-19 antibodies in which good candidates must contain charged amino acids in the area of contact with the virus.

General Counteraction Exerted by Sugars against Denaturants

The conformational stability of globular proteins is strongly influenced by the addition to water of different co-solutes. Some of the latter destabilize the native state, while others stabilize it. It is emerging that stabilizing agents are able to counteract the action of destabilizing agents. We have already provided experimental evidence that this counteraction is a general phenomenon and offered a rationalization. In the present work, we show that four different sugars, namely fructose, glucose, sucrose, and trehalose, counteract the effect of urea, tetramethylurea, sodium perchlorate, guanidinium chloride, and guanidinium thiocyanate despite the chemical and structural differences of those destabilizing agents. The rationalization we provide is as follows: (a) the solvent-excluded volume effect, a purely entropic effect, stabilizes the native state, whose solvent-accessible surface area is smaller than the one of denatured conformations; (b) the magnitude of the solvent-excluded volume effect increases markedly in ternary solutions because the experimental density of such solutions is larger than that of pure water.

Small-Angle X-ray Scattering Signatures of Conformational Heterogeneity and Homogeneity of Disordered Protein Ensembles

An accurate account of disordered protein conformations is of central importance to deciphering the physicochemical basis of biological functions of intrinsically disordered proteins and the folding-unfolding energetics of globular proteins. Physically, disordered ensembles of nonhomopolymeric polypeptides are expected to be heterogeneous, i.e., they should differ from those homogeneous ensembles of homopolymers that harbor an essentially unique relationship between average values of end-to-end distance R_EE and radius of gyration R_g. It was posited recently, however, that small-angle X-ray scattering (SAXS) data on conformational dimensions of disordered proteins can be rationalized almost exclusively by homopolymer ensembles. Assessing this perspective, chain-model simulations are used to evaluate the discriminatory power of SAXS-determined molecular form factors (MFFs) with regard to homogeneous versus heterogeneous ensembles. The general approach adopted here is not bound by any assumption about ensemble encodability, in that the postulated heterogeneous ensembles we evaluated are not restricted to those entailed by simple interaction schemes. Our analysis of MFFs for certain heterogeneous ensembles with more narrowly distributed R_EE and R_g indicates that while they deviate from MFFs of homogeneous ensembles, the differences can be rather small. Remarkably, some heterogeneous ensembles with asphericity and R_EE drastically different from those of homogeneous ensembles can nonetheless exhibit practically identical MFFs, demonstrating that SAXS MFFs do not afford unique characterizations of basic properties of conformational ensembles in general. In other words, the ensemble to MFF mapping is practically many-to-one and likely nonsmooth. Heteropolymeric variations of the R_EE-R_g relationship were further showcased using an analytical perturbation theory developed here for flexible heteropolymers. Ramifications of our findings for interpretation of experimental data are discussed.

Uncovering Differences in Hydration Free Energies and Structures for Model Compound Mimics of Charged Side Chains of Amino Acids

Free energies of hydration are of fundamental interest for modeling and understanding conformational and phase equilibria of macromolecular solutes in aqueous phases. Of particular relevance to systems such as intrinsically disordered proteins are the free energies of hydration and hydration structures of model compounds that mimic charged side chains of Arg, Lys, Asp, and Glu. Here, we deploy a Thermodynamic Cycle-based Proton Dissociation (TCPD) approach in conjunction with data from direct measurements to obtain estimates for the free energies of hydration for model compounds that mimic the side chains of Arg+, Lys+, Asp-, and Glu-. Irrespective of the choice made for the hydration free energy of the proton, the TCPD approach reveals clear trends regarding the free energies of hydration for Arg+, Lys+, Asp-, and Glu-. These trends include asymmetries between the hydration free energies of acidic (Asp- and Glu-) and basic (Arg+ and Lys+) residues. Further, the TCPD analysis, which relies on a combination of experimental data, shows that the free energy of hydration of Arg+ is less favorable than that of Lys+. We sought a physical explanation for the TCPD-derived trends in free energies of hydration. To this end, we performed temperature-dependent calculations of free energies of hydration and analyzed hydration structures from simulations that use the polarizable Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) force field and water model. At 298 K, the AMOEBA model generates estimates of free energies of hydration that are consistent with TCPD values with a free energy of hydration for the proton of ca. -259 kcal/mol. Analysis of temperature-dependent simulations leads to a structural explanation for the observed differences in free energies of hydration of ionizable residues and reveals that the heat capacity of hydration is positive for Arg+ and Lys+ and negative for Asp- and Glu-.

Modern Kinetics and Mechanism of Protein Folding: A Retrospective

Modern experimental kinetics of protein folding began in the early 1990s with the introduction of nanosecond laser pulses to trigger the folding reaction, providing an almost 10⁶-fold improvement in time resolution over the stopped-flow method being employed at the time. These experiments marked the beginning of the "fast-folding" subfield that enabled investigation of the kinetics of formation of secondary structural elements and disordered loops for the first time, as well as the fastest folding proteins. When I started to work on this subject, a fast folding protein was one that folded in milliseconds. There were, moreover, no analytical theoretical models and no atomistic or coarse-grained molecular dynamics simulations to describe the mechanism. Two of the most important discoveries from my lab since then are a protein that folds in hundreds of nanoseconds, as determined from nanosecond laser temperature experiments, and the discovery that the theoretically predicted barrier crossing time is about the same for proteins that differ in folding rates by 10⁴-fold, as determined from single molecule fluorescence measurements. We also developed what has been called the "Hückel model" of protein folding, which quantitatively explains a wide range of equilibrium and kinetic measurements. This retrospective traces the history of contributions to the "fast folding" subfield from my lab until about 3 years ago, when I left protein folding to spend the rest of my research career trying to discover an inexpensive drug for treating sickle cell disease.

Molecular Transfer Model for pH Effects on Intrinsically Disordered Proteins: Theory and Applications

We present a theoretical method to study how changes in pH shape the heterogeneous conformational ensemble explored by intrinsically disordered proteins (IDPs). The theory is developed in the context of coarse-grained models, which enable a fast, accurate, and extensive exploration of conformational space at a given protonation state. In order to account for pH effects, we generalize the molecular transfer model (MTM), in which conformations are re-weighted using the transfer free energy, which is the free energy necessary for bringing to equilibrium in a new environment a "frozen" conformation of the system. Using the semi-grand ensemble, we derive an exact expression of the transfer free energy, which amounts to the appropriate summation over all the protonation states. Because the exact result is computationally too demanding to be useful for large polyelectrolytes or IDPs, we introduce a mean-field (MF) approximation of the transfer free energy. Using a lattice model, we compare the exact and MF results for the transfer free energy and a variety of observables associated with the model IDP. We find that the precise location of the charged groups (the sequence), and not merely the net charge, determines the structural properties. We demonstrate that some of the limitations previously noted for MF theory in the context of globular proteins are mitigated when disordered polymers are studied. The excellent agreement between the exact and MF results poises us to use the method presented here as a computational tool to study the properties of IDPs and other biological systems as a function of pH.

Revealing the Hidden Sensitivity of Intrinsically Disordered Proteins to their Chemical Environment

Intrinsically disordered protein-regions (IDRs) make up roughly 30% of the human proteome and are central to a wide range of biological processes. Given a lack of persistent tertiary structure, all residues in IDRs are, to some extent, solvent exposed. This extensive surface area, coupled with the absence of strong intramolecular contacts, makes IDRs inherently sensitive to their chemical environment. We report a combined experimental, computational, and analytical framework for high-throughput characterization of IDR sensitivity. Our framework reveals that IDRs can expand or compact in response to changes in their solution environment. Importantly, the direction and magnitude of conformational change depend on both protein sequence and cosolute identity. For example, some solutes such as short polyethylene glycol chains exert an expanding effect on some IDRs and a compacting effect on others. Despite this complex behavior, we can rationally interpret IDR responsiveness to solution composition changes using relatively simple polymer models. Our results imply that solution-responsive IDRs are ubiquitous and can provide an additional layer of regulation to biological systems.

Frustration and inhomogeneous environments in relaxation of open chains with Ising-type interactions

Frustration can contribute to very slow relaxation times in large open chains, as in spin glasses and in biopolymers. However, frustration may not be sufficient to produce broken ergodicity in finite systems. Here we employ a system-plus-reservoir approach to investigate how strongly inhomogeneous environments and frustration compete in the relaxation of finite open chains. We find a sufficient condition for our inhomogeneous environments to break ergodicity. We use the microscopic model to derive a Markovian quantum master equation for a generic chain with ultrastrong intrachain couplings. We show that this microscopic model avoids a spurious broken ergodicity we find in the phenomenological model. We work out an explicit example of broken ergodicity due to the inhomogeneous environment of an unfrustrated spin chain as far as simulating a recent experiment on protein denaturation (where environment inhomogeneity is especially relevant). We finally show that an inhomogeneous environment can mitigate the effects of frustration-induced degeneracies.

Characterizing osmolyte chemical class hierarchies and functional group requirements for thermal stabilization of proteins

Osmolytes are naturally occurring organic compounds that protect cellular proteins and other macromolecules against various forms of stress including temperature extremes. While biological studies have correlated the accumulation of certain classes of osmolytes with specific forms of stress, including thermal stress, it remains unclear whether or not these observations reflect an intrinsic chemical class hierarchy amongst the osmolytes with respect to effects on protein stability. In addition, very little is known in regards to the molecular elements of the osmolytes themselves that are essential for their functions. In this study, we use differential scanning fluorimetry to quantify the thermal stabilizing effects of members from each of the three main classes of protecting osmolytes on two model protein systems, C-reactive protein and tumor necrosis factor alpha. Our data reveals the absence of a strict chemical class hierarchy amongst the osmolytes with respect to protein thermal stabilization, and indicates differential responses of these proteins to certain osmolytes. In the second part of this investigation we dissected the molecular elements of amino acid osmolytes required for thermal stabilization of myoglobin and C-reactive protein. We show that the complete amino acid zwitterion is required for thermal stabilization of myoglobin, whereas removal of the osmolyte amino group does not diminish stabilizing effects on C-reactive protein. These disparate responses of proteins to osmolytes and other small molecules are consistent with previous observations that osmolyte effects on protein stability are protein-specific. Moreover, the data reported in this study support the view that osmolyte effects cannot be fully explained by considering only the solvent accessibility of the polypeptide backbone in the native and denatured states, and corroborate the need for more complex models that take into account the entire protein fabric.

Electrostatic Interactions Govern Extreme Nascent Protein Ejection Times from Ribosomes and Can Delay Ribosome Recycling

The ejection of nascent proteins out of the ribosome exit tunnel, after their covalent bond to transfer-RNA has been broken, has not been experimentally studied due to challenges in sample preparation. Here, we investigate this process using a combination of multiscale modeling, ribosome profiling, and gene ontology analyses. Simulating the ejection of a representative set of 122 E. coli proteins we find a greater than 1000-fold variation in ejection times. Nascent proteins enriched in negatively charged residues near their C-terminus eject the fastest, while nascent chains enriched in positively charged residues tend to eject much more slowly. More work is required to pull slowly ejecting proteins out of the exit tunnel than quickly ejecting proteins, according to all-atom simulations. An energetic decomposition reveals, for slowly ejecting proteins, that this is due to the strong attractive electrostatic interactions between the nascent chain and the negatively charged ribosomal-RNA lining the exit tunnel, and for quickly ejecting proteins, it is due to their repulsive electrostatic interactions with the exit tunnel. Ribosome profiling data from E. coli reveals that the presence of slowly ejecting sequences correlates with ribosomes spending more time at stop codons, indicating that the ejection process might delay ribosome recycling. Proteins that have the highest positive charge density at their C-terminus are overwhelmingly ribosomal proteins, suggesting the possibility that this sequence feature may aid in the cotranslational assembly of ribosomes by delaying the release of nascent ribosomal proteins into the cytosol. Thus, nascent chain ejection times from the ribosome can vary greatly between proteins due to differential electrostatic interactions, can influence ribosome recycling, and could be particularly relevant to the synthesis and cotranslational behavior of some proteins.

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.

Folding perspectives of an intrinsically disordered transactivation domain and its single mutation breaking the folding propensity

Transcriptional regulation is a critical facet of cellular development controlled by numerous transcription factors, among which are E-proteins (E2A, HEB, and E2-2) that play important roles in lymphopoiesis. For example, primary hematopoietic cells immortalisation is promoted by interaction of the conserved PCET motif consisting of the Leu-X-X-Leu-Leu (LXXLL) and Leu-Asp-Phe-Ser (LDFS) sequences of the transactivation domains (AD1) of E-proteins with the KIX domain of CBP/p300 transcriptional co-activators. Earlier, it was shown that the LXXLL motif is essential for the PCET-KIX interaction driven by the PCET helical transition. In this study, we analyzed the dehydration-driven gain of helicity in the conserved region (residues 11-28) of the AD1 domain of E-protein. Particularly, we showed that AD1 structure was dramatically affected by alcohols, but was insensitive to changes in pH or the presence of osmolytes sarcosine and taurine, or high polyethylene glycol (PEG) concentrations and DOPC Liposomes. These structure-forming effects of solvents were almost completely absent in the case of L21P AD1 mutant characterized by weakened interaction with KIX. This indicates that KIX interaction-induced AD1 ordering is driven by PCET motif dehydration. The L21P mutation-caused loss of molecular recognition function of AD1 is due to the mutation-induced disruption of the AD1 helical propensity.

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.

Structural Origins of FRET-Observed Nascent Chain Compaction on the Ribosome

A fluorescence signal arising from a Förster resonance energy transfer process was used to monitor conformational changes of a domain within the E. coli protein HemK during its synthesis by the ribosome. An increase in fluorescence was observed to begin 10 s after translation was initiated, indicating the domain became more compact in size. Since fluorescence only reports a single value at each time point it contains very little information about the structural ensemble that gives rise to it. Here, we supplement this experimental information with coarse-grained simulations that describe protein conformations and transitions at a spatial resolution of 3.8 Å. We use these simulations to test three hypotheses that might explain the cause of domain compaction: (1) that poor solvent quality conditions drive the unfolded state to compact, (2) that a change in the dimension of the space the domain occupies upon moving outside the exit tunnel causes compaction, or (3) that domain folding causes compaction. We find that domain folding and dimensional collapse are both consistent with the experimental data, while poor-solvent collapse is inconsistent. We identify alternative dye labeling positions on HemK that upon fluorescence can differentiate between the domain folding and dimensional collapse mechanisms. Partial folding of domains has been observed in C-terminally truncated forms of proteins. Therefore, it is likely that the experimentally observed compact state is a partially folded intermediate consisting, according to our simulations, of the first three helices of the HemK N-terminal domain adopting a native, tertiary configuration. With these simulations we also identify the possible cotranslational folding pathways of HemK.

Salt Effects on Protein Folding Thermodynamics

Salts differ in their ability to stabilize protein conformations, thereby affecting the thermodynamics and kinetics of protein folding. We developed a coarse-grained protein model that can predict salt-induced changes in protein properties by using the transfer free-energy data of various chemical groups from water to salt solutions. Using this model and molecular dynamics simulations, we probed the effect of seven different salts on the folding thermodynamics of the DNA binding domain of lac repressor protein ( lac-DBD) and N-terminal domain of ribosomal protein (NTL9). We show that a salt can act as a protein stabilizing or destabilizing agent depending on the protein sequence and folded state topology. The computed thermodynamic properties, especially the m values for various salts, which reveal the relative ability of a salt to stabilize the protein folded state, are in quantitative agreement with the experimentally measured values. The computations show that the degree of protein compaction in the denatured ensemble strongly depends on the salt identity, and for the same variation in salt concentration, the compaction in the protein dimensions varies from ∼4% to ∼30% depending on the salt. The transition-state ensemble (TSE) of lac-DBD is homogeneous and polarized, while the TSE of NTL9 is heterogeneous and diffusive. Salts induce subtle structural changes in the TSE that are in agreement with Hammond's postulate. The barrier to protein folding tends to disappear in the presence of moderate concentrations (∼3-4 m) of strongly stabilizing salts.

Fourier Transform-Ion Mobility-Orbitrap Mass Spectrometer: A Next-Generation Instrument for Native Mass Spectrometry

A new instrument configuration for native ion mobility-mass spectrometry (IM-MS) is described. Macromolecule ions are generated by using a static ESI source coupled to an RF ion funnel, and these ions are then mobility and mass analyzed using a periodic focusing drift tube IM analyzer and an Orbitrap mass spectrometer. The instrument design retains the capabilities for first-principles determination of rotationally averaged ion-neutral collision cross sections and high-resolution measurements in both mobility and mass analysis modes for intact protein complexes. Operation in the IM mode utilizes FT-IMS modes (originally described by Knorr ( Knorr , F. J. Anal. Chem . 1985 , 57 ( 2 ), 402 - 406 )), which provides a means to overcome the inherent duty cycle mismatch for drift tube (DT)-IM and Orbitrap mass analysis. The performance of the native ESI-FT-DT-IM-Orbitrap MS instrument was evaluated using the protein complexes Gln K (MW 44 kDa) and streptavidin (MW 53 kDa) bound to small molecules (ADP and biotin, respectively) and transthyretin (MW 56 kDa) bound to thyroxine and zinc.

Dynamics of Disordered Proteins under Confinement: Memory Effects and Internal Friction

Many proteins are disordered under physiological conditions. How efficiently they can search for their cellular targets and how fast they can fold upon target binding is determined by their intrinsic dynamics, which have thus attracted much recent attention. Experiments and molecular simulations show that the inherent reconfiguration timescale for unfolded proteins has a solvent friction component and an internal friction component, and the microscopic origin of the latter, along with its proper mathematical description, has been a topic of considerable debate. Internal friction varies across different proteins of comparable length and increases with decreasing denaturant concentration, showing that it depends on how compact the protein is. Here we report on a systematic atomistic simulation study of how confinement, which induces a more compact unfolded state, affects dynamics and friction in disordered peptides. We find that the average reconfiguration timescales increase exponentially as the peptide's spatial dimensions are reduced; at the same time, confinement broadens the spectrum of relaxation timescales exhibited by the peptide. There are two important implications of this broadening: First, it limits applicability of the common Rouse and Zimm models with internal friction, as those models attempt to capture internal friction effects by introducing a single internal friction timescale. Second, the long-tailed distribution of relaxation times leads to anomalous diffusion effects in the dynamics of intramolecular distances. Analysis and interpretation of experimental signals from various measurements that probe intramolecular protein dynamics (such as single-molecule fluorescence correlation spectroscopy and single-molecule force spectroscopy) rely on the assumption of diffusive dynamics along the distances being probed; hence, our results suggest the need for more general models allowing for anomalous diffusion effects.

Perspective: Chain dynamics of unfolded and intrinsically disordered proteins from nanosecond fluorescence correlation spectroscopy combined with single-molecule FRET

The dynamics of unfolded proteins are important both for the process of protein folding and for the behavior of intrinsically disordered proteins. However, methods for investigating the global chain dynamics of these structurally diverse systems have been limited. A versatile experimental approach is single-molecule spectroscopy in combination with Förster resonance energy transfer and nanosecond fluorescence correlation spectroscopy. The concepts of polymer physics offer a powerful framework both for interpreting the results and for understanding and classifying the properties of unfolded and intrinsically disordered proteins. This information on long-range chain dynamics can be complemented with spectroscopic techniques that probe different length scales and time scales, and integration of these results greatly benefits from recent advances in molecular simulations. This increasing convergence between the experiment, theory, and simulation is thus starting to enable an increasingly detailed view of the dynamics of disordered proteins.

Crowding Stabilizes DMSO-Water Hydrogen-Bonding Interactions

Up to 40% of intracellular water is confined due to the dense packing of macromolecules, ions, and osmolytes. Despite the large body of work concerning the effect of additives on the biomolecular structure and stability, the role of crowding and heterogeneity in these interactions is not well understood. Here, infrared spectroscopy and molecular dynamics simulations are used to describe the mechanisms by which crowding modulates hydrogen bonding interactions between water and dimethyl sulfoxide (DMSO). Specifically, we use formamide and dimethylformamide (DMF) as molecular crowders and show that the S═O hydrogen bond populations in aqueous mixtures are increased by both amides. These additives increase the amount of water within the DMSO first solvation shell through two mechanisms: (a) directly stabilizing water-DMSO hydrogen bonds; (b) increasing water exposure by destabilizing DMSO-DMSO self-interactions. Further, we quantified the hydrogen bond enthalpies between the different components: DMSO-water (61 kJ/mol) > DMSO-formamide (32 kJ/mol) > water-water (23 kJ/mol) ≫ formamide-water (4.7 kJ/mol). Spectra of carbonyl stretching vibrations show that DMSO induces the dehydration of amides as a result of strong DMSO-water interactions, which has been suggested as the main mechanism of protein destabilization.

Denaturants Alter the Flux through Multiple Pathways in the Folding of PDZ Domain

Although we understand many aspects of how small proteins (number of residues less than about hundred) fold, it is a major challenge to quantitatively describe how large proteins self-assemble. To partially overcome this challenge, we performed simulations using the self-organized polymer model with side chains (SOP-SC) in guanidinium chloride (GdmCl), using the molecular transfer model (MTM), to describe the folding of the 110-residue PDZ3 domain. The simulations reproduce the folding thermodynamics accurately including the melting temperature (T_m), the stability of the folded state with respect to the unfolded state. We show that the calculated dependence of ln k_obs (k_obs is the relaxation rate) has the characteristic chevron shape. The slopes of the chevron plots are in good agreement with experiments. We show that PDZ3 folds by four major pathways populating two metastable intermediates, in accord with the kinetic partitioning mechanism. The structure of one of the intermediates, populated after polypeptide chain collapse, is structurally similar to an equilibrium intermediate. Surprisingly, the connectivities between the intermediates and hence, the fluxes through the pathways depend on the concentration of GdmCl. The results are used to predict possible outcomes for unfolding of PDZ domain subject to mechanical forces. Our study demonstrates that, irrespective of the size or topology, simulations based on MTM and SOP-SC offer a theoretical framework for describing the folding of proteins, mimicking precisely the conditions used in experiments.

Changing relations between proteins and osmolytes: a choice of nature

The stabilization and destabilization of the protein in the presence of any additive is mainly attributed to its preferential exclusion from protein surface and its preferential binding to the protein surface, respectively.

Globular Protein Folding In Vitro and In Vivo

In vitro, computational, and theoretical studies of protein folding have converged to paint a rich and complex energy landscape. This landscape is sensitively modulated by environmental conditions and subject to evolutionary pressure on protein function. Of these environments, none is more complex than the cell itself, where proteins function in the cytosol, in membranes, and in different compartments. A wide variety of kinetic and thermodynamics experiments, ranging from single-molecule studies to jump kinetics and from nuclear magnetic resonance to imaging on the microscope, have elucidated how protein energy landscapes facilitate folding and how they are subject to evolutionary constraints and environmental perturbation. Here we review some recent developments in the field and refer the reader to some original work and additional reviews that cover this broad topic in protein science.

Fast Protein Translation Can Promote Co- and Posttranslational Folding of Misfolding-Prone Proteins

Chemical kinetic modeling has previously been used to predict that fast-translating codons can enhance cotranslational protein folding by helping to avoid misfolded intermediates. Consistent with this prediction, protein aggregation in yeast and worms was observed to increase when translation was globally slowed down, possibly due to increased cotranslational misfolding. Observation of similar behavior in molecular simulations would confirm predictions from the simpler chemical kinetic model and provide a molecular perspective on cotranslational folding, misfolding, and the impact of translation speed on these processes. All-atom simulations cannot reach the timescales relevant to protein synthesis, and most conventional structure-based coarse-grained models do not allow for nonnative structure formation. Here, we introduce a protocol to incorporate misfolding using the functional forms of publicly available force fields. With this model we create two artificial proteins that are capable of undergoing structural transitions between a native and a misfolded conformation and simulate their synthesis by the ribosome. Consistent with the chemical kinetic predictions, we find that rapid synthesis of misfolding-prone nascent-chain segments increases the fraction of folded proteins by kinetically partitioning more molecules through on-pathway intermediates, decreasing the likelihood of sampling misfolded conformations. Novel to this study, to our knowledge, we observe that differences in protein dynamics, arising from different translation-elongation schedules, can persist long after the nascent protein has been released from the ribosome, and that a sufficient level of energetic frustration is needed for fast-translating codons to be beneficial for folding. These results provide further evidence that fast-translating codons can be as biologically important as pause sites in coordinating cotranslational folding.

Exploring the Denatured State Ensemble by Single-Molecule Chemo-Mechanical Unfolding: The Effect of Force, Temperature, and Urea

While it is widely appreciated that the denatured state of a protein is a heterogeneous conformational ensemble, there is still debate over how this ensemble changes with environmental conditions. Here, we use single-molecule chemo-mechanical unfolding, which combines force and urea using the optical tweezers, together with traditional protein unfolding studies to explore how perturbants commonly used to unfold proteins (urea, force, and temperature) affect the denatured-state ensemble. We compare the urea m-values, which report on the change in solvent accessible surface area for unfolding, to probe the denatured state as a function of force, temperature, and urea. We find that while the urea- and force-induced denatured states expose similar amounts of surface area, the denatured state at high temperature and low urea concentration is more compact. To disentangle these two effects, we use destabilizing mutations that shift the T_m and C_m. We find that the compaction of the denatured state is related to changing temperature as the different variants of acyl-coenzyme A binding protein have similar m-values when they are at the same temperature but different urea concentration. These results have important implications for protein folding and stability under different environmental conditions.

Dimension conversion and scaling of disordered protein chains

To extract protein dimension and energetics information from single-molecule fluorescence resonance energy transfer spectroscopy (smFRET) data, it is essential to establish the relationship between the distributions of the radius of gyration (Rg) and the end-to-end (donor-to-acceptor) distance (Ree). Here, we performed a coarse-grained molecular dynamics simulation to obtain a conformational ensemble of denatured proteins and intrinsically disordered proteins. For any disordered chain with fixed length, there is an excellent linear correlation between the average values of Rg and Ree under various solvent conditions, but the relationship deviates from the prediction of a Gaussian chain. A modified conversion formula was proposed to analyze smFRET data. The formula reduces the discrepancy between the results obtained from FRET and small-angle X-ray scattering (SAXS). The scaling law in a coil-globule transition process was examined where a significant finite-size effect was revealed, i.e., the scaling exponent may exceed the theoretical critical boundary [1/3, 3/5] and the prefactor changes notably during the transition. The Sanchez chain model was also tested and it was shown that the mean-field approximation works well for expanded chains.

Decoupling of size and shape fluctuations in heteropolymeric sequences reconciles discrepancies in SAXS vs. FRET measurements

Unfolded states of proteins and native states of intrinsically disordered proteins (IDPs) populate heterogeneous conformational ensembles in solution. The average sizes of these heterogeneous systems, quantified by the radius of gyration (R_G ), can be measured by small-angle X-ray scattering (SAXS). Another parameter, the mean dye-to-dye distance (R_E ) for proteins with fluorescently labeled termini, can be estimated using single-molecule Förster resonance energy transfer (smFRET). A number of studies have reported inconsistencies in inferences drawn from the two sets of measurements for the dimensions of unfolded proteins and IDPs in the absence of chemical denaturants. These differences are typically attributed to the influence of fluorescent labels used in smFRET and to the impact of high concentrations and averaging features of SAXS. By measuring the dimensions of a collection of labeled and unlabeled polypeptides using smFRET and SAXS, we directly assessed the contributions of dyes to the experimental values R_G and R_E For chemically denatured proteins we obtain mutual consistency in our inferences based on R_G and R_E , whereas for IDPs under native conditions, we find substantial deviations. Using computations, we show that discrepant inferences are neither due to methodological shortcomings of specific measurements nor due to artifacts of dyes. Instead, our analysis suggests that chemical heterogeneity in heteropolymeric systems leads to a decoupling between R_E and R_G that is amplified in the absence of denaturants. Therefore, joint assessments of R_G and R_E combined with measurements of polymer shapes should provide a consistent and complete picture of the underlying ensembles.

Preferential binding effects on protein structure and dynamics revealed by coarse-grained Monte Carlo simulation

The effect of preferential binding of solute molecules within an aqueous solution on the structure and dynamics of the histone H3.1 protein is examined by a coarse-grained Monte Carlo simulation. The knowledge-based residue-residue and hydropathy-index-based residue-solvent interactions are used as input to analyze a number of local and global physical quantities as a function of the residue-solvent interaction strength (f). Results from simulations that treat the aqueous solution as a homogeneous effective solvent medium are compared to when positional fluctuations of the solute molecules are explicitly considered. While the radius of gyration (Rg) of the protein exhibits a non-monotonic dependence on solvent interaction over a wide range of f within an effective medium, an abrupt collapse in Rg occurs in a narrow range of f when solute molecules rapidly bind to a preferential set of sites on the protein. The structure factor S(q) of the protein with wave vector (q) becomes oscillatory in the collapsed state, which reflects segmental correlations caused by spatial fluctuations in solute-protein binding. Spatial fluctuations in solute binding also modify the effective dimension (D) of the protein in fibrous (D ∼ 1.3), random-coil (D ∼ 1.75), and globular (D ∼ 3) conformational ensembles as the interaction strength increases, which differ from an effective medium with respect to the magnitude of D and the length scale.

Evidence for the principle of minimal frustration in the evolution of protein folding landscapes

Theoretical and experimental studies have firmly established that protein folding can be described by a funneled energy landscape. This funneled energy landscape is the result of foldable protein sequences evolving following the principle of minimal frustration, which allows proteins to rapidly fold to their native biologically functional conformations. For a protein family with a given functional fold, the principle of minimal frustration suggests that, independent of sequence, all proteins within this family should fold with similar rates. However, depending on the optimal living temperature of the organism, proteins also need to modulate their thermodynamic stability. Consequently, the difference in thermodynamic stability should be primarily caused by differences in the unfolding rates. To test this hypothesis experimentally, we performed comprehensive thermodynamic and kinetic analyses of 15 different proteins from the thioredoxin family. Eight of these thioredoxins were extant proteins from psychrophilic, mesophilic, or thermophilic organisms. The other seven protein sequences were obtained using ancestral sequence reconstruction and can be dated back over 4 billion years. We found that all studied proteins fold with very similar rates but unfold with rates that differ up to three orders of magnitude. The unfolding rates correlate well with the thermodynamic stability of the proteins. Moreover, proteins that unfold slower are more resistant to proteolysis. These results provide direct experimental support to the principle of minimal frustration hypothesis.